
[Sponsors] 
December 14, 2010, 06:07 
Mach number Vs Convergence trend

#1 
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
I was trying out a ANSYS CFX workshop  Transonic flow over a NACA 0012 Airfoil and the Mach number used in the workshop was 0.7. Everything went well when the same Mach number was used but while running the simulation for higher mach numbers it was found that the solver is struggling more to reach the convergence target (= 1E6 R.M.S). The Reynolds number is 9e6.The SST turbulence model is used. Dynamic viscosity is 1.82e5. The temperature for the simulation is 288 K. The image attached shows the calculations for the pressure ( which is applied at the outlet boundary) and the velocity ( which is applied at the inlet boundary). Can some one suggest the ways by which I can have a more smoother convergence trend at Mach numbers higher than 0.7 in this problem.
__________________
Best regards, Santhosh. 

December 14, 2010, 06:49 

#2 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
Convergence is often harder as the Ma number increases, this is not unusual. Sometimes you have to significantly change the numerical approach when the shock waves start appearing so a setup which ran fine subsonic does not work with a small increase in speed.


December 14, 2010, 08:41 

#3 
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Can you recommend some changes in a focussed way so that I can try them out in CFX for acheiving convergence at the supersonic speeds. Infact I've tried these settings for Mach numbers 1.2 and 1.6 but then the solver terminated abruptly and I recieved some Notices in the third iteration some thing like:
For Mach number = 1.6: Notice: The maximum Mach number is 2.431E+02. For Mach number = 1.2: Notice: The maximum Mach number is 1.208E+02. So after giving out such notices at the end of the third iteration the solver abruptly terminated. Can you please suggest something out of this.It would be better even if you provide some generic steps to be followed in making the settings to achieve proper convergence for the supersonic flows.
__________________
Best regards, Santhosh. 

December 14, 2010, 17:33 

#4 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
I find local timescale factor very useful in getting these sort of flows to start converging. Once you are on the road to convergence then swap back to physical time scales.
More details here: http://www.cfdonline.com/Wiki/Ansys...gence_criteria 

December 15, 2010, 03:46 

#5 
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Thanks for the suggestion Glenn, it works. I could indeed obtain the convergence with the local time scale control itself but can you please tell me why is it not recommended to go with the local time scale alone all the way to convergence (or) why is it stressed upon to run the final few iterations necessarily with a physical time scale instead.
And then after observing a convergence trend with the local time scale control, how do we shift the solver settings to the physical time scale control. Is it done by editing the run in progress? I don’t really have a clarity upon this as this is the first time that I’m working with such options. Can you please guide me upon this.
__________________
Best regards, Santhosh. 

December 15, 2010, 05:57 

#6 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
Do a search of the forum, this issue has been discussed a few times. To be honest, I can't remember the reason, I just remember you should not run local timescale factor all the way to convergence.
Yes, you can do a edit run in progress to switch over. No need to stop and restart. 

December 18, 2010, 14:16 
Hypersonic flow  reg...

#7 
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Thanks Glenn, yep I could find such information about the Local Timescale control on these Links of the forum:
Convergence on Local Time Scale Acceptable Local Timescale Factor So now everything is fine till a Mach number = 2.5 but there was a problem in convergence when I was trying out the things at Mach number = 3. My objective is indeed to acheive convergence at Hypersonic speeds and as such I was proceeding in a step by step manner trying out different Machnumbers. I couldn't acheive convergence at Mach 3. Inspite of using a Local Timescale Factor = 10 and even by enabling the advanced options in the solver settings like: Global Dynamic Model control Velocity Pressure coupling>RhieChow option>High resolution Compressibility control> High speed Numerics The problem is especially seen with the convergence of the continuity(PMass) equation. I'm using the Total energy Heat transfer model and the SST turbulence model. Can you please help me out in resolving this issue.
__________________
Best regards, Santhosh. Last edited by saisanthoshm88; December 18, 2010 at 21:37. 

December 19, 2010, 17:05 

#8 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
Use a lower time scale factor as the Mach number increases.


December 22, 2010, 01:33 

#9 
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Hello Glenn,
I tried out reducing the Local Timescale Factor for the Mach 3 speed but it didn’t work. The PMass equation seems to be highly non linear. In the ANSYS CFX documentation on compressible flows it was mentioned that such non linearity in the PMass equation is natural at Mach numbers greater than 2. It was also mentioned that the CFXSolver will take action to stabilize the solution with an inner iteration that relinearizes the continuity equation. But in my case no such action is taken by the solver and in spite of that I was able to achieve a convergence at Mach 2.5 speed while I had a problem at Mach 3.
__________________
Best regards, Santhosh. 

December 22, 2010, 07:15 

#10  
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
Quote:
I think you can manually set additional PMass equation solves with an expert parameter. 

December 28, 2010, 12:26 

#11  
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Quote:
However, I found an option called tef numerics the description of which was quoted above from the documentation but my problem wasn't solved even after changing this option value to 1. It is indeed mentioned in the documentation for compressible flows that the solver takes a step by itself to invoke such internal PMass equation when the local mach number exceeds 2 but that isn't happening for my case.
__________________
Best regards, Santhosh. 

December 28, 2010, 16:39 

#12 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
No, that' snot the option I was talking about. I do not have access to the software right now so cannot give you the exact name, it is something like high speed numerics, extra pmass conservation loop or something like that. There should be a comment in the output file suggesting you activate it for the high Ma number runs.


December 29, 2010, 02:49 

#13  
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Thanks Glenn I found that option.
Quote:
And A Happy New Year In Advance!
__________________
Best regards, Santhosh. 

December 29, 2010, 06:46 

#14 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
Excellent, yes that was the parameter I was thinking about. Good work for finding it from my vague description.


December 31, 2010, 02:03 

#15 
Senior Member

Join Date: Oct 2010
Posts: 303
Rep Power: 10 
Hi Glenn
The solution converged well till Mach 3 when the max continuity loops option was used but there seemed to be a problem when I've attempted to work at Mach 6. As I've mentioned it earlier I was indeed working out these things upon a Mesh for NACA Airfoil from a Ansys workshop when I've checked the mesh quality by importing it to ICEM CFD (it was a fluent mesh file indeed) I found it to be 0.16 So to have a convergence in case of a Mach 6 speed do you suggest me to try using any relaxation factors in the solver that help in resolving some issues with a poor mesh quality and also to have more continuity loops. By the way I'm running a steady state simulation
__________________
Best regards, Santhosh. Last edited by saisanthoshm88; January 12, 2011 at 23:59. 

January 13, 2011, 01:59 

#16 
Super Moderator
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 12,717
Rep Power: 99 
Mach 6 is getting very fast. Are you sure you don't have dissociation effects? If yes then CFX cannot model this flow.
I would try to improve mesh quality. Mesh quality becomes exponentially more important the faster you go. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Mesh Refinement  Luiz Eduardo Bittencourt Sampaio (Sampaio)  OpenFOAM Mesh Utilities  41  January 17, 2013 03:43 
Unaligned accesses on IA64  andre  OpenFOAM  5  June 23, 2008 10:37 
Abt: Mach number.. Help! Plz!  jinwon park  Main CFD Forum  0  February 5, 2008 14:57 
The Farfield Pressure Boundary & Low Mach Number  SaifDeen Akanni  FLUENT  1  March 5, 2007 06:05 
Trimmed cell and embedded refinement mesh conversion issues  michele  OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ...  2  July 15, 2005 04:15 