Pair dispersion in inhomogeneous turbulent thermal convection
Olivier Liot David-Martin-Calle Amélie Gay Julien Salort Francesca Chillà Mickaël Bourgoin























 []olivier.liot@imft.fr
 [Presently at ]Institut de Mécanique des Fluides de Toulouse, Toulouse, France
 
 [Presently at ]Institut Lumière Matière, CNRS, 8 rue Ada Byron, Villeurbanne, France

[Presently at ]IRPHE - UMR 7342, Technopôle de Château-Gombert, 49, rue Joliot Curie - B.P. 146, 13384 Marseille Cedex 13, France



Univ Lyon, ENS de Lyon, Univ Claude Bernard, CNRS, Laboratoire de Physique, F-69342 Lyon Cedex 7, France
















Due to large scale flow inhomogeneities and the effects of temperature, turbulence small-scale structure in thermal convection is still an active field of investigation, especially considering sophisticated Lagrangian statistics. 
Here we experimentally study Lagrangian pair dispersion (one of the canonical problems of Lagrangian turbulence) in a Rayleigh-Bénard convection cell. 
A sufficiently high temperature difference is imposed on a horizontal layer of fluid to observe a turbulent flow. 
We perform Lagrangian tracking of sub-millimetric particles on a large measurement volume including part of the Large Scale Circulation (LSC) revealing some large inhomogeneities. 
Our study brings to light several new insights regarding our understanding of turbulent thermal convection: 
(i) by decomposing particle Lagrangian dynamics into the LSC contribution and the turbulent fluctuations, we highlight the relative impact of both contributions on pair dispersion; 
(ii) using the same decomposition, we estimate the Eulerian second-order velocity structure functions from pair statistics and show that after removing the LSC contribution, the remaining statistics recover usual homogeneous and isotropic behaviours which are governed by a local energy dissipation rate to be distinguished from the global dissipation rate classically used to characterise turbulence in thermal convection; 
and (iii) we revisit the super-diffusive Richardson-Obukhov regime of particle dispersion and propose a refined estimate of the Richardson constant.


















Pair dispersion in inhomogeneous turbulent thermal convection
    Mickaël Bourgoin
    December 30, 2023
=============================================================








§ INTRODUCTION
 

Various natural and industrial flows are buoyancy-driven. 
Atmospheric or oceanic dynamics, Earth's mantle flows and heat exchangers are some examples where the natural thermal convection – without mechanical forcing – has a dramatic influence on the flow. 

The main model system is the Rayleigh-Bénard cell. 
A horizontal layer of fluid, confined in a cell is cooled from above and heated from below. 
If the thermal forcing is large enough, the flow is turbulent and the fluid is well mixed. 
A consequence is that the temperature is nearly homogeneous on average, and thermal gradients are confined in boundary layers, close to the horizontal plates. 
Boundary layers are thermally unstable, which leads to the emission of coherent structures, named thermal plumes. 
Despite the significant progress of the modelling of turbulent thermal convection in the past decade <cit.>, the link between these structures and the global heat flux is still not fully understood. 
Examples of open questions are the dynamics of the turbulent plumes <cit.>, the interactions between the mean flow and the boundary layers, and the influence of the inhomogeneity on the heat transport.

The system is controlled by two non dimensional parameters. 
(i) The Rayleigh number, Ra, is the balance of the buoyancy effects and the dissipative ones. 
It accounts for the thermal forcing,

    Ra=gαΔ T H^3/νκ,

where g is the acceleration due to gravity, Δ T is the temperature difference between the cooling and the heating plates, H is their separation distance, α represents the thermal expansion coefficient of the fluid, ν the fluid's kinematic viscosity and κ its thermal diffusivity. 
(ii) The Prandtl number, Pr, compares the two dissipative processes (thermal and viscous diffusion): 

    Pr=ν/κ.

The response of the system is estimated by the Nusselt number, Nu, which compares the heat flux through the convection cell to the purely diffusive one,

    Nu=Q H/λΔ T,

