Walid Tarraf,Diogo Queiros-Condé,Patrick Ribeiro,Rafik Absi
Fractal Geometric Model for Statistical Intermittency Phenomenon
scale entropy,entropic skin theory,turbulence,intermittency,nonlinear dynamics,statistical physics,complex systems,multiscale fractal geometry,clustering
The phenomenon of intermittency has remained a theoretical concept without any attempts to approach it geometrically with the use of a simple visualization. In this paper, a particular geometric model of point clustering approaching the Cantor shape in 2D, with a symmetry scale θ being an intermittency parameter, is proposed. To verify its ability to describe intermittency, to this model, we applied the entropic skin theory concept. This allowed us to obtain a conceptual validation. We observed that the intermittency phenomenon in our model was adequately described with the multiscale dynamics proposed by the entropic skin theory, coupling the fluctuation levels that extended between two extremes: the bulk and the crest. We calculated the reversibility efficiency γ with two different methods: statistical and geometrical analyses. Both efficiency values, and , showed equality with a low relative error margin, which actually validated our suggested fractal model for intermittency. In addition, we applied the extended self-similarity (E.S.S.) to the model. This highlighted the intermittency phenomenon as a deviation from the homogeneity assumed by Kolmogorov in turbulence.
Fractal Geometric Model for Statistical Intermittency Phenomenon

The phenomenon of intermittency has remained a theoretical concept without any attempts to approach it geometrically with the use of a simple visualization. In this paper, a particular geometric model of point clustering approaching the Cantor shape in 2D, with a symmetry scale θ being an intermittency parameter, is proposed. To verify its ability to describe intermittency, to this model, we applied the entropic skin theory concept. This allowed us to obtain a conceptual validation. We observed that the intermittency phenomenon in our model was adequately described with the multiscale dynamics proposed by the entropic skin theory, coupling the fluctuation levels that extended between two extremes: the bulk and the crest. We calculated the reversibility efficiency γ with two different methods: statistical and geometrical analyses. Both efficiency values, and , showed equality with a low relative error margin, which actually validated our suggested fractal model for intermittency. In addition, we applied the extended self-similarity (E.S.S.) to the model. This highlighted the intermittency phenomenon as a deviation from the homogeneity assumed by Kolmogorov in turbulence.

Kolmogorov’s theory in developed turbulence (K41) assumes that turbulence is statistically homogeneous (invariant due to spatial translation, far from limits), isotropic (invariant due to rotation), self-similar (invariant due to scale expansion) and stationary (invariant due to time translation) [1,2,3]. According to Kolmogorov’s theory, the developed turbulence can be described using the structure functions of dissipated energy and velocity , with ε(r) the rate of energy dissipation per unit mass averaging over an r-dimensional ball and δV(r) the velocity fluctuation V between two points separated by a distance of r. This theory, based on the hypothesis of the homogeneity of energy dissipation, predicts universal laws of scales and for high Reynolds numbers and scales in the inertial range (η , where η and are, respectively, the Kolmogorov scale and the integral scale), such as and being the scaling exponents. Considering the assumption of the refined similarity , this means that [4,5,6,7]. Indeed, Kolmogorov predicted that followed a linear law, as the energy dissipation rate was assumed to be scale-independent according to the homogeneity’s hypothesis (), hence, . However, the Kolmogorov’s theory does not agree with the results of many numerical and experimental studies carried out later, showing significant deviations from the linearity of scaling exponents assumed by Kolmogorov [1,8,9,10,11,12,13,14]. This deviation from linearity (homogeneity’s assumption) is due to the intermittent characteristics of energy dissipation, where the gradient of turbulent flows fluctuates significantly over the volume occupied by the fluid [15,16,17]. This phenomenon is known as “intermittency”. Intermittency in turbulence is a phenomenon of nonlinear dynamics resulting from complex and nonlinear interactions between different variables of a system [18,19,20]. This phenomenon has been studied with the use of various geometrical theories, such as fractal geometry. One of the most important was the theory of Procaccia and Grassberger on the strange attractor method [21]. Therefore, several models have been developed to describe it and determine the exponents of intermittency, such as the She and Lévêque’s model [22,23], the multifractal model [24,25], the β-model [14,26] and the log-normal model of Kolmogorov, Obukhov and others [13,24,27]. Turbulent flow has been presented as a hierarchical structure of fluctuation levels distinguished using some thresholding processes. Thus, these fluctuation levels are characterized through structure functions, such as the difference between the velocity or energy dissipation rate with average values, respectively. Therefore, the hierarchy of intermittent turbulence extends between disordered structures corresponding to low fluctuations and ordered structures, where the highest fluctuations are generated. This theoretical concept has been validated through experimental observations, such as for high Reynolds numbers, as well as for large thresholds, with the strongest fluctuations taking place on a filamentary structure of a fractal dimension equal to 1 [10,28,29], while, for a low threshold, the fluctuations take place on a structure of fractal dimension equal to 2.36. Thus, this variation in fractal dimensions with fluctuation levels emphasizes the hierarchical multiscale characteristics of intermittent turbulence [24,25]. Note that these multiscale characteristics have been discussed in several studies that proposed that turbulent interfaces, such as those in clouds, turbulent jets, boundary layers and turbulent combustion, exhibit fractal behavior for scales belonging to the range of inertial scales [30,31,32].

