CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Multiphase simulation with RPG table

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 20, 2018, 12:05
Default Multiphase simulation with RPG table
  #1
New Member
 
Join Date: Nov 2016
Posts: 4
Rep Power: 9
Lucy089 is on a distinguished road
Hi everyone,


I am trying to simulate with Ansys CFX a very simple quasi 2D case with liquid nitrogen in order to investigate cavitation.
In order to better match the numerical results with the experimental ones I want to use real gas properties. I have Refprop so I used the RPG generator recommended by Ansys to generate a table for nitrogen (https://support.ansys.com/portal/sit...extfmt=default).

For the table generation I use following ranges:
Tmin-Tmax: 63.5K – 200K; Pmin-Pmax: 6.9Pa – 50e5Pa with respectively 150 points.
In Ansys CFX I create two materials using the same RPG table: for liquid nitrogen I specify N2 and for gaseous nitrogen N2VAP.
First, I just want to simulate both phases but without cavitation in order to get an initial solution. I use the homogeneous model Total Energy for the heat transfer specifying the static temperature of 77.64 K at inlet. The boundary conditions are:
Inlet --> velocity 10 m/s and static temperature
Outlet --> static pressure 7e5 Pa

The simulation starts, however I get this warning from the first iteration:
+--------------------------------------------------------------------+
| Table bounds warnings at: END OF TIME STEP |
+--------------------------------------------------------------------+

+--------------------------------------------------------------------+
| |
| Independent variables went out of bounds while computing the |
| variables listed below using table interpolation. In each case |
| the bounds error was handled by clipping or extrapolation. |
| If this situation persists, consider increasing the table range. |
| |
+--------------------------------------------------------------------+
| |
| Location Name : Fluid |
| Mesh location : VERTICES |
| Routine : UPD_LOCALE_PROP1 |
| Partition : 32 |
| Variable Name : N2Liquid.Temperature |
| Ind. Variable : N2Liquid.Static Enthalpy |
| Bound : Upper |
| Max Value : -1.2073E+05 |
| Handled By : Extrapolation |
| |
+--------------------------------------------------------------------+


It seems that the ranges for the table generation have to be increased. Therefore, I have changed it to the possible max and min but the problem persists.
Furthermore, the results that I am getting do not make any sense.
The enthalpy should be around -121 kJ and I get -2.5 kJ. Consequently, the temperature is not correct: it should be 77.64 K and I get over the entire domain 123.6 K, which is by the way the maximal saturation temperature. And of course all other fluid properties are wrong.

It’s four days that I am trying to fix the problem, I have read several threads here in the forum, but I really don’t know how to fix the problem.
I don’t understand if CFX is wrongly reading the RPG table or if the RPG table itself is wrong.

Any idea or suggestions? Please I really need your help!!!
Lucy089 is offline   Reply With Quote

Old   September 20, 2018, 19:37
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,665
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
You should not define your properties table so tightly about your expected result. As CFX converges the solver is likely to use some way-off values of properties before it homes in as it converges. Also, numerical artefacts like overshoots and wiggles means that even a converged solution will probably exceed expected bounds on properties.

So make the property table much wider for temperature and pressure and that will help.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   September 21, 2018, 05:28
Default
  #3
New Member
 
Join Date: Nov 2016
Posts: 4
Rep Power: 9
Lucy089 is on a distinguished road
Thank you very much for the reply!

Ok, now I have increased the bounds to:
Tmin-Tmax: 22K – 5000K
Pmin-Pmax: 6Pa – 2e9 Pa

But it doesn’t help, I get the same problem.
The interesting thing is that in the warning it is written that the max value of the static enthalpy is -121 kJ/kg which is in fact the enthalpy that I should get with 77 K. But:
1) Looking in the table this is not the max value
2) Looking into the results the static enthalpy is about 5.8 kJ/kg

I don’t understand what CFX is doing because the static enthalpy is defined in the Theory guide as: h_stat = u_stat+p_stat/rho_stat. I suppose u_stat is the enthalpy taken from the table, p_stat is calculated but what about rho_stat?

I have attached also the CCL of my simulation. Perhaps there is a mistake that I can’t see…
Attached Files
File Type: txt CCL_CryoN2_rpg.txt (9.2 KB, 60 views)
Lucy089 is offline   Reply With Quote

Old   September 21, 2018, 07:32
Default
  #4
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 17,665
Rep Power: 143
ghorrocks is just really niceghorrocks is just really niceghorrocks is just really niceghorrocks is just really nice
By the way, does this simulation converge properly if you use constant properties? Define your nitrogen with representative constant properties and make sure you can converge with that simple model. If the simple does not converge you have no hope with the RGP model.

