Simulation Analysis of Air Flow and Turbulence Statistics in a Rib Grit Roughened Duct

The implementation of variable artificial roughness patterns on a surface is an effective technique to enhance the rate of heat transfer to fluid flow in the ducts of solar air heaters. Different geometries of roughness elements investigated have demonstrated the pivotal role that vortices and associated turbulence have on the heat transfer characteristics of solar air heater ducts by increasing the convective heat transfer coefficient. In this paper we investigate the two-dimensional, turbulent, unsteady flow around rectangular ribs of variable aspect ratios by directly solving the transient Navier-Stokes and continuity equations using the finite elements method. Flow characteristics and several aspects of turbulent flow are presented and discussed including velocity components and statistics of turbulence. The results reveal the impact that different rib lengths have on the computed mean quantities and turbulence statistics of the flow. The computed turbulence parameters show a clear tendency to diminish downstream with increasing rib length. Furthermore, the applied numerical method is capable of capturing small-scale flow structures resulting from the direct solution of Navier-Stokes and continuity equations.


Introduction
Turbulent flow is a complex phenomenon even for simple geometries like square obstacles mounted on a surface which finds application in many physical and engineering problems. Such flows occur across a wide range of applications from aerodynamics, heat exchange systems to solar air heaters. For the flow investigation, researchers usually conduct experimental and numerical analyses in both two and three dimensions. The analyses focus on flow either around the sides of a surface-mounted obstacle, which is characterized as three-dimensional, or around a surface-mounted rib where the flow is then characterized as two-dimensional provided that the rib extends to the walls of the duct/simulation domain in the spanwise direction.
Because of their diverse and fundamental significance, numerous investigations both experimental as well as numerical have been conducted to particularly examine turbu-lent flow around two-dimensional surface-mounted ribs. The geometry that is commonly used when studying twodimensional turbulent flows consists of square ribs. Such twodimensional obstacles have been employed from a number of researchers to deduce conclusions about the different aspects of the respective turbulent flows. For example, Acharya et al. [1], Hwang et al. [2], and Panigrahi and Acharya [3] examined different aspects of flow movement and how it is modulated when using a surface-mounted square rib. More recently, Panigrahi et al. [4] experimentally investigated the reattaching shear layer development behind a surface-mounted square rib. The same geometry has been also used by Liu et al. [5] and Shah and Tachie [6] to investigate the shedding and flapping of vortical flows and the dynamics of separated and reattached flows, respectively, caused by a surface-mounted rib. An extension to the above experimental studies was the investigation of turbulent flow modulation when the flow passes over a series of square rib elements to draw useful conclusions concerning turbulence characteristics [7] and the wall roughness influence on flow configuration and heat transfer [8][9][10].
From the above it is evident that the study of twodimensional, turbulent, separated, and reattached flow around ribs possesses a big portion of the researchers' interest due to the fact that such flow is frequently encountered in many technological and engineering applications. In particular, turbulent flow around ribs has practical importance when studying how fluid-related devices, such as gas turbines, turbo machines, and heat exchange systems, behave under certain flow conditions. Additionally, ribs serve as artificial roughness elements in a fluid flowing channel to improve the convective heat transfer by creating turbulence in the flow. This technique is of particular interest in renewable energyrelated applications, especially when investigating the surface roughness effect and heat transfer in solar air heater devices where roughness element geometry is the most important parameter influencing the turbulent flow characteristics [11][12][13]. However, the experimental investigation of such phenomenon in order to find the optimum roughness pattern among different geometries is often a time consuming and tiresome process [14]. In this regard, the scope of the present study is to assess the flow behavior over surface-mounted ribs of varying lengths by using a dynamic numerical approach based on computational aspects and technique which, in many aspects, saves time and reduces cost. More specifically, the current work is investigating how different rib lengths affect the two-dimensional turbulent flow and turbulence characteristics such as turbulent intensity and turbulent kinetic energy for the case of Re = 1000. Moreover, the statistical processing of the simulation results depicts the impact that different rib lengths have on the separation and The Scientific World Journal   reattachment points of the flow, which directly affect the rate of heat transfer ( [15], and references therein).