Clustering is a data analysis technique widely used in the field of dynamical systems, especially turbulence, to identify spatial and temporal structures in turbulent flows [33,34,35,36]. Clustering can be used to divide a flow into fields with similar characteristics in terms of turbulent dynamics. The identified fields can correspond to vortices, laminar flows and recirculation zones [37,38,39,40]. This technique tends to be used in conjunction with a fractal analysis to study the phenomenon of developed turbulence. The fractal analysis is a geometric analysis method that quantifies the irregularity and complexity of structures in dynamic systems and their variability at different spatial scales [41,42,43]. Since turbulence is a complex and chaotic phenomenon, characterized by rapid and random fluctuations in flow variables such as velocity and pressure, the fractal analysis has been widely used to quantify the geometric complexity of turbulent structures and their variability at different spatial scales [44,45]. The complexity can be measured well with the fractal dimension, having been measured using several methods. One of the well-known methods is box counting. Based on the concept of Hausdorff dimensions, this method consists of covering the object or structure with boxes of different sizes and counting the number of boxes needed to cover the object. In this way, the fractal dimension is equal to the ratio of the logarithm of the number of boxes to the logarithm of the size of the boxes [46]. This has been used to identify the most important spatial scales for turbulent energy dissipation, such as the Kolmogorov scale [47]. The combined use of these two techniques has led to a better understanding of turbulence mechanisms.

One of these models, the entropic skin geometry (E.S.G.), shown to be in agreement with other experiments [17,48], is an interpretation of this intermittency phenomenon achieved through a coupling between the multiscale fractal geometry and statistical intermittent analysis within the turbulence by assuming the existence of a certain hierarchy of correlated fractal structures, characterizing the dynamics of fluctuation levels (skins). Each fluctuation level of the order p has been associated with a fluctuating field , representing a structure of dimension , included in the field of dimension , such as All fluctuating fields are included in the initial field . Therefore, according to this theory, turbulence is considered a superposition of structures ranging from the extremely disordered one, showing less fluctuations (at ) with dimension equal to the embedding dimension (), to the extremely ordered one (filamentary form), where the largest fluctuations (crest) in regard to the mean occur (when ). For , there is dimension , designating the classical fractal dimension , and corresponding to the structure called the “bulk emergence”. Likewise, for , there is the extremely ordered structure that we called the “ crest convergence “ with crest dimension [17,49]. Therefore, the bulk and the crest represent the two extremities of the hierarchy of skin structures. Furthermore, at the geometric level, the concept of scale entropy for each fluctuation level p in a scale r has been introduced to quantify the spatial extension of each skin in the intermittency hierarchy. It is defined as the logarithm of the space occupation rate [36,50]: where and is the volume occupied by at scale r in a total volume corresponding to the integral scale . Since is the minimum number of coverages with thickness covering the set , the volume of the set can be determined using , with showing that the set is fractal. In addition, the entropic skin theory proposes the notion of geometric reversibility efficiency , defined as the scale entropy flux ratio, characterizing its transmission through two consecutive skins, and : This assumes that the reversibility efficiencies are constant throughout the whole skin hierarchy. Thus, the relationship in the reversibility efficiency across all structures of the overall hierarchy obtained with the recurrence method can be given by: such as . For , this means that the loss of entropy through the skins is zero and, therefore, the phenomenon is nonintermittent (homogeneous). For the phenomenon becomes infinitely intermittent. In the statistical approach, the value of a structure function is given by its “active part”, offering two contributions; the first one is the average energy dissipation rate over the actual active part of the field, and the second is the volume fraction, which indicates an intermittency factor and allows for the “lacunary aspect” of energy dissipation [3,20]. Thus, we would have: for , . It should be noted that follows a power law [49]: Inspired by the idea of coupling between the fractal geometry and the statistics, these quantities allow to determine experimentally, proposing the relative structure functions: such as and , through which we obtain the relationship of the relative structure function [22]: Since = , for we could obtain . The same applies to the velocity fluctuation, introducing relative structure functions: Similarly, for the velocity, we could still write based on the relationship and the refined similarity hypothesis [3]: with and (r) = . Note that this equation was also obtained by She and Lévêque in their model [22]. Thus, could be obtained statistically with Equations (3) and (4). The reversibility factor calculated with the geometric analysis (using the entropic scale) would be the same when calculated with the statistical analysis (using the structure functions). For this purpose, in order to distinguish them, we noted them, respectively, as and , whose equality can be expressed through . This equality has an intrinsic and fundamental consequence in the theory of E.S.G, allowing to link the statistical and geometric aspects of a spatial distribution [51,52]. Note that the expression of the exponents allows, via the Legendre transformation, to obtain the multifractal spectrum. Therefore, a link between the entropic skin intermittency theory and the concept of multifractals was established [49,53].