where Q is the heat flux and λ the thermal conductivity of the fluid.

High-resolution spatio-temporal Lagrangian measurements in turbulent flows are now possible thanks to improved digital imaging techniques and computing resources <cit.>. 
Flows with important mixing and transport properties deserve to be studied with a Lagrangian approach <cit.>.
Moreover, the Lagrangian point of view is relevant to stochastic models which describe some aspects of the turbulence, such as finite-Reynolds effects <cit.>, and intermittency <cit.>. 
Important transport properties of turbulent flows are related to multi-particle dispersion <cit.>, among which pair dispersion is the most fundamental and has been pioneered by Sir Richardson <cit.>. 
The way two particles go away from each other has been largely studied in turbulent flows simultaneously experimentally, numerically and theoretically (see e.g. the review articles by Sawford <cit.> and Salazar et al. <cit.>). 
It has important applications including oceanic plankton or atmospheric pollutant dispersion <cit.> for which thermal convection is a key ingredient. 
Three regimes can be distinguished. 
(i) For time scale shorter than t^*=(Δ_0^2/ϵ)^2/3, where Δ_0 is the initial separation, as represented in Fig. <ref>, and ϵ the kinetic energy dissipation rate, a ballistic regime is observed. 
It is called the Batchelor regime, and the squared pair dispersion D^2_Δ_0(t) (as defined in Eq. <ref>) scales as t^2 <cit.>. 
(ii) For t^*<t<T_L, where T_L is the Lagrangian correlation time, the super-diffusive or Richardson-Obukhov regime appears and D^2_Δ_0(t)∝ t^3 <cit.>. 
(iii) For t>T_L, a diffusive regime is reached and D^2_Δ_0(t)∝ t <cit.>.

Nevertheless these three regimes are poorly investigated in turbulent thermal convection, contrary to isothermal turbulence. 
Indeed, the presence of thermal plumes and the strong inhomogeneities of the flow should lead to turbulent statistics far from the usual Homogeneous and Isotropic Turbulence (HIT) framework. 
To our knowledge, few experimental works were reported.
The first one used a single large particle with embedded thermistors to combine Lagrangian temperature and velocity measurements <cit.>. 
The second one <cit.> used sub-millimeter particles, focusing on the very center of convection cells where the flow is quite homogeneous and isotropic. 
Recently we reported on a Lagrangian study of the velocity and acceleration statistics in a large measurement volume of a Rayleigh-Bénard cell where the mean flow is highly inhomogeneous <cit.>. 
Very recently, a study focused on the aspect ratio influence on Lagrangian statistics <cit.>. 
Some experimental studies have been performed in rotating Rayleigh-Bénard turbulence <cit.>. 
Moreover, heat flux and particle dispersion were numerically studied by Schumacher and co-workers all along the decade. 
Among these previous studies, the pair dispersion was explored numerically in thermally-driven flows <cit.>. 
The Batchelor and Richardson-Obukhov regimes were observed. 
In experimental studies, mostly the case of initial separations close to the dissipative length-scale was investigated <cit.>, with some evidence of the Richardson regime.

The purpose of this paper is to explore how flow inhomogeneity observed in a wide central zone of a thermal flow affects the turbulence. 
We deliberately place this study out of the HIT scope and address the role of the Large Scale Circulation (LSC) on Lagrangian transport. 
The Lagrangian tracking of particle pairs in a turbulent Rayleigh-Bénard cell is presented. 
A significant part of the LSC is observed thanks to the size of the measurement volume. 
The growth of the mean-square separation of particles is analysed. 
The main aspect of this article concerns the influence and the modelling of the inhomogeneous mean flow on pair dispersion. 
We also discuss the estimation of the Richardson constant in the Richardson-Obukhov regime.




§ EXPERIMENTAL SETUP AND MEASUREMENT TECHNIQUES
 


