do not necessarily reflect the views of UKDiss.com.
The unsteady flow features behind a heliostat in a narrow channel at a high Reynolds number: Experiment and Large Eddy Simulation
ABSTRACT
Several heliostats are arranged on the field for concentrating beam solar radiation onto a receiver. The flowinduced dust deposition and vibration of heliostats will lead to failure of the receiver and the entire system in eventuality. In view of this, fluid flow analysis behind a heliostat model in a narrow channel of limited length is presented depicting a fieldtype condition. For this purpose, experiment and large eddy simulation (LES) are performed. Laser Doppler velocimetry (LDV) technique is used for horizontal velocity measurement around a heliostat with 25° angle of inclination at a Reynolds number of 60000 basing on the length of an inclined plate. The wallresolved coarse and fine grids are utilized for LES with dynamic SmagorinskyLilly and subgridscale kinetic energy transport equation based models. Comparison between local instantaneous and statistical quantities, like horizontal velocity and its rootmeansquare values are presented. The flow development in the wake region of heliostat model is described using the timedependent and the statistical horizontal velocity. Based on the adopted LES approach the coherent structures are identified with the λ_{ci}criterion. The analysis of power spectral density and swirling strength isosurface reveal the responsible coherent structures in the inertial subrange for transport of energy to smaller structures. Interestingly, the presence of hairpintype vortices is observed extending from the two sides of model heliostat using the λ_{ci}criterion. Further, instantaneous spanwise vorticity contour reveals their alignment with flow represented by streamlines. The provided insight to these flow features and their organization is expected to frame strategies for mitigating the unsteady flowinduced vibrations and dust depositions on heliostats.
 INTRODUCTION