According to the extended self-similarity principle proposed in developed turbulence (for high and low Reynolds numbers), any function of order p (either of velocity or of the rate of energy dissipation) is a power law of any other function of order q, whose logarithmic plot slope is [54]. In the case of homogeneous and isotropic turbulent flows, according to the Kolmogorov assumption, the refined similarity relation () implies that . However, in an intermittent case, a deviation from p/q can arise due to fluctuations produced through the energy dissipation rate. In this paper, we propose a multiscale model to describe and represent the phenomenon of statistical intermittency. This phenomenon has never been defined in a geometrical pattern, and has remained a theoretical concept without any attempts to approach it through the use of a simple visualization. Indeed, the hierarchical structure of fluctuation levels has been assumed with several models and theories (such as the She and Leveque model [22]), but this remained within the limits of a hypothesis. Note that in the search for a geometrical explanation of intermittency, Queiros-Condé (1997) showed that diffusion-limited aggregates (DLAs) displayed analogous intermittent features to fully developed turbulence [55]. Therefore, we introduce a cluster model that generates intermittency across scales in a fractal manner, modeling this phenomenon in the visualization step in order to better understand it. We take up tools and statistics already used to describe intermittency in developed turbulence, mentioned above, such as the entropic skin theory and the extended self-similarity hypothesis, to verify the ability of this model to describe statistical intermittency [55]. The structure of this paper is presented as follows: After Section 1, Section 2 details the proposed cluster model to describe intermittency. Section 3 introduces the measures applied to perform the required geometric and statistical analyses, as well as the parameters that apparently influence them. In Section 4, we verify the validity of the model to describe intermittency by applying the entropic skin theory. The Section 5 and Section 6 include multiscale and statistical analyses, respectively. In Section 7, we discuss the results of the data analyses and validate the ability of the proposed geometric model to describe intermittency. Finally, Section 8 synthesizes the paper and discusses the next steps in the research. In the following, we present our model.

