November 14, 2019 at 5:56 amBartSubscriber
I am conducting a simulation of cavitation in a simple geometry. However, i didn't use the cavitation models which had been implemented in Fluent. I am trying to implement a single-phase model (a barotropic models, pressure is only a function of density) into pressure based solver in Fluent by using UDFs. However, when the simulation was run, i got the error information about "Divergence detected in AMG solver".
The simple geometry is meshed with all-quadrilateral element.
I am using a Realizable k-epsilon turbulence model with scalable wall functions. Pressure-inlet of 30MPa and pressure-outlet of 1Mpa are adopt.
Validation of UDFs which is a piecewise function by an excel spreadsheet.
The error is shown as follows:
When i checked the contour of density or pressure, there was a black cavity in the model.
November 14, 2019 at 10:05 amRobAnsys Employee
Turn the UDF off, does the model run? If so, focus on what the UDF is doing in the highlighted region.
November 14, 2019 at 11:07 amBartSubscriber
Thanks for your reply. If turn the UDFs off, the simulation is changed into a simple single-phase fluid. The model can run well, though the negative pressure appears due to the absence of cavitation model.
The content of UDF is a piecewise function (density = function (pressure)). And the funciton is continuous in terms of the excel spreadsheet. Therefore i don't understand what you mean.
November 14, 2019 at 11:33 amRobAnsys Employee
So, what is the density in the negative pressure region if you plot on the graph?
November 14, 2019 at 11:55 amBartSubscriber
Actually, the occurance of negative pressure is wrong because the operation pressure is set to 0. Generally, in this case, the cavitation model can prevent pressure from falling below saturation vapor pressure.
In the case of UDF using, a pressure value corresponds to a density value. In this simulation, the error appears as before the pressure falls below the saturation vapor pressure. For example, as the min pressure is 100000 Pa, the calculation diverges and the corresponding density is changed into a black cavity as the graph shown. In this case, the negative pressure is forbidden. And sometimes AMG devergence occurs as the presence of negative pressure.
November 14, 2019 at 1:38 pmRobAnsys Employee
Not exactly. Typically the solver sorts itself out within the time step as the pressure, density, volume etc are resolved. By adding in a UDF you may be hindering this process. Try running with a much lower dP and slowly increase this to test the UDF and models.
November 15, 2019 at 7:36 amBartSubscriber
Thanks for your suggestion.
I tried to fix the pressure-inlet as 50MPa and decrease gradually the pressure-outlet from 40MPa. As the pressure-out is 22MPa, the model can run well. While as the pressure-out is set to 21MPa, the devergence occurs.
Strangely, the Min density is 602 kg/m3, and the corresponding min pressure is 662482Pa. I can't understand this result because the density should be 998kg/m3 under the pressure of 662482Pa. A decrement of 1 MPa induces such big variation.
When i change the pressure-out from 25MPa to 20MPa, there is another error. The min pressure is 754074Pa and the density is 999kg/s. There is also a black cavity in the graph of density contour.
I'm not sure if I need to make any changes in the multigrid panel.
Moreover, i read the HELP document. Updating properties (user-defined properties in my UDFs) is behind the equation solving. In this UDFs, the value of density is obtained from the value of pressure. Therefore, is there something wrong with the pressure solving or the AMG setup?
November 15, 2019 at 11:45 amRobAnsys Employee
Turn off node values & replot. You're looking at a smoothed result, but the solver will be using the cell values.
November 15, 2019 at 12:13 pmBartSubscriber
I don't think it's a replot problem. The calculation diverges, then the graph obtained is wrong.
November 15, 2019 at 2:38 pmRobAnsys Employee
Turn off node values & replot. You're looking at a smoothed result, but the solver will be using the cell values. This may well give you a different range.
November 16, 2019 at 1:42 pm
November 28, 2019 at 3:04 amBartSubscriber
I have tried many times. The simulation diverges as the first grid adjacent to the wall at the corner is below 0 Pa. Both a small URF and time step cannot solve this problem. I think the transition of pressure is too fast, but the grids at the corner are dense enough and y+ is ~0.5
Actually, i got a successful calculation of cavitation phenomenon at a small pressure difference in a bigger channel. However, the same setup failed under the condition of a larger pressure difference in a tiny channel.
I will be very appreciate if you can give further advices.
November 28, 2019 at 5:50 amDrAmineAnsys EmployeeWhat about starting by smaller pressure differences than ramping it.
November 28, 2019 at 2:37 pmBartSubscriber
Thank you for your reply.
I tried and it didn't work. Small pressure differences can make the calculation smooth. But, an increment of 0.1MPa can lead to the divergence as the injection pressure is 7MPa ans the outlet pressure is 3MPa. There is also a grid with negative pressure in the contour of pressure. Is there a special setup in multigrid panel that can solve my problem? Because i think the problem can be attributed to the processing mesh.
November 28, 2019 at 4:33 pmDrAmineAnsys EmployeeCan you please add all information regarding the case here?
November 30, 2019 at 3:54 amBartSubscriber
The model is 2d step nozzle meshed by ICEM.
The turbulent model is shown as follow. The correction is applied during the calculation turbulent viscosity, where density is replaced with a function.
The cavitation mode is referred in a paper titled "Performance of turbulence and cavitation models in prediction of incipient and developed cavitation". One needs an appropriate equation of state (EOS) that corresponds to the phase change of the liquid to the liquid–vapour mixture.
The solution method is shown.
November 30, 2019 at 11:05 amDrAmineAnsys EmployeeAs you are using a barotropic model instead of the cavitation models please use udrgm for eos.
November 30, 2019 at 11:06 amDrAmineAnsys EmployeeMoreover switch on pseudo transient and start conservative probably you will need to be transient. Again for stability reasons I recommend using udrgm real gas.
November 30, 2019 at 11:55 amBartSubscriber
Thank you for your reply.
Firstly, what do the words "start conservative" mean?
Pseido transient will make the calculation diverge more fast. I tried to make the case transient, but it didn't work.
Moreover, UDRGM makes the energy model on and needs more equations, such as enthalpy, entropy, thermal_conductivity, which aren't provided in the published paper. And there is no example of UDRGM in users' guide that involved the phase transition.
I hope to use the model without energy equation. I think the large pressure difference at the corner makes the pressure term diverge. Is there any option to correct pressure term? Because I have a successful case in the model of large channel under small pressure difference based on the similar Reynolds number, I don't want to give up this.
November 30, 2019 at 1:11 pmDrAmineAnsys EmployeeJust switch off the energy equation keeping only the density dependent on pressure. I cannot comment on your density UDF as I do not know the paper (what is the ref). Try with unsteady solver and start with small time steps. Now I wish happy weekend
November 30, 2019 at 2:49 pmBartSubscriber
Thank you very much! I will try to use UDRGM.
Finally, have a nice weekend!
- 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.
- Suppress Fluent to open with GUI while performing in journal file
- Floating point exception in Fluent
- What are the differences between CFX and Fluent?
- Heat transfer coefficient
- Getting graph and tabular data from result in workbench mechanical
- Difference between K-epsilon and K-omega Turbulence Model
- The solver failed with a non-zero exit code of : 2
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error: Received signal SIGSEGV
© 2023 Copyright ANSYS, Inc. All rights reserved.