The convection cell consists of an octagonal-shaped setup with eight
transparent polymethylmetacrylate vertical walls. 
The two horizontal plates are made of anodised aluminium. 
They are 40cm in diameter, tangent to the octagon made by the walls and vertically spaced 30cm apart. 
A  custom-made 40cm in diameter spiral electrical resistance imposes the heat flux from the bottom while the top plate temperature is regulated by a glycol circulation pump. 
The temperature of the plates is monitored using eight PT100 temperature sensors – four inside each plate. 
The working fluid is deionised and degased water, with a density ρ_f = 0.99. 
The imposed heat power is 800W, the mean temperature is fixed to 40 and the temperature difference reaches Δ T=19.2K. 
The corresponding control parameter values are Ra=2.0×10^10 and Pr=4.3. 
The Nusselt number is consistent with the Grossmann-Lohse theory <cit.>, and is Nu=139 (after removing 15 % of thermal losses). 
In these conditions the flow is turbulent. 
The subsequent Large Scale Circulation (LSC) consists in a convection roll confined between two diametrically opposite sides of the octagonal cell. 
We observe using shadowgraphs that the convection roll orientation can change spontaneously with a typical time of few hours. 
Nevertheless, we do not observe such reversals on the time scale of our acquisitions (180s, which is confirmed by the shape-similar velocity distributions of each run <cit.>. 
However, LSC sloshing in the x-direction is observed during the acquisitions <cit.>.





The flow is seeded with thousands of polystyrene particles. 
Their diameter is d=250 and their density is ρ_p=1.03. 
The Stokes number in a turbulent flow can be estimated as <cit.>:

    St=ρ_p/18ρ_f(d/η)^2.
 
In our case we get St≈ 7× 10^-3, which means that particles can be considered as tracers.
Nevertheless we observe a decrease of suspended particle concentration on a typical time of few dozen minutes, probably due to the entrapment by the thermal boundary layers. 
This problem is solved by splitting the acquisitions into six runs with re-seeding between each one. 

The fluid volume is illuminated by two sets of four vertical LED bars of 864lumens, put in front of two different vertical walls, as shown in Fig. <ref>.
Moreover, to improve the particle visibility, the vertical wall on the far side of the tank from where the camera is situated is covered with black paper. 
Thus only three faces are still available to place the cameras. 
To avoid any angle between the camera lens and the wall and to have the largest volume captured by all the cameras, we choose to use only one camera per free face. 
Consequently, three cameras are positioned in the horizontal plane situated at mid-height of the walls, at polar angles ψ=0, 135 and 270, as illustrated in Fig. <ref>).
The measurement volume where we perform 3D Lagrangian tracking is about 11cm per side, 17cm high and is centered in the cell (see Fig. <ref> and <ref>). 
The cameras have a resolution of 1088× 2048 pixels^2. 
The maximum resolution is 6pixels per Kolmogorov length, defined as η = (ν^3/⟨ϵ⟩)^1/4 where ⟨ϵ⟩ is the average mass rate of kinetic energy dissipation, estimated as <cit.>:

    ⟨ϵ⟩ = ν^3/H^4RaPr^-2(Nu-1).




The acquisitions are split in six independant runs of 180 s. 
The total measurement time is about 2880 times the Kolmogorov time – τ_η = √(ν/⟨ϵ⟩)≈0.36 s – and 575 times the free-fall time – τ_f=√(H/gαΔ T).
We can therefore consider that the flow is fairly averaged as long as no LSC orientation change occurs.
The sampling frequency is fixed to 200Hz (about 70 times the dissipative time scale), which ensures a good resolution of small scale Lagrangian dynamics with sufficient oversampling to filter noise on reconnected tracks <cit.>. 
A Gaussian kernel is used to filter the trajectories and their temporal derivatives; its width (0.3 τ_η) does not affect the smallest turbulent scales.