Concentrated solar thermal (CST) based technologies are versatile, in view of their applications. One such novel application is solar convective furnace system, which was recently proposed by Patidar et al. (2015). In this system, heliostat based central receiver is used for generating hot air. For this purpose, an open volumetric air receiver was specially designed by Sharma et al. (2015a, 2015b). Such a system will be useful in desert regions, like Rajasthan having direct normal irradiance of 56 kWh/m^{2}/day. However, a high windspeed of nearly 47m/s, see e.g. Lakshmanan et al. (2009), may cause vibration of heliostat and even initiate the saltation process. This will lead to deposition of dust on heliostats and reduction in their reflectivity, see e.g. Niknia et al. (2012). The soil composition of Rajasthan has high fraction of fine sand particles (< 60μm) Gupta (1986). These are prone to saltation and deposition on the surface of heliostat. This was verified by SEM based analysis in Yadav et al. (2014a). The deposition of dust on an inclined square finite plate was investigated by Yadav et al. (2014a, 2014b) in a wind tunnel. They demonstrated using experiment and simulation that the deposition is affected by location and separation between the inclined finite flat plates at a given angle of attack. Thus, for a better understanding of the deposition process, both experimental and numerical studies on fluid flow around an inclined flat plate are reviewed in this section.
Nature of fluid flow behind a bluff body, such as flat plate or a circular cylinder, is very complex including separation of flow, unsteady vortex shedding and presence of small and large –scale turbulent structures. The pioneering experimental investigation by Fage and Johansen (1927) was conducted at 18 different angles of inclination ranging from 0.15^{0} to 90^{0}over a wide range of horizontal velocity. They reported a relation between angle of inclination and vortex shedding frequency for infinite flat plate. The Strouhal number of 0.148 for an angle of attack between 30° to 90° was calculated based on the projected length of the inclined flat plate normal to the flow direction. They also explained the mechanism of vortex shedding from the leading and trailing edges at the same rate. Roshko (1954b) extended this study to cylinder and wedge in addition to a flat plate. A wake based Strouhal number was defined and observed to be ubiquitous compared to the one based on the projected length. The freestreamline theory based notchedhodograph method by Roshko (1954a) was used to estimate the drag from vortex shedding frequency. This method was found by Wu (1961) to yield results in agreement with that of Fage and Johansen (1927) at angles of inclination greater than 30^{0}. A new model was proposed to understand the twodimensional wake region behind an obstacle. This model accounted for coefficients of lift, drag and normal force as stated by Fage and Johansen (1927). Abernathy (1962) extended the freestreamline theory of Roshko (1954b) for an inclined flat plate. Separation between freevortex shear layers was found to be independent of flow constriction. This distance was used to devise a wake based Strouhal number. These wake based Strouhal numbers were mentioned besides the customary Strouhal number in Roshko (1954b) and Abernathy (1962). Perry and Seiner (1987) experimentally investigated turbulent wakes past nominal twodimensional bluff bodies along their plane of symmetry. They concluded that the flow was still twodimensional in view of the observed differences with twodimension based theory. This was attributed to the waviness in vortex cores along the spanwise direction. Kiya and Matsumura (1988) showed that the shearing stresses in the turbulent wake behind a normal thin beveled flat plate is mainly due to the random velocity fluctuations with half of the vortex shedding frequency. Chein and Chung (1988) determined the strength and location of vortices shed from an inclined flat plate using the Kutta condition. Their discrete vortex method predicted both macroscopic and microscopic features better than previous numerical practices. An experimental study for estimating Strouhal numbers with rectangular bluff bodies with at various aspect ratios was conducted by Knisely (1990). The authors investigated drag and lift characteristics with varying angles of attack. A rapid increase in Strouhal number was observed at low inclinations, which was associated with reattachment of separated shear layers. The dynamics of flow in the nearwake region behind a vertical flat plate was experimentally investigated by Leder (1991). Fluid mixing was reported in the separation region. The deficiency of using isotropic eddyviscosity approaches in some of the modeling approaches was concluded. Systematic two and three dimensional numerical simulations of flow over a flat plate are performed to investigate drag and vortex shedding. For instance, Najjar and Vanka (1995a) used a range of Reynolds number (Re) flows in twodimension to study the drag and wake characteristics including shedding frequencies and compare the results with previous studies. They stated that the flow become threedimensional at a Re ~200.Later on it was established with threedimensional direct numerical simulation (DNS) by Najjar and Vanka (1995b). They showed that the calculated drag agreed well with experiments. A study focusing on transition from two to threedimensional flow behind various bluff bodies was made by Thompson et al. (2006).They identified two modes of instabilities having onset at Re ~107 and 125 for a normal flat plate. Thus, at a higher Re, say beyond 10000, extremely complex flow structures and the related instabilities are expected behind a flat plate. Najjar and Balanchandar (1998) employed threedimensional computation to correlate fluctuating aerodynamic coefficients with vortex shedding frequencies.
Laser Doppler anemometry (LDA) investigation by Lam (1996) showed that the trailing edge vortices, around a plate inclined at an angle of 30^{0} and Re ~30000, are more dominant than that of leading edge. Lam and Leung (2005) performed similar experimental investigations using particle image velocimetry (PIV) at 20^{0}, 25^{0} and 30^{0}angle of inclinations at a Re ~5300. They reported that the vortex formation at the leading edge involves more complicated mechanism than the KevinHelmholtz instability of trailing edge. The trailing edge vortices showed to retain more energy than their leading edge counterparts in wakeregion at a point. Lam and Wei (2010) conducted unsteady Reynoldsaveraged NavierStokes (URANS) simulations for analyzing these observations. Their analysis provided an insight to the differences in the vortex shedding processes at the leading and trailing edges. They also explored how vortex shedding induces fluctuating lift and drag on the plate as in Najjar and Balachandar (1998). Narasimhamurthy and Andersson (2009) performed DNS of the wake behind a normal flat plate to investigate turbulence characteristics. Statistics like mean flow and Reynolds stresses exhibited characteristics similar to that of a cylinder. A mean recirculation zone was found to extend upto twotimes the height of flat plate in the wake region. They also explored the relation between base pressure and vortex formation. They revealed that the reduction in base pressure increases the vortex shedding frequency with a shorter recirculation zone. Yang et al. (2012) conducted two and three dimensional DNS at Re ~1000 and mentioned the reliability of the threedimensional approach. They found unequal vortices shedding from leading and trailing edges behind the inclined plate. Saha (2007) reported a Hopf bifurcation using DNS in the separated nearwake behind a normal flat plate that instills unsteadiness at a Re ~ 30 – 35. Further, the steadiness of farwake for a Re ~75 – 140 was reported to become unsteady at a Re = 145. This was attributed to largescale vortices with low frequency. This is similar to that of the reported flow characteristics in Thompson et al. (2006). A tertiary vortex street was observed to appear in farwake region at a Re = 150. Their frequency decreases with increasing Re. The drag and Strouhal numbers were reported to be nearly unaffected by these transitions. Further investigation by Zhang et al. (2009) revealed various transition zones between steady and chaotic flow through Hopf and subsequent bifurcations. This was verified by calculating the largest Lyapunov exponent value, which is nearly zero for nonchaotic flows, else distinctly positive. Drabble at al. (1990) studied the streamwise fluctuating forces on a flat plate normal to the flow direction with inlet velocity fluctuations. The response was defined through an admittance function, which increased with fluctuations in coherent flows and decreased for turbulent flows. This is relevant for further studies with varying inlet velocity to mimic wind gusts as in Lakshmanan et al. (2009). Indeed, DNS is limited only to low Re in view of computational requirements, see e.g. Chandra and Grötzbach (2007, 2008).
URANS, Detached Eddy Simulation (DES) and Large Eddy Simulation (LES) are some of the available options, especially, for high Re. Their degree of simplicity leads to uncertainties, which are being evaluated. Breuer et al. (1998, 2000) used various subgrid stress models for LES analysis of flow past a circular cylinder at Re~3900 and 140000. It was found that the importance of models increases with Re. Breuer and Jovicic (2001) conducted LES using the Smagorinsky (1963) model of flow around a flat plate inclined at 18^{0} with Re ~20000. These results were consistent with the experimental findings by Lam (1996). They found that the Smagorinsky model has drawbacks in transitional flows in the nearwall region and for high Re. At high Re turbulent flows Smagorinsky coefficient should be local time dependent. In order to calculate this parameter as a function of space and time, Germano et al. (1991) proposed a dynamic approach. The phenomenon of backscattering, which is found in transition flows, was addressed in this model. This serves as an extension to Smagorinsky model. Thus, the dynamic subgrid stress model proposed by Germano with Lilly’s modification Lilly (1992) and even with kinetic energy transport model as in Kim and Menon (1997) can be used to capture the physics of separated flows behind the heliostat plate. Breuer et al. (2003) employed DES, URANS in addition to LES on a coarser grid to evaluate their relative capabilities and found LES that adequately describes the flow characteristics. Afgan et al. (2013) carried out extensive LES based fluidflow analysis over a flat plate using the Smagorinsky model with Van Driest walldamping to represent the flow in nearwall region as in Moin and Kim (1982). These are summarized in Table I.
It may be noted that most of the experiments are performed for Re ≤ 20000. Indeed, a few wallresolved LES based flow analysis is reported at a higher Re, mostly, with an inclinedisolated flat plate or around a cylinder. The current work aims at understanding flow around a heliostat that comprise of a finite square inclined flat plate, stand and base at a higher value of Re. Heliostats are radially located on the field in a staggered layout. The distance between these heliostats varies depending on their radial positions. Thus, the reported observations using an isolated infinite or finite flat plate only serves as a reference point. The flow features behind a heliostat model is expected to be more complex considering the presence of surrounding heliostats or walls. For analyzing flow features behind a heliostat, based on recommendations wallresolved LES approach is adopted. In addition, experiments are performed for validating the adopted numerical approach. These are described subsequently.
 EXPERIMENT
