The Solar Tower Chimney Geometry Impact on Its Performance
Abstract:
In this study, the effect of the collector–chimney junction geometry on the performance of a Solar Chimney Power Plant (SCPP) is investigated using numerical simulations based on a two-dimensional axisymmetric model with the standard k–ε turbulence model, coupled with the Discrete Ordinates (DO) radiation model. An SCPP with a collector diameter of 3 m and a chimney height of 3 m is considered. Several junction configurations are analyzed by varying the junction inclination angle, arc radius, and fillet radius. The results show that the junction configuration with a collector inclination angle of 20° and a fillet radius of 14 cm provides the best performance among the inclined cases, enhancing the mass flow rate by approximately 3.2% and increasing the power density by about 15% compared to the same inclination angle without a fillet. In addition, the configuration incorporating an 18 cm arc radius combined with a 14 cm fillet radius increases the mass flow rate by around 1.7% and improves the power density by approximately 7.5% compared to the corresponding arc configuration without a fillet Furthermore, a comparative analysis between configurations with and without a fillet radius reveals that the introduction of a fillet significantly improves the overall system performance, yielding increases of up to 8% in mass flow rate and about 5% in power density. These enhancements are attributed to smoother airflow transitions and more efficient inlet acceleration at the collector–chimney junction, leading primarily to improved aerodynamic performance of the solar chimney system.
1. Introduction
Solar energy has emerged as one of the most widespread and reliable forms of renewable energy due to its availability and environmental compatibility. Numerous technologies have been developed to harness this resource and convert it into electricity. Among these technologies, Solar Chimney Power Plant (SCPP) represent a promising solution, particularly in regions characterized by high levels of high solar irradiation throughout the day. These systems enable the efficient utilization of large land areas to generate electricity in a safe, sustainable, and emission-free manner. Since the demonstration of the first solar chimney prototype in Spain by engineers Haaf et al. [1], this concept has attracted increasing interest from researchers community.
Numerous studies have been conducted to assess the performance of the solar chimneys and evaluate their energy generation capacity and operational efficiency. Other investigations have focused on identifying system limitations and challenges related to practical implementation. To enhance the energy efficiency of SCPP, various optimization strategies have been proposed. In particular, the Manzanares prototype has been widely used as a reference to assess the viability of the concept in arid desert climates. Rabehi et al. [2] employed CFD simulations to evaluate SCPP performance across four Algerian regions, incorporating meteorological parameters such as solar radiation and ambient temperature. Their results indicated an average monthly output of 68 to 73 kW, with an hourly peak of 113 kW in June in Tamanrasset. Sellami et al. [3] analyzed the Manzanares prototype applied to Tamanrasset and Constantine, demonstrating that secondary or tertiary collector roofs could achieve annual outputs of 51 to 79 kW, with hourly peaks of 80 kW in June. Collectively, these studies highlight not only the technological feasibility of SCPP in arid regions such as southern Algeria but also their potential for optimization and dynamic adaptation to local energy demand.
Recent research on SCPP demonstrates a growing interest in hybrid configurations, aiming to overcome the inherent limitations of purely solar systems. Indeed, this technology still suffers from low energy efficiencies and intermittent operation due to its direct dependence on solar radiation. Faced with these constraints, several authors, including Al-Kayiem et al. [4] performed a CFD-based aerodynamic and thermodynamic analysis of a hybrid solar chimney integrated with a gas turbine to recover exhaust heat for continuous operation. The results showed a notable rise in power output (49.04 kW to 67 kW) and collector efficiency (34.3% to 78.5%), demonstrating that exhaust gas heat recovery greatly enhances the energy efficiency and stability of the hybrid system. Habibollahzadea et al. [5] proposed combining a SCPP with a waste-to-energy (WTE) unit to improve electricity production at night. In this system, hot air from the WTE unit is injected into the chimney, which increases the power output. Singh et al. [6] investigated the integration of photovoltaic (PV) panels into a hybrid solar chimney power plant (HSCPP). Their results showed that using a converging collector and a diverging chimney enhanced PV efficiency by approximately 7%, with passive PV cooling proving particularly effective over 80% of the collector surface. More recently, Pandey et al. [7] analyzed the performance of a solar chimney integrated with a subsurface heat exchanger (HE) designed to recover residual ground heat and enable nighttime power generation. Their results showed that this concept can extend chimney operation after sunset and improve the stability of electricity supply. Mohammed et al. [8] studied a HSCPP incorporating phase-change materials (PCM) in the collector and modified with a cavity structure to improve heat storage and distribution. Using three-dimensional CFD modeling, the authors showed that adding PCM to the collector increases the air temperature by nearly 20% and the power output by approximately 25% compared to the conventional model. Furthermore, the cavity structure optimizes heat flow and extends system operation after sunset. Azad et al. [9] conducted a multi-objective optimization of a solar chimney system integrated with electricity-generation and water desalination, by combining CFD simulations with artificial neural networks (ANN), they predicted collector temperature and turbine velocity based on geometric parameters. Results showed that the optimized geometry increased system efficiency by over 20% and improved desalination output. Fadaei et al. [10] studied a HSCPP incorporating PCM in the collector to store and release thermal energy. Using ANN, they accurately predicted key performance indicators such as air velocity, outlet temperature, and power output, with correlation coefficients greater than 0.99. The integration of PCM improved thermal stability and enabled extended operation into the night.
Due to the high cost of building and operating hybrid systems, the exploration of numerous configurations is limited. Consequently, researchers have concentrated on geometric optimization to maximize energy performance while keeping costs manageable. To date, most studies have primarily addressed the technical aspects of SCPP. Mandal et al. [11] used ANN to analyze the effect of geometric parameters (chimney height, air inlet diameter, collector area) on the performance of a solar power plant with a chimney. Using CFD data, their model predicted performance with over 98% accuracy. The results indicate that a greater chimney height and a larger collector area improve mass flow rate and power density, while an excessively large inlet diameter reduces efficiency. The study confirms the effectiveness of AI for the geometric optimization of solar power plants with chimneys. Abdul Hussein et al. [12] investigated a solar chimney system incorporating geometric modifications to the nozzle and canopy to optimize performance. Using CFD modeling, they demonstrated that adding a convergent nozzle and an optimized canopy increased the updraft velocity to 17.8 m/s (compared to 16.0 m/s for the conventional SCPP) and the power output to 67 kW (from 49 kW), significantly enhancing the system’s aerodynamic and energy efficiency. Cuce et al. [13] studied the effect of chimney height on the therm-aerodynamic performance of the Manzanares prototype, showing that greater height improves power output up to a threshold where pressure losses and structural stresses reduce additional gains. Subsequently, in their more recent work [14], the authors combined CFD modeling and ANN to analyze the impact of system geometry on the efficiency of SCPP. They demonstrated that a chimney divergence close to 2$^{\circ}$ and a collector angle between 5$^{\circ}$ and 7$^{\circ}$ promote an increase in air velocity, pressure gradient, and power generated, surpassing conventional cylindrical designs.
Ghernaout et al. [15] conducted a parametric CFD study on the airflow structure in a SCPP to evaluate the effect of geometric and thermal parameters on its performance. The authors showed that increasing the chimney radius and the temperature of the absorbing soil increases the upward air velocity by approximately 20% and the power density by nearly 18%. Conversely, certain configurations cause recirculation zones that reduce efficiency. Kasaeian et al. [16] and Ghalamchi et al. [17] conducted experimental and numerical investigations using a scaled prototype at the University of Tehran to optimize the geometric parameters of solar chimneys. Their results demonstrated that reducing the collector inlet height, increasing chimney height, and employing high thermal conductivity materials such as aluminum significantly enhanced thermal performance. Furthermore, Mandal et al. [18] introduced an innovative absorber surface design in which the conventional flat surface was replaced with a wavy triangular profile; their analysis confirmed improvements in both energy output and efficiency. Singh et al. [19] performed a numerical analysis of a divergent SCPP configuration that included an arc and fillet radius at the chimney base, resulting in enhanced airflow uniformity and increased power output. Mandal et al. [20] investigated the effects of chimney divergence and collector inclination on the performance of solar chimney power plants. Their findings indicate that these parameters significantly enhance thermal efficiency and energy output by optimizing airflow and solar radiation absorption. Singh and Kumar [21] introduced a technical innovation by incorporating a bell-mouth inlet into a convergent-divergent solar chimney system. CFD simulations revealed an approximately 270% increase in air velocity and an electricity output of up to 1,738 kW compared to conventional designs.
The performance of solar chimney power plants depends on several key design parameters, including the chimney height, collector geometry, and collector area. According to the existing literature, most previous studies have primarily focused on these global geometric parameters because of their direct influence on airflow behavior, heat transfer, and overall energy efficiency. However, the geometry of the collector-chimney junction has received relatively little attention, even though it represents a critical transition zone where complex flow phenomena occur. This research gap highlights the need for a more comprehensive investigation into how junction configuration affects the overall performance of solar chimney systems.
In this context, the present study focuses on the structural optimization of the collector-chimney junction. Targeted geometric modifications to the collector-chimney connection are proposed (Table 1) with the aim of enhancing flow stability, aerodynamic performance, and overall energy conversion efficiency.
Parameter | Index | Size |
|---|---|---|
Collector with inclination angle | $\theta_{\text{Col}}$ | 10–30$^{\circ}$ |
Collector with arc | $\text{R}_{\text{Ac}}$ | 6–18 cm |
Chimney base with fillet | $\text{R}_{\text{fillet}}$ | 8–14 cm |
2. Physical Model
Figure 1 shows schematic illustration of the operating principle of a SCPP. Ambient air is heated beneath the solar collector (transparent cover), creating a density gradient that induces an upward airflow through the vertical chimney. This airflow is exploited by a wind turbine installed at the chimney base to convert kinetic energy into electricity [22].
In this study, the geometric dimensions of the solar chimney were adopted from the experimental setup of Ghalamchi et al. [17] to ensure accurate validation of the proposed numerical model. A two-dimensional CFD model is developed using ANSYS Design Modeler, as illustrated in Figure 1, which presents the configuration of the conventional solar chimney system together with the applied boundary conditions (Table 2).
The proposed geometric configurations are illustrated, highlighting the modifications at the junction between the collector outlet and the chimney inlet
Figure 2 illustrate variations in the collector inclination angle (10$^{\circ}$–30$^{\circ}$) and the arc-shaped junction with radius ranging from 6 to 18 cm, respectively. Figures 3 present the combined configurations, incorporating a curved chimney base (fillet) with radius varying from 8 to 14 cm.