The measurement volume is wide enough to capture a significant part of the LSC. 
Figure <ref> presents the part of the convection roll captured by the measurement volume within the convection cell. 
To visualise the inhomogeneity of the flow we compute pseudo-Eulerian maps from the Lagrangian data. 
They are averaged over a spatial grid, and the resulting velocity field is interpolated. 
We call v^E_k (k=x,y,z) the resulting average velocity. 
Figures <ref> (a) and (b) show three-dimensional representations of the average velocity vectors field, with two distinct viewing angles. 
In the y-direction, we observe a very inhomogeneous flow with large positive velocities at the top of the measurement volume, large negative ones at the bottom and null velocities at the centre. 
The other horizontal average velocity component (v_x^E) and the vertical one (v_z^E) are both nearly null, as highlighted in Fig. <ref> (b) and (a) respectively. 










§ PAIR DISPERSION AND INHOMOGENEITY
 

 We address the question of the impact of large scale inhomogeneities on pair dispersion. 
 Figure <ref> shows the principle of the pair dispersion: we study the evolution of the separation between pairs of particles with time by defining

    D^2_Δ_0(t)     =    ⟨(Δ(t) - Δ_0)^2 ⟩_0, 
    
    R^2_Δ_0,k(t)     =    ⟨(Δ_k(t) - Δ_0,k)^2 ⟩_0,

where Δ(t) and Δ_0 are respectively the vectors connecting two particles at times t and t_0 and Δ_k and Δ_0,k their projections along the k axis (k = x, y, z); ⟨·⟩_0 represents the statistical mean over pairs with an initial separation Δ_0 (see Fig. <ref>). 



The leading term of the Taylor development at short time of Eq. <ref> leads to a ballisitic (or Batchelor <cit.>) regime: 

    D^2_Δ_0 =  S^2_v(Δ_0) t^2 + 𝒪(t^3)

where S^2_v(Δ_0) is the second-order Eulerian velocity structure function for a separation Δ_0. 
This regime is expected to last for times t≲ t^*=(Δ_0^2/⟨ϵ⟩)^1/3. 
For time larger than t^* the super-diffusive (or Richardson-Obukhov <cit.>) regime appears <cit.>. 
It corresponds to a t^3 dependence of D^2_Δ_0(t) and was observed in turbulent convection both numerically <cit.> and experimentally <cit.>. 

At a given ⟨ϵ⟩, the observation of the super-diffusive regime requires several conditions. 
(i) t^* must be significantly smaller than the Lagrangian correlation time T_L, otherwise the ballistic regime transitions directly to the diffusive regime. 
(ii) Experimental tracks of particle pairs must be longer than t^*. 
Both conditions imply that the observation of the Richardson-Obukhov regime is most favourable for small initial separations Δ_0, as confirmed in numerical simulations for both isothermal <cit.> and thermal <cit.> turbulence.  
Experimentally, in high Reynolds turbulence, this regime is hard to observe <cit.>. 
This is mostly due to the difficulty to access long tracks and to have good statistical convergence for small initial separations (which require high seeding densities making the tracking more complex and noisy). 
In thermal convection, Ni & Xia <cit.> used very small initial separations (0.9–1.3 η) in a small measurement volume and observed a fleeting super-diffusive regime. 
In our case we choose to track long trajectories over a large measurement volume (up to 30 s corresponding to 85 τ_η) but with initial separations starting from 2.7 η (η≈0.7 mm).

We present in Fig. <ref> (a) and (b) the mean-square separation D^2_Δ_0 versus t/t^* for short (range 1.9–2.9 mm) and larger (range 4.3 mm–80.3 mm) initial separations respectively. 
In terms of Kolmogorov scale, the ranges are respectively [2.7 η–4.0 η] and [6.0 η–114 η].
For all initial separations we observe the early ballistic (t^2) regime. 
For the shortest separations (Fig. <ref> (a)), the trajectories are long enough compared to t^* to observe the transition towards the super-diffusive (t^3) regime. 
This regime will be discussed further in this paper. 
For the larger initial separations (Fig. <ref> (b)), the Richardson-Obukhov regime could not be reached. 









