hugo CFD


I'm simulating a flat plate heat exchanger, almost counter flow.

In the study of knit independence, I used the richardson method. The ratio between the dimension of the mesh element is 2. Three mesh sizes 0.5mm, 1mm, 2mm were studied. The shape of the elements used is hexagonal.

According to the variable in study, heat flux, the three meshes studied converge to the value extrapolated by the richardson method. However, what is the reason for the p-value, order of convergence is higher than the stimulated ratio? Is it supposed to be smaller due to the associated errors, or am I wrong?