The objective of our work was to find a multiscale geometric model for the intermittency phenomenon. For this purpose, we proposed a simple model, comprising a clustering of points taking the form of a two-dimensional Cantor with a symmetry scale θ that represented the intensity of the clustering intermittency. The model started by distributing points randomly on a plane in a homogeneous way using MATLAB software. Then, we introduced a differentiation for the density by defining zones where the density of the points would be higher. In these areas, the points would be distributed randomly, but with a higher density than in the first step, following the 2D Cantor geometry. The differentiation process was then reproduced across scales. We, thus, associated the homogeneous distribution of points with a fractal differentiation mechanism that reproduced many features of the intermittent phenomenon. This mechanism allowed for creating a clustering structure showing high fluctuations in the point number and widely differentiated local densities. The aim was to simulate the bulk–crest dynamic supposed to describe the hierarchy of intermittency in developed turbulence. This model was iteratively developed following the steps detailed below: Step 1: Consider a domain of size pixels, throwing in a few points uniformly at random: iteration 0 (Figure 1a). Step 2: Remove the areas at the four corners from the initial draw (known in normal 2D Cantor), having sides of size ; θ is a symmetry scale and is the entire scale. Then, throw the same number of points as in the previous iteration into each of these areas (also uniformly random) (Figure 1b). Step 3: Consider each generated “sub-block” square and apply this process. Step 4: Iterate to infinity. Note that scale represents the size of the sides of the generated squares at the ith iteration, such as . The model is shown in Figure 1c. Generally, we worked in a domain with 4096 × 4096-sized pixels, except when we wanted to know the effect of the size variation on the measurements. The scaling factor () showed a fundamental factor reflecting the degree of intermittency varying in the range [0, 0.5], such that if was relatively low, the distribution would display strong fluctuations, and, thus, considerable intermittency. If was high (close to 0.5), the distribution would become close to homogeneous.

The interest in the 2D Cantor clustering model is that it is very simple and efficient to describe, and allows for a clear visualization of the geometry of an intermittent structure regarding the bulk and crest dynamic. To perform the required analysis, we covered the square box with a grid of the same size, having an elementary mesh length r representing the scale. We counted the different parameters at each scale r (see Figure 2). For the scale analysis, we measured the number of meshes affected by the distributed points for each scale r in order to represent the shape of in a logarithmic scale as a function of scale , according to the box-counting method mentioned in Section 1.3. The fractal dimension presented the opposite of the slope of the curve. At the statistical analysis level, we counted the number of points in each cell in order to measure the local fluctuations in point numbers. We determined the mean number of points at the given scale as follows: where is the number of blocks at scale and the number of points in mesh . We, therefore, defined the structure functions (mean fluctuation rate of point numbers) at fluctuation level p: According to the principle of cascading fluctuating structures through the scales of entropic skins, we assumed that the structure functions followed a scaling law: Due to of the random jets (uniform distribution), some deviations in the number of affected meshes could exist. This produced uncertainty problems. For this reason, we remedied this problem by repeating the jets of points several times in order to reduce these uncertainties and by averaging the measurements on all the jets. Several parameters seemed to influence the geometric and statistical analyses, including the symmetry scale and the topic scale , as well as the population factor . We defined as the topic scale, the mean distance between particles, at a certain iteration : where represents the number of points per block associated with iteration , is the size of the block, such that with is the symmetry scale (intermittency parameter), the population factor represents the ratio of the number of points per block and and , are associated, respectively, to two certain consecutive iterations I and i+1, such that (see Figure 3).

Similar to the dynamics of entropic skins applied in developed turbulence [51,52], we assumed that the intermittent statistic in the fractal clustering model was driven by two specific fractal objects: the bulk and the crest. As such, these two objects were associated, respectively, with the fluctuating structures , with dimension = (at ), such as as the fractal dimension of the phenomena, and , with dimension (). Therefore, using analogy with the assumption that giving the mean fluctuation rate of energy dissipation would be associated with the p-level of fluctuation through (in Equation (1)), we proposed the mean fluctuation rate of point numbers in scale r, associated with the p-level, given by: where is the mean fluctuation rate of point numbers on the affected blocks, corresponding to the fluctuation level p. Likewise: Dividing Equation (9) by Equation (8) gave: Let us define the relative structure function as follows: Obviously: Therefore: This gave us: Therefore, we could obtain: When generalizing, we obtained: where . In our study, we took as a thresholding path.

For the scaling analysis, we used the box-counting method, counting the grid cells affected by the studied structure. This analysis allowed us to measure the fractal dimension () of the 2D Cantor clustering of the points, for which we calculated the reversibility efficiency according to the geometrical method: In order to perform this calculation, first, the crest dimension had to also be calculated. We could then obtain the geometric efficiency of reversibility .