Also, I have never seen anybody use real gas properties with a cavitation model before. Getting cavitation models to converge in constant properties is hard enough without adding the extra complexity of real gas properties. You may have to simplify your gas model to get convergence.
__________________
Note: I do not answer CFD questions by PM. CFD questions should be posted on the forum.
ghorrocks is offline   Reply With Quote

Old   September 21, 2018, 08:30
Default
  #5
New Member
 
Join Date: Nov 2016
Posts: 4
Rep Power: 9
Lucy089 is on a distinguished road
Several researchers numerically investigate cavitation using real gas properties. Here some examples of the same case/ geometry that I am trying to simulate:
https://www.sciencedirect.com/scienc...11227517303521
https://www.sciencedirect.com/scienc...11227515000570
http://fluidsengineering.asmedigital...icleid=1430135

The first one even use CFX and the NIST databank, so I know that it is feasible. I have tried to contact the author but I didn’t get an answer yet.

The simulation with constant properties converge very well. However, the numerical results do not match properly the experimental one and I am quite sure that the cause is the missing temperature dependency of the fluid properties. Therefore, I have started to use RPG tables but something is wrong and I can’t figure out what…
I have the impression that the table that I am using is somehow defect. For example, there is a saturation table after TABLE_9, but the saturation table is not listed in the header.
However, since I define in the material tab the Thermodynamic State Liquid or Gas (see pictures below) it shouldn’t make any difference…
Attached Images
File Type: png N2Liquid_1.png (9.8 KB, 58 views)
File Type: png N2Liquid_2.png (11.5 KB, 48 views)
File Type: png N2Vapor_1.png (9.7 KB, 42 views)
File Type: png N2Vapor_2.png (11.9 KB, 38 views)
Lucy089 is offline   Reply With Quote

Old   October 15, 2018, 10:13
Default
  #6
New Member
 
luo dan
Join Date: Sep 2018
Posts: 27
Rep Power: 7
LUO DAN is on a distinguished road
Why you chose "Gas Phase Combustion" in material group rather than "user"??
LUO DAN is offline   Reply With Quote

Old   October 15, 2018, 10:38
Default
  #7
New Member
 
Join Date: Nov 2016
Posts: 4
Rep Power: 9
Lucy089 is on a distinguished road
@ Luo Dan: Which benefit would bring me switching to user?

By the way, I have spoken with Ansys and fixed the probelm!
Since my fluid is cryogenic I have to "tell" to CFX how to deal with the RPG table. First, you have to enable the beta features Edit->Options->General->Beta Options. The selection of the RPG Liquid Properties is then visible in the GUI and you have to switch from "saturated" to "subcooled".
Now the RPG table is working! However, when using the RPG it is better to do a sensitivity study of the parameters Tmin Tmax Pmin and Pmax for the table generation. I have noticed that depending on these parameters the fluid properties such as static enthalpy, static entropy, density, etc vary a lot and I have tried different ranges in order to match the values given in NIST.
Lucy089 is offline   Reply With Quote

Old   October 16, 2018, 02:54
Default
  #8
New Member
 
luo dan
Join Date: Sep 2018
Posts: 27
Rep Power: 7
LUO DAN is on a distinguished road
But liquid Phase Combustion is a material group, while RGP is generated by external procedure, so we choose user. As for RGP liquid properities, we can edit in "material properties".
Attached Images
File Type: png QQ??20181024204216.png (12.6 KB, 46 views)

Last edited by LUO DAN; October 24, 2018 at 09:43.
LUO DAN is offline   Reply With Quote

Old   October 16, 2018, 02:57
Smile
  #9
New Member
 
luo dan
Join Date: Sep 2018
Posts: 27
Rep Power: 7
LUO DAN is on a distinguished road
hi lucy,
Can you share the rgp generator with me by the dropbox link? I donot have an ANSYS account. I will appreciate it. Thank you!

Last edited by LUO DAN; October 24, 2018 at 09:34.
LUO DAN is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with 2-Way FSI Multiphase Simulation Juli Fluent Multiphase 4 March 21, 2017 09:27
2-Way FSI MultiPhase Simulation Juli ANSYS 0 March 17, 2017 05:52
Multiphase heat exchange simulation efer2109 STAR-CCM+ 1 September 15, 2016 14:22
Single phase result file for multiphase simulation Kushagra CFX 2 July 8, 2008 22:14
multiphase simulation... 2D flow through an elbow Tim CFX 10 April 3, 2008 19:13


All times are GMT -4. The time now is 09:06.