Collector Roof | Mixed, h = 10 w / m$^2$ | Wall (Semi Transparant) |
|---|---|---|
Ground | Mixed, h = 10 w / m$^2$ | Wall (opaque) |
Chimney | Adiabatic (Heat flux $ = 0$ ) | Wall (opaque) |
Collector inlet | Tamb = 300 k, Pgauge = 0 Pa | Pressure inlet |
Chimney outlet | Pgauge = 0 Pa | Pressure outlet |


Following the assumptions reported by Kaplan [23], the flow within the solar chimney system is modeled as a two-dimensional, axisymmetric problem. The governing equations include the continuity equation, the momentum equations, and the energy equation, as detailed by Mirzamohammad et al. [24].
where, $u$ and $v$ (m/s), are the velocity components in the $x$ and $y$ directions, respectively, and $\Delta T$ ($K$) represents the temperature difference between the air inside the chimney and the ambient air.
In this study, the standard k-$\varepsilon$ turbulence model is adopted due to its proven robustness and computational efficiency in mixed convection flows, as well as its extensive validation in solar chimney configurations [16], [25], [26]. Although the RNG and realizable k-$\varepsilon$ variants may provide improved predictions for strongly swirling or highly an-isotropic flows, the standard model offers an optimal balance between accuracy, numerical stability, and computational cost. This choice is consistent with recent numerical investigations on solar chimneys, including the work of Ayadi et al. [26].
where, $G_k$ represents the generation of turbulent kinetic energy due to gradients of the mean velocity. $G_b$ denotes the generation of turbulent kinetic energy due to buoyancy effects. The model constants have the following values: $C_1=1.44, C_2=1.92, \sigma_k=1.0, \sigma_{\varepsilon}=1.3$.
Regarding irradiation heat transfer, the Discrete Ordinates (DO) model available in ANSYS Fluent is employed. This model is based on the numerical solution of the radiative transfer equation in participating and semi-transparent media, such as the glass cover of the solar collector considered in the present study. The DO approach is particularly suitable due to its ability to accurately capture directional radiation effects and radiation–matter interactions [11].
The Boussinesq approximation was used to calculate the density variation [27].
In order to objectively compare the efficiency of the investigated configurations, the potential aerodynamic power at the system outlet is evaluated without considering the presence of a wind turbine. This assumption ensures a fair assessment of the impact of geometric modifications. Under this framework, the air mass flow rate is determined using the following expression [28].
The chimney efficiency [29] is given by Eq. (10)
In the absence of a wind turbine, the airflow reaches its maximum velocity, and the entire pressure difference across the system is utilized to accelerate the air, thereby converting it into kinetic energy [30]. Under these conditions, the theoretical power output of the system can be reasonably estimated using the updraft power density concept [19].
Inside the chimney, the maximum air velocity [31].
The convective heat transfer coefficient ($h$) represents the rate of heat exchange between the air inside the collector and the transparent cover. In solar chimney studies, $h$ typically ranges from 5 to 12 W/m$^2$.K, based on experimental and numerical results [32].
In the present work, a value of W/m$^2$.K is adopted, consistent with the literature Rabehi et al. [2]. A constant solar irradiation of 850 W/m$^2$ is assumed, following previous studies on solar chimney power plants [3], [19], [21]. This benchmark value enables a direct comparison of results and facilitates the analysis of peak system performance.
The computational domain is divided into three distinct regions: The collector, the junction, and the chimney. A structured mesh conforming to the geometry of each region was generated using ANSYS. The pressure-velocity coupling is achieved through the SIMPLE method, while gradients are computed using the Green-Gauss cell-based approach. Pressure discretization is performed with a second-order scheme, and the momentum and energy equations are solved using a second-order upwind scheme. For the turbulence and the DO equations, a first-order discretization scheme is employed.
To assess the influence of mesh resolution on the numerical results, three different meshes were compared: a coarse mesh with 13,300 elements, a medium mesh with 19,300 elements, and a fine mesh with 26,200 elements. As shown in Figure 4, the differences in chimney exit velocity between the medium and fine meshes was less than 2%. This analysis confirmed that the medium mesh offered the best balance between accuracy and computational cost.

