^{1}

^{1}

^{1}

^{2}

Analysis of the entrance and wall dynamics of a high-flux gas-solid riser was conducted using embedded solid concentration time series collected from a 76 mm internal diameter and 10 m high riser of a circulating fluidized bed (CFB) system. The riser was operated at 4.0 to 10.0 m/s air velocity and 50 to 550 kg/m
^{2}s solids flux of spent fluid catalytic cracking (FCC) catalyst particles with 67 μm mean diameter and density of 1500 kg/m
^{3}. Data were analyzed using prepared FORTRAN 2008 code to get correlation integral followed by determination of correlation dimensions with respect to the hyperspherical radius and their profiles, plots of which were studied. It was found that correlation dimension profiles at the centre have single peak with higher values than the wall region profiles. Towards the wall, these profiles have double or multiple peaks showing bifractal or multifractal flow behaviors. As the velocity increases the wall region profiles become random and irregular. Further it was found that, as the height increases the correlation dimension profiles shift towards higher hyperspherical radius at the centre and towards lower hyperspherical radius in the wall region at
* r/R *= 0.81. The established method of mapping correlation dimension profiles in this study forms a suitable tool for analysis of high-flux riser dynamics compared to other analyses approaches. However, further analysis is recommended to other gas-solid CFB riser of different dimensions operated at high-flux conditions using the established method.

