Program crashes if turbulence model in control file is none or dns

[mhadziabdic]

Dear guys,

If we specify none in control file the program crashes due to this line in Main.f90
169 ! Compute deltas for Spalart-Allmaras models
170 call Turb_Mod_Calculate_Deltas(turb)

I suggest to move these two lines to Turb_Mod/Init.f90 like this
22 ! Compute deltas for Spalart-Allmaras and DES models
23 if(turbulence_model .eq. SPALART_ALLMARAS .or. &
24 turbulence_model .eq. DES_SPALART) then
25 call Turb_Mod_Calculate_Deltas(turb)
26 end if

If we specify dns in control file it crashes due to this line in Save_Results.f90

420 kin_vis_t(:) = 0.0
421 if(turbulence_model .ne. NO_TURBULENCE) then
422 kin_vis_t (-grid % n_bnd_cells:grid % n_cells) = &
423 turb % vis_t(-grid % n_bnd_cells:grid % n_cells) / &
424 flow % viscosity(-grid % n_bnd_cells:grid % n_cells)
425 call Save_Scalar(grid, IN_4, IN_5, "EddyOverMolecularViscosity", &
426 plot_inside, &
427 kin_vis_t(-grid % n_bnd_cells), f8, f9)
428 end if

Program also crashes if we use none or dns and try to collect statistics.

In addition, user function Force has grid argument and there should be flow.

[Niceno]

I will work on this in the afternoon.

[Niceno]

I can't reproduce these issues. I cloned the T-Flows, stayed in development_branch. I took the case from Tests/Laminar/Backstep/Nonorthogonal and all works fine.

Furthermore, if I check User_Mod/Force.f90, it starts as:

subroutine User_Mod_Force(grid, ui, a_matrix, b_vector)

It is called from Compute_Momentum.f90 as:

call User_Mod_Force(grid, ui, a, b)

which looks the same to me.

Could you tell exactly which case did you run when you encountered problems? Are you sure you were in development_branch?