First, we computed the crest dimension for the simplest cases that corresponded to a population factor . This meant that the number of points per block was constant across the iterations. Since fluctuation and intermittency develop on the spatial support or “effective zone” between two extremities of the hierarchy of structures, we could deduce that in our 2D Cantor clustering, the crest presented with the densest and most clustered areas compared to the remainder. Given that the clustering of points occurred across iterations in the form of blocks, exhibiting scale invariance, the most densely clustered zones in points corresponded to the last iteration. Moreover, since the areas among the dense blocks were neglected across the fluctuation levels due to the adopted thresholding in the entropic skin model, we could compute the crest dimension (case a = 1) with the classical box-counting method at by counting the number of blocks at the last iteration i (the highest), as follows: where is the number of blocks at the last iteration i. On the other hand, the phenomenon of the clustering of points in the blocks highlighted the crucial role of the topical interparticle scale. The scaling analysis fully reached the points generated until the last iteration i (the highest), such that the number of affected meshes was equal to the number of these points at small variable scales r relatively close to the topical scale . From this perspective, we noticed that the scale analysis remained homogeneous up to the topic scale of the homogeneous distribution (iteration 0) ; thus, the curve began to enter a transitory phase before taking on an asymptotic appearance, such that its slope at smaller scales tended towards (see Figure 4). To this end, when generalizing , we undertook the analytical study of the shape of the scaling analysis curve through the use of iterations among the topical scales and ∀ i, particularly at small scales. We, therefore, introduced the local fractal dimension , corresponding to iteration i, among the topical scales and . Indeed, our adopted geometric attractor generated self-similar distributed blocks containing homogeneously distributed points. There were obviously crossovers among these blocks through the iterations, in which there was also a homogeneous distribution of points. During the counting of points, we had to distinguish between the two zones that were called the block zone and cross zone, as shown in Figure 5. The two areas formed the whole of the generated geometry. Knowing that at the topic scale the number of affected meshes was theoretically equal to the number of points accumulated until this iteration i. Therefore, we defined the equation of the total number of points thrown into the 2D Cantor clustering up to iteration : such that and represent, respectively, the number of points having accumulated at the cross zone and at the block zone at iteration i. Moreover, is the number of points that accumulated at the cross zone at the first iteration and is the total number of points at iteration 0, where represents the populating factor (see Figure 3). We understood that: For a very important number of iterations, we obtained (for more details see Appendix A): Using , and Equation (7) we obtained: Therefore, we found: For and , which was the case of a normal one-third 2D Cantor: This verified the shape of the curve, such that its slope at small scales tended towards the crest dimension , as seen in Figure 4. We verified this result with the use of both analytical and graphical methods for different values of ϵ [1, 1.4] and ϵ [0.25; 0.43]. This motivated us to search for the topical scale from which the fractal analysis curve adopted an asymptotic behavior of slope in order to define the range of scales allowing for an approximate measurement of the crest’s dimensions. Indeed, since the local fractal dimension completely left the value from up through the scales as shown in Figure 4, this highlighted the role of this scale. We, therefore, defined the crest scale range between the global scale and the topical scale at the first iteration . For this purpose, we called the crest scale .

Let us go back to the definition of the bulk, since it was the lowest extremity in the hierarchy of skins showing the least fluctuation, which included all the highest fluctuation levels. Therefore, the bulk extension presented the spatial support, where the whole studied distribution occurred. Furthermore, we noticed in our 2D Cantor clustering that the density of points increased through the iterations, indicating that the least-fluctuating structure corresponded to the first iteration. Therefore, we could say that from , as we moved down the scale, the analysis curve began to leave the bulk domain towards the crest domain (by descending with the scale) when the bulk was in the scale range between and . To test our hypothesis, we measured the local fractal dimension at variable steps of scales by fixing the smallest scale end on the topical scale of the last iteration and varying the other end of the scale r or by performing the measurement for each ∆r. We noted that when ∆r was close to the value ; i.e., for ; (in the logarithmic representation), the local dimension began to touch values very close to as it descended with the scale range ∆r and, hence, the scale r (see Figure 6). This verified our observation. Thus, the topical scale at the first iteration exhibited a characteristic scale to determine the two scale ranges of the crest and body. Therefore, to calculate the fractal dimension or () of the distribution, we measured the linear regression slope of the scaling curve over the range between what we called the crest scale and the integral scale as the distribution extended over the entire space (see Figure 7). To sum up, the 2D Cantor clustering scale analysis was explicitly dependent on topical scales through successive iterations. It was, therefore, divided into three phases. First, it began with the homogeneous phase extending between and . In this phase, the analysis had a dimension equal to the Euclidean dimension (d = 2). Then, it continued in a transition phase between and , where the local fractal dimension changed from the Euclidean dimension value to the one of the crest dimension . Finally, the scaling analysis entered the crest phase, occupying the scale range , for which the curve exhibited an oblique asymptote with a slope equal to the crest dimension (see Figure 7). The scale analysis in the homogenous distribution gave us the fractal dimension , equal to the Euclidean dimension of plan .