Experiments are performed with one heliostat model to analyze fluid flow around this object. For instantaneous local measurement a twodimensional laser Doppler velocimetry (LDV) technique is used, see Figure 1. LDV is a wellknown nonintrusive in which the Doppler frequency of the seed particles crossing the formed interference fringe by unshifted and shifted laser beam, as shown in this figure, is measured. The product of spacing between two adjacent fringes and the Doppler frequency provides the velocity component. This is in line with the plane formed due to this dual laser beam crossing.
(a) 
(b) 
FIG. 1. (a) Image of beam crossing, showing fringes, idealized seed particle and intensity signal at photo detector and (b) twodimensional LDV setup at IIT Jodhpur 
 Geometry and Conditions
In the reported flow measurement experiment a single heliostat is considered. Heliostat model comprise of three distinct structures, namely, inclined square flat plate (P), stand (S) and base (B). The model is placed at an inclination of 25^{0} with respect to the horizontal plane. The plate is located at 7 mm above the stand inside the test section of a subsonic windtunnel as shown in Figure 2. The angle of inclination is comparable to the latitude of Jodhpur and serves as the reference point for further investigations. Indeed the angle of inclination will vary throughout the day in order to concentrate solar radiation to a fixed absorber. Preliminary analysis revealed its variation from 25 to 65°. The dimensions of experiment heliostat model are given in Table II.
TABLE II. Dimensions (in mm) of heliostat model
Components  Dimensions (length × width× height/thickness) 
Plate (P)  60 × 60 × 3 
Stand (S)  18 × 6 × 50 
Base (B)  74 × 68 × 20 
(a)  (b) 
FIG. 2.(a) Experimental heliostat model with (b) twodimensional LDV system 
A schematic of the measurement section at the symmetry plane is shown in Figure 3(a). In the performed experiment a diffuser with 5° divergence angle is attached to this section having a length of about 10 times the length of inclined plate. The constant freestream average inlet air velocity at the wind tunnel is set to U^{0 }= 16 m/s. The required length scale (D) for calculating Re_{D}, length (L) and height (H) of the windtunnel are given in Table III. The heliostat model is placed inside the windtunnel with its center at a height H_{c} = 82.5 mm. In the measurement section, the domain size along horizontal (z), vertical (y) and spanwise (x) direction is 8.3D, 2.5D and 2.5D, respectively. The domain length along the horizontal direction is comparable to Breuer and Jovicic (2001). Moreover, the selected length of test section represents separation between heliostats along radial direction at some azimuthal angle in a fieldlayout. The narrow channel, especially in the spanwise direction, depicts the presence of neighboring heliostats with stand as wall. Thus, the selected domain size is an attempt to simulate flow around a single heliostat in the simplest possible way. The outcome will serve as the starting point for further detailed investigations with multiple heliostats. Time dependent horizontal velocity component is measured along the vertical direction at a distance of 5 mm, upstream of the leading edge and 5, 10, 30 and 200 mm downstream of the trailing edge of heliostat model. The Reynolds number
, over a time step Δt with spatial discretized length of Δx. Subscript “i” means center of element with Δx/2 length on each side, “r,” “l,” “t” and “b” mean right, left, top and bottom. Least squares cellbased gradient evaluation assumes that the solution varies linearly between adjoining grid elements. Suppose a vector r_{i} joins the cell c_{0} to adjoining c_{i}, then the gradient of a variable ϕ can be represented as:
(23) 
To obtain the cell gradient, let . Each of these components is calculated as; here,W_{i0} is a weight factor and is determined by GramSchmidt orthogonalization, Anderson and Bonhus (1994). The simulations have been performed for a total of 17 flow through times considering Afgan et al. (2013). Out of this 13 flow through times based LES analysed data is used for instantaneous and statistical analysis.
 DISCUSSION AND RESULTS
 Validation
