
[Sponsors] 
November 19, 2013, 08:57 
high speed flow, phase change problem

#1 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
I am solving a problem in which the inlet is a high speed cryofluid in vapor state and there is a possibility of phase change. Please help
OUTER LOOP ITERATION = 8000 ( 720) CPU SECONDS = 9.357E+05 (1.895E+05)   Equation  Rate  RMS Res  Max Res  Linear Solution  ++++++  UMomBulk  0.96  2.9E04  4.7E02  1.1E02 OK  VMomBulk  0.83  1.6E04  1.8E02  4.5E02 OK  WMomBulk  0.85  4.7E06  7.8E04  4.2E02 OK  PVol  1.00  2.2E05  2.4E03  6.2 5.3E02 OK ++++++  MassVapor  1.00  9.4E11  1.2E08  5.9 6.4E07 OK ++++++  HEnergyVapor  0.98  1.4E03  1.2E01  1.0E03 OK  HEnergyLiquid  1.01  6.5E05  7.5E03  1.4E07 OK  TEnergy  0.90  1.3E05  1.2E03  7.1 1.0E03 OK ++++++  KTurbKEBulk  1.01  1.3E04  1.7E02  6.8 9.9E04 OK  OTurbFreqBulk  0.75  2.5E04  9.2E02  7.4 2.1E04 OK ++++++ ++  Writing backup file 8000_full.bak   Name : Backup Results 1   Type : Standard   Option : Iteration Interval  ++ Slave: 3 Slave: 3 Fatal bounds error detected Slave: 3  Slave: 3 Variable: Vapor.Density Slave: 3 Locale : Domain Interface 1 Side 2 1 Parallel run: Received message from slave  Slave partition : 3 Slave routine : ErrAction Master location : RCVBUF,MSGTAG=1082 Message label : 001100279 Message follows below  : ++  ERROR #001100279 has occurred in subroutine ErrAction.   Message:   Stopped in routine ENFORCE_BOUNDS            ++ ++  An error has occurred in cfx5solve:     The ANSYS CFX solver exited with return code 1. No results file   has been created.  ++ End of solution stage. 

November 19, 2013, 18:54 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
Assuming your physics is correctly set up then this problem is similar to the linear solver failure error. This FAQ tells you the key things to look at: http://www.cfdonline.com/Wiki/Ansys...do_about_it.3F


November 20, 2013, 15:40 

#3 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Yes, I refined the mesh, tried running the model with double precision, viscous dissipation, reducing time steps, the model does not work. The model when setup to run thermal energy equation works just fine.


November 20, 2013, 17:13 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
You are modelling high speed compressible flow with phase change. You have to expect convergence difficulties with a tough model like that. If you are not familiar with how to get convergence with tricky simulations then I suggest you simplify it and practise on more straight forward simulations  like try to model compressible ideal gas flow with some shock waves. You will learn a lot about obtaining convergence when you get those sort of simulations to converge.


November 21, 2013, 11:58 

#5 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Thanks Glenn. I appreciate your help. I am running a cases as per your suggestion.
Last edited by PYJG; November 21, 2013 at 13:43. 

November 21, 2013, 13:42 

#6 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Update
(a)ideal gas compressible cases has not issues converging, I changed the material properties from ideal gas to real gas and then it started giving the same density error. Mach number calculated is high(~3.5) in the second iteration. (b)I am using Argon, there is Souve Redlich Kwong, peng robinson, and redlich kwong model how to ensure which of them would be the correct one for real gas properties. I have been using ArRK model. (c)Will increasing residual target from 1e6 to 1e4 help? Thanks for your suggestions, they have been very valuable 

November 21, 2013, 17:48 

#7 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
So it sounds like you have a problem with introducing real gas properties. This is very common, it is difficult to get these to converge at the best of times.
Loosening the convergence tolerance is OK for development, but once it is working you need to do a convergence sensitivity check to ensure you are tight enough for accuracy. 

November 21, 2013, 18:18 

#8 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Yes, I used 104, 105 and 106 same issue persists. Any suggestions on getting convergence when using real gas properties?


November 22, 2013, 05:43 

#9  
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
Let me quote post #2:
Quote:


November 26, 2013, 16:58 

