Question on residual convergence
I ran a solution and monitored some physical parameters (mass flow rate at inlet and outlet and linear pressure drop) as well as the residuals. My continuity and z-momentum residuals only drop about an order of magnitude (the z-direction is the dominant flow direction), but they are normalized to the initial values, which were set pretty close to the converged values. So it makes sense to me if they don't drop that far from the beginning, however, they start to oscillate as the iterations continue on and I'm not sure if this should be expected or not. See the attached picture. Is this oscillation acceptable as long as my inlet/outlet mass balance is good and the pressure distribution makes sense?
Okay, so I inspected some additional physical parameters, namely average flow area velocity and it turns out that the solution wasn't converged. The area averaged velocity never settled on a number, but just kept changing through the iteration. I played with many variables for a while before I found something that allowed my solution to converge on realistic physical parameters AND for the residuals to look reasonable. My near wall prism layer was very thin in comparison to the 2nd prism layer and surrounding polyhedral cells, so I doubled its thickness. I also increased the gap fill percentage to 48% to try and get rid of some of the very small cells in the tight crevices of the models. Running the model without a prism layer exhibited ugly residuals, but using either the polyhedral mesher or trim mesher with the thicker near wall prism layer gave me better results. My residuals now look as follows:
|All times are GMT -4. The time now is 23:55.|