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

Supersonic flow crashing in CFX Solver

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

Like Tree4Likes
  • 1 Post By tomson199
  • 1 Post By tomson199
  • 1 Post By tomson199
  • 1 Post By turbo

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 7, 2016, 21:12
Default Supersonic flow crashing in CFX Solver
  #1
Member
 
Faizan
Join Date: Mar 2014
Posts: 76
Rep Power: 12
Mfaizan is on a distinguished road
Hi all,

I am modeling a supersonic flow through convergent-divergent nozzle with throat dia 2.6mm. I am using a physical timescale of 1e-6[s] and sometimes 1e-7[s]. In both situations the CFX solver crashing by causing an huge increase in mach no value of greater than 8 or 9 after 50 iterations with 1e-6[s] physical timescale and after 200 iterations in case of 1e-7[s] timescale. I plotted the velocity profile of failed solver file and found that the highest no of mach no was found in divergent section of the nozzle. Means in convergent section and upto the throat the solver was stable. I have changed the mesh nodes by making it more dense but still CFX solver is crashing. Please suggest what should I do to make the supersonic flow up and running. I used total pressure and total temperature inlet conditions at nozzle inlet and the values are 550 deg centigrade and 1.4MPa.

Any help and support will be much appreciated.

Thanks in advance

Faizan
Mfaizan is offline   Reply With Quote

Old   March 8, 2016, 04:22
Default
  #2
Member
 
Thomas
Join Date: Dec 2014
Location: Poland
Posts: 49
Rep Power: 11
tomson199 is on a distinguished road
Maybe your divergent section has quite large angle and expand too fast, so it may cause a big separation vortex? If yes, you should simulate this case as a transient? Did you perform hand calculation of this nozzle before?
Mfaizan likes this.
tomson199 is offline   Reply With Quote

Old   March 8, 2016, 04:41
Default
  #3
Member
 
Faizan
Join Date: Mar 2014
Posts: 76
Rep Power: 12
Mfaizan is on a distinguished road
Hi Tomson,

Thanks for your response. The nozzle which I used in the simulation is an exact replica of actual nozzle used in the lab equipment. So I am very much sure that nozzle geometry is absolutely correct. The flow I am modeling is steady state flow.

Please suggest.

Thanks-faizan
Mfaizan is offline   Reply With Quote

Old   March 8, 2016, 04:44
Default
  #4
Member
 
Thomas
Join Date: Dec 2014
Location: Poland
Posts: 49
Rep Power: 11
tomson199 is on a distinguished road
Could you past a few photos and print CEL from solver and show? It will shows me in more details with what you are stuggling.
Mfaizan likes this.
tomson199 is offline   Reply With Quote

Old   March 8, 2016, 04:54
Default
  #5
Member
 
Faizan
Join Date: Mar 2014
Posts: 76
Rep Power: 12
Mfaizan is on a distinguished road
Hi Tomson,

I have just uploaded the pics of mach no, velocity and mesh from a nbackup of crashed solver file.

The message in solver before crashing is given below:

================================================== ====================
OUTER LOOP ITERATION = 233 ( 17) CPU SECONDS = 4.458E+05 (2.782E+04)
----------------------------------------------------------------------
| Equation | Rate | RMS Res | Max Res | Linear Solution |
+----------------------+------+---------+---------+------------------+
| U-Mom | 0.83 | 6.0E-04 | 1.4E-01 | 1.9E-03 OK|
| V-Mom | 0.88 | 6.0E-04 | 8.8E-02 | 1.7E-03 OK|
| W-Mom | 1.04 | 1.2E-03 | 5.1E-01 | 1.6E-03 OK|
| P-Mass | 1.02 | 4.8E-04 | 3.8E-01 | 15.1 5.0E-01 ok|
+----------------------+------+---------+---------+------------------+
| H-Energy | 0.96 | 6.3E-04 | 8.2E-01 | 1.1E-02 OK|
| T-Energy | 0.95 | 2.2E-04 | 7.0E-03 | 6.6 1.1E-02 OK|
+----------------------+------+---------+---------+------------------+
| K-TurbKE | 2.76 | 5.0E-04 | 1.7E-01 | 6.7 1.5E-03 OK|
| E-Diss.K | 1.75 | 4.9E-04 | 1.7E-01 | 12.4 7.9E-05 OK|
+----------------------+------+---------+---------+------------------+
| HELIUM HPGas Materia | 0.00 | 0.0E+00 | 0.0E+00 | 6.7 0.0E+00 OK|
| NITROGEN HPGas Mater | 0.86 | 1.2E-03 | 5.0E-01 | 6.7 1.6E-02 OK|
+----------------------+------+---------+---------+------------------+
----------------------------------------------------------------------
| Equation | Rate | RMS Res | Max Res | Linear Solution |
+----------------------+------+---------+---------+------------------+
| P-Mass |99.99 | 1.0E-01 | 8.2E-02 | 15.1 5.3E-02 OK|
+----------------------+------+---------+---------+------------------+
+--------------------------------------------------------------------+
| Notice: The maximum Mach number is 4.589E+01. |
+--------------------------------------------------------------------+