As described by Eq <ref>, the short term pair dispersion is dominated by the t^2-(ballistic) regime. 
The pre-factor is given by the second-order Eulerian structure function S_v which in non-homogeneous thermal flows, embeds both the statistics of the turbulent spatial fluctuations and the spatial inhomogeneities of the LSC. 
We compute the structure functions with a Δ_0 separation for each velocity component from: 

    R^2_Δ_0,k(t) = S^2_v_k(Δ_0) t^2 + 𝒪(t^3),

where S^2_v_k(Δ_0) corresponds to the Eulerian second-order longitudinal structure function of the k-component of the velocity for an initial separation Δ_0. 

Figure <ref> shows these structure functions compensated by Δ_0^2/3 following classic HIT scaling. 
We observe a strong anisotropy. 
S^2_v_x(Δ_0) and S^2_v_z(Δ_0) present the same plateau which is in agreement with the Kolmogorov scaling, S^2_v_k ∝ Δ_0^2/3 <cit.>. 
From this plateau (1.8± 0.1mm^4/3/s^2) we can estimate the local average kinetic energy dissipation rate ⟨ϵ⟩_loc in the sub-domain delimited by the measurement volume – instead of ⟨ϵ⟩ for the whole cell computed using Eq <ref>. 
Indeed, in the Kolmogorov theory we have S^2_v_k=(11/9)C_2⟨ϵ⟩_loc^2/3 with C_2 the Kolmogorov constant <cit.>. 
In our case of weak turbulence – the Reynolds number based on the Taylor microscale reaches R_λ≈ 75, we take C_2=1.7 <cit.>. 
We obtain ⟨ϵ⟩_loc=0.8±0.05×10^-6 m^2/s^3. 
We can compare this value to the estimation of ⟨ϵ⟩ in the whole cell using Eq <ref>, corrected by inhomogeneity effects. 
Kunnen et al. <cit.> performed numerical simulations in a cylindrical cell for Pr=6.4 and Ra=1×10^9. 
From this work we estimate that in the center of the cell ⟨ϵ⟩_loc is about 20 % to 30 % of its mean value calculated with Eq. <ref>.
Assuming these corrections we choose ⟨ϵ⟩_loc = (0.25±0.05) ⟨ϵ⟩ and we finally have ⟨ϵ⟩_loc = 1.3±0.3×10^-6 m^2/s^3, which is quite consistent with the experimental deductions. 
Note that the kinetic energy dissipation rate value used to estimate t^* is computed as 25% of the value given by Eq <ref> – ⟨ϵ⟩_loc≈ 1.3×10^-6m^2/s^3.

S^2_v_y(Δ_0) is dramatically different. 
It does not match with the Kolmogorov theory at all, but can be explained by the influence of the mean flow. 
First we decompose the velocity components, in a manner similar to previous Lagrangian works in von Kármán and thermal turbulent flows <cit.>. 
The flow velocity at a given time t and position (x,y,z) can be seen as the superposition of the time average velocity at this position and a time-dependent fluctuation. 
Since we use tracers, their velocity at a given time and position can be decomposed in the same way. 
Practically, the local time average velocity corresponds to the pseudo-Eulerian velocity v_k^E whose computation is explained at the end of section <ref>. 
The Lagrangian fluctuation v_k'(t) is the difference between actual particle velocity and pseudo-Eulerian velocity at the particle position. 
This can be formulated as:

    v_k(t)=v^E_k(x(t),y(t),z(t))+v'_k(t),