Since the structure functions exhibited a rate of fluctuation across scales and fluctuation levels following the scaling law , and in order to perform the statistical analysis of the studied geometry, we measured the scaling exponents. For this purpose, we calculated the structure functions according to Equation (5), and represented them in logarithmic coordinates showing the scaling exponent as the slope. In the homogeneous case, we obtained straight lines with a low dispersion of points, having slopes equal to , from which we found that varied as () (see Figure 8). In the case of the intermittent Cantor, a certain log-periodicity appeared in the moments (at the different scale ) as a function of the normalized scale among the successive values of where (see Figure 9). Note that was different from p, unlike the homogeneous case. Some deviation from homogeneity occurred [50,56]. Note that to define the inertial range, we represented as a function of , in the homogenous case; we observed that there was a scale range where was constant, which corresponded to the inertial range (see Appendix B).

In order to highlight the phenomenon of intermittence and its characteristics, we used the hypothesis of extended self-similarity previously seen in studies of intermittence in developed turbulence. This hypothesis assumed the existence of the extended self-similarity of structure functions through different levels of fluctuations [3,54]. Therefore, we assumed that any structure function of point numbers at fluctuation level p was a power law of any other structure function at fluctuation level q. Through representing as a function of and measuring a slope value , we obtained the slope equal to in the case of homogeneity. A deviation in had to occur in the case of intermittency. Indeed, we undertook the verification of this law for both the homogeneous and intermittent cases (see Figure 10).

To compute the statistical reversibility efficiency, we used Equation (10). We presented and in logarithmic coordinates for all scales on the same graph. Assuming that was constant for all the scales, we obtained a linear curve; its slope was (see Figure 11). To ensure that we worked in a range of scales where intermittency developed, we avoided scales smaller than as the last iteration topic scale; thus, the measure scale range of was []. In the homogeneous case, the statistical analysis gave a factor .

As we mentioned at the beginning, our objective was to propose a simple fractal geometric model to explain the statistical intermittency phenomenon. We applied the entropic skin theory concept and the extended self-similarity hypothesis to verify the ability of the proposed model to describe this phenomenon. Indeed, we accomplished the statistical and geometrical analyses by performing all the measurements required for different degrees of intermittency through varying the intermittency intensity chosen as the symmetry scale . It was the most effective factor that explicitly and proportionally presented the density of intermittency. The parameter varied within a range of [0; 0.5], such that if tended to 0, the 2D Cantor’s clustering became highly intermittent. In addition, if neared the value of 0.5, the clustering became less intermittent and affected the homogeneity in such a way that the fluctuation to the average strongly decreased. We measured the values of as a function of p by applying the self-similarity assumption. We found that the curve showed a linear shape very close to p, identical to the homogeneous case, such that signified the existence of self-similarity; the structure function was a power law of . This was analogous to the Kolmogorov theory in developed turbulence. A certain deviation from the power law, beginning from p = 1.2, was observed in the cases of 2D Cantor clustering for different intermittency parameters θ (symmetry scale). This highlighted the intermittent characteristics of the clustering as a structure linking the levels of fluctuations (see Figure 12). The deviation from homogeneity increased inversely, proportional to the intermittency parameter . Moreover, when increased the ratio , self-similarity approached the homogeneous case. This was compatible in the case of intermittency in fully developed turbulence.

The measurement of the bulk and crest dimensions for various intermittence parameters revealed the correspondence and harmonization between and the two dimensions, explicitly highlighting the density of intermittence. We represented as a function of for , according to Equation (11). Figure 13 shows the variation in as a function of , which seemed to be parabolic, so that could grow with . The same applied for , while was linear with the slope equal to . The interpolation of the two-dimensions curves and gave two equations: The two curves intersected at an abscissa ; hence, and 2. This in fact corresponded to homogeneity.

