May 18, 2021 at 1:52 pmsmarx2Subscriber
I am using R2021 VOF-to-DPM for simulating a nozzle, and am having my simulation error out due to high courant number. I tried smaller time steps, coarser mesh, etc. but found that the problem continues. I decided to plot the courant number and it showed that my courant number was exceeding 1000. So I then plotted my velocity magnitude to see what was happening. In a configuration where I wouldn't expect more than 100 m/s velocity at any given point, the simulation is showing up to 65,000 m/s (image attached is the plane parallel to the vee cut of the nozzle and the iso surface of multiphase interface). Any ideas on what I'm missing or what's happening?May 18, 2021 at 3:33 pmRobAnsys EmployeeThat suggests the model was diverging before that time step. If you go back some steps are you seeing any convergence issues? Also check the mesh in that region, look for cell quality, cell size change and whether you're resolving the flow.
May 18, 2021 at 3:41 pmsmarx2SubscriberRob I am using the Adaptation Controls similar to the VOF-to-DPM Crosswind tutorial, with a minimum cell size of 1e-13 and refinement level of 5. My initial cell size in that region is quite large though, so maybe that has something to do with it? How do I check cell quality for a specific region? (I am a self taught Fluent user, so don't know all the ins and outs). Thanks.
May 18, 2021 at 3:57 pmRobAnsys EmployeeIf it's a hex mesh quality should be fine as the adaption should turn it into smaller hex cells. If you look in the reports or contours there are skew metrics under mesh: check what the ranges are.
May 19, 2021 at 2:57 pmMay 19, 2021 at 3:17 pmRobAnsys EmployeePossibly, turn off node values and alter the range to see what is going on in that area.
May 19, 2021 at 3:29 pmMay 20, 2021 at 7:41 amaitor.amatriainSubscriberCould you please post some mesh pictures? For example, the initial mesh, the mesh just before Co<1 and the mesh just after Co>1. As Rob wrote, it is interesting to detect that transition to see what is happening. Just some data with Co>>1 does not help to detect the problem.
May 21, 2021 at 2:44 pmsmarx2SubscriberAitor and Rob It appears to have to do with my adaptive grid criteria. I attempted to re-mesh the domain and re-run the model. The Courant number again spiked, and did so as the interface reached the larger cells.
For the model, my minimum droplet size is 50 microns, so I calculated that to have two cells across with 8 levels of refinement, my maximum cell size is 6.4mm (.05mm/2=.025mm... .025mm*2^8=6.4mm... hopefully I did that correctly).
Could my issue be a misunderstanding of the adaptive grid criteria?
Here's an image of the mesh during adaptation (Courant#~1):
Here's a zoomed in image of the interface mesh after the Courant number went over 250:
Here's images of the interface when Courant#~1:
I'm happy to post more pictures. Not sure what's helpful... thank you for the quick responses and aid in this issue.
May 21, 2021 at 3:18 pmRobAnsys EmployeeYes, that's the refinement+coarsening to many levels messing with the Poly mesh: it's because it's refined a lot then coarsened repeatedly so the mesh is no longer what you created. We've seen it before but that's about the worst I've seen, even in testing. Easiest fix is to use poly+hexcore or reduce the max refinement level: the former works as the hex mesh doesn't distort.
May 21, 2021 at 6:43 pmsmarx2SubscriberRob I attempted to use the Poly-Hexcore meshing method. I reduced my max cell size to 1.6mm so that I could use 6 levels of refinement. After everything was set up and initialized, my volume fraction plot appears to have gaps in a couple cells (below). When I tried to run calculations, I received a floating point exception (which I assume is due to the missing cells). Do I need to try to re-mesh the domain? Also, please feel free to let me know if I misunderstood any of your recommendations.
May 24, 2021 at 7:22 amaitor.amatriainSubscriberI created a Service Request some time ago related to the white cells. I was told that this bug was also previously reported, and a provisional solution was to use CFD post. The problem is just visual, and does not depend on cell quality. The cells are there; you can try to simulate a simple single-phase flow and you will see that flow variables are smooth in that area.
In regards to your specific problem, this is related to mesh refinement and I have not enough experience to give you a solid advice. However, I would like to mention that I have faced a similar problem in the past related to Poly-Hexcore mesh. In detail, during some transient LES (or hybrid RANS/LES) simulations, there was a (poly) cell at the outlet (or a couple of cells) where Co increased gradually until tending to infinite, and this caused divergence. The cells in the neighbourhood had Co<<1, so this was not supposed to be a problematic area.
I checked everything (mesh quality, domain size, maximum Co number and so on), and after some discussion with ANSYS members, the solution was to patch Co number in that cell, and set v = 0 there. In principle, this worked, even though after some iterations the problem appeared with temperature. I did not made more progress because I changed the approach to my specific problem.
May 24, 2021 at 7:23 amRobAnsys EmployeeThe gaps are just a surface display glitch, the facet isn't in the surface so you see white "holes". The solution is fine, and it's been an occasional issue since iso surfaces and tet mesh first came out." target="blank"> We do have a best practice but it's not in the public area so I can't share the document.
I do suggest that you patch in the whole nozzle with liquid, and watch this
Viewing 12 reply threads
Ansys Innovation Space
- You must be logged in to reply to this topic.
Simulation World 2022
Check out more than 70 different sessions now available on demand. Get inspired as you hear from visionary companies, leading researchers and educators from around the globe on a variety of topics from life-saving improvements in healthcare, to bold new realities of space travel. Take a leap of certainty and check out a session today here.
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.Trending discussions
- Suppress Fluent to open with GUI while performing in journal file
- Heat transfer coefficient
- What are the differences between CFX and Fluent?
- Floating point exception in Fluent
- The solver failed with a non-zero exit code of : 2
- Difference between K-epsilon and K-omega Turbulence Model
- Getting graph and tabular data from result in workbench mechanical
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error in cfd post
Top Rated Tags
© 2022 Copyright ANSYS, Inc. All rights reserved.Ansys does not support the usage of unauthorized Ansys software. Please visit www.ansys.com to obtain an official distribution.