To validate the numerical model, the simulation results were compared with reference experimental data [17]. The dimensions and boundary conditions were kept identical to those of the experimental setup. The prototype geometry has a horizontal solar radiation absorption surface with a radius of 3 m, a solar collector inlet height of 6 cm, a chimney height of 3 m, and a chimney diameter of 10 cm. The standard k-$\varepsilon$ turbulence model with wall-function treatment was employed, and the near-wall mesh resolution was verified. The dimensionless wall distance falls within the range 0.54 < Y + < 11.94 which is consistent with the recommended values for the standard k-$\varepsilon$ model and ensures an appropriate representation of the near-wall boundary-layer behavior.

Figure 5 presents the velocity and temperature profiles, showing that the numerical simulation closely follows the experimental trends. The peak velocity 1.55 m/s occurs near the chimney axis (0.1 m from the inlet) and decreases toward the collector entrance. A minor deviation is observed near the inlet, where the simulated peak velocity is 1.40 m/s compared to the experimental value of 1.55 m/s, corresponding to a relative error of 9.6%. Similarly, the temperature profile exhibits a maximum of 329 K at the chimney center, decreasing radially toward the ambient temperature of 305 K. The simulation successfully reproduces this behavior, with a slight deviation at 0.1 m from the inlet, where the predicted temperature is 326 K compared to the experimental value of 329 K, resulting in a relative error of 0.9%. Overall, despite localized differences, the Standard k-$\varepsilon$ model accurately predicts the dominant flow and thermal characteristics, confirming its suitability for the present study.
The thermal and physical properties of the key components of the SCPP studied here are summarized in Table 3 and evaluated at ambient conditions. All values are obtained from experimental measurements and validated literature sources, ensuring consistency and accuracy in the CFD simulations [30].
Property / Material | Collector | Ground | Chimney | Units |
|---|---|---|---|---|
Density | 2500 | 2160 | 2719 | (kg/m$^3$) |
Thermal conductivity | 1.15 | 1.83 | 202.4 | (W/mK) |
Specific heat | 750 | 710 | 871 | (J/kgK) |
Absorption coefficient | 0.03 | 0.9 | 0 | - |
Transmissivity | 0.9 | Opaque | Opaque | - |
Emissivity | 0.1 | 0.9 | 1 | - |
Refractive index | 1.526 | 1 | 1 | - |
3. Results and Discussion
This section presents in detail how the shape of the chimney base geometry influences the performance and behaviour of airflow and heat transfer within the SCPP. Four distinct geometric configurations are considered in the analysis. The simulation was conducted for various collector inclination angles (10°, 15°, 20°, 25°, and 30°) and arc radius (6 cm, 8 cm, 12 cm, and 18 cm).
First, a comparison was made between two approaches: gradually increasing the collector inclination angle and widening the radius of curvature of the junction between the collector outlet and the chimney inlet, while maintaining the conventional geometry of the chimney base.
Second, the same comparison was repeated, varying the fillet radius of the chimney base using different fillet radii. The results obtained for the various configurations are presented and discussed in the following sections.
In a conventional solar chimney system, the chimney is designed as a vertical cylinder positioned perpendicular to the horizontal solar collector in order to generate an upward airflow, The objective of this study is to compare the effects of the collector inclination angle and the arc radius while maintaining a conventional chimney base geometry as shown in Figure 2.
Figure 6 illustrate the velocity distribution along the chimney height for different collector–chimney junction geometries. For the inclined collector configuration, relatively high inlet velocities are observed, reaching a maximum value of 1.23 m/s at the chimney inlet. This maximum velocity fully justifies the choice of locating the turbine at this position to take advantage of the more active airflow. In comparison, the arc-shaped collector configuration promotes a smoother acceleration of the airflow starting from the collector outlet. This geometry leads to a more uniform velocity distribution along the chimney height. For an arc radius of 18 cm, the inlet velocity reaches a maximum value of approximately 1.21 m/s. Although this peak velocity is slightly lower than that obtained with the inclined configuration, the arc-shaped configuration provides a more stable and evenly distributed airflow, which is favorable for system efficiency.

