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

Rocket Ma=2.089 - void CSolver::SetResidual_RMS(const CGeometry*, const CConfig*)

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 19, 2021, 18:07
Default Rocket Ma=2.089 - void CSolver::SetResidual_RMS(const CGeometry*, const CConfig*)
  #1
Member
 
Alain Islas
Join Date: Nov 2019
Location: Mexico
Posts: 43
Rep Power: 4
alainislas is on a distinguished road
Hello SU2 community,

I am modeling a rocket (fuselage only) at supersonic flow (Ma = 2.089) conditions.

- Models:

Turbulence: kw-SST
Viscosity: SUTHERLAND
Conductivity: CONSTANT_PRANDTL
EOS : IDEAL_GAS

- Fluid properties:

Ref. Viscosity = 1.716e-05 N.s/m^2
Sutherland Temp. = 273.15 K
Sutherland Const = 110.4 K

Prandtl (Lam.) = 0.72
Prandtl (Turb.) = 0.9

Gas Constant = 287.058 N.m/kg.K
Spec. Heat Ratio = 1.4

Initial and free-stream conditions:


Static Pressure = 84257.3 Pa
Density = 1.04824 kg/m^3
Temperature = 280.011 K
Total Energy = 453854 m^2/s^2
Velocity-X = 700.769 m/s
Velocity-Y = 0 m/s
Velocity-Z = 0 m/s
Velocity Magnitude = 700.769 m/s

Viscosity = 1.74976e-05 N.s/m^2
Conductivity = 85310.9 W/m^2.K
Turb. Kin. Energy = 7366.16 m^2/s^2
Spec. Dissipation = 4.41291e+07 1/s

Mach Number = 2.089
Reynolds Number = 1.8262e+06



My geometry consists of a 3D structured grid generated via ICEM (.cgns formati) with minimum orthogonality 21.81.

+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 21.1814| 89.9999|
| CV Face Area Aspect Ratio| 1.05976| 1.35169e+06|
| CV Sub-Volume Ratio| 1.00003| 2930.34|
+--------------------------------------------------------------+


I only have FAR_FIELD, SYMMETRY and WALL markers. The walls are assumed as adiabatic, and my turbulent properties are: FREESTREAM_TURBULENCEINTENSITY= 0.10 & FREESTREAM_TURB2LAMVISCRATIO= 10.0

I am solving the flow with 1st order upwind schemes, linear solver: BCGSTAB with ILU preconditioner. My CFL is 0.2. I Use 3 multigrid V cycles.


Can anyone help me to figure out what´s going wrong with the setup? I get the following error:

void CSolver::SetResidual_RMS(const CGeometry*, const CConfig*). My residual for omega is always positive during the calculation. So I guess is related to the turbulence initialization...
Attached Images
File Type: jpg Layout.jpg (91.4 KB, 19 views)
alainislas is offline   Reply With Quote

Old   October 20, 2021, 06:36
Default
  #2
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 352
Rep Power: 8
pcg is on a distinguished road
Try starting at a lower mach number.
It's a bit challenging to start a turbulent simulation Mach 2 from scratch.
pcg is offline   Reply With Quote

Old   October 20, 2021, 09:13
Default
  #3
Member
 
Alain Islas
Join Date: Nov 2019
Location: Mexico
Posts: 43
Rep Power: 4
alainislas is on a distinguished road
Dear Pedro,

Thank you for your reply. I started a new simulation with Ma = 1. I could run the case for more iterations, this time the Omega Residual reached negative values, but I had the same error.

I checked the flow in paraview, and I found some weird cells with excessive temperature (clipped by >2000 K). Do you think the divergence is attributed to this? Are there any temperature limiters in SU2?
Attached Images
File Type: png su2_0.png (130.6 KB, 13 views)
File Type: png su2_1.png (26.6 KB, 11 views)
alainislas is offline   Reply With Quote

Old   October 21, 2021, 12:20
Default
  #4
pcg
Senior Member
 
Pedro Gomes
Join Date: Dec 2017
Posts: 352
Rep Power: 8
pcg is on a distinguished road
Try turning off the multigrid and using fgmres instead of bcgstab.
Our multigrid does not work very well with the type of mesh you have (high aspect ratio in the Fairfield).
pcg 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
[snappyHexMesh] snappyHexMesh error "Cannot determine normal vector from patches." lethu OpenFOAM Meshing & Mesh Conversion 1 June 3, 2020 08:49
Fatal error: invalid wall function specification kcc49 OpenFOAM Running, Solving & CFD 13 September 26, 2018 05:07
[mesh manipulation] RefineMesh Error and Foam warning jiahui_93 OpenFOAM Meshing & Mesh Conversion 4 March 3, 2018 12:32
[mesh manipulation] refineMesh Error mohsen.boojari OpenFOAM Meshing & Mesh Conversion 3 March 1, 2018 23:07
chtMultiRegionSimpleFoam: crash on parallel run student666 OpenFOAM Running, Solving & CFD 3 April 20, 2017 12:05


All times are GMT -4. The time now is 02:57.