Simulation Setup Aspects
The results presented here are derived by executing a modified version of a noncommercial Fortran code initially developed by Fragos et al. [17], which employs the standard Galerkin finite elements method [18,19] to solve the governing equations along with the appropriate initial and boundary conditions. The code has been experimentally verified by investigating the airflow patterns around obstacles with arched and pitched roofs [20]. The fluid is considered to have constant viscosity ] and density . A uniform stream flow with velocity 0 is used as boundary condition at the entrance of the computational domain. The fluid moves uniformly through the duct until it hits the surface-mounted rib of height ℎ with a 90 ∘ angle of attack. No-slip boundary conditions are imposed along the walls of the duct and the rib where the fluid is decelerated. The outlet boundary condition is a free boundary condition which lets the fluid exit the computational domain freely without any distortion. The two-dimensional turbulent flow is assumed to have a constant temperature and the fluid to be incompressible. Under these assumptions the flow behavior can be accurately estimated by directly solving the dimensionless Navier-Stokes and continuity equations which represent conservation of momentum and mass, respectively. Consider where tensor notation has been used and , are dummy indices which go from 1 to 2 and stand for and coordinates, respectively. are the instantaneous components of the velocity vector, is the pressure, and Re is the Reynolds number.
is the total time, which in the simulation runs is 150 dimensionless time units. The above equations have been rendered dimensionless by choosing 0 and ℎ as the characteristic velocity and length, respectively. Re = 0 ℎ/] is the Reynolds number, where ] is the kinematic viscosity of the fluid and ℎ is the rib height. Additionally, the pressure and time have been normalized by using the terms 2 0 and , where is defined by the ratio ℎ/ 0 . The values of the input parameters yield a Reynolds number of 1000. However, it should be borne in mind that the flow is fully turbulent when using the duct height for the calculation of the Reynolds number. The computational meshes used in the present study for the variable rib lengths are shown in Figure 1. They consist of rectangular finite elements of different sizes with nine nodes in each of them. The unknown velocities and the pressure of the governing equations (1) are expanded in terms of Galerkin basis functions as where , are the quadratic and linear basis functions in each element, respectively.
The governing equations, weighted integrally with the basis functions, result in the following continuity, , and momentum, , residuals, where I is the unit matrix and T = ∇U + (∇U) is the stress tensor of the Newtonian fluid. By applying the divergence theorem in order to decrease the order of differentiation, (4) reduces to The volume integral in (5) along with (3) is evaluated at all the interior nodes of the computational domains. The inclusion of the surface integral in (5) at all exit nodes of the domains essentially results in the imposition of the free boundary condition. By not imposing any arbitrary boundary condition at the exits of the computational domains the outflow is given by the exact solution of the governing equations minimizing in this way any outflow distortion that may propagate towards the interior of the domains. The nonlinear algebraic equations system resulting from (3) and (5) is solved with the Newton-Raphson iterative scheme. In Table 1 flow domain tessellations are shown in relation to the number of finite elements, nodes, and unknowns.

Results and Discussion
The results are presented here in the following sequence: (1) instantaneous and mean streamlines, (2) mean velocity The Scientific World Journal where, in our case, the instantaneous flow quantities, , are obtained from the direct solution of the incompressible, dimensionless Navier-Stokes and continuity equations. The streamwise velocity 1 and transverse velocity 2 fluctuations are computed by decomposing the instantaneous velocities into their mean and fluctuating parts, Figure 2, for five rib aspect ratios (1 : 0.5, 1 : 1, 1 : 2, 1 : 3, 1 : 4) at = 70. At this particular time instant, the phenomenon of vortex shedding in its full development downstream of the rib is clearly demonstrated. In all cases, the isolated vortical structures are produced at the right corner of the rib and their shape changes in a random way as they travel downstream. Their continuous number and shape variation in the shear layer region lead to a periodical motion of the reattachment point, as has been previously demonstrated by Fragos et al. [21]. Vortices are also produced upstream and over the rib due to the separation of the flow with the over-the-rib vortical structures being more pronounced as the rib aspect ratio increases. A very interesting feature that one notices is the way vortices behave just downstream of the rib. While at small rib ratios, particularly for 1 : 0.5 and 1 : 1 ratios, merged-like vortical structures are observed that are rather difficult to be distinguished, greater rib ratios have the effect of severing the vortical flows resulting in well-formed, individual eddies.