As shown in Figure 7, the velocity profiles exhibit a generally similar trend across all configurations, characterized by a gradual increase in velocity from the collector inlet. This similarity persists up to the point where the collector geometry begins to vary (at approximately 130 cm). These findings highlight the significant influence of the collector-chimney junction geometry on the velocity distribution within the system.
The inclination of the collector has a significant impact on the aerodynamic performance of the system. For low inclination angles (0$^{\circ}$–10$^{\circ}$), the airflow acceleration is enhanced, indicating a more effective thermal updraft. However, when the inclination angle exceeds 20$^{\circ}$, a noticeable reduction in velocity occurs, suggesting increased flow disturbance and a degradation of aerodynamic performance.
Configurations with smaller arc radius (6–8 cm) promote greater heat concentration and consequently higher airflow velocity. In contrast, increasing the arc height (12–18 cm) results in flow expansion and a reduction in the average velocity, indicating a loss of available kinetic energy.
Excluding the conventional configuration ($\theta$ = 0$^{\circ}$), the collector geometry with an arc radius of 6 cm exhibits the highest average chimney outlet velocity, reaching approximately 0.55 m/s, compared to 0.36 m/s for the configuration with a 10$^{\circ}$ inclined collector. This represents an improvement of about 53%, calculated according to:

Figure 8 illustrates the static pressure contours, indicates that the inclined collector ($\theta$ = 20$^{\circ}$) presents sharp pressure gradients and localized losses at the chimney inlet, whereas the arc radius configuration ($R_{Arc}$ = 18 cm) achieves a smoother pressure distribution with reduced drops. This improved pressure field enhances suction, leading to higher mass flow rates and increased power density. The static pressure distribution along the collector radius for different configurations is presented in Figure 7. Introducing an arc radius at the collector-chimney junction significantly enhances flow acceleration beneath the collector.The configuration with a 6 cm arc radius exhibits a pressure drop of -0.181 Pa, which is approximately 16.8% higher in magnitude than that of the 10$^{\circ}$ inclined collector configuration (-0.155 Pa). However, when compared to the conventional chimney geometry (-1.834 Pa), this value corresponds to a reduction in pressure losses of about 90%. This significant decrease confirms that a smoother collector-chimney junction effectively reduces flow resistance and enhances momentum transport along the collector.
Figure 9 show that, although the optimized inclined collector configuration ($\theta$ = 20$^{\circ}$) slightly improves the temperature distribution near the chimney base, pronounced thermal gradients persist at the collector-chimney junction. In contrast, the arc radius configuration ($R_{Arc}$ = 18 cm) provides a smoother and more uniform thermal field, contributing to reduced local thermal gradients and improved aerodynamic conditions. As illustrated in Figure 7, the differences in outlet temperature among the investigated configurations remain limited. The conventional configuration (0$^{\circ}$) reaches the highest outlet temperature of 314 K, which is only about 0.5% higher than that of the 6 cm arc configuration (312.5 K) and slightly higher than the 10$^{\circ}$ inclined case (313.4 K, $\approx$ 0.3%). These small variations indicate that geometric modifications primarily influence the aerodynamic performance of the system, while their impact on the overall thermal behavior remains marginal.