where k=x,y,z. Using this decomposition, we can develop the structure function S^2_v_y as:


    S_v_y^2(Δ_0)     =    ⟨(v_y^E(r+Δ_0)-v_y^E(r))^2⟩_r+⟨(v'_y(r+Δ_0)-v'_y(r))^2⟩_r
           +2 ⟨(v_y^E(r+Δ_0-v_y^E(r))(v'_y(r+Δ_0)-v'_y(r))⟩_r,


where r is a position in space. 
Since S^2_v_y(Δ_0) is an Eulerian quantity, ⟨·⟩_r represents the spatial average over all the accessible positions. 
The third term of Eq. <ref> is a cross-correlation term between the local mean flow and the fluctuations. 
In a previous study <cit.>, a sensor-embedded particle was used to explore the flow in a parallelepipedic cell for similar Pr and Ra. 
We observed that the correlations between the mean flow and the Lagrangian fluctuations are very small compared to the auto-correlations of the mean flow and the fluctuations. 
Consequently we neglect the third term. 
The second term corresponds to the second order Eulerian structure function of the fluctuations. 
We have seen (section II), for k=x, z that v_k^E≈ 0, which means that v_k≈ v'_k according to Eq <ref>.
Consequently, S^2_v_k≈⟨(v'_k(r+Δ_0)-v'_k(r))^2⟩_r. 
We also assume isotropy of turbulent fluctuations as supported by our recent investigation of single particle statistics <cit.>. 
The large similarity of S^2_v_x and S^2_v_z in Fig <ref> reinforces the hypothesis for two-particle statistics. 
From this observation, we assume that the second term of S^2_v_y from Eq. <ref> is similar to S^2_v_x:


    ⟨(v'_y(r+Δ_0)-v'_y(r))^2⟩_r≈⟨(v'_x(r+Δ_0)-v'_x(r))^2⟩_r≈ S^2_v_x.


 Finally, the first term of Eq. <ref> is related to the mean flow structure. 
 As we observe in Fig <ref>, inside our measurement volume the mean flow is mostly a shear flow in the {y⃗,z⃗} planes for every x, with zero-velocity in the center of the cell. 
 Thus we have v_y^E(r)=β z (within a constant) where β is a shear rate defined from the y-component velocity gradient in the z direction (which is about uniform). 
 We estimate it to β≈ 0.14 s^-1 from mean vertical velocity profile. 
 The mean flow structure is sketched in Fig. <ref>. 
 Then the first term in Eq. <ref> can be written as: 

    ⟨(v_y^E(r+Δ_0))-v_y^E(r))^2⟩_r = β^2 Δ_0,z^2.

Furthermore, the initial pair separations are assumed to have three similar Δ_0,k in order to avoid considering particles with a very large separation in one direction and a very short separation in an other one. 
Consequently we have Δ_0≈√(3)Δ_0,z. 
Finally we can write: 

    S_v_y^2(Δ_0)≈β^2(Δ_0^2/3) + S_v_x^2(Δ_0).




The black line in Fig. <ref> compares this no-free-parameter model to experimental data. 
It is valid while the three components of the initial separation are similar. 
Since the volume is larger in the z⃗ direction than in the horizontal ones, this condition is violated for Δ_0 ≳ 50 mm. 
The solid part of the line corresponds to the zone where the model is expected to be valid, whereas the dashed part, which departs from experimental data, corresponds to a range of initial separations where the condition Δ_0≈√(3)Δ_0,k is not satisfied anymore. 
Fig. <ref> and its insert show the same experimental results, except for S^2_v_y.
In the insert, the mean-flow component in Eq. <ref>, β^2Δ_0^2/3 is subtracted from  experimental data S^2_v_y. 
This leads to a collapse of the curve on the S^2_v_x and S^2_v_z on the model validity domain. 
This collapse shows that the Eulerian second-order structure function S^2_v_y recovers the HIT behaviour if the mean flow influence is removed.



§ ABOUT THE RICHARDSON-OBUKHOV REGIME


The observation of the Richardson-Obukhov regime is experimentally subtle in turbulent convection. 
Moreover, the super-diffusive regime could include extra diffusion due to the LSC-generated shear rate. 
For short initial separations, we observe in Fig. <ref> (a) a transition towards a t^3 regime similar to that of Richardson-Obukhov. 
In this regime, the expected pair separation expression is: 

    D^2_Δ_0=g⟨ϵ⟩_loc t^3,