The measured and LES based timedependent horizontal velocity is compared at 5 mm downstream of the trailing edge of the inclined plate. This serves as the first test of the adopted LES approach. The measured horizontal velocity consists of more than 1000 LDV bursts over 90 seconds following Lam (1996).
(a) Experiment – full data set  (b) Experiment – partial data set 
(c) Coarse grid –k equation model  (d) Coarse grid – DSGS model 
(e) Fine grid – k equation model  (f) Fine grid – DSGS model 
FIG. 9. (a), (b) Measured and (c), (d), (e) (f) LES analyzed horizontal velocity at 5mm downstream of the trailing edge 
Figure 9(a) clearly depicts that the flow is, practically, independent of time. For an elaborated view the measured horizontal velocity for 10 seconds is shown in Figure 9(b). The random oscillations are, mostly, within ±69 m/s. Furthermore, LES analyzed timedependent horizontal velocity at the same location is shown for the selected models and grids in Figures 9(c) to 9(f). These are analyzed over 13 flowthrough times. As in the experiment, the LES analyzed oscillations are bounded within ±610 m/s. For a reference, dotted lines are provided at ±6 m/s. Thus the applicability of LES approach to capture flow features behind a heliostat model is inferred. The time resolution in LES is much better than that of the experiment. This is evident even from visual inspection of Figures 9(a) to 9(f). Therefore, LES is expected to provide an acceptable estimate of St, if any, for shedding of vortices having the same or different characteristic frequencies. Analysis of LES data shows repeated occurrence of maxima and minima in horizontal velocity as indicated in Figure 9(d) by vertical dotted lines. However, periodicity is not very obvious. This may be associated with transport and interaction of multiple largescale structures at the measurement location having different frequencies. This is expected in view of different structures in the modeled geometry viz. P, S and B as in Figure 2. In view of lack of any apriori estimation phaseaveraged values are not adopted for analysis. More detailed analysis is provided for analyzing the underlying phenomena. The range of fluctuations is summarized in Table VI. A comparative qualitative and quantitative assessment reveals that LES analyzed values are well within the measured data and thus suitable for further analysis.
Further qualitative and quantitative evaluation of the adopted LES approaches is performed using the statistical horizontal velocity and the corresponding RMS values at various horizontal positions along the centerline of heliostat model. These locations are depicted in Figure 10(a) and numbered 1 to 7. Their unequal spacing can be inferred from the indicated relative positions. The first point indicated by 1 is located at 5 mm and the last point indicated by 7 is located at 200 mm downstream of the trailing edge of inclined plate. As explained, the locations 1 and 7 are in near and far wake region, respectively. For statistical analysis of LES analyzed data thirteen flowthrough times are considered. A comparison between LES and experimentally analyzed statistical horizontal velocity is shown in Figures 10(b) and (c), which depict a qualitative agreement. The coarse and finegrid based LES analyzed data overpredicts and underpredicts the measured statistical horizontal velocity values, respectively. The increasing trend of statistical horizontal velocity is predicted by LES. The lowest value of horizontal velocity is a signature of highest level of fluctuations and vice versa. Mean flowfield development downstream of the heliostat model is inferred from increasing statistical horizontal velocity. This behavior is captured by LES and experiment. The extent of fluctuations indicate turbulent nature of the flow past heliostat and is quantified using statistical RMS values of horizontal velocity at the considered horizontal positions. The standard deviation in the measured horizontal velocity and its RMS values are estimated using eight temporal datasets including 100, 200, 300 upto 800 measurements. The relative standard deviation of mean and RMS values allowed us underlining the measurement uncertainty. The maximum standard deviation or uncertainty is about 9% for the mean horizontal velocity at 5 mm and about 10% for the RMS of the horizontal velocity at 200mm downstream of the trailing edge of inclined plate.
(a) 

(b)  (c) 
FIG. 10.(a) Measurement locations along centerline downstream of model; measured and LES analyzed statistical (b)horizontal velocity and (c)RMS of horizontal velocity at these locations 
The comparison between measured and LES analyzed statistical RMS values within 3.54 m/s of horizontal velocity clearly depict both the qualitative and quantitative predictive capability of the modeling approaches. High values of statistical RMS horizontal velocity are observed upto 30 mm and the smallest value is obtained at 200 mm from the trailing edge of inclined plate. Thus, it is inferred that indeed the strong flow disturbance dissipates from near to farwake region. The underprediction of RMS values with LES beyond 30mm, within uncertainty bound, may be attributed to the applied atmospheric pressure outlet condition considering a realistic separation between heliostats in a fieldlayout. The LES analyzed statistical and RMS horizontalvelocities are practically independent of the selected sgsmodels for the fine grid. This observation is interesting in view of DNS in which models are not expected to influence the flow features. Furthermore, the LES index of quality (LES_IQ) as explained in e.g. Geurts and Froehlich (2002), Celik et al. (2005) is computed using the following expression:
here, is the Kolmogoroff microscale of length, which is obtained using the characteristic lengthscale of the testsection (H=150 mm) and h indicates the smallest gridsize near the wall. The computed value of and is acceptable for a LES. Using the suggested values of results in and is an indicator of a good LES. For DNS this parameter should be greater than 0.95. However, the value is beyond 0.8 for LES, which suggests that ca. 80% of the energy is resolved see e.g. Pope (2009) and even a step closer to DNS, as observed.
 Discussion