#10 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Thank you for your kind help Glenn. It is very very valuable. Yes a better mesh did the job. Also, with real gas I cannot get completely converged properties, but I can monitor some important variables and see if they have reached equillibrium and the model shows that.
The problem so far has adiabatic boundary condition. Next step is to look at heat flux boundary condition. As soon as I apply a heat flux on the model that uses total energy and real gas properties the model fails with errors. So I decided to follow the earlier route of getting a converged solution with ideal gas properties+total energy model and using that as a starting solution for the real gas model. In doing that, if I monitor the temperature on the wall where I apply heat flux, the temperature does not rise very much even after I keep increasing the heat flux applied. Any suggestions? ( Mesh along wall has good inflation) 

November 26, 2013, 17:51 

#11 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
As your material properties get more complex then convergence gets harder. So that makes a high quality mesh, double precision, good initial conditions and those factors even more important.
I bet if you simplified your geometry down so you could apply a high quality hex mesh on it that you would have no problems with convergence at all. So the challenge is to generate a mesh of adequate quality on the real geometry. 

November 26, 2013, 18:11 

#12 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
My geometry is 2D( one element in thickness direction) axisymmetric 2 degree wedge. The mesh is refined to over 3.5million nodes for approximate rectangluar dimensions of 25mmx0.65mm. Inflation on walls. Y+ for adiabatic case was less than 5.Any suggestions on getting the heat transfer case to work?
The temperature on the wall where I am applying heat flux is not increasing 

November 27, 2013, 05:32 

#13 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
If you have tight restrictions on mesh quality it is going to be tough to do a 2D axisymmetric simulation. It means mesh the elements along the axis are wedges with a 2° angle and that is a terrible mesh quality.
I would try a larger wedge angle if possible. But even better, if you want to do serious 2D simulations I would change to another solver which supports proper 2D simulations (eg Fluent). 

November 30, 2013, 09:59 

#14 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Thank you for your help it is very valuable. The same problem exists even if the geometry is rotated 360 degrees and mesh refined, inflation along wall . I used the high speed compressible wall heat transfer model and following notice/error appears. I don't see the term in expert parameter. Please help.
****** Notice ******   The nondimensional near wall temperature (T+) has been clipped   for calculation of Wall Heat Transfer Coefficient.     Boundary Condition : front wall   T+ clip value = 1.0000E10     If this situation persists and you are using the High Speed Model,   consider enabling Mach number based blending between low speed and   high speed wall functions. You can do so by specifying a Mach   number threshold as follows:     EXPERT PARAMETERS:   highspeed wf mach threshold = 0.1 # default=0.0 (off)   END  Slave: 3 Slave: 3 Fatal bounds error detected Slave: 3  Slave: 3 Variable: Density Slave: 3 Locale : Default Domain Parallel run: Received message from slave  Slave partition : 3 Slave routine : ErrAction Master location : BRCBUF,MSGTAG=1062 Message label : 001100279 Message follows below  : ++  ERROR #001100279 has occurred in subroutine ErrAction.   Message:   Stopped in routine ENFORCE_BOUNDS  

November 30, 2013, 12:32 

#15 
New Member
Join Date: Oct 2013
Posts: 22
Rep Power: 5 
Expert Parameter does not have this option to select, however, using command prompt type the following and click process.
FLOW: Flow Analysis 1 EXPERT PARAMETERS: highspeed wf mach threshold = 0.1 END END I still get the same density error. Last edited by PYJG; November 30, 2013 at 18:14. 

December 1, 2013, 05:29 

#16 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 13,304
Rep Power: 103 
The root cause of this problem is much more likely to be a physical or numerical instability which is leading to temperatures less than absolute zero. This is quite common as you develop simulations using the real gas meterial property models.
This all points to the same thing I have been posting about all along  you need better mesh quality, a better initial condition, smaller timesteps, etc. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Thermal phase change set up problem  Danial Q  Main CFD Forum  1  July 10, 2014 14:41 
Flow in porous media + phase change  styleworker  FLOW3D  1  June 27, 2012 14:51 
Two Phase Flow Problem  miner15kick  CFX  5  October 28, 2010 18:03 
compressible two phase flow in CFX4.4  youngan  CFX  0  July 1, 2003 23:32 
Geometry setup for phase change problem  Fred G. Kang  Main CFD Forum  1  October 14, 1998 11:41 