In this section, the analysis is extended to examine the effect of fillet radius variations at the chimney base, considering values of 8 cm, 12 cm, and 14 cm, as illustrated in Figure 3. To ensure a consistent and systematic evaluation, the fillet radius is kept fixed for each set of simulations, while the effects of the arc radius and collector inclination angle are investigated. Figure 10 present the variations of these parameters for each fixed fillet radius. This stepwise approach enables a clear comparison among configurations and facilitates the identification of the best-performing geometry for the present study.

For all configurations, the air velocity increases rapidly near the collector center, reaches a maximum in the inner region (r $\approx$ 0.15–0.30 m), and then gradually decreases toward the chimney.
By comparing Figure 10 a, b, and c, it is observed that increasing the chimney base fillet radius leads to a systematic increase in the maximum velocity. The velocity profiles become smoother and more stable, indicating a more efficient aerodynamic transition at the collector–chimney junction.
Since Figure 10 c records the highest velocity values, the configuration with a fillet radius of 14 cm is selected as the reference case to evaluate the effects of the collector arc radius and collector inclination angle. It is clearly observed that increasing the collector arc radius results in a progressive rise in the maximum velocity, along with more uniform and radially extended velocity profiles. In particular, the maximum velocity increases from 0.18 m/s for a collector arc radius of 6 cm to 0.57 m/s for an arc radius of 12 cm, corresponding to a percentage increase of 216.7%.
In contrast, variations in the collector inclination angle induce only limited changes in the airflow characteristics. The configuration with a 10$^{\circ}$ inclination exhibits the highest peak velocity, reaching approximately 0.39 m/s.
As the inclination angle increases from 10$^{\circ}$ to 30$^{\circ}$, the peak velocity decreases and tends to stabilize for inclination angles of 20$^{\circ}$, 25$^{\circ}$, and 30$^{\circ}$, with values of about 0.28 m/s. This behavior indicates that higher inclination angles do not significantly enhance flow acceleration and that, beyond moderate values, the influence of the inclination angle on the velocity field becomes marginal when compared to the effect of the collector arc radius.
Overall, the configuration combining a chimney base fillet radius of 14 cm with a collector arc radius of 12 cm achieves the highest airflow velocity, confirming that the collector inclination angle plays a secondary role compared to the fillet and arc radius in governing the aerodynamic performance of the solar chimney power plant. Consistent with the velocity contour distributions shown in Figure 11 a and b, the airflow within the collector exhibits low velocities near the inlet, followed by a gradual acceleration along the collector length and a pronounced acceleration region near the outlet.
The maximum velocity of 1.24 m/s observed for the 12 cm arc-radius configuration highlights the effectiveness of the optimized junction geometry in enhancing airflow characteristics.

