July 13, 2023 at 1:24 pmLucySubscriber
I'm relatively new in CFD. I'm trying to model the multiphase flow through a thin (5mm) porous internal using a porous media model. The simulation converges nicely, but then suddenly crashes. Looking at the gas volume fraction contour I can see that in some regions the fluid flows freely (the porous media model does not create any resistance), but in another, gas gets stuck under the object and there is no flow through it. This creates a kind of "dead zone" after the object. I think that causes deconvergence, but I don't know how to solve this problem. I couldn't find any similar problems anywhere. I tried different meshing strategies including local sizing in porous zone regions. The time step is small (0.001). The mesh there is very fine (max.1 mm cells) The quality of mesh is good. The model is as simple as possible.
Does anyone has an idea what can cause such issue? I truly appreciate any ideas.
July 13, 2023 at 2:46 pmRobAnsys Employee
Please can you post some images?
July 14, 2023 at 6:57 am
July 14, 2023 at 6:58 amLucySubscriber
July 14, 2023 at 6:59 amLucySubscriber
Thank you Rob for such a quick answer. More or less it looks like this. I have a similar problem with different meshes so I think that it could be more set up. Generally, the gas phase does not cross the porous media model in some regions.
July 14, 2023 at 8:07 amRobAnsys Employee
Do you have any interface zones? Also, can you zoom into the porous region and show the mesh? The checks will tell you the cell quality, but if you've got a huge jump in cell size between zones it will cause problems.
Please can you plot volume fraction?
July 14, 2023 at 8:44 am
July 14, 2023 at 9:28 amRobAnsys Employee
The cells are OK but the mesh isn't. One of the requirements is aspect ratio and another is growth rate. I suspect you've got a change in flow (ie a gradient) in the left-right direction and it's then causing the solver to fail.
Have a look at porous jump (it may not be suitable here) and also some of the meshing videos/tutorials to see what's been done there. As an aside, for multiphase I tend to avoid poly-hexcore and just use a poly mesh to avoid the jump in cell size in the free stream. That's fine for aero models where they can refine wakes ahead of time but not so good for multiphase models where we generally aren't as sure where stuff is going to go.
July 27, 2023 at 11:17 amLucySubscriber
Hi Rob! Indeed polyhedral mesh seems to provide more stable simulations. Thanks for this recommendation. I also decreased the growth rate which solved the problem that I described above. However, now I see that there is another issue with no difference in a gas fraction before and after the porous media model. As porous zone causes resistance I expect to see such a difference. The model converges nicely now and I don't see any discrepancies that would cause such inaccuracy in results.
Do you have any advice regarding this one? Thanks a lot for helping!
July 14, 2023 at 9:37 amLucySubscriber
Thank you for your valuable feedback. I will look into those videos but as far as I know porous jump unfortunately won't be suitable for my case. I'll also change the mesh to poly, decrease the growth rate and let you know how it works. Thanks again for your help!
July 27, 2023 at 11:22 amRobAnsys Employee
What are you seeing now? Why would gas fraction be different?
July 27, 2023 at 11:32 amLucySubscriber
Thanks for the quick reply. Generally, gas fraction during simulation is the highest in the porous zone and higher above the porous zone than below. As the simulation runs, gas fractions after and before the porous zone become equal.
I expect gas fraction to be different based on other research and as the porous media in the model is supposed to cause a physical obstacle hindering the flow of gas.
July 27, 2023 at 12:45 pmRobAnsys Employee
If you set phase specific resistance coefficients and the gas has somewhere to go along/away from the porous media then, yes it will. Otherwise the gas has to go somewhere!
August 14, 2023 at 4:29 pmLucySubscriber
Hi Rob! As you said, I decreased the growth rate and changed the mesh to polyhedral. In the beginning, everything was fine and I thought that it solved the problem. However, then I lengthened the simulation time and the error with gas phases stocking before the porous zone came back. The growth rate is only 1.05 and the aspect ratio is 8,42.
Furthermore, I am getting the upward flow near the walls and downward in the centrum of geometry, while from experimental data it should be opposite. Before the porous zone, the flow is correct.
Do you think that it's still the mesh? If I decrease the growth rate even further then I'm getting mesh over a million. Now, I'm using 250k.
August 15, 2023 at 10:08 amRobAnsys Employee
Please can you post some images, phase, pressure and velocity would be helpful. Depending on volume fraction etc you can get liquids collecting and then dropping down in a porous zone in the experiment and gas then diverts up due to zone pressure effects. If the porous zone is very thin with low resolution in the CFD model you may not see that.
I assume gravity is pointing in the correct direction?
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.
- Floating point exception in Fluent
- What are the differences between CFX and Fluent?
- Heat transfer coefficient
- Difference between K-epsilon and K-omega Turbulence Model
- Getting graph and tabular data from result in workbench mechanical
- The solver failed with a non-zero exit code of : 2
- Suppress Fluent to open with GUI while performing in journal file
- Mesh Interfaces in ANSYS FLUENT
- Time Step Size and Courant Number
- error: Received signal SIGSEGV
© 2023 Copyright ANSYS, Inc. All rights reserved.