|
[Sponsors] |
Problem with rotating testcase "caradonna_tung" |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 10, 2014, 09:32 |
Problem with rotating testcase "caradonna_tung"
|
#1 |
New Member
JinZhiyi
Join Date: Feb 2014
Posts: 14
Rep Power: 12 |
Dear developers and users,
I am using SU2 version 3.0 "eagle", when I am running the rotating testcase of "caradonna_tung", it seems can not be convergence, and when it runs to 22 iterations, some errors occur like follows: CsysSolve::modGramSchmidt:: w[i+1] = NaN libc++abi.dylib: terminating with uncaught exception of type int Abort trap: 6 Can any nice guy, tell me what happened?? Or figure out a solution, I want to compare the numerical results from SU2 with the classic experiments. Thanks in advance Jin |
|
April 10, 2014, 11:08 |
|
#2 |
New Member
emily
Join Date: Mar 2014
Posts: 15
Rep Power: 12 |
If you ran in parallel, try changing the np (e.g. from -p 14 to -p 16) and see if the error disappear.
|
|
April 10, 2014, 13:10 |
|
#3 |
New Member
JinZhiyi
Join Date: Feb 2014
Posts: 14
Rep Power: 12 |
||
April 10, 2014, 13:25 |
|
#4 |
New Member
emily
Join Date: Mar 2014
Posts: 15
Rep Power: 12 |
Sorry, I don't know the exact solution to this error. It seems that the error occurred frequently in both serial and parallel run. So far, I haven't figured out the reason of the error...
|
|
April 10, 2014, 16:26 |
|
#5 |
New Member
Santiago Padron
Join Date: May 2013
Posts: 17
Rep Power: 13 |
Hi Jin,
We are aware of this problem, this is due to the SU2 version 3.0 "eagle" having new options that are not compatible with the test case you are trying to run. We will be releasing an updated version of the code and test cases shortly. In the mean time you can try using an older version of the code that should be compatible with your configuration file. Also, a comparison to experiment of this case in SU2 can be found here http://adl.stanford.edu/papers/AIAA-2012-3018.pdf Santiago |
|
April 16, 2014, 04:51 |
|
#6 |
New Member
Join Date: Oct 2013
Posts: 12
Rep Power: 13 |
Dear all,
I have download and setup the new version of your very good CFD solver SU2: the 3.1 version. Moreover, I have seen that the Caradonna & Tung case has been presented in the 2014 AIAA "Open-source Analysis and Design Technology for Turbulent Flows". I would like to know if the tutorial case (which is located in the new archive test cases of the 3.1 version) worked? In fact, I have tried and the solver seems to diverge. In a first step, I have generated the mesh with the cyclic boundaries using SU2_PBC which correctly generated the periodic mesh. Then, I have launch CFD simulation using SU2_CFD and the calculation seems to diverge. Have I done a mistake? I would really like to make comparison between SU2 and OpenFOAM on this case, my results using OpenFOAM seems to be not too bad. Best regards, Martin |
|
April 16, 2014, 05:09 |
|
#7 | |
New Member
JinZhiyi
Join Date: Feb 2014
Posts: 14
Rep Power: 12 |
Quote:
Can you please tell me how to use SU2_PBC to generate the periodic mesh with cyclic boundaries. I would like to try your method in my machine. Best regard Jin |
||
April 16, 2014, 06:35 |
|
#8 |
New Member
Join Date: Oct 2013
Posts: 12
Rep Power: 13 |
Hi JinZhiyi,
As it doesn't work, I am not sure at 100% that I am not doing mistake. - In a first step, open the rot_caradonna_tung.cfg and set the MESH_FILENAME=mesh_caradonna_tung.su2 - Execute the command: SU2_PBC rot_caradonna_tung.cfg. It will create you the periodic mesh: mesh_caradonna_tung_periodic.su2. - Change the MESH_FILENAME=mesh_caradonna_tung.su2 and execute the CFD flow by: SU2_CFD rot_caradonna_tung.cfg Best regards, Martin |
|
April 16, 2014, 12:53 |
|
#9 | |
New Member
JinZhiyi
Join Date: Feb 2014
Posts: 14
Rep Power: 12 |
Quote:
However, when checking the surface-flow.dat, compared with the AIAA paper the developer provided, it seems the pressure coefficient don't match the data showed on the paper well. Have a try Jin |
||
April 24, 2014, 06:14 |
|
#10 |
Super Moderator
Thomas D. Economon
Join Date: Jan 2013
Location: Stanford, CA
Posts: 271
Rep Power: 14 |
Hi,
For the time being (we will likely be overhauling the periodic boundary condition soon in the code), I would recommend creating a mesh of the full computational domain for the two-bladed Caradonna & Tung case without periodic boundaries. Also, note that the results given in the paper here http://adl.stanford.edu/papers/AIAA-2012-3018.pdf are based on the Euler equations. In this case, due to the inviscid assumption, the numerically computed pressures and shock location do not perfectly agree with experiment (although they are similar to inviscid results from other codes). For a more accurate validation of the C & T rotor with SU2 using viscous flow, please see the case mentioned above that was presented in this recent paper http://su2.stanford.edu/documents/SU...ciTech2014.pdf Hope this helps, Tom |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
VOF method- Fluid in a rotating Cylinder Problem | Krishna Sandeep | FLUENT | 1 | July 3, 2012 17:44 |
Problem with a rotating problem | samolcue | FLUENT | 4 | June 25, 2012 23:27 |
Rotating problem / transient / 2D | Logan86 | CFX | 6 | July 8, 2011 07:57 |
Problem with rotating SOLID domain | Roland R | CFX | 0 | May 8, 2009 04:38 |
Rotating Vessel: Problem Set-Up | Harmeet | FLUENT | 0 | February 4, 2005 06:41 |