Instantaneous and Mean Streamlines. A time instant of simulated streamlines is shown in
The separation point upstream of the rib appears to be unaffected by the change of the rib length. However, it should be mentioned that the upstream separation point is moving with time, changing the shape of the recirculation region, as it has been previously demonstrated by Fragos et al. [21]. Additionally, there is no apparent relationship between the ribs' length and the vortices' shape; they rather seem to change their shape in a random fashion without following any particular formation pattern. The calculated mean streamlines for the five different rib aspect ratios are shown in Figure 3(a). It can be observed that the shape of the recirculated flow upstream of the rib stays rather stable as rib ratio increases with the only exception of 1 : 2 rib ratio where recirculation flow demonstrates a considerable shrinkage. Concerning the recirculation zone downstream of the rib, it is generally decreasing in both extent and intensity; the longer the rib is the less extensive and less intense the recirculation zone is.
In Figure 3(b) zoomed-in averaged flow profiles for all five aspect ratios are shown. In this figure the capability of our numerical scheme of resolving the tertiary eddies being present in the lower corners of the ribs (thin open red circles) is demonstrated. The capturing of these small-scale vortical flows under time-averaged conditions has an additional value since it implies that these structures are a permanent (at least for the time scale we use) and not a transient feature of the adjacent-to-the-rib flows.
By closer inspection of Figure 3(a) we were able to define with high accuracy the length and thus the point of reattachment downstream of the ribs. The reattachment point was found by determining the -location where the dividing streamline reattaches on the floor. Figure 4 shows the linear fitting of the five reattachment points which are defined in relation to the rib lengths. The coefficient of determination, 2 , is estimated to be 0.9518 which is a relatively high value.

6
The Scientific World Journal This indicates that the linear fitting defined can be used as a prediction tool to estimate the lengths of reattachment under variable rib lengths. This is very useful especially in solar air heater applications where the recirculation region extent is the factor that mainly affects the heat transfer coefficient [11]. In order to strengthen the efficiency of our numerical method, in Figure 5, a qualitative comparison is made with the experimental results shown in the upper panel of Figure 3 in Agelinchaab and Tachie [16] for the square rib case. The comparison of the computed flow streamlines of the present study with the streamlines of the above laboratory experiment, for Re = 1920, and the sense of vortex rotation, with the downstream primary and secondary eddies demonstrating a counterclockwise and clockwise sense of recirculation, respectively, confirms that the present numerical model is simulating satisfactorily the studied airflow patterns. The difference in the recirculation length of the flow between the simulation and the experiment is due to the different upstream boundary conditions, the thickness of the boundary layer, and the turbulence level. More specifically, Agelinchaab and Tachie [16] are using a turbulence level of 5% which is very close to the turbulence level that is considered in the present study. Therefore the differences in flow conditions are mainly due to the different initial conditions. In both panels, two pairs of separations and reattachments are found with the exception of a small separation bubble evident on top of the rib for the numerical case (Figure 5(a)). Similar direct numerical simulations have been conducted using the authors' code for different obstacles geometries in wind tunnel/duct and the validity of the code and its respective numerical results were verified by laboratory experiments carried out in a wind tunnel [20]. Figures 6(a) and 6(b) show the streamwise and transverse mean velocity color-coded profiles, respectively, for the whole computational domain. As it is evident, the velocity profiles in the recirculation region downstream of the ribs deviate significantly from the upstream profiles, in qualitative agreement with the experimental results of Agelinchaab and Tachie [16]. The streamwise mean flow close to the base is almost equally distorted for all rib aspect ratios, although the dividing streamlines, as defined in Agelinchaab and Tachie [16], have significantly variable points of reattachment relative to the rib lengths, as it is verified by the relative extent of the recirculation zones when looking at Figure 3(a). By close comparison of left and right panels we can note that the streamwise mean velocity ⟨ 1 ⟩ is quite comparable to the transverse mean velocity ⟨ 2 ⟩ inside the recirculation region downstream of the ribs. This is particularly true at 1 distance unit downstream of the ribs where the velocity ratio is almost 1 providing a numerical confirmation for the experimental results of Agelinchaab and Tachie [16] who have demonstrated that for a square rib the classical thin shear-layer assumption where ⟨ 2 ⟩ is small compared to ⟨ 1 ⟩ is invalid in the recirculation region close to the base wall.