The LES analyzed instantaneous swirling strength isosurface or λ_{ci}cretiron as in Zhou et al. (1999) corresponding to 900Hz using the DSGS model for coarse and fine grids is shown in Figure 11 at 0.34 second. This criterion identifies similar structures as Q or λ_{2}, Green et al. (2007). The relationship between these criteria is explained by Charaborty et al. (2005). Swirling strength, given by λ_{ci}, is a measure of the local swirl rate inside the vortex. In other words, the time period for completing one revolution of the streamline is given by 2π/λ_{ci}. The identified coherent structures appear to be in the inertial subrange, which is inferred by analyzing power spectral density in Figure 14. Such structures are responsible for transport of turbulent energy to smaller scales, as investigated by Terashima et al. (2015) for a planer jet. The presence of classical horseshoe vortex is clearly visible in front of base of the heliostat model, see Figure 11(a). It is evident that owing to a higher spatial resolution much smaller gridresolved structures are captured by the fine grid as shown in Figure 11(b). This is attributed to lower values of gridbased local Courant number for the finegrid in comparison to the coarsegrid. In the case of finegrid, except in the walllayer, this is found to be less than 1. The aligned coherent structures with the flow direction are visible, especially near the lowerwall of test section, for LES with fine grid. Both the coarse and fine grids based swirling strength isosurface show the emanating coherent structures from the two sides of finite inclined flat plate. Their nonlinear interaction is leading to eventual formation of smaller structures. The smaller structures are clearly visible in Figure 11(b). It may be safely concluded that the overall flowfield is captured by both the gridbased LES simulations. Thus, for a clearer insight to the flowfield and on vortex shedding the coarse grid based wallresolved LES is used.
(a)  (b) 
FIG. 11. Swirling strength isosurface (λ_{ci}criterion)for (a) coarse and (b) fine grids with DSGS at 0.34 second. 
In Figure 12, the vortices shed from the trailing edge of the plate are shown for the DSGS model at an interval of 0.01 second. The selected timeperiod corresponds to the estimated vortexshedding frequency of about 109Hz from an inclined plate. Figure 12(b) and 12(c) indeed shows the detachment of coherent structures from the trailing edge during 0.32 to 0.33 second. Finally, at a later instant the attached structure on the trailing edge of plate reappears in Figure 12(c).Thus the entire series of event leads to continuous vortex shedding from the trailing edge. As expected, the flowfield around heliostat is asymmetric along vertical direction, which is also inferred from the measured mean horizontal velocity, say, at a distance of 5 mm downstream of the leading edge in Figure 4. For completeness the fine grid based LES analyzed instantaneous swirling strength isosurface is shown in Figure 12(d). As in the coarse grid, at this instant of time, attached structures on the trailing edge of the finite inclined flat plate is observed. Further, the vortex path is depicted from the leading edge into the wake region and is explained with more details. As a result of higher grid resolutions, the resolved smaller structures are captured. These vortices will exert dynamic loading on a structure placed in the wake of heliostat, Lam et al. (2005). It was found in Lam (1996) that the trailing edge vortices are stronger in sense that they carry more energy than that of leading edge. However, the current study involves a square plate, stand and base leading to even more complex interactions and resulting flow structures.
Interestingly, amidst these multiple vortices, hairpintype structures are visible behind the square inclined finite flat plate in the LES analyzed data at 0.34secondin Figure 13(a) and 13(b). Similar structure with flow past a sphere was reported by Johnson and Patel (1999). The hairpintype vortex extends from the twosides of the inclined square finite plate rolls upto the center forming a closed structure as depicted with swirling strength isosurface corresponding to 650 Hz in Figures 13(a) and 13(b). Indeed the λ_{ci}criterion identifies such structures as expected, Zhou et al. (1999).The velocity vector with horizontalvorticity contour is shown in Figure 13(c) and 13(d) at a distance of z/D = 1.0 and 2.0 from the trailing edge of inclined plate. The observed counterrotating vortices on these planes are regarded as a signature of hairpin vortex as explained by Adrian et al. (2000). It appears that these structures ensemble forming larger structures as inferred in Pope (2009). It is further observed that their angle of inclination with horizontal direction increases with distance from the trailing edge. Apparently, these structures are transported with the flow enclosing smaller structures and are extending upto 3 times the length of plate i.e. z/D=3.0. Nonlinear interaction between the vortex structures at the trailing edge as in Figure 11 and hairpintype structures as in Figure 13 may lead to (a) generation of smaller structure, (b) excludes the presence of a unique Strouhal number and (c) the energy redistribution among structures with different frequencies. This aspect will be analyzed using frequency domain analysis of timedependent velocity.
To understand the transport of coherent structures and the related energy dissipation process the power spectral density (PSD) is used. It employs the fast Fourier transform (FFT) to analyze the temporal data on a frequency domain. This may provide any prevailing dominant frequency at a particular location, which is essentially an indicator of vortex shedding. This is defined as,
.  (24) 
Where, N is the sample size, u(t_{k}) is the horizontal velocity at time t_{k} and f_{s} as the sampling frequency.The PSD of LES analyzed horizontal velocity using oneequation model and finegrid is shown in Figure 14.These are performed at different horizontal positions, viz. 5, 10, 30 and 200 mm downstream of the inclined square flatplate along its centre line as shown in Figure 10(a). The energy is distributed over a widerange of frequencies upto 200 Hz. The LES analyzed energy spectrum of horizontal velocity with a slope of about f^{1/2}in this region is attributed to the presence of verylarge scale motions Katul et al. (2012), which is unlike f^{1} scaling as in Perry and Marusic (1995). The experimental specific kinetic energy upto 100Hz based on PSD is underpredicted by about 15% than that of LES analyzed values. As expected, f^{5/3} depicts that the inertial subrange and is captured with LES as in Calaf et al. (2013) in between 200 – 1000 Hz. Furthermore, the selected gridresolution resolves structures beyond this region. Thus, this is an extremely well resolved LES; at this position the spatial resolution is near to a DNS. Unusual is that the inertial subrange is rather short in terms of frequency range. This poses a question that whether such a fluid flow can be termed as turbulent. The second argument would be that such a scale, say beyond 1000Hz, is affected by the diffusive nature of these sgsmodels and numerical approach. Indeed energy levels decrease with increasing distance along the flow direction, which is inferred from PSD at the selected locations along the centreline of model. This indicates nonlinear interactions leading to generation of smaller structures, which eventually dissipate the energy, Tennekes and Lumley (1972). As expected, PSD of LES analyzed turbulent kinetic energy provided the same slopes and is not included avoiding repetitions. Apparently, at the considered locations along the centreline dominant energy containing structure is not observed and thus, a unique St is not expected.
f ^{1/2}
f ^{5/3} 
FIG. 14. PSD of LES analyzed horizontal velocity with finegrid and kequation model at 5, 10, 30 & 200mm downstream of heliostat model. 
In order to understand the vortex dynamics in the wakeregion LES analyzed positive and negative spanwise vorticity contours are shown in Figure 15 at certain instants of time. These are shown at the midplane of heliostat model. The considered instants account for the one flowthrough time. The flow is indicated by instantaneous streamlines colored by vorticity. The vortices are found to be aligned with the streamlines and are thus inferred as convected downstream in the wake region of the model as shown for all instants of time. This is indicated by dotted arrows and is consistent with findings using PIV based measurement by Lam and Leung (2005). Such interaction between shear and eddies lead to generation of turbulent stresses and are consistent with observations by Breuer and Jovicic (2001), Lam (1996).The negative values of spanwise vorticity correspond to clockwise rotating vortices. The positive values of spanwise vorticity indicate counter clockwise rotating vortices. Figures 15(a) to (f) show the attached vortex on the trailing edge is rotating in the clockwise direction and that of leading edge along the counter clockwise direction. Careful observations at the trailing edge show the stretching of clockwise rotating vortices, detachment, formation and transport of localized vortices from 0.31 to 0.34 second along the streamline. Behind the stand on this plane the nonlinear interaction of vortices leading to formation of smaller structures is inferred. The presence of clockwise rotating vortices is observed behind the base. Furthermore, the analysis of positive spanwise vorticity at 0.32 second indicates the presence of horseshoe vortex in front of the base. Dust particlevortex interaction will be analysed at the next step and reported.
The LES analyzed Reynolds stresses viz. crosscorrelation of velocity components are computed using twelve flowthrough times. These are depicted in Figures 16(a) to (c). These clearly depict the influence of vortices, as described in Figure 15, especially in the component. The other two components are localized in view of the affected region behind the heliostat model. The positive values of Reynolds stresses depicts that the velocity components are, more frequently, in the same direction. The negative values depict the opposite effect. Analyzing these quantities allows us to conclude that indeed the dominating will lead to shear stress on the wakeaffected heliostat or mirror surface and its erosion with wind and dust. Thus, mitigating the same would be beneficial for a durable coating or mirror surface.
(a)  (b)  
(c)  
FIG. 16. LES analyzed specific Reynolds stresses; (a)(b) and (c) 
 CONCLUSION