The gas-solid circulating fluidized bed (CFB) systems have been shown to be complex and largely governed by the non-linear relationships. The solid concentration or voidage and pressure time series from these fluidized beds systems exhibit the characteristics of low dimensional deterministic chaos [

Correlation dimension is one of the chaos analysis parameters which is widely used to account for the complexity of the gas-solid flow systems in riser reactors [

This study investigated the flow dynamics behaviors of the gas-solids suspensions flow in the riser reactor using the measured solid concentration time series signals. In particular, this paper presents chaos analysis of the temporal variations of solid concentration in the entrance section, flow development section, developed flow sections and in the wall region of a high flux gas-solid riser using correlation dimension parameter.

This paper makes the use of the fact that the correlation dimension varies along the correlation function as the hyperspherical radius (r) is changed across the attractor. While only the maximum value is normally reported by different researchers [

Chaos analysis applies to non-linear and deterministic systems whose behavior can in principle be predicted for a while and then appears to become random. Studies show that, dilute gas-solid suspension flow is predictable, until when clusters form and pass across a specific point, events which appear at short intervals as depicted by cluster frequency [

Literatures show that non-linear time series analyses techniques are founded on various theories of dynamical system of which its time of evolution is expressed in some phase space [

Thus, being one of the non-linear techniques, chaos analysis is achieved by first reconstructing the attractor in the phase space or state space [

The state of the gas-solid dynamical system in a CFB riser can be defined by a time series of one of the system’s operating conditions or parameters such as the voidage or solid concentrations and pressure fluctuations. This is referred to as time series of single variable. It is from such a time series of a single measured parameter where reconstruction of a phase space or state-space is done through embedding method which utilizes time delay and embedding dimension [

The reconstructed phase space by using a single variable signal to represent other true variables, through time lag which is applied to the single variable signal, mimics the true phase space which is created from the true process variables. The CFB riser, for instance, is operated under various true variables such as gas velocity and solids flux; the measured solid concentrations, ε_{s}(t), can be used to reconstruct the true attractor by using the time delay. The number of variables that affect the system can be estimated from the reconstructed phase space and the determined system invariants such as correlation dimensions. However, in most studies once the correlation dimension is established, and the variation of such a dimension with operating conditions and spatial location is known, the outcome is taken to be a sufficient description of the complex system.

Correlation dimension is a measure of the fractal dimension of the time series, which express the complexity that quantifies or gives the characteristics of the geometry and shape of the attractor. It measures the scaling of the points along trajectories on the attractor, and tells to which exponential law these points follow. The higher the correlation dimension, the more the complex microscopic flow dynamics and vice versa in relation to the time scale. Microscopic level is depicted by smaller volume of measurement point used, especially when optic fibers are used, as compared to pressure drops across a riser section which can be referred to as macro scale. The values of the correlation dimension are commonly obtained by first plotting the graph of log(Cr) versus log(r) which shows the linear dependence in some ranges of r values. Correlation dimension estimator of the attractor is the slope of the scaling region of the correlation integral, C(r), versus hypersphere radius, r, [

The process of determining correlation dimension involves first computing the correlation sum which counts the vector pairs (X_{i}, X_{j}) whose distance is smaller than the selected hyperspherical radius, r, as given in Equation (1). The counting is started by letting the count i constant, and then traversing all counts j to the total number of data points, N-1, in the time series before i is also changed up to N in this case. When the total data points approaches to infinity, i.e., N → ∞ , then the value of N(N-1) approaches to N^{2}. The correlation sum, C(r), is determined as shown in Equation (1) [

C m ( r ) = lim N → ∞ 1 N 2 ∑ i = 1 N ∑ j = i + 1 N Θ ( r − [ ∑ k = 1 m ( ε s ( i , k ) − ε s ( j , k ) ) 2 ] 1 2 ) (1)

where ε s ( i , k ) and ε s ( j , k ) represents the vector points of the solid concentration data signals and Θ ( u ) is the Heaviside step function, for which:

u = r − [ ∑ k = 1 m ( ε s ( i , k ) − ε s ( j , k ) ) 2 ] 1 2 and Θ ( u ) = { 1 , i f u ≥ 0 0 , i f u < 0

Thus, the correlation sum, C(r), is a cumulative distribution function (CDF) which reflects the distribution probability of the vector distances which are less than r between pairs of points on the chaotic attractor [_{min} < r < r_{max}, then C(r) scale like a power law. Further, if the attractor characterises the time series of a single dynamical variable, in the positive ranges of the hyperspherical radius, r, the correlation sum, C(r), is related to the hyperspherical radius, r, by the power law as per Equation (2) [

C ( r ) ∝ r D (2)

where D is the characteristic parameter of an attractor referred as correlation dimension which is the. Thus, when there is an attractor in the phase space, the relationships can be presented as C(r) ~ r^{D}^{'}, where D' is the topological dimension of the attractor. The topological dimension can be regarded as local correlation dimension determined at every short distances, r. The dimension, D, is also equal to the topological dimension D' of the attractor if the system dynamics is classified as periodic or quasi-periodic. Consequently, the fractal dimension D is not an integer if the system is chaotic and the attractor is a strange attractor. Thus according to the power law property, basing on the behavior of a correlation sum, a dimension value D can be defined as [

D ( r ) = lim r → 0 lim N → ∞ d [ log C ( r ) ] d [ log ( r ) ] (3)

Since the number of points is fixed during computations of correlation sum, C(r), then Equation (4) simplifies to Equation (5):

D ( r ) = lim r → 0 d [ log C ( r ) ] d [ log ( r ) ] (4)

Using the scaling region, D is the slope of ln(Cr) versus ln(r) plot. The slope D of the linear part of the log-log curve provides all necessary information for characterizing the attractor [

Chaos analysis begins by reconstruction of the phase space through embedding techniques [

Data were collected from a CFB system shown in ^{2}s solids flux and 4.0 to 10.0 m/s gas velocity. In particular U_{g} = 5.5 m/s, 8.0 m/s and 10 m/s were used in this study. Fluid catalytic cracking catalyst particles with 67 μm mean diameter and density of 1500 kg/m^{3} were used. A 70% to 80% humid air was used for transporting the solid particles. Signals were sampled from eight (8) axial levels (i.e., Z = 0.98, 1.52, 2.73, 3.96, 5.13, 6.34, 8.74, and 9.42 m) and 11 radial points (i.e., r/R = 0.00, 0.16, 0.38, 0.50, 0.59, 0.67, 0.74, 0.81, 0.87, 0.92, and 0.98) at each level where r/R is the normalized radial distances from the center to the wall of the riser. To each point, 29,100 data points of solid concentration were sampled in 30 seconds using optical fiber probe at 970 Hz.

The embedding dimension was determined by the method of saturation of correlation dimension. The maximum correlation dimension was determined at

different embedding dimensions and the profile drawn. The minimum embedding dimension was determined at the point where maximum correlation dimension saturate. The time delay was determined from the maximum correlation dimension versus time delay profile. The time delay (τ) can be any value as it is not going to change the attractor of the system. But it decides the independency of the coordinates. The time delay that keeps the coordinates more time independent is desirable. Too large values cause loss of information contained in the data and too temporally close vectors become rather remote, giving rise to uncertainties i.e. large delays lead to vectors whose components are (almost) uncorrelated and the data are thus randomly distributed in the embedding space. On the other hand, small time delay yields correlated vector (crossing among the trajectories) elements. _{3}, against the embedding dimension, m, and the time delay, τ, at Z = 1.52 m for U_{g} = 5.5 m/s and G_{s} = 300 kg/m^{2}s. From

The delay time, τ, was estimated from point where the curve gives the higher correlation dimension after the first minimum value. From

Determination of correlation dimensions (D’s) were achieved by first computing correlation integrals, C(r), at various hyperspherical radius, r, by applying Equation (1) incodes prepared in FORTRAN 2008. Correlation dimensions were then determined as the slope of the correlation integral, ln(Cr), versus hypersphere radius, ln(r), which were estimated between point intervals of three to nine steps. Slope was determined as the ratio of change in ln(Cr) to the change in ln(r), to get D_{3}, D_{4}, …, D_{9} using Equation (4). Determined correlation dimensions were plotted to get D’s versus ln(r) plots which were assessed. The D_{3} versus ln(r) plots gave the maximum value at each signal and hence these plots were used for further analysis, Maximum values were also determined, plotted and analysed. This method is established in this study to analyse the high-flux gas-solid CFB riser.

Given the correlation function, C(r), in discrete form, values of which were established at several values of hyperspherical distance, r, the first derivative of the ln(Cr) function with respect to ln(r) yields the correlation dimension function of ln(r). Thus, the D_{3} values were obtained as per equation:

D 3 ( ln ( r ) ) = d [ ( ln C ( r ) ) ] d [ ( ln ( r ) ) ] | n n + 2 (5)

This is the slope determined at intervals of consecutive three (3) points. The values of D_{3} corresponding to the n^{th} value of ln(r), a plot of which is defined in this study as D_{3} profile. Profile of D from different intervals, 4 to 9, along ln(r) did not yield the desired maximum values, hence the study focused on D_{3} values or D_{3} profiles. The use of these profiles to study the flow dynamics in the high-flux gas-solid riser is also established in this study.

_{g} = 8 m/s and G_{s} = 300 kg/m^{2}s. Correlation dimension can be determined from the

slope of ln(Cr) and ln(r) plot as per Equation (5). The question of how many steps to use in the determination of slope is answered in _{3} to D_{9} (i.e., slopes obtained from three to nine steps) at the upper part and below there is a correlation integral profile from which the values of D_{3} to D_{9} were estimated. From D_{3} to D_{9}, i.e., increasing the number of steps for determination of slope implies determination of the optimum distance, r, which brings the highest value of topological dimension, referred to as correlation dimension. Higher number of steps implies larger local r, r → ∞. Comparing the values of D_{3} to D_{9}, a three (3) steps correlation dimension (i.e., D_{3}) profile was found to give a maximum correlation dimension and thus used for further analysis of the correlation dimension. Since the single step of r, the smallest unit of r is known setting D_{3} implies 3Δr, which is well known. Thus, this study used D_{3} values to estimate the correlation dimension for the given signals. It is noted also that across the whole correlation sum function, denoted as ln(Cr) in _{3} ranging from zero to D_{3}, as ln(r) changes. Plots of D_{3} versus ln(r) form different profiles of which characteristics are extensively studied in this work.

The correlation dimension profiles, i.e., D_{3} at the centre region (r/R = 0.0) have only one peak value on the D_{3} profile in all axial levels as shown in the bottom (Z = 1.52 m), middle (Z = 5.9 m) and top (Z = 9.42 m) sections. However, correlation dimension profiles for signals near the wall (r/R = 0.98) in all axial levels have multiple peaks for D_{3} for the entrance (Z = 1.52 m), flow development (Z = 5.9 m) and fully developed flow (Z = 9.42 m) sections.

According to Equation (4), the correlation dimension value is defined when r → 0 and N → ∞. That is, at small values of r obtained from correlation sum using large number of data points in the original time series.

The single peak values of the correlation dimension profiles at the centre in all axial levels indicate the uniform flow structure or low tendency to the formation of segregated flow or solid distribution fluctuations compared to the wall region, where the flow is characterized by the formation of dense phase such as particle aggregates or clusters in the dilute phase. The multiple peak values of correlation dimensions with varying magnitudes towards the wall as shown at r/R = 0.98 in the entrance section (Z = 1.52 m), flow development (Z = 5.9 m) and top or fully developed (Z = 9.42 m) sections indicates the formation of multifractal flow structures which is attributed to the segregated flows, i.e., particle aggregates with varying solid concentrations. From

Further observation of _{3} profile. A maximum correlation dimension is a good indicator that can reveal the dynamical differences between various points or states of the system. Thus, since D_{3} profile gives a maximum correlation dimension it is therefore, a good and strong indicator of the flow complexity, variations of which can provide details of the micro-flow structure of the gas-solid flow suspensions compared to the rest of D values, that is, D_{4} to D_{9}. Moreover, a three-step slope show good linear correlation compared to D_{4} to D_{9} due to small range/width of the scaling region observed when using solid concentration signals. It should also be noted that chaos analysis is about effect of small changes in the gas-solid flow on the dynamics of the system such as noted in instantaneous solid concentrations, ε s ( t ) , thus D_{3} being dependent on small changes in distances across the attractor, forms a good estimator of the dynamics of the riser. In this case, the correlation dimension, D_{3}, was selected to be used as a measure of the gas-solid flow behaviors in this study.

Based on discussion above, and also from literature, correlation dimension is a strong tool which can be used to indicate the complexity of the gas-solid flow behavior in the CFB riser reactors. A higher correlation dimension is related to a more complex microscopic flow dynamics while a lower value (flatter peak) shows a less complex flow dynamics with low particle-particle, gas-particle and particle-wall interactions based on shorter time scales depicted by sampling frequency. Figures 4-6 show the plots of the correlation dimension, D_{3}, (referred to D in figures) versus ln(r), at U_{g} = 5.5 m/s, 8.0 m/s and 10 m/s respectively at G_{s} = 300 kg/m^{2}s.

The D_{3} plots show a peaked shape at the centre (r/R = 0.0) with highest values compared to the profiles in the wall region (0.74 ≤ r/R ≤ 0.98). The Correlation dimension decreases from the centre of the riser towards the wall for all operating conditions investigated, i.e., U_{g} = 5.5, 8 and 10 m/s and G_{s} = 300 kg/m^{2}s. Similar observation have been reported in the literature [

Contrary to the wall region where there is low degree of freedom for particles, there is less interactions due to high solid concentration , presence of solid clusters and the boundary layer or wall effects, all of which limits particle’s motions in the gas-solid suspensions leading to low values of the correlation dimension. Studies show that the dilute phase and the dense or cluster phase have different dynamical properties and therefore shows different dynamic behaviors [

Lower peak values correlation dimension are generally found for conditions with high solids concentration, due to the fact that the presence of high solids concentration in the gas phase near the wall reduces the amount of turbulence and complexity of flow [

Towards the wall both phases are partially coupled system with cluster or dense phase dominating leading to low correlation dimension [

Results show a tendency of change in location of the peak value, i.e., the maximum value of D_{3} along the ln(r) axis (horizontal position). The maximum D_{3} shifts from left to right along the ln(r) axis towards the origin of the D_{3} versus ln(r) plane as one move towards the wall. Also for all elevations, towards the wall of the riser reactor (i.e., 0.74 ≤ r/R ≤ 0.98) correlation dimension plots show the presence of the multi-fractal behavior which manifests by the presence of more than one local maximum value indicating two or more peak values of D_{3}. Presence of two local maximum values indicates the bifractal behavior of the gas-solid flow. Such behaviors were observed and become more dominant towards the wall when the gas velocity, U_{g}= 10 m/s and solid flux, G_{s}, =300 kg/m^{2}s, as shown in _{g} = 8 m/s and solid flux, G_{s} = 300 kg/m^{2}s this behavior is found in the bottom section at Z = 1.52 m and r/R = 0.0, Z = 2.73 m and 0.74 ≤ r/R ≤ 0.87; in the flow developing section at Z = 3.96 m and 0.81 ≤ r/R ≤ 0.98 and also in the top section of the riser at Z = 9.42 m, 0.74 ≤ r/R ≤ 0.98.

According to _{g} = 10 m/s, G_{s} = 300 kg/m^{2}s (which is presented by

The single peak values of the correlation dimension profiles at the centre in all axial levels indicates the uniform flow structure or low tendency to the formation of segregated flow. The multiple peaks of correlation dimension profiles with varying magnitudes towards the wall indicates the presence of multifractal flow structures which can be attributed to the presence of dense phases dominated by particle aggregates with varying solid concentrations leading to the mixed flow behaviors, as stated above.

The analysis of the shapes of correlation dimension profiles reveals three different types, i.e., A, B and C as shown in

Type B profiles are double peaked profiles, which emerge for all operating conditions in various spatial locations of the riser. They are mostly found in the wall region, for instance, in the entrance section, Z = 1.52 m, close to the wall for which r/R = 0.92 and at r/R = 0.74, Z = 2.73 m. They are also found in the flow development section at r/R = 0.81 and 0.98 for Z = 3.96 m and Z = 3.16 m for U_{g} = 5.5 m/s and G_{s} = 300 kg/m^{2}s. At U_{g} = 8 m/s and G_{s} = 300 kg/m^{2}s such profiles are observed in the flow development section at r/R = 0.87 and Z = 3.16 m. When U_{g} was increased to 10 m/s at G_{s} = 300 kg/m^{2}s, the double peaked correlation dimension profiles appeared, for instance, in the entrance section at r/R = 0.81 for Z = 1.52 and 2.73 m, in the developed flow section at r/R = 0.81 and 0.87 for Z = 5.90 m and in the fully developed flow section r/R = 0.87 and Z = 8.74 m. Type B correlation dimension profiles indicate two dominant flow types in the wall region of the riser, attributable to relatively dilute phase and dense or cluster phase.

Type C correlation dimension profiles are multiple peaked profiles. They are found in some locations of the wall region for all operating conditions and in all sections of the riser. For instance, they are found in the entrance section at r/R = 0.87 and Z = 2.73 m; in the flow developing section at r/R = 0.87 and Z = 3.96 m; and in the fully developed flow section towards the riser exit at r/R = 0.87 and Z = 8.74 m and r/R = 0.74 and 0.87 for Z = 9.42 m for U_{g} = 5.5 m/s and G_{s} = 300 kg/m^{2}s. For U_{g} = 8 m/s and G_{s} = 300 kg/m^{2}s these profiles were observed in the entrance section at r/R = 0.87 and 0.92 for Z = 2.73 m. When U_{g} was increased to 10 m/s at G_{s} = 300 kg/m^{2}s these profiles were observed in the entrance section for instance at r/R = 0.81 and Z = 1.52 m; in the developing flow section at r/R = 0.74 and Z = 3.96 m; in the transition at r/R = 0.87 and Z = 6.34 m and in the fully developed section towards the exit section they are found at r/R = 0.81 and Z = 9.42 m. Multiple peak correlation dimension profiles indicate more than two dominant flow patterns, attributable to dilute phase, clusters and intensive interactions between these two phases with the wall.

Analysis of mapped correlation dimension profiles was used to investigate the change of gas-solid flow dynamics in the riser by observing shapes of the profiles, their positions on the ln(r) axis and the maximum value of D_{3}. Mapping of correlation dimension, D_{3}, on the D_{3} − ln(r) plane is shown in _{g} at three axial elevations. This approach provides a new method of studying the flow dynamics in a high-flux riser.

_{3} profiles along the ln(r) axis at U_{g} = 5.5, 8 and 10 m/s for G_{s}= 300 kg/m^{2}s in the entrance, developing flow and developed flow sections, that is, Z = 1.52, 3.96 and 9.42 m respectively. It was further observed that the D_{3} profiles at the center (r/R = 0.0) forms a higher peak which relatively increases towards the top of the riser in the fully developed flow section. The rest of the profiles in the wall regions are flatter at the top of the riser, compared to those in the flow developing and entrance sections. That is, in the entrance section the profiles close to the wall (r/R = 0.98) have higher peaks or higher correlation dimension than the profiles at the top of the riser. This is attributable to intensive mixing in the entrance section after the gas passes the distributor plate.

The position of the correlation dimension profiles at r/R = 0.0 along the ln(r) axis is relatively the same for all sections. However, the profiles in the wall region have the tendency to shift their positions along the ln(r) axis depending on gas velocity, U_{g}, and axial elevation, Z. This indicates that the scaling regions of the reconstructed attractors change their position and size. When the profiles shift to the right, it implies that attractors expand in size and the distances between trajectories also increases. It can be observed that the profile at r/R = 0.81 changes position towards the negative direction along the ln(r) axis from the entrance (0.5 to −2.5), developing (−0.5 to −3.2) and the developed (−2.5 to −5.2) sections at the top of the riser. The profiles also show a clear double peak in the flow development section (Z = 3.96 m). The profile close to the wall (r/R = 0.98) in the entrance section (Z = 1.52 m) is located in the same position (around −2.5 to −5.5) with the centre profile with its peak being the highest. This suggests resemblance in the gas-solid flow behavior. The peak of this profiles decrease along the high flux riser from the entrance section through the developing up to the developed section. This behavior is also shown by the profiles at r/R = 0.92.

At U_{g} = 8 m/s, the centre profile (at r/R = 0.0) shows a small shift towards the right side or positive direction along the ln(r) axis from the entrance section where it is located around −3.7 to −6.5 towards the fully developed flow section around −2.5 to −4.9. The correlation dimension profiles at the wall show the tendency to shift towards the negative direction along the ln(r) axis as Z increases as observed also at U_{g} = 8 m/s. In the entrance section, the profile at r/R = 0.81 is located from 0 to −3 and shifts to a location between −1 and −4 in the developing flow section and then from −3.4 to −6 in the fully developed flow section. The profiles close to the wall, at r/R = 0.92 and 0.98, are located in the same span in the flow developing section (1 to −1.5) and fully-developed flow section (−2.4 to −5.1) while in the entrance section the profile at r/R = 0.98 is located on the left (−1 to −4) of the profile at r/R = 0.92 (that is from 0.2 to −2.8). Also these profiles have low peaks compared to the other profiles.

Further _{3} profiles at U_{g} = 10 m/s and 300 kg/m^{2}s in the three flow development sections, that is, entrance (Z = 1.52 m), flow development (Z = 3.96 m) and fully developed (Z = 9.42 m) sections.

The tendency of D_{3} profiles to form multiple peaks closer to the wall at r/R = 0.92 and 0.98, increases especially in the entrance section at U_{g} = 10 m/s compared to the center profiles. These observations can be associated to the increasing solid concentration and formation of dense phases which reduces the degree of freedom of particles. Also the influence of the wall on the gas phase and solid particles becomes significant. All these reduce particles momentum and interactions with gas phase and the wall resulting into less complex flow behavior and hence low values of correlation dimensions towards the wall. Multiple peaks with lower varying magnitudes of correlation dimensions towards the wall are attributed to the presence of low dimensional multifractal flow behavior which results from the presence of dense phases which are dominated by clusters or particle aggregates with varying solid concentrations which leads to the mixed phases flow behaviors.

The span of the D_{3} profiles from left to right indicate changes in the size of the reconstructed attractor in phase space. For example, in _{g} = 5.5 m/s, at r/R = 0.92 and Z = 9.42 m, the profile spans within −1 and 0, while at Z = 1.52 m, the profile span from −2 to 1 along the ln(r) axis. A profile of D_{3} on the right hand side indicates that the scaling region(s) of the attractor is within longer trajectory distances, implying expanded attractor. The radial position r/R = 0.81 exhibits a wide range of the changes in the flow dynamics in the gas-solids flow along the axial elevations of the riser. This has been explained in literature using core-annulus model, which was based on changes in radial solids concentration and particle velocities [_{3} vary for the profiles from different locations and at different operating conditions. Also there are wide variations in the location of the profiles along ln(r) axis from 0 to −5.5 which forms on other property of gas-solid flow dynamics developed in this study. This indicates that the microscopic flow scales (originally captured by single solids concentration signal) diminishes from largest at the entrance to the smallest at the top of the riser. That is, at r/R = 0.81, changes occur from densest to the most dilute flow behavior along the riser.

Variation of the maximum correlation dimension, D_{max}, determined from the correlation dimension profiles of D_{3}, along the radial position at a given axial elevations and operating conditions are shown in _{max} along the axial direction are relatively higher at U_{g} = 5.5 m/s and G_{s} = 300 kg/m^{2}s except at at Z = 2.73 m and 8.74 m. Lower values were observed in the bottom section (Z = 1.52 m), transition (Z = 3.96 m) and top (Z = 9.42) sections. Also away from the wall towards the centre, D_{max} profiles are irregular revealing the non-uniform time-variations in the flow behavior of the gas-solid flow depicted by solids concentrations signals in the riser reactor. For r/R = 0.81, again, using normal radial profiles of D_{max}, one cannot notice the inherent dynamical changes along the riser in this radial position. Thus, the current study used D_{3} profiles to show that this radial position exhibit strongest variations in the gas-solids flow dynamics. Using the values of D_{max} could not reveal this phenomenon (

The effects of changing gas velocity on the correlation dimension and its profiles were studied. _{g}) on the Correlation Dimension (D_{3}) at different radial position (r/R) in selected sections of the riser, using profiles generated from local solid concentration.

_{g} = 8 and 10 m/s, the correlation dimension profiles are located within a higher hyperspherical radius of the phase space especially at r/R = 0.92 and 0.98.

In the flow development section, Z = 3.96 m, (as per _{g} = 5.5, 8 and 10 m/s are located in the higher hyperspherical radius along ln(r) axis, which signifies expanded attractors with longer trajectory distances near the wall.

As shown in _{g} = 10 m/s, the correlation dimension (D_{3}) profiles are located within a higher hyperspherical radius of the phase space especially at r/R = 0.81, 0.92 and 0.98. Such observation of changes in the dynamics have not been reported using average values of solids concentration signals or single value chaos invariants in the past, when r/R, U_{g} and Z were changed. This study focused on the dynamics of gas-solid flow in the wall region of the high-flux gas-solid riser while other studies considered the whole cross-section of the risers. The shift of maximum correlation dimension or scaling region along the ln(r) axis is being reported in this study for the first time.

The shifting tendency of the correlation dimension profiles is further summarized in _{3} is located along the ln(r) axis. The analysis is presented for different radial positions (the centre and the wall region) along the riser height from the entrance section (Z = 1.52 m), flow development section (Z = 3.96 m) and in the fully developed flow section (Z = 9.42 m). For instance, in the entrance section (Z = 1.52 m) close to the wall (r/R = 0.98), when U_{g} = 5.5 m/s the maximum D3 is located at ln(r) = −4.05, while for U_{g}= 8.0 m/s it is located at ln(r) = −3.38 and for U_{g} = 10 m/s it is located at ln(r) = −0.26.

The effect of the axial elevations was investigated in three sections of the riser, i.e., the entrance, developing and the fully developed sections. _{3}) at different radial position (r/R) for the selected operating conditions.

In _{3}) profiles at different axial elevation (Z)at U_{g} = 5.5 m/s and G_{s} = 300 kg/m^{2}s are shown. The peak values of the correlation dimension in the fully developed section (Z = 9.42 m) are lower compared to other elevations at the wall. Also at r/R = 0.81 as the height increases from Z = 1.52 m to 9.42 m the profiles shifts along the horizontal axis towards the lower hyperspherical radius. Close to the wall at r/R = 0.98 the correlation dimension profiles in the entrance section (Z = 1.52 m) has higher peak values at U_{g} = 5.5 and 8 m/s but are located at the lowest value of the hyperspherical radius which indicates contraction of attractor which implies higher solid concentration.

Correlation dimension profiles at different axial elevation (Z) for U_{g} = 8 m/s and G_{s} = 300 kg/m^{2}s are also shown in

Further, the correlation dimension profiles at different axial elevation (Z) for U_{g} = 10 m/s and G_{s} = 300 kg/m^{2}s is also shown in _{g} = 8 m/s and G_{s} = 300 kg/m^{2}s except at r/R = 0.92 and 0.98. The correlation dimension profiles at the centre of the riser (r/R = 0.0) in the entrance sections (Z = 1.52 m), flow

Z (m) | r/R | U_{g} (m/s) | Value of maximum D_{3} | Location of maximum D_{3} along ln(r) |
---|---|---|---|---|

9.42 | 0.0 | 5.5 | 7.36 | −3.31 |

8 | 5.50 | −3.47 | ||

10 | 5.60 | −4.04 | ||

0.81 | 5.5 | 3.35 | −4.64 | |

8 | 3.73 | −5.56 | ||

10 | 3.90 | −4.46 | ||

0.92 | 5.5 | 1.11 | −0.77 | |

8 | 2.15 | −3.86 | ||

10 | 3.13 | −0.35 | ||

0.98 | 5.5 | 1.77 | −0.30 | |

8 | 1.36 | −3.55 | ||

10 | 1.00 | −1.38 | ||

3.96 | 0.0 | 5.5 | 7.09 | −3.26 |

8 | 5.50 | −4.14 | ||

10 | 4.93 | −4.61 | ||

0.81 | 5.5 | 3.25 | −1.46 | |

8 | 3.97 | −3.08 | ||

10 | 5.00 | −2.28 | ||

0.92 | 5.5 | 2.79 | −1.58 | |

8 | 0.97 | 0.36 | ||

10 | 3.43 | 0.14 | ||

0.98 | 5.5 | 4.26 | −0.56 | |

8 | 2.14 | −0.11 | ||

10 | 1.07 | −0.26 | ||

1.52 | 0.0 | 5.5 | 5.05 | −3.98 |

8 | 5.68 | −5.24 | ||

10 | 4.94 | −4.31 | ||

0.81 | 5.5 | 4.35 | −1.21 | |

8 | 4.13 | −2.46 | ||

10 | 3.72 | −0.93 | ||

0.92 | 5.5 | 4.58 | −0.08 | |

8 | 6.16 | −2.83 | ||

10 | 2.31 | −0.05 | ||

0.98 | 5.5 | 6.80 | −4.05 | |

8 | 4.20 | −3.38 | ||

10 | 1.35 | −0.26 |

developing section (Z = 3.96 m) and in the fully developed section (Z = 9.42 m) shifts towards the higher hyperspherical radius as the height increases. Also in the wall region at r/R = 0.81 the trend is vice versa. The profile at Z = 1.52 m is located at higher hyperspherical radius and as the height increases the location of the profiles shifts to the lower hyperspherical radius. Profiles at r/R = 0.92 and 0.98 for all axial elevations studied are located in the same locations of the hyperspherical radius with higher peak value at r/R = 0.92 than at r/R = 0.98.

The effect of riser height can be observed further through the shifting tendency of the correlation dimension profiles as shown also in _{3} is located along the ln(r) axis in the centre and the wall region for all axial elevations of the riser. For instance, close to the wall (r/R = 0.98) when U_{g} = 5.5 m/s, in the entrance section, Z = 1.52 m, the maximum D_{3} is located at ln(r) = −4.05, while in the flow developing section, Z = 3.96 m, it is located at ln(r) = −0.56 and in the fully developed section, Z = 9.42 m, it is located at ln(r) = −0.30.

Gas-solid flow dynamics in the high-flux CFB riser was studied by mapping the correlation dimension, D_{3}, at different spatial locations for different operating conditions (U_{g}) and locations along the ln(r) axis. The number of peaks for each profile was also assessed. Based on the discussion, it can be concluded that the three-step correlation dimension (D_{3}) profiles are good indicator of the differences in the flow dynamics over various states of the gas-solid flow in the riser. These correlation dimension profiles at the centre are mostly single peaked with higher peak values than the wall region showing uniform flow dynamics. These peak values decrease towards the wall. Towards the wall these profiles exhibit double or multiple peaks profiles showing bifractal or multifractal flow behavior. As the velocity increases the correlation dimension profiles close to the wall become more random and irregular affirming chaotic nature in these locations. Further, towards the wall these profiles change location towards higher hyperspherical radius of the attractor at Z = 9.42 m. As the height increases the correlation dimension plot shifts towards higher hyperspherical radius at the centre and towards lower hyperspherical radius in the wall at r/R = 0.81. Mapping of the established correlation dimension (D_{3}) profiles forms a good tool of analysis of the flow dynamics in the high-flux gas-solid riser compared to other tools of analysis. The study recommends further analysis to other risers of different dimensions operated at high-flux conditions using the established method.

The authors declare no conflicts of interest regarding the publication of this paper.

Jeremiah, J.M., Manyele, S.V., Temu, A.K. and Zhu, J.-X. (2018) Mapping Correlation Dimension along the Wall Region of a High-Flux Gas-Solid Riser Using Embedded Solid Concentration Time Series. Engineering, 10, 655-679. https://doi.org/10.4236/eng.2018.1010048