The obtained results showed equality between the statistical reversibility factor and the other geometric , with a maximum relative error margin of 1.7%. The two intermittency coefficients and , as functions of factor , showed almost the same trends, with negligible errors for both curves. They increased with up to a value of 1 for θ = 0.5 and then decreased, showing the increase in the intermittency intensity. Thus, this equality remained validated with a narrow margin of error. Indeed, the plot of and showed a clear correlation between these two quantities (see Figure 14). To verify this equality, we performed the experiment by varying all the parameters that could influence the distribution, such as the symmetry scale , the number of points thrown per block , the box size and the population factor a. However, we presented the ratio as a function of intermittency intensity θ. The values were scattered around the right-hand side: y = 1 (see Figure 15).

In this paper, we proposed a simple fractal geometric model to explain the statistical intermittency phenomenon. For this aim, we were inspired by the E.S.G entropic skin theory, characterizing the phenomenon of intermittency in fully developed turbulence using coupling fractal geometric and statistical intermittency properties. For this aim, a particular geometric model of a cluster of points approaching the Cantor shape in 2D, with the symmetry scale θ being an intermittency parameter, was proposed. This model involved the introduction of point density differentiation across scales in a multiscale mechanism that reproduced several features of the intermittency phenomenon. To verify its ability to describe intermittency, we applied this model to the concept of the entropic skin theory, obtaining a conceptual validation. This validation was confirmed through the intrinsic equality of the two reversibility factors and , calculated using two different analysis methods: geometric and statistical. This equality was accompanied well by a correlation as a function of the symmetric scale θ for different parameter values that could influence the geometry. We found that the intermittency phenomenon in the 2D Cantor cluster was adequately described with the multiscale dynamics, proposed by the entropic skin theory and coupling the fluctuation levels that extended between two extremes: the bulk and the crest. The dimensions and , as well as γ, explicitly expressed the intensity of cluster intermittency. The crest displayed the densest and most aggregated areas with respect to the remainder of clusters, while the bulk displayed the spatial extension where all the studied clusters took place. Then, at the statistical analysis level, we applied the extended self-similarity hypothesis (E.S.S.) to the model; a deviation from the homogeneity assumed by Kolmogorov in the fully developed turbulence was noticed. This highlighted the intermittency dynamics as a multiscale structure linking the levels of fluctuations. One of the potential applications of our work is related to wall turbulence, which needs better understanding of the involved phenomena. Our approach is related to experimental near-wall turbulence. The study of the phenomenon of intermittency in wall turbulence based on ESG and the notion of scale entropy used to validate the clustering fractal model was applied to an experimental database of boundary layer flows [3]. This geometrical approach was validated with experimental data [57] considering a set of 50 cases corresponding to 50 values of y+ ranging from 1 to 2722. These boundary layer experiments showed the interest in our approach for a wide range of wall distances (1 ≤ y+ ≤ 2722). Using extended self-similarity, the scaling exponents were found to vary, with y+ displaying a tendency toward saturation in the vicinity of the wall, a behavior captured with the ESG. The theoretical and experimental results showed that the entropic skin geometry based on conceptual yet simple geometrical arguments could describe intermittency features of wall turbulence. The aim of our study is to find potential applications in computational fluid dynamics (CFDs). For RANS and hybrid RANS/LES modeling of wall-bound turbulent flows, the demand for a fine grid resolution is great near walls, particularly for high Reynolds number flows. Hybrid RANS/LES, used in the outer region large Eddy simulation (LES) approach and in the near-wall region low-Re RANS model, remains computationally intensive and time consuming due to the high demand for wall-normal grid resolutions. It is possible to decrease the demand for fine grid resolutions at the wall region by using wall functions. In some cases, these methods have shown that they allowed to obviate the need for solving the turbulent kinetic energy (TKE) equation near walls, resulting in quicker and more accurate engineering computations [58,59,60]. Our proposed fractal model could present an interesting framework for the development of new tools for modeling near wall turbulence, especially since it highlighted the fractal aspect of statistical intermittency. This study also provides a basis for various potential future applications for complex systems and for ecological and social phenomenon.