================================================== ====================
OUTER LOOP ITERATION = 234 ( 18) CPU SECONDS = 4.475E+05 (2.944E+04)
----------------------------------------------------------------------
| Equation | Rate | RMS Res | Max Res | Linear Solution |
+----------------------+------+---------+---------+------------------+

+--------------------------------------------------------------------+
| ERROR #004100018 has occurred in subroutine FINMES. |
| Message: |
| Fatal overflow in linear solver. |
+--------------------------------------------------------------------+

+--------------------------------------------------------------------+
| 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.

+--------------------------------------------------------------------+
| The following transient and backup files written by the ANSYS CFX |
| solver have been saved in the directory |
| M:\Multicomponent_Models\T550P14_MC_PIV_N2_MESH2_B _070316\T550P14- |
| _MC_PIV_N2_MESH2_B_070316_003: |
| |
| 232_full.bak |
+--------------------------------------------------------------------+


+--------------------------------------------------------------------+
| The following user files have been saved in the directory |
| M:\Multicomponent_Models\T550P14_MC_PIV_N2_MESH2_B _070316\T550P14- |
| _MC_PIV_N2_MESH2_B_070316_003: |
| |
| PT1.trk, solver.setup, mon |
+--------------------------------------------------------------------+


This run of the ANSYS CFX Solver has finished.


Please suggest.

thanks
Attached Images
File Type: jpg Mach No.jpg (30.5 KB, 45 views)
File Type: jpg Mesh.jpg (97.6 KB, 42 views)
File Type: jpg Velocity1.jpg (44.9 KB, 33 views)
File Type: jpg Vel2.jpg (34.2 KB, 34 views)
Mfaizan is offline   Reply With Quote

Old   March 8, 2016, 05:09
Default
  #6
Member
 
Thomas
Join Date: Dec 2014
Location: Poland
Posts: 49
Rep Power: 11
tomson199 is on a distinguished road
Ok, I saw that your grid looks good. You've got the problem with "Overflow". I had this problem few times. If it appeared, I change my turbulence model to SST and in the most cases it works.
Mfaizan likes this.
tomson199 is offline   Reply With Quote

Old   March 8, 2016, 16:22
Default
  #7
Senior Member
 
Join Date: Jun 2009
Posts: 174
Rep Power: 16
turbo is on a distinguished road
It looks you have multiple fluids in the domain. Try to run a simple air case at first to see if the crash came from fluid material modelings. Building a nice initialization with multiple fluids would not be that easy.
Mfaizan likes this.
turbo is offline   Reply With Quote

Old   March 8, 2016, 21:28
Default
  #8
Member
 
Faizan
Join Date: Mar 2014
Posts: 76
Rep Power: 12
Mfaizan is on a distinguished road
Hi turbo,

Thanks for your suggestion. Actually I have already done the case with nitrogen only and this is my new project to deal with multicomponent fluid (Nitrogen and Air). The CFX set up is OK. I reckon the mesh I am using is not good for supersonic flow. Can you assist me to overcome this problem.

Thanks

faizan
Mfaizan is offline   Reply With Quote

Old   March 9, 2016, 08:13
Default
  #9
Senior Member
 
Join Date: Jun 2009
Posts: 174
Rep Power: 16
turbo is on a distinguished road
If the nitrogen case was with a different mesh, you will need to run it again to check if the crash was mesh-related or material-related. Whenever you have unknown errors, the best way is to narrow down the cause for debugging.
turbo 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
compressible flow calculation error using rhoSimpleFoam solver student4326 OpenFOAM Running, Solving & CFD 7 November 2, 2015 12:34
how to simulate orifice in supersonic flow? YING FLUENT 2 May 16, 2013 11:41
CFX fails to calculate a diffuser pipe flow shenying0710 CFX 7 March 26, 2013 05:13
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 16:45
Supersonic flow past a wedge with counter flow Mahesh Bailakanavar FLUENT 0 February 14, 2008 01:21


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