{"action":"opened","issue":{"url":"https://api.github.com/repos/HiFiLES/HiFiLES-solver/issues/121","repository_url":"https://api.github.com/repos/HiFiLES/HiFiLES-solver","labels_url":"https://api.github.com/repos/HiFiLES/HiFiLES-solver/issues/121/labels{/name}","comments_url":"https://api.github.com/repos/HiFiLES/HiFiLES-solver/issues/121/comments","events_url":"https://api.github.com/repos/HiFiLES/HiFiLES-solver/issues/121/events","html_url":"https://github.com/HiFiLES/HiFiLES-solver/issues/121","id":152710263,"number":121,"title":"Suggest to me for calculation channel flow, please ","user":{"login":"Natky","id":18625187,"avatar_url":"https://avatars.githubusercontent.com/u/18625187?v=3","gravatar_id":"","url":"https://api.github.com/users/Natky","html_url":"https://github.com/Natky","followers_url":"https://api.github.com/users/Natky/followers","following_url":"https://api.github.com/users/Natky/following{/other_user}","gists_url":"https://api.github.com/users/Natky/gists{/gist_id}","starred_url":"https://api.github.com/users/Natky/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/Natky/subscriptions","organizations_url":"https://api.github.com/users/Natky/orgs","repos_url":"https://api.github.com/users/Natky/repos","events_url":"https://api.github.com/users/Natky/events{/privacy}","received_events_url":"https://api.github.com/users/Natky/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":0,"created_at":"2016-05-03T06:42:50Z","updated_at":"2016-05-03T06:42:50Z","closed_at":null,"body":"Hi, \r\nI'm trying to test channel flow case on HiFiLES-code, when I calculated periodic channel flow 3D for lamina flow follow set up my input file, eles.cpp and my domain , at calculation round massflux = -60.7985xxx , Ubulk = -9.67758xxx and bodyforce(1) = -1, the code have show “Fatal error” \r\nThere is a Fatal error: \r\n \r\n**Fatal error 'NaN in residual, exiting.' at ../src/eles.cpp:1638**\r\n\r\nFor set up and edit “eles.cpp” at line\r\n\r\n6401    // case-specific parameters\r\n6402    // periodic channel:\r\n6403    area = 2*pi; \r\n6404    vol = 4*pi*pi;\r\n\r\n6406    // periodic hill (HIOCFD3 Case 3.4):\r\n6407    //area = 9.162;\r\n6408    //vol = 114.34;\r\n6409    //mdot0 = 9.126;  initial mass flux   ********** i dont know that 's set up comment \"mdot0\" , you know? or suggest to me, ****************\r\n\r\nand \r\n\r\n6438  // set body force for streamwise momentum and energy\r\n6439    body_force(0) = 0.;\r\n6440    //body_force(1) = alpha/area/dt*(mdot0 - mass_flux); // modified SD3D version\r\n6441    //body_force(1) = 1.0/area/dt*(mdot0 - 2.0*mass_flux + mdot_old); // HIOCFD C3.4 version\r\n6442    body_force(1) = -1.0;    //Hack for bodyforce constant\r\n6443    body_force(2) = 0.;\r\n6444    body_force(3) = 0.;\r\n6445    body_force(4) = body_force(1)*ubulk; // energy forcing\r\n\r\nHere is my set up domain periodic channel flow and grid element \r\nLx = 0 : 2*pi  \r\nLy = 0 : pi \r\nLz = -1 : 1 \r\nand cell element : nonuniform grid = ( 24 * 12 * 28 ) \r\n\r\nand here is my input file \r\n----------------------------\r\n----------------------------\r\n Solver parameters\r\n----------------------------\r\n// 0: Euler/Navier-Stokes, 1:Advection/Adv-Diffusion\r\nequation  0\r\nviscous   1\r\n\r\n// 0: Rusanov, 1: Lax-Friedrich, 2: Roe\r\nriemann_solve_type       0\r\nvis_riemann_solve_type   0\r\n\r\n// 0: Isentropic Vortex, 1: Uniform flow, 2: Sine Wave\r\nic_form    6\r\n\r\n// 0: Normal (doesn't have an analytical solution), 1:Isentropic Vortex, 2: Advection-Equation\r\ntest_case  0\r\norder      5          // Order of basis polynomials\r\ndt_type    0          // 0: User-supplied, 1: Global, 2: Local\r\ndt         0.0001\r\nCFL        1\r\nn_steps    1000000\r\nadv_type   3          // 0: Forward Euler, 3: RK45\r\ntau        1.0\r\npen_fact   0.5\r\n\r\n-----------------------\r\nMesh options\r\n-----------------------\r\nmesh_file  ch3d_pi.msh\r\ndx_cyclic  6.283185307179586\r\ndy_cyclic  3.141592653589793\r\ndz_cyclic  2.0\r\n\r\n-----------------------------------\r\nMonitoring, plotting parameters\r\n-----------------------------------\r\np_res        4            // Plotting resolution, # of nodes per direction\r\nwrite_type   0            // 0: Paraview, 1: Tecplot\r\n// Choose extra fields to be written to file: u v w energy pressure mach vorticity q_criterion. \r\n// Set to 0 or comment out for no diagnostic fields\r\nn_diagnostic_fields    5 u v w pressure mach\r\ninters_cub_order   4     // Order of cubature rule for integrating over element interfaces\r\nvolume_cub_order   4     // Order of cubature rule for integrating over element volumes\r\n\r\nmonitor_res_freq   100\r\nplot_freq          500\r\ndata_file_name     ch3d_v01zy_\r\nrestart_dump_freq  10000\r\nres_norm_type      2     // 1:L1 norm, 2:L2 norm, 0:infinity norm <not implemented>\r\nerror_norm_type    2     // 1:L1 norm, 2:L2 norm, 0:infinity norm <not implemented>\r\nres_norm_field     0     // Density\r\n\r\n---------------------------\r\nElement parameters\r\n---------------------------\r\n==== Tris ====\r\nupts_type_tri      0\r\nfpts_type_tri      1\r\nvcjh_scheme_tri    3\r\nc_tri              0.0\r\nsparse_tri         0\r\n\r\n==== Quads ====\r\nupts_type_quad     0    // 0: Gauss, 1: Gauss-Lobatto\r\nvcjh_scheme_quad   1    // 0: VCJH, 1: DG, 2: SD, 3: Hu, 4: c_+\r\neta_quad           0.0\r\nsparse_quad        0\r\n\r\n==== Hexs ====\r\nupts_type_hexa     0\r\nvcjh_scheme_hexa   1\r\neta_hexa           0.\r\nsparse_hexa        0\r\n\r\n==== Tets ====\r\nupts_type_tet      1\r\nfpts_type_tet      0\r\nvcjh_scheme_tet    1\r\neta_tet            0.0\r\nsparse_tet         0\r\n\r\n==== Prisms ====\r\nupts_type_pri_tri  0\r\nupts_type_pri_1d   0\r\nvcjh_scheme_pri_1d 1\r\neta_pri            0.0\r\nsparse_pri         0\r\n\r\n------------------------------------\r\nFluid Parameters\r\n------------------------------------\r\ngamma         1.4\r\nprandtl       0.72\r\nS_gas         120.\r\nT_gas         291.15\r\nR_gas         286.9\r\nmu_gas        1.827E-05\r\n\r\n-----------------------------------\r\nBoundary conditions\r\n-----------------------------------\r\n==== Viscous ====\r\nfix_vis           0                   // 0: Sutherland's law, 1: Constant viscosity\r\nMach_free_stream  0.2\r\nnx_free_stream    1.\r\nny_free_stream    0.\r\nnz_free_stream    0.\r\nRe_free_stream    180.\r\nL_free_stream     1.\r\nT_free_stream     300 \r\nT_wall            300\r\n//Mach_wall         0.\r\n//nx_wall           1.\r\n//ny_wall           0.\r\n//nz_wall           0.\r\n//T_wall            300\r\n\r\n==== Inviscid ====\r\nrho_bound         1.\r\nu_bound           1.0\r\nv_bound           1.0\r\nw_bound           0.0\r\np_bound           17.857142857142854098\r\n\r\n-----------\r\nForcing\r\n-----------\r\nbody_forcing  1  // Used for the periodic channel test case. 0 = no forcing, 1 = forcing\r\n-----------\r\nInitial Profile: If ic_form=6, use these parameters to describe polynomial initial condition. Useful for cyclic problems. Specify coeffs of poly of form a + bx + cx^2 + dx^3 + ex^4 + by + cy^2 + dy^3 + ey^4 + bz + cz^2 + dz^3 + ez^4\r\n-----------\r\n// Coeffs of x-momentum  [need the 13 so reader know how many quantities]\r\nx_coeffs  13  1.875 0 0 0 0 0 0 0 0 0 -3.75 0 1.875   \r\n// Coeffs of y-momentum\r\ny_coeffs  13  0 0 0 0 0 0 0 0 0 0 0 0 0\r\n// Coeffs of y-momentum\r\nz_coeffs  13  0 0 0 0 0 0 0 0 0 0 0 0 0\r\nperturb_ic   1  // If 1, add velocity perturbation to kickstart turbulence in periodic channel. 0 otherwise.\r\n\r\n\r\nPlease, suggest to me \r\nthank you, \r\n\r\nsuggest to me or contact to me : nattakit_ch@cmu.ac.th"}}