[dzaky.naufal]

I’m performing a transient simulation of a 3D solar chimney system with 24 hours time range (starting from midnight until the next midnight) for my thesis with ANSYS Fluent to see the natural convection phenomenon. Right now, I’m only simulating the fluid domain (the solid domain isn’t simulated) and the fluid that I'm simulating is air. I'm using pressure-based solver, PISO coupling scheme, and Boussinesq density function. The computational domain and boundary conditions are shown below.




The result variables that I want to see are temperature, velocity, and volumetric flow rate. I wanted to make sure that the mesh size and time step size are small enough and the number of iterations per time step is big enough to give me valid results, so I did 'convergence tests' for each of them separately. Unfortunately, I failed to reach convergencies. Here are the details:

1) Convergence test of mesh size: 270 mm, 90 mm, 70 mm, 50 mm (time step size = 300 s, number of iterations per time step = 15)
2) Convergence test of time step size: 300 s, 200 s, 100 s (mesh size = 50 mm, number of iterations per time step = 15)
3) Convergence test of the number of iterations per time step: 5, 10, 15, 20, 40 (mesh size = 50 mm, time step size = 300 s)

Here are the results of the volume-averaged temperature and volume-averaged velocity for each of the convergency tests above.








I’m not so sure about the rest but aren’t the temperature plot for both time step and iteration convergence test show that those plots diverge?

Did I make mistakes in the convergence tests? Or is it possible that I went wrong on other parameters, such as density function, coupling scheme, or even in the meshing process? Is it also possible that the cause of the problem is the absence of the solid domain?

I'm relatively new to CFD, so I rely heavily on trial-and-error method but it seemed ineffective and time-consuming, so I'm really looking forward to your help.

If you need more details, please let me know.
Thank you.