where g is called the Richardson constant <cit.>. 
In homogeneous isotropic turbulence (HIT), the expected value is g∈[0.65-0.7] in the range of Reynolds number R_λ corresponding to our experiment <cit.>. 
Figure <ref> shows the pair separation compensated by ⟨ϵ⟩_loc t^3 for the smallest initial separations (corresponding to Fig. <ref> (a)). 
A plateau is well defined for t between 0.7 t^* and 1.1 t^* with g∈ [2-4]. 
This range is consistent although a bit larger than HIT value. 
However it is significantly larger than values reported in turbulent convection by Ni & Xia <cit.> (g≈ 0.1 for Δ_0∈ [0.9 η-1.3 η]). 
Moreover, the plateau is observed for lower t/t^* than in Ni & Xia. 
Aditionally we find a systematic dependency on the plateau in Fig. <ref> with the initial separation, Δ_0.
 


 Both observations can be explained from the peculiar behaviour of pair dispersion for initial separation close to the dissipative scale, and by the value of the average kinetic energy dissipation rate used. 
 The shift in the time value for the occurrence of the plateau between our data and that of Ni & Xia can be explained by our use of ⟨ϵ⟩_loc to compute t^* which leads to a higher t^* than when using ⟨ϵ⟩. 
 Concerning the plateau value, numerical simulations performed by Boffetta & Sokolov <cit.> and Sawford et al. <cit.> show that for initial separations of the order η, the t^3 Richardson-Obukhov regime is preceded by a local minimum, leading to an apparent short t^3 lower plateau. 
 This phenomenon disappears for higher Δ_0. 
 This is attributed to a contamination of the initial range by dissipation effects <cit.>. 
 Based on the ballistic cascade model proposed by Bourgoin <cit.>, we plot in Fig. <ref> the pair separation compensated by (ϵ t^3) in HIT for initial separations between η and 10 η. 
 The local minimum zone before the super-diffusive plateau is clearly visible for an initial separation Δ_0∼η. 
 A progressive disappearance of the minimum zone is observed as Δ_0 increases. 
 This highlights that tracks too short can lead to a biased (misleadingly too small) estimation of the Richardson constant, misleadingly taken as the apparent plateau observed for small Δ_0 near the local minimum between ballistic and super-diffusive zone in the D^2_Δ_0/(⟨ϵ⟩_loct^3) plot. 
 In this zone the apparent plateau also leads to an erroneous dependency of g on Δ_0. 
 Figure <ref> suggests that the plateau due to the local minimum zone is nearly level with the actual Richardson plateau for Δ_0≳ 4 η. 
 These observations very likely explain the apparent initial separation dependency observed in Fig. <ref> and also reported in previous work by Ni & Xia <cit.>. 
 As in both studies, due to the limited track length, it is likely that the local minimum zone is explored, rather than the actual plateau of the Richardson regime. 
 Ni & Xia proposed an estimation of the Richardson constant based on the smallest initial separation they had (Δ_0≲η), for which the plateau related to the local minimum zone significantly under-estimates the actual value of g. 
 In this previous pioneering study, the tracks for an initial separation larger than η were too short and not usable to estimate g.
 In our study, the tracks of pairs with initial separations in the range [2.7 η-4 η] are marginally longer and allow for computation of the plateau due to the local minimum zone.
 Since we have larger initial separations, the plateau related to the local minimum zone is naturally higher (see Fig. <ref>) than for Ni & Xia. 
 Based on the previous discussion, the plateau around the minimum zone roughly levels to the actual Richardson plateau for initial separations close to 4η.
 Considering the range of initial separations in Fig. <ref> ([2.7 η-4 η]) we can therefore expect that the observed plateau leads to a reasonable estimate of the actual Richardson constant, hence g ∈ [2-4].
 Note that the use of ⟨ϵ⟩ instead of ⟨ϵ⟩_loc also further under-estimates the plateau value in the work by Ni & Xia <cit.>.



 

