Investigating the Thermal Behavior and Productivity of a
Double Slope Solar Still System Integrated with a Basin of
Phase Change Material Supported by an Experimental Data
AuthorsHazim Jassim Jaber, Qais A. Reshack, Qahtan A. Abed, Saif W.
Mohammed Ali
JournalAIP Conference Proceedings
AbstractA numerical model was formulated in the present study to investigate the performance of a double slope
solar still system integrated with a PCM (Phase Change Material) basin. The thermal behavior and productivity of the
solar still was investigated under different solar irradiation fluxes corresponding to different operating time hours.
OpenFOAM V.8 code was used to solve the model using finite volume method for solution and scheme discretization. In
addition, the solidity phase fraction gradient was studied under similar conditions. The numerical results were supported
by an experimental data from literature conducted to the same project of the current work to validate the formulated
numerical model. The numerical test of the double slope solar still was investigated under two different solar fluxes
applied on both the south and north sides of the glass cover. Results show that temperature gradient on the south part
side of the solar still is higher than the gradient values in the north part side. Additionally, it was found that temperature
increases toward water surface because of the heat stored in the water and PCM. Regarding the solidity phase fraction of
the PCM, it was found that the PCM is fused partially as the melting temperature can’t be reached. However, because
the upper layers in the PCM region are closer to the heat source comes from the solar still, solidity fraction of these
layers is higher than the fraction values in the layers below. Indeed, the phase fraction gradient is dependent on the solar
flux, higher solar flux, higher phase fraction gradient. Moreover, interesting finding is that the productivity increases
with the operating time hours despite the reduction in the solar irradiation where the maximum productivity reached is
0.450 (l/m2.hr) at 4:00 PM and solar irradiation equals to 466 (W/m2) while at 12:00 PM where the irradiation reaches its
maximum value (1208 W/m2), the productivity is lower (0.19 l/m2.hr). The numerical behavior of the accumulated and
instantaneous productivity is agreeable and acceptable in compare with the experimental behavior where the average
error of all instantaneous productivity values does not pass 20%.
Investigating the thermal performance of air-cooling system
integrated with TEC modules: Transient numerical simulation
AuthorsTaif M. Mansoor, Saif W. Mohammed Ali, Hussam H. Jabbar
JournalElsevier
AbstractThe current study represents a developed 3-D numerical transient model to a complete indirect
non-evaporative cooling system integrated with an assisted TEC unit using OpenFOAM V.9 source
code. The TEC model was developed using a hand-written code as a source added to the energy
equation. The primary air temperature was a criterion to compare between three cases of environment’s
temperature: low mode (299 K), medium mode (309 K), and high mode (319 K). A new
dimensionless factor called by Heat Flux Ratio (HFR) is derived to evaluate the thermal performance.
In addition, wall heat flux and temperature difference along the vertical air channel were
the parameters used for comparison. The complete system was constructed and validated according
to number of experimental and theoretical studies from literature. Furthermore, all cases
were tested under five different air inlet velocities: 0.75 m/s, 1 m/s, 1.25 m/s, 1.5 m/s, and 1.75
m/s. Results show that HFR values increase as time increases. In addition, it was found that HFR
increases with the reduction in the velocity. Moreover, it was found that higher temperature
difference can be reached at the lowest velocity (0.75 m/s) while the minimum difference would
be at the highest velocity (1.75 m/s) with the error between modes doesn’t pass (2.5%). On the
other hand, the highest wall heat flux approached at the highest velocity (1.75 m/s) and vice
versa where the error doesn’t exceed (2.85%) between the modes. Finally, the validation shows
good and acceptable agreement in COP values of the system with both experimental and theoretical
works form literature where the error is ranged between (20.3% and 3.64%).
Assessment of Thermal Performance of Non-Evaporative Cooling
System Assisted with TEC Models at High Temperature Climate
Authors
Hussam H. Jabbar, Mohammed A. Al-Tayyar, Taif M. Mansoor, Saif W. Mohammed Ali
JournalApplied Mechanics and Materials
AbstractIn this work, a non-evaporative cooling system is used with an assisted thermoelectric cooler (TEC) devices module. The system was proposed as an alternative cooling system in the high temperature climate to overcome the high energy consumption of traditional air-conditioning compression cycle. The open source Open FOAM V.9 was used to solve the transient effect of 3D model of indirect non-evaporative cooling system. The primary air temperature was set to 319 . While, the air flow was tested under four different air inlet velocities: 0.75 m/s, 1 m/s, 1.25 m/s, and 1.5 m/s. the validation shows good and acceptable agreement in COP values of the system with both experimental and theoretical works from literature within an error between (12.9 % and 9.5 %). Results show that the temperate difference value on a slice through the length of the air channel starts to decrease as velocity increasing. For example, at the last timesteps of each velocity, the temperature difference reaches about (~10 oK) when velocity is (0.75 m/s) starting from the first quarter of the channel, while the same difference in temperature not reached until the half way of the channel from the channel inlet when velocity is (1.5 m/s). Revealing that even though the percentage increase in the velocity is about 50%, the change in the temperature difference value between the inlet and outlet of the channel is about 1.2%. The local Nusselt number shows that steady state heat transfer reached very quickly as the velocity increased (i.e., at 0.75 m/s at 12s while for 1.5 m/s at 4s). Notwithstanding, as the time processed the ( increases for all cases but becomes lower as the velocity increased. .
Performance of Kaplan Turbine Operating at Design Condition
AuthorsMuhannad Altimemy, Saif Watheq, Justin Caspar, Alparslan Oztekin
JournalASME Collection
AbstractDesign and optimization using computational fluid dynamics to enhance the hydro turbine’s performance are becoming gradually more common because of its flexibility, minor detailed flow description, and cost-effectiveness. These features are not easily achievable in model testing. k–ω simulations conducted in OpenFOAM 7 characterize the flow structure inside an industrial-sized Kaplan turbine module operating at the peak design flowrate. The power signal, velocity, vorticity, and pressure field are presented over the blades and throughout the draft tube. Additionally, pressure fluctuations were probed along the draft tube wall. The simulation shows a tip vortex rope in the narrow gap between the blade tip and turbine casing. The strong influence of the swirl leaving the runner had a negative impact on the flow pressure fluctuation. Also, high vortical activity was presented near the draft tube wall, leading to turbine instability. It was demonstrated that the turbine generates 14.923 MW of average power. The power signal showed minor fluctuations induced by the vortical activity close to the runner region and the corresponding pressure fluctuations. The Fast Fourier Transform showed the system is dominated by low frequency, high amplitude fluctuations.
Flow Characterization of an Industrial Size Francis Turbine Operating at Ultra-Low Load – The Effect of Water Injection
AuthorsMuhannad Altimemy, Justin Caspar, Saif Watheq, Alparslan Oztekin
JournalASME Collection
AbstractLarge Eddy Simulations (LES) are carried out for a Francis turbine operating at an ultra-low load with and without injection. The flow rate of the turbine is 40% of the design value. The injection aims to improve turbine operation for the already unstable base case away from the design flow rate. Tangential water injection was introduced through the draft tube wall in the same and opposite runner rotation direction. The injection angle was varied (15°, 30°, 45° and 60°). Two water injection rates were applied at 4% and 8 % of the optimal design flowrate. While injection with the 4% rate and 30° in the opposite runner rotation direction helped reduce pressure fluctuations downstream of the injection inlets; no injection configuration could completely mitigate the power and pressure fluctuations. The injection was found to increase the amplitude of pressure fluctuations close to the injection inlets by 2 to 20 times the magnitude of fluctuations without injection. There was a slight reduction in mean power production (4–10% loss) by injection. The high amplitude fluctuations were observed in power signals with and without the injection.
Large Eddy Simulations of Francis Turbine Operating at Ultra-Low Loads
AuthorsMuhannad Altimemy, Justin Caspar , Saif Watheq , Alparslan Oztekin
JournalASME
AbstractHigh-fidelity Large Eddy Simulations (LES) were conducted to characterize the spatial and temporal structure of turbulent flows in an industrial-sized Francis turbine running at 120% of the design flow rate. Injection at a 4% and 8% flow rate is applied and investigated as a mitigation method for pressure-induced fluctuations along the draft tube. Contours of velocity and vorticity in the draft tube are presented to examine the effects of water injection. Probes placed alongside the draft tube measure the pressure signal and compare both operational regimes to characterize the pressure fluctuations. The intensity of pressure fluctuations along the draft tube wall is an order of magnitude smaller compared to that at the center. As the injection is applied, the intensity of the pressure fluctuations along the draft tube wall is increased while the intensity of pressure fluctuations in the center of the draft tube is reduced. Pressure probes in the center of the draft tube measure an 86% to 57% reduction in amplitude for 4% to 8% flow rate injection, respectively. There is a 30% to 40% increase in fluctuations along the wall for 4% to 8% flow rate injection, respectively. These changes in flow structure are due to the dissipation of the vortex rope as the injection is applied.