Non-orthogonality and pressure jumps
2 Attachment(s)
Hello everyone,
I am facing a problem with my simulation. I try to simulate a particle within a tubeflow. It is a steady state problem. The velocity profile works perfectly, but all my simulations lead to pressure jumps around the particle and on the surface of the particle. They seem to be caused by the non-orthogonality of the mesh, because they only appear near the particle where the mesh changes. Attachment 22088 Attachment 22089 I tried different things to solve this problem. First I raised nNonOrthogonalCorrectors to 2, 3, 5 and even 15, but the problem still appeared. I also changed in fvschemes sngradschemes and laplacianschemes from “corrected” to “limited 0,333”, but still nothing changed. Last I changed the solver for p (GAMG-->PCG), but there were still pressure jumps. I use simpleFoam as a solver and I attached the fvschemes, fvsoluten, controldict and the result of the checkMesh to this post. Thank you in advance. controlDict Code:
/*--------------------------------*- C++ -*----------------------------------*\ Code:
/*--------------------------------*- C++ -*----------------------------------*\ Code:
/*--------------------------------*- C++ -*----------------------------------*\ Code:
/*---------------------------------------------------------------------------*\ |
Hello,
I solved the case now with icoFoam (with 3 nNonOrthogonalCorrectors) instead of simpleFoam, but the results stayed the same and there where still pressure jumps. Does really noone has any idea what i could try to solve the problem? Greetings, Florian |
Hi Florian,
You are using Gauss linear for the grad(p) and grad(U), which is not so accurate for non-orthogonal faces. Try leastSquares, which should maintain second order accuracy on non-orthogonal faces. Best regards, Philip |
Hey Philip,
thank you for your reply. I tried it with leastSquares, but there were still pressure jumps around the sphere. Do you know any other option to solve this problem? Best regards, Florian |
Quote:
try limited schemes for pressure gradient, like Code:
grad(p) faceMDLimited Gauss linear 1; M |
Hey,
I tried it with Quote:
Best regards, Florian |
Hi Florian,
Do these peaks appear in the internal domain or just near the boundary? i.e. if you clip the domain in ParaView through the internal domain, do you see these jumps at the non-orthogonal faces? Philip |
1 Attachment(s)
Quote:
they also appear in the internal domain Attachment 22451 Florian |
Hi Florian,
It doesn't seem to be a non-orthogonality problem of the mesh. Max. skewness ratio is well within the limits. could u post a pic of your mesh ? Maybe, try coarsening the mesh and see what happens to pressure simulation. Also try putting 10e-08 as residual control for pressure if not less. Check your pressure boundary condition on 'ÇYLINDER' based on case setup. regards, achyutan |
1 Attachment(s)
Quote:
first a pic of the mesh (only cylinder and symmetryplane): Attachment 22496 I ran the simulation with a coarsed mesh, but the results were the same. Around the sphere were pressure jumps. I also ran a simulation with the residual condition e-08, but this did not change anything. The pressure boundary condition for the sphere is "zeroGradient". I think this should be okay and the pressure around the sphere looks okay except for the pressure jumps. It seems like these jumps are at the positions where the mesh changes. Florian |
Hi Florian,
Hmmnn this is clearly an issue with mesh non-orthgonality and/or mesh skewness. Did you try leastSquares for grad(U) and grad(p)? Also, you could try skewCorrected for some of the divergence/laplacian terms e.g. Code:
div((nuEff*dev(T(grad(U))))) Gauss skewCorrected linear; |
Quote:
I tried it a few simulations with your suggestions, but the results stayed the same. Florian |
Hey Florian, i had the same problems while simulating a micro channel flow because there were areas, where circular mesh parts and rectangular parts were connected.
I think the problem occurs because there are points where more than 4 cells are in contact with each other (two dimensional, more than 8 in 3d). I am pretty sure that this is the reason why you get the peaks in the pressure field (i had them in my velocity fields, maybe resulting from pressure jumps). Do you get the same problems when you try to simulate the flow around a non spherical particle, maybe a cube? That could eliminate other possible influences. I could not find a solution when i was working on that case, so i switched to another mesh because of time issues, which then worked fine (unstructured, gmsh). I hope that this might help you to track down the reason for this |
Hello Lutz,
I am facing this problem for my Bachelor thesis and at the moment I don't have the time to try out the mesh with a non-sphericle particle. When there is some time left, I will try it. Florian |
All times are GMT -4. The time now is 08:44. |