§ DISCUSSION AND CONCLUSION
  
 
 To explore the influence of inhomogeneity and ani­sotropy on turbulent statistics, an experimental study of pair dispersion in a turbulent thermal flow was performed. 
 Our experimental setup and analysis/post-processing tools allow us to obtain long trajectories, compared to the Kolmogorov time, and exceeding Batchelor time t^* in some cases.
 
 The quantitative analysis of the ballistic regime of pair dispersion highlights the influence of the Large Scale Circulation on the turbulent transport and gives a way to measure the mean kinetic energy dissipation rate in the considered volume. 
 Given its spatial inhomogeneity, a measurement of this rate is necessary to compare thermal convection turbulence to HIT. 
 This assessment is in good agreement with an estimation from global and Eulerian approaches. 
 We also used pair dispersion to access Eulerian velocity structure functions using only particle displacement (without needing to derive their velocity). 
 This is a way to study the influence of the inhomogeneity on thermal flow turbulence. 
 The particle dispersion in the y⃗ direction is highly influenced by the convection roll. 
 This is visible on the Eulerian velocity structure function which departs from the Δ_0^2/3 Kolmogorov scaling. 
 We proposed a model to describe the shape of this structure function which mixes a phenomenological approach of the mean flow, and experimental data from the other horizontal velocity structure functions, not affected by the mean flow. 
 This last point is useful to take into account the shape of the velocity structure function at low initial separations before reaching the Kolmogorov scaling. 
 With this choice, the non-negligible viscous dissipation for initial separations close to the Kolmogorov scale is considered. 
 This model is in good agreement with experimental data, except for large initial separations, due to the loss of the hypothesis of equipartition for initial separation components. 
 While turbulent convection is intrinsically inhomogeneous and anisotropic, this approach demonstrates that we can recover statistics from usual HIT simply by removing the mean flow.
 This is not trivial, especially for two-particle statistics. 
 This also shows that there is no temperature influence on the turbulence organisation, which is important for understanding the temperature role in thermal turbulent flows. 
 Our experiments validate this point in the center part of the convection cell, where plumes are scarce. 
 
 For the smallest initial separations, we obtain trajectories longer than t^* and we are able to observe the transition from ballistic to super-diffusive regime with a Richardson constant comparable but larger than HIT. 
 For short initial separations, this non-trivial transition reveals some complex behaviours even in HIT. 
 Moreover, the extra diffusion due to LSC-generated shear rate could affect the observations of the super-diffusive regime. 
 The difference with other experimental results could be explained from a well-predicted short plateau due to a local minimum zone which appears between ballistic and super-diffusive zones at short initial separations on compensated plots.
 
To summarise, our study shines a light on three new physical insights about turbulent thermal convection.  Using a decomposition between Large Scale Circulation and turbulent fluctuations, we can compute the relative impact of each contribution to pair dispersion. Then, we compute Eulerian second-order velocity structure functions from pair separations. Using the same decomposition to remove LSC contribution, we reveal that the remaining statistics recover usual HIT behaviours, being careful to use a local estimation of the kinetic energy dissipation rate. Finally, we propose a revisited and more precise estimation of the Richardson constant. 
 
In addition to these points, the pair separation is a good statistical tool to study the transport properties of turbulent thermal flows, especially in the presence of an inhomogeneous mean flow. 
Because of the large field of view in our experiment, we are able to study the effects of inhomogeneity due to Large Scale Circulation specific to the Rayleigh-Bénard convection. 
In future studies we aim to explore an even wider measurement volume and be able to study the jets where a high concentration of thermal plumes is observed. 
These coherent structures should affect the different components of the pair dispersion <cit.>. 
Furthermore, measurements of particle dispersion in the jets could give information on thermal transfer in the Rayleigh-Bénard convection. 
Possible analysis would look at the next order of the Taylor expression of Eq. <ref>. 
This can be done by computing the difference between the forward-in-time and the backward-in-time pair dispersions <cit.>, and would open new perspectives to study the energy cascade in turbulent convection.



  We gratefully thank Marc Moulin for manufacturing of the cell and Miguel López-Caballero for his help with the camera calibration. 
  The data processing and analysis were made possible with the PSMN computer resources. 
  We warmly acknowledge Theresa Oehmke for the manuscript proofreading.

 