Figure 12 a, b, and c illustrate the effects of the collector inclination angle and arc radius on the static pressure distribution for different chimney base configurations with fillet radius of 8 cm, 12 cm, and 14 cm. The pressure drop within the chimney is mainly attributed to flow acceleration and losses caused by abrupt flow transitions [21]. Increasing the fillet radius clearly improves the static pressure distribution, indicating smoother flow paths and reduced local pressure losses at the collector-chimney junction.

For a fillet radius of 8 cm, the pressure profile exhibits a pronounced negative peak near the collector center, with a minimum static pressure of approximately -0.17 Pa for an arc radius of 6 cm, indicating significant local pressure losses at the collector-chimney junction.
Increasing the fillet radius to 12 cm reduces the magnitude of the negative pressure to about -0.12 Pa and results in a smoother pressure gradient along the collector radius. For the largest fillet radius of 14 cm, the pressure distribution becomes more uniform and stable, with a minimum static pressure of approximately -0.19 Pa for an arc radius of 12 cm. Despite this local minimum, the overall pressure profile is smoother, suggesting reduced flow separation and improved aerodynamic behavior at the junction.
A comparison of the minimum pressure values for fillet radius of 8 cm (-0.17 Pa) and 12 cm (-0.12 Pa) reveals an approximate 29% reduction in the maximum pressure drop. This result confirms that increasing the fillet radius at the chimney base significantly improves aerodynamic performance by mitigating local pressure losses at the collector-chimney junction.
The comparison between arc-shaped geometries and inclined collector configurations, presented in Figure 12 c, is conducted by fixing a 14 cm fillet radius at the chimney base, as this configuration provides the smoothest and most uniform pressure profiles, indicating reduced local pressure losses. Arc configurations generate stronger depressions near the chimney inlet, with minimum pressure values of approximately -0.13 Pa and -0.19 Pa for arc radius of 18 cm and 12 cm, respectively, reflecting enhanced flow acceleration.
In contrast, inclined configurations exhibit more moderate pressure variations, with minimum values of about -0.07 Pa and -0.12 Pa for inclination angles of 15$^{\circ}$ and 10$^{\circ}$, confirming that, for a large fillet radius, the use of an arc is more effective than inclination alone in improving the aerodynamic performance of the SCPP. These trends are further confirmed by the pressure fields shown in Figure 11 c and d.
Figure 13 and Figure 14 present the variation of the mass flow rate and the theoretical power density as functions of the fillet radius for different arc radii and collector inclination angles.


