Fatal error 'NaN' at ../src/eles.cpp:4521

[paipaizhu]

Recently,I calculate 3D mesh for LES ,there is a Fatal error:
Writing Tecplot file number 0....
In calculation of plot_quantutiy u ：Fatal error 'NaN' at ../src/eles.cpp:4521
In calculation of plot_quantutiy u ：Fatal error 'NaN' at ../src/eles.cpp:4521
In calculation of plot_quantutiy u ：Fatal error 'NaN' at ../src/eles.cpp:4521
In calculation of plot_quantutiy u ：Fatal error 'NaN' at ../src/eles.cpp:4521
at first,I think my 'dt' is too large ,but not matter how little it is ,the fatal error still happen！

Here is my input file

equation 0
viscous 1

riemann_solve_type 2 // 0: Rusanov, 1: Lax-Friedrich, 2: Roe
vis_riemann_solve_type 2

ic_form 1 // 0: Isentropic Vortex, 1: Uniform flow, 2: Sine Wave
test_case 0 // 0: Normal (doesn't have an analytical solution), 1:Isentropic Vortex, 2: Advection-Equation
order 1 // Order of basis polynomials
dt_type 0 // 0: User-supplied, 1: Global, 2: Local
dt 0.0000000000000000001
CFL 0.01
n_steps 10000000
adv_type 3 // 0: Forward Euler, 3: RK45
tau 1.0
pen_fact 0.5
LES 1
filter_type 0
filter_ratio 2
SGS_model 1
wall_model 0
wall_layer_thickness 3

---

Restart options

restart_flag 0
restart_iter 25000
n_restart_files 8

---

Mesh options

mesh_file naca0012_les_4.neu
dx_cyclic 100000
dy_cyclic 100000
dz_cyclic 0.075

---

Monitoring, plotting parameters

p_res 1 // Plotting resolution, # of nodes per direction
write_type 1 // 0: Paraview, 1: Tecplot
inters_cub_order 2 // Order of cubature rule for integrating over element interfaces
volume_cub_order 2 // Order of cubature rule for integrating over element volumes
plot_freq 1000
data_file_name NACA0012
diagnostic_fields 6 u v w energy pressure mach
restart_dump_freq 10000
monitor_res_freq 10
monitor_cp_freq 100 // specify output frequency of file containing cp point values on surface. 0: no output.
res_norm_type 2 // 0:infinity norm, 1:L1 norm, 2:L2 norm
error_norm_type 2 // 0:infinity norm, 1:L1 norm, 2:L2 norm
res_norm_field 0 // Density

---

Element parameters

==== Tris ====
upts_type_tri 0
fpts_type_tri 0
vcjh_scheme_tri 1
c_tri 0.0
sparse_tri 0

==== Quads ====
upts_type_quad 0 // 0: Gauss, 1: Gauss-Lobatto
vcjh_scheme_quad 4 // 0: VCJH, 1: DG, 2: SD, 3: Hu, 4: c_+
eta_quad 0.0
sparse_quad 0

==== Hexs ====
upts_type_hexa 0
vcjh_scheme_hexa 2
eta_hexa 0.
sparse_hexa 0

==== Tets ====
upts_type_tet 1
fpts_type_tet 0
vcjh_scheme_tet 0
eta_tet 0.0
sparse_tet 0

==== Prisms ====
upts_type_pri_tri 0
upts_type_pri_1d 0
vcjh_scheme_pri_1d 0
eta_pri 0.0
sparse_pri 0

---

Fluid Parameters

gamma 1.4
prandtl 0.72
S_gas 120.
T_gas 291.15
R_gas 286.9
mu_gas 1.827E-05

---

Boundary conditions

==== Viscous ====
fix_vis 0 // 0: Sutherland's law, 1: Constant viscosity
Mach_free_stream 0.3
nx_free_stream 1.
ny_free_stream 0.
nz_free_stream 0.
Re_free_stream 10000.
L_free_stream 1.
T_free_stream 300
T_wall 300
Mach_wall 0.
nx_wall 1.
ny_wall 0.
nz_wall 0.
T_wall 300
nx_c_ic 1
ny_c_ic 0
nz_c_ic 0
Re_c_ic 10000.
T_c_ic 300

[mlopez14]

Can you check combinations of dt = 0; LES = 0; order > 1 ?

[paipaizhu]

Hi:
I check what you say. But some error happen .I do not know why.
What you say is HiFiLES canont use LES in 3D mesh?
best regards
zhu

[paipaizhu]

Hi
I update HIFILES, there are no such error,while 'dt' is too small,it is 10^-12.and I found the new version cannot use parallel.when i 'make' it there are some mistake. In previous version,there is not such error.

CC libmetis/libmetis_a-wspace.o
AR libmetis.a
make[3]: Leaving directory /home/zhuzhu/Application/temp/HiFiLES-solver-0.1.0/lib/parmetis-4.0.2/metis' make[3]: Entering directory/home/zhuzhu/Application/temp/HiFiLES-solver-0.1.0/lib/parmetis-4.0.2'
CC libparmetis/libparmetis_a-akwayfm.o
In file included from ./libparmetis/parmetislib.h:21,
from libparmetis/akwayfm.c:14:
./include/parmetis.h:18:19: error: metis.h: No such file or directory
./include/parmetis.h:33:4: error: #error "Incorrect user-supplied value fo IDXTYPEWIDTH"
./include/parmetis.h:42:4: error: #error "Incorrect user-supplied value fo REALTYPEWIDTH"
In file included from ./libparmetis/parmetislib.h:21,
from libparmetis/akwayfm.c:14:

best regards
zhu