i dont understand: grid validation - laminar flow
1 Attachment(s)
in order to validate a grid to simulate the laminar flow in a pipe (2D, diameter=5 cm, lenght=50 cm, inlet velocity=0.1 m/s), I have used five different grids:
Case 1: 40 elements Case 2: 168 elements Case 3: 1000 elements Case 4: 25021 elements Case 5: 156250 elements I expected that for the denser grid (case 5), I would get the most accurate solution. However, by plotting the velocity profile (at a distance x = 50 cm, having verified that the flow is already developed) I saw the result in Figure 1. According to these results, case 3 would be the most appropriate since the parabolic profile. Why the best result corresponds to an intermediate grid? Did not I should get the best result for denser grid? |
What is the Reynolds number in your simulation?
Are you sure that your solutions are converged? What does your mesh look like? I mean in which direction did you refine the mesh? |
3 Attachment(s)
First of all, thank you very much for answering.
The reynolds number is around 340. Respect to convergence, I first used residual values = 0.001, and then 0.0001. Convergence is always achieved. Figures 2 and 3 show respectively mesh parameters and a printing of the mesh. Figure 4 shows the average velocity at x = 50 cm, variable chosen to validate the mesh. According to Figure 4, the mesh is not yet validated (although I think it is sufficiently dense), but according to Figure 1, the Case 3 is the one that coincides with the expected solution. |
Quote:
With a Reynolds number this high, the flow is not fully developed after a length/height ratio of 10. The velocity profile is still influenced by the profile at the inlet, which is a uniform profile I guess. The reason why the solution at intermediate mesh resolutions appears to be the best one is the numerical diffusion caused by the coarse mesh. You can either decrease the Reynolds number, apply a velocity profile at the inlet or use periodic boundaries between inlet and outlet to obtain a solution which meets your expectations. Another issue: if you want to compare double to single precision results, you should run both simulations until the round-off accuracy of the machine is reached, indicated by the residuals leveling out. The way you did it you compare the influence of the residual criterion (10e-3 against 10e-4). |
Thank you for your help! I thought that the Reynolds number was low enough to consider that the flow is not fully developed. So I have no reason to expect a parabolic velocity profile.
|
You're right! thank you very much again.
|
All times are GMT -4. The time now is 04:30. |