Turbulent Intensities and Reynolds Stresses.
The magnitudes of streamwise turbulent intensity √ ⟨ 1 2 ⟩ and transverse turbulent intensity √ ⟨ 2 2 ⟩ for the five different rib lengths are shown in Figure 7, while Reynolds shear stress −⟨ 1 2 ⟩ is shown in Figure 8(a). Streamwise turbulent intensity maximizes downstream of the ribs in the recirculation region close to the base wall, while its maximum value diminishes as we proceed to higher rib lengths. Similar to the streamwise, transverse turbulent intensity also peaks downstream of the ribs; however, the extent of the maximum values is much higher in both streamwise and transverse directions. It is also worth noting that the region of peak values tends upwards away from the base wall, while their magnitudes gradually diminish as we proceed to greater rib lengths, similar to the streamwise case.
The computed values of Reynolds stresses (Figure 8(a)) maximize downstream of the ribs and by reexamining Figure 3(a) we do observe that the regions of maximization coincide with the cores of the primary flow eddies. This is something expectable since Reynolds stress represents the average momentum flux due to the velocity fluctuations, characterizing the transfer of momentum by turbulence [22].

Turbulent Kinetic Energy. The turbulent kinetic energy (TKE) is given by the equation
which is characterized by the measured root-mean-square (RMS) velocity fluctuations, that is, the turbulent intensities already plotted in Figure 7. Hence, as it is evident in Figure 8(b), turbulent kinetic energy can be considered as the superposition of the corresponding panels in Figure 7. By close inspection we do see that TKE follows the same dissipation pattern as we proceed to higher rib lengths. Concerning the production of TKE, −⟨ ⟩( ⟨ ⟩/ ⟨ ⟩), the discrepancy between Reynolds stresses and turbulent kinetic energy is apparently attributed to the fact that the production of TKE is also affected by the strain rate (gradient of mean velocity) which in the case of incompressible flows is linearly dependent on viscous stress.

Conclusions
The present paper investigates the two-dimensional flow around a rib with five different aspect ratios using the Direct Numerical Simulation method. The study concentrates on the characteristics of the fluid flow passing through a duct simulating a solar air heater with single rectangular ribs as roughness elements being attached on one wall. The analysis clearly shows how the rib length can affect the extent of the recirculation region and thus the point of reattachment downstream of the rib. The extent and intensity of turbulence are more pronounced for smaller rib lengths, something which is confirmed by both the lengths of reattachment and the turbulent intensity profiles. This fact clearly suggests that The Scientific World Journal  The Scientific World Journal  The Scientific World Journal as we proceed to smaller rib element lengths a simultaneous increase in heat transfer efficiency will occur.
The analysis although carried out for single rectangular ribs of different aspect ratios can be extended to include cases of multiple ribs of different geometries, for example, pitchedor arched-shaped ribs. Such an analysis can contribute significantly to the investigation of how the pitch ratio and roughness element geometries would affect the flow pattern and consequently the heat transfer coefficient.