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_schemequad 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