Figure 13 illustrates the influence of the arc radius on the mass flow rate for various fillet radius. For all investigated cases, increasing the arc radius from 6 cm to 18 cm leads to a monotonic increase in the mass flow rate. For example, at a fillet radius of 14 cm, the mass flow rate increases from approximately 0.0343 kg/s for an arc radius of 6 cm to about 0.0349 kg/s for an arc radius of 18 cm, corresponding to an improvement of nearly 1.7%. This trend is consistently observed for smaller fillet radius, confirming that larger arc radius promote smoother flow guidance at the collector-chimney junction.
Moreover, the introduction of a fillet radius systematically enhances the mass flow rate compared to the corresponding configurations without a fillet. The optimal configuration among the arc radius cases is obtained for an 18 cm arc radius combined with a 14 cm fillet radius, which achieves a maximum mass flow rate of approximately 0.0349 kg/s, representing an improvement of about 8% relative to the same arc-radius configuration without a fillet. This improvement can be attributed to reduced flow separation and lower aerodynamic losses at the chimney base.
The corresponding effect on the power density is shown in Figure 14. For a given fillet radius, increasing the arc radius from 6 cm to 18 cm results in a systematic increase in power density. At a fillet radius of 14 cm, the power density rises from approximately 1.07 W/m$^2$ to about 1.15 W/m$^2$, corresponding to an improvement of nearly 7.5%. This confirms that the aerodynamic enhancement induced by a larger arc radius directly translates into improved power output.
Figure 13 presents the variation of the mass flow rate with the collector inclination angle for different fillet radius. For a given fillet radius, increasing the inclination angle from 10$^{\circ}$ to 30$^{\circ}$ leads to a noticeable enhancement in the mass flow rate. At a fillet radius of 14 cm, the mass flow rate increases from approximately 0.0321 kg/s at 10$^{\circ}$ to about 0.0326 kg/s at 20$^{\circ}$, corresponding to an improvement of nearly 1.5%.
Consequently, the optimal inclination angle is identified as 20$^{\circ}$, when combined with a 14 cm fillet radius, yielding the highest mass flow rate among all inclination angle cases. Compared to the same inclination angle without a fillet radius (0.0316 kg/s), this configuration exhibits an increase of approximately 3.2%, highlighting the significant role of the fillet in improving flow smoothness and momentum transfer.
The influence of the collector inclination angle on the power density is depicted in Figure 14 b. For all fillet radius, the power density increases as the inclination angle rises from 10$^{\circ}$ to 20$^{\circ}$. At a fillet radius of 14 cm, the power density increases from approximately 1.12 W/m$^2$ at 10$^{\circ}$ to about 1.29W/m$^2$ at 20$^{\circ}$, corresponding to an improvement of nearly 15%. Beyond this angle, the rate of improvement diminishes, confirming that excessive inclination adversely affects the aerodynamic performance.
4. Conclusions
This paper presented a numerical investigation of a SCPP to analyze the influence of the collector–chimney junction geometry on airflow and thermal performance. Using the standard k-$\varepsilon$ turbulence model coupled with the DO radiation model, the study evaluated various geometric configurations under a constant solar irradiation of 850 W/m². Based on the simulation results, the following conclusions are drawn:
$\bullet$ Impact of Junction Geometry: The geometry of the collector chimney interface is critical for system performance. Introducing a curved junction (fillet) significantly suppresses flow recirculation zones and reduces pressure drops compared to the conventional sharp-edged design.
$\bullet$ Increasing the arc radius from 6 cm to 18 cm consistently enhances the mass flow rate and power density, with the most favorable performance achieved for an 18 cm arc combined with a 14 cm fillet radius. increases the mass flow rate by approximately 1.7% and the power density by about 7.5% compared to a smaller arc radius of 6 cm. The improvement is mainly attributed to smoother airflow and reduced aerodynamic losses at the collector-chimney junction.
$\bullet$ The optimal inclination of 20${}^{\circ}$ combined with a 14 cm fillet radius results in an increase in mass flow rate of approximately 3.2% and an improvement in power density of nearly 15% compared to the same angle without fillet. due to smoother airflow transition and reduced friction losses at the junction. Deviations from this optimal angle, either higher or lower, lead to reduced performance due to less favorable flow acceleration and increased aerodynamic losses.
$\bullet$ The comparison between configurations with and without fillet radius reveals that introducing a rounded junction significantly enhances the performance. These results confirm that the curved transition contributes to smoother airflow and improved aerodynamic efficiency at the collector-chimney interface.
$\bullet$ Overall, the study confirms that geometric optimization of the collector-chimney interface is a simple yet effective strategy to enhance both aerodynamic and thermal efficiencies of SCPP, offering valuable insights for the design of future large-scale solar chimney systems.
Limitations and Future Work:
While this study demonstrates the benefits of geometric optimization, certain limitations must be acknowledged. First, the reliance on a 2D axisymmetric model assumes perfect flow symmetry and may not fully capture complex three-dimensional vortex shedding at the chimney inlet. Second, the simulations were conducted under constant solar irradiance, neglecting the transient thermal inertia effects caused by diurnal solar cycles. Finally, the model validation relied on small-scale experimental data. Future research will focus on full 3D transient simulations to account for environmental variations and will seek validation against larger-scale pilot facilities to ensure scalability of these findings.
The data used to support the research findings are available from the corresponding author upon request.
The authors declare no conflicts of interest.
$A_{chimney}$ | Chimney area, $\mathrm{m}^2$ |
$c_p$ | Specific heat of air, kJ/kg K |
$C_1, C_2$ | Constants for the turbulence model |
$D_{coll}$ | Collector diameter, m |
$D_{chimney}$ | Chimney diameter, m |
$g$ | Acceleration due to gravity, $\mathrm{m} / \mathrm{s}^2$ |
$H_{chimney}$ | Chimney height, m |
$I$ | Solar irradiation, W/m${}^2$ |
$\dot{m}$ | Mass flow rate, kg/s |
$P$ | Theoretical power produced, $\mathrm{w} / \mathrm{m}^2$ |
$P_{tot}$ | Power output, W |
$R_{Arc}$ | Arc radius, cm |
$R_{fillet}$ | Fillet radius, cm |
$T_a$ | Ambient temperature, K |
$v$ | Maximum air velocity, m/s |
$x$ | Horizontal coordinate, m |
$y$ | Vertical coordinate, m |
Greek symbols | |
$\beta$ | Thermal expansion coefficient |
$\mu$ | Dynamic viscosity of air, Ns/m$^2$ |
$\eta_{chimney}$ | Chimney efficiency, % |
$\rho$ | Ambient air density, kg/m$^3$ |
${\theta}_{{Coll}}$ | Inclination angle, ${}^{\circ}$ |