Unsteady flow features behind a single heliostat in view of a fieldlayout is analyzed using the laser Doppler velocimetry and the large eddy simulation (LES). For this purpose a narrow channel of a limited length is considered indicating the presence of neighboring heliostats in a fieldlayout. The local measurement of the horizontal velocity depicts strong fluctuations at 5mm from the trailing edge of heliostat model. The decreasing level of oscillations from neartofar wake region is inferred from the measured and LES analyzed rootmeansquare values of horizontal velocity. Thus, dissipation of energy is inferred downstream of heliostat in the wake region. Analysis of swirling strength isosurface or λ_{ci}criterion reveals that LES with finegrid resolves smaller structures owing to lower values of Courant number with the same timeresolution. The observed hairpinlike structure in the wake behind heliostat is consistent with the presence of velocity vector at two different vertical planes showing the presence of counter rotating vortices. These are transported along the wake and may interact with heliostat in this region. Such coherent structures are regarded as responsible for energy transfer from large to smaller structures. The intensity of which may lead to its vibration. The PSD of LES analyzed horizontal velocity reveals a slope of about f^{1/2}in the low frequency region up to 200Hz infers the presence of very large scale motions. An indication of which is found in local timedependent horizontal velocity with timescale of the order of flowthrough time. Further analysis of spanwise vorticity reveals the alignment and transport of counter and clockwise rotating vortices along the streamlines in the wake region. Thus, a heliostat in the wake will experience such unsteady flow features. Further analysis with multiple heliostats in a fieldlayout will be performed to verify the observed flowfeatures and to establish their influence on heliostats. However, it may be safely concluded, while designing a heliostat field, flow based parameters are to be established. Such parameters will hold the key to mitigate flowinduced dust deposition and vibration effect on heliostat. This will, in turn, enhance the performance of concentrated solar thermal system based on heliostat in near future.
REFERENCES
Abernathy, F. H., “Flow Over an Inclined Plate”, J. Basic Eng. 84(3), 380388 (1962).
Adrian, R. J., Meinhart, C., D., and Tomkins, C., D., “Vortex organization in the outer region of the turbulent boundary layer”, J. Fluid Mech. 422, 154 (2000).
Afgan, I., Benhamadouche, S., Han, X., Sagaut, P. and Laurence, D., “Flow over a flat plate with uniform inlet and incident coherent gusts”, J. Fluid Mech. 720, 457485 (2013).
Alexandrova, O., Saur, J., Lacombe, C., Mangeney, A., Mitchell, J., Schwartz, S. J., and Robert, P., “Universality of SolarWind Turbulent Spectrum from MHD to Electron Scales”, PRL 103, 165003 (2009).
Batchelor, G. K. and Townsend, A. A., “The nature of turbulent motion at large wavenumbers,” Proc. R. Soc. London, Ser.A 199, 238–255 (1949).
Boddupalli, N., Yadav, N. K.and Chandra, L., “Experimental Setup for Investigation of AirFlow and Dust Deposition in Heliostat Field: Design and Evaluation“,Proc. Solar World Congress 2015. (Accepted)
Breuer, M., “Numerical and modeling influences on large eddy simulations for the flow past a circular cylinder”, Int. J. Heat and Fluid Flow 19, 512521, (1998).
Breuer, M., “A challenge test case for large eddy simulation: high Reynolds number circular cylinder flow,” Int. J. Heat and Fluid Flow 21, 648654 (2000).
Breuer, M. and Jovicic, N., “Separated flow around a flat plate at high incidence: an LES investigation”, J. Turbulence 2, 393398 (2001).
Breuer, M., Jovicic, N. and Mazaev, K., “Comparison of DES, RANS and LES for the separated flow around a flat plate at high inclination”, Int. J. Num. Meth. Fluids 41(4), 357388 (2003).
Calaf, M., Hultmark, M., Oldroyd, H. J., Simeonov, V. and Parlange, M. B., “Coherent structures and the k 1 spectral behavior”, Physics of Fluids 25, 125107 (2013).
Celik, I. B., Cehreli, Z. N., and Yavuz, I., “Index of resolution quality for large eddy simulations”, J. Fluid Engineering 127, 949958 (2005).
Chandra, L., and Grötzbach, G., “Analysis and Modeling of the turbulent diffusion of turbulent kinetic energy in natural convection,” Flow, Turbulence and Combustion 79(2), 133154 (2007).
Chandra, L. and Grötzbach, G., “Analysis and modeling of the turbulent diffusion of turbulent heat fluxes in natural convection,” Special TSFP5 issue of Int. J. Heat and Fluid Flow 29 (3), 743751 (2008).
Chakraborty, P., Balachandar, S. and Adrian, R. J., “On the relationships between local vortex identification schemes”, J. Fluid Mech. 535, 189–214 (2005).
Chein, R. and Chung, J. N., “Discretevortex simulation of flow over inclined and normal plates”,Computers & Fluids 16(4), 405427 (1988).
Chen, J. M. and Fang, Y.C.,“Strouhal numbers of inclined flat plates”, J . Wind Eng. Ind. Aero. 61, 99112 (1996).
Choi, H. and Moin, P., “Gridpoint requirements for large eddy simulation: Chapman’s estimates revisited”, Annual Brief, CTR (2011).
Drabble, M. J., Grant, I., Armstrong, B. J.and Barnes, F. H., “The Aerodynamic admittance of a square plate in a flow with a fully coherent fluctuation”, Phys. Fluids 2(6), 10051013 (1990).
Fage, A. and Johansen, F. C., “On the flow of air behind an inclined flat plate of infinite span”, Proc. R. Soc. LondonSer. A 116, 170197 (1927).
Gaskell, P. H. and Lau, A. K. C., “Curvaturecompensated convective transport: SMART, a new boundednesspreserving transport algorithm”, Int. J. Num. Methods Fluids 8, 617–641 (1988).
Georgiadis, N. J., Rizetta, D. P. and Fureby, C., “LargeEddy Simulation: Current Capabilities, Recommended Practices, and Future Research”, AIAA–2009–948 (2009).
Germano, M., Piomelli, U., Moin, P. and Cabot, W. H., “A dynamic subgridscale eddy viscosity model”, Phys. Fluids A3(7), 17601765 (1991).
Green, M. A., Rowley, C. W. and Haller, G., “Detection of Lagrangian coherent structures in threedimensional turbulence”, J. Fluid Mech. 572, pp. 111–120 (2007).
Grötzbach, G., “Revisiting the resolution requirements for turbulence simulations in nuclear heat transfer”, Nucl. Eng. Des. 241, pp. 4379 – 4390 (2011).
Geurts, B. J., and Froehlich, J., “A Framework for predicting accuracy limitations in large eddy simulation”, Phys. Fluids, 14, 41–44 (2002)
Gupta, J. P., “Moisture and thermal regimes of the desert soils of Rajasthan, India, and their management for higher plant production”, Des Hydrologiques 31 (3), 9 (1986).
Johnson, T. A. and Patel, V. C., “Flow past a sphere up to a Reynolds number of 300”, J. Fluid Mech. 378, 1970 (1999).
Katul, G., Porporato, A. and Nikora, V., “Existence of k1 powerlaw scaling in the equilibrium regions of wallbounded turbulence explained by Heisenberg’s eddy viscosity”, Phys. Rev. E 86, 066311 (2012).
Kim, WW., and Menon, S., “Application of the localized dynamic subgridscale model to turbulent wall bounded flows”,Technical Report AIAA970210, 35th Aerospace Sciences Meeting, RenoNV, (1997).
Kim, SE., “Large Eddy Simulation Using an Unstructured Mesh Based FiniteVolume Solver”, 34th AIAA Fluid Dynamics Conference and Exhibit, AIAA 20042548(2004).
Kiya, M.and Matsumura, M., “Incoherent turbulence structure in the near wake of a normal plate”, J. Fluid Mech. 190, 343–356 (1988).
Knisely, C. W., “Strouhal numbers of rectangular cylinders at incidence: a review and new data”,J. Fluids and Structures4, 371393 (1990).
Kravchenko, A. G., Moin,P. and Moser, D. R., “Zonal Embedded Grids for Numerical Simulations of WallBounded Turbulent Flows”, J. Comp. Phys. 127 (2), 412–423 (1996).
Lakshmanan, N., Gomathinayagam, S., Harikrishna, P., Abraham, A. and Ganapathi, S. C., “Basic wind speed map of India with longterm hourly wind data”, Current Science, 96 (7), 911922 (2009).
Lam, K. M., “PhaseLocked education of vortex shedding in flow past an inclined flat plate”, Physics of Fluids 8(5), 11591168 (1996).
Lam, K. M., and Leung, M. Y. H., “Asymmetric vortex shedding flow past an inclined flat plate at high incidence”, European Journal of Mechanics B/Fluids 24, 33–48 (2005).
Lam, K. M. and Wei, C. T., “Numerical Simulation of Vortex Shedding from an Inclined Flat Plate”, Engineering Applications of Computational Fluid Mechanics 4(4), 569579 (2010).
Launder, B. E.and Spalding, D. B., “Lectures in Mathematical Models of Turbulence”,(Academic Press London, 1972).
Leder, A., “Dynamics of Fluid mixing in Separated flow”, Phys. Fluids 3(7), 1741 1748 (1991).
Leonard, B., P., “The ultimate conservative difference scheme applied to unsteady onedimensional advection”, Comp. Methods Appl. Mech. Eng. 88, 17–74 (1991).
Lilly, D. K., “A proposed modification of the Germano subgridscale closure method”, Phys. Fluids A4, 633635 (1992).
Mimura, Y., “The Flow with Wake past an Oblique Plate”, J. Phys. Soc. Japan 13(9), (1958).
Moin P. and Kim, J., “Numerical investigation of turbulent channel flow”, J. Fluid Mech. 118,341 (1982).
Najjar, F. M. and Vanka, S. P., “Effects of intrinsic three dimensionality on the drag characteristics of a normal flat plate” Phys. Fluids A7(10), 2516 – 2518 (1995).
Najjar, F. M. and Vanka, S. P., “Simulations of the unsteady separated flow past a normal flat plate”, Int. J. Num. Methods in Fluids 21, 525 – 547 (1995).
Najjar, F. M.and Balachandar, S., “Lowfrequency unsteadiness in the wake of a normal flat plate”, J. Fluid Mech. 370, 101–147 (1998).
Narasimhamurthy, V. D. and Andersson, H. I., “Numerical simulation of the turbulent wake behind a normal flat plate”, Int. J. Heat and Fluid Flow 30(6), 1037–1043 (2009).
Niknia, I., Yaghoubi, M. and Hessami, R., “A novel experimental method to find dust deposition effect on the performance of parabolic trough solar collectors”, Int. J. Env. Studies 69, 233252 (2012).
Patankar, S. V. and Spalding, D. B. “A calculation procedure for heat, mass and momentum transfer in threedimensional parabolic flows”, Int. J. of Heat and Mass Transfer 15(10), 17871806 (1972).
Patidar, D., Tiwari, S., Pardeshi, R., Chandra, L. and Shekhar, R., ““Solar Convective Furnace for Metals Processing”, JOM 67, 26962704 (2015).
Perry, A. E. and Steiner, T. R., “Large scale vortex structures in turbulent wakes behind bluff bodies. Part 1: Vortex formation process”, J. Fluid Mech. 174, 233270 (1987).
Perry, A. E. and Marusic, I., “A wallwake model for the turbulence structure of boundary layers. Part 1. Extension of the attached eddy hypothesis,” J. Fluid Mech. 298, 361–388 (1995).
Piomelli, U. and Balaras, E., “WallLayer Models for LargeEddy Simulations”, Annual Review of Fluid Mechanics 34, 349–374 (2002).
Piomelli, U., Cabot, W. H., Moin, P. and Lee, S., “Subgrid scale backscatter in turbulent and transitional flows”, Phys. Fluids A3, 17661771 (1991).
Pope, S. B., “Turbulent Flows”, (Cambridge University Press, 2009)
Rodi, W., “On the simulation of turbulent flow past bluff bodies”, J. Wind Eng. Ind. Aero. 4647, 319 (1993).
Rodi, W., “Comparison of LES and RANS calculations of the flow around bluff bodies”,J. Wind Eng. Ind. Aero. 6971, 5575 (1997).
Roshko, A., “On the Drag and Shedding Frequency of TwoDimensional Bluff Bodies”, NACA TN 3169(1954).
Roshko, A., “A New Hodograph for FreeStreamline Theory”, NACA TN 3168 (1954).
Saha, A. K., “Farwake characteristics of twodimensional flow past a normal flat plate”, Phys. Fluids 19, 128110 (2007).
Sharma, P., Sarma, R., Chandra, L., Shekhar, R. and Ghoshdastidar, P. S., “Solar Tower Based Aluminum Heat Treatment System: Part I. Design and Evaluation of an Open Volumetric Air Receiver”, Solar Energy111, 135150 (2015a).
Sharma, P., Sarma, R., Chandra, L., Shekhar, R. and Ghoshdastidar, P. S., “On the design and evaluation of open volumetric air receiver for process heat applications”, Solar Energy 121, 4155 (2015b).
Sharma, R. P., Yadav, N., Kumari, A., “Coherent structures and turbulent spectrum in solar wind plasmas”, Phys. Plasmas 20, 082308 (2013).
Singh, G., Saini, D., Yadav, N., Sharma, P., Chandra and L., Shekhar, R., “Dust deposition mechanism and cleaning strategy for open volumetric air receiver based solar tower subsystems”, Energy Proc. 69, 2081 – 2089 (2015).
Smagorinsky, J., “General circulation experiments with the primitive equations”, Monthly Weather Review, 91(3), (1963).
Tennekes, H. and Lumley, J. L.,“A first course in turbulence”, (MIT Press, Cambridge, 1972).
Terashima, O., Sakai, Y., Goto, Y., Onishi, K., Nagata, K. and Ito, Y., On the turbulent energy transport related to the coherent structures in a planar jet, Exp. Thermal Fluid Science 68, 697710 (2015).
Thompson, M. C., Hourigan, K., Ryan, K.and Sheard, G. J., “Wake transition of twodimensional cylinders and axisymmetric bluff bodies”, J. Fluids Struct. 22, 793–806 (2006).
Wu, T. Y., “A wake model for freestreamline flow theory part I. fully and partially developed wake flows and cavity flows past an oblique flat plate”, Engineering Division, California Institute of Technology, Pasadena, California. Report No. 972 (1961).
Xie, ZT., and Castro, I. P., “LES for streetscale environments and its prospects”, In:Quality and Reliability of LargeEddy Simulations, 271282 (2008).
Yadav, N. K, Pala, D.and Chandra, L., “On the understanding and analyses of dust deposition on heliostat”, Energy Proc. 57, 3004 – 3013 (2014a).
Yadav, N. K., Sharma D. K.and Chandra, L., “Dust deposition on heliostat: understanding and evaluation,” Proc. FMFP 2014, In: FMFP Contemporary Research, Springer (2016).
Yang, D., Pettersen, B., Andersson, H. I., and Narasimhamurthy, V. D., “Vortex shedding in flow past an inclined flat plate at high incidence”, Phys. Fluids 24, 084103 (2012).
Zhang, J., Liu,NanSheng and Lu, XiYun, “Route to a chaotic state in fluid flow past an inclined flat plate”,Physical Review E 79, 045306 (2009).
Zhou, J., Adrian, R. J., Balachandar, S. and Kendall, T. M. Mechanisms for generating coherent packets of hairpin vortices in channel flow. J. Fluid Mech. 387, 353–396 (1999).