>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Setting up parameters... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> problem_type = steady Project name = steady >>> This is a steady problem -> One time step with dt=infinity. >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Reading the grid and boundary condition files... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Reading the grid file....steady.grid Total numbers: nodes = 4221 triangles = 8000 quads = 0 Boundary nodes: segments = 4 boundary 1 bnodes = 201 bfaces = 200 boundary 2 bnodes = 21 bfaces = 20 boundary 3 bnodes = 201 bfaces = 200 boundary 4 bnodes = 21 bfaces = 20 Reading the boundary condition file....steady.bcmap Boundary conditions: boundary 1 bc_type = dirichlet_all boundary 2 bc_type = dirichlet_u boundary 3 bc_type = dirichlet_all boundary 4 bc_type = dirichlet_all >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Constructing grid data... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Constructing grid data.... --- Node neighbor data: --- 2 neighbors for the node = 21 --- 2 neighbors for the node = 4201 ave_nghbr = 5 min_nghbr = 2 at node 21 max_nghbr = 6 at node 23 Generating CC scheme data...... --- Vertex-neighbor data: ave_nghbr = 11 min_nghbr = 3 elm = 39 max_nghbr = 12 elm = 43 >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Allocating arrays... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Checking the grid data... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Checking grid data.... --- Max sum of directed area vector around a node: max(sum_dav_i_x) = 1.1188966420050406E-016 max(sum_dav_i_y) = 1.1102230246251565E-016 --- Global sum of the directed area vector: sum of dav_x = -1.518E-18 sum of dav_y = -1.735E-17 --- Global sum of the boundary face vector: sum of bfn_x = 2.220E-16 sum of bfn_y = 0.000E+00 minimum element volume = 2.179992900428395E-05 maximum element volume = 2.372721461713976E-04 average element volume = 1.250000000000004E-04 minimum dual volume = 4.166666666666656E-05 maximum dual volume = 3.420239302384017E-04 average dual volume = 2.369106846718791E-04 ------ Skewness check (NC control volume) ---------- L1(e_dot_n) = 0.35955265304501349 Min(e_dot_n) = 5.5029205525528514E-002 Max(e_dot_n) = 0.99999990611670886 ---------------------------------------------------- ------ Aspect ratio check (NC control volume) ---------- Interior nodes only L1(AR) = 9.3588743186632364 Min(AR) = 5.6966963770782213 Max(AR) = 12.873785396870369 Boundary nodes only L1(AR) = 9.7217967799298055 Min(AR) = 5.1815119314117908 Max(AR) = 12.564040943770658 -------------------------------------------------------- Grid data look good! >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Computing LSQ coefficients (at nodes)... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Constructing LSQ coefficients... ---(1) Constructing Linear LSQ coefficients... ---(2) Constructing Quadratic LSQ coefficients... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Checking the LSQ coefficients... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> --------------------------------------------------------- --------------------------------------------------------- - Checking Linear LSQ gradients... - Storing a linear function values... - Computing linear LSQ gradients.. - Computing the relative errors (L_infinity).. Max relative error in ux = 9.7699626167013776E-015 Max relative error in uy = 3.0642155479654321E-014 --------------------------------------------------------- --------------------------------------------------------- --------------------------------------------------------- --------------------------------------------------------- - Checking Quadratic LSQ gradients... - Storing a quadratic function values... - Computing quadratic LSQ gradients.. - Computing the relative errors (L_infinity).. Max relative error in ux = 1.588E-13 at (x,y) = 0.00000E+00 0.00000E+00 At this location, LSQ ux = 1.0000000000E+00: Exact ux = 1.0000000000E+00 Max relative error in uy = 1.467E-11 at (x,y) = 1.56267E-01 9.79354E-01 At this location, LSQ uy = 1.6272548876E-01: Exact uy = 1.6272548876E-01 --------------------------------------------------------- --------------------------------------------------------- >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Setting initial and exact solutions... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Writing an initial solution Tecplot filesteady_nc_tecplot_0.dat ... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Starting physical time stepping... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> *********************************************************************************** Time step = 1 : dt = 999999986991104.00 We're now going to compute the solution at t = 999999986991104.00 -------------------------------------------------------------------- Beginning of Jacobian-Free Newton-Krylov Solver u-equation p-equation q-equation CFLp= 0.00E+00 steps= 0 L1(res)= 8.61E-07 4.53E-06 1.12E-05 ---------------------------------------------------------------------------------- - preconditioner: GS(sweeps:cr) 26:9.95E-02 - preconditioner: GS(sweeps:cr) 48:9.74E-02 - preconditioner: GS(sweeps:cr) 48:9.48E-02 - preconditioner: GS(sweeps:cr) 35:9.86E-02 CFLp= 1.00E+02 steps= 1 L1(res)= 3.35E-09 7.91E-08 1.23E-07 GCR(projections:sweeps:cr)= 4: 157:1.59E-02 Note: We take these residuals (after the 1st iteration) as the initial residual to measure the reduction. ---------------------------------------------------------------------------------- - preconditioner: GS(sweeps:cr) 44:9.89E-02 - preconditioner: GS(sweeps:cr) 40:9.70E-02 - preconditioner: GS(sweeps:cr) 34:9.85E-02 - preconditioner: GS(sweeps:cr) 33:9.63E-02 CFLp= 2.53E+06 steps= 2 L1(res)= 3.49E-10 5.37E-09 1.11E-08 c.r.: 0.1044 0.0678 0.0904 GCR(projections:sweeps:cr)= 4: 151:7.10E-02 ---------------------------------------------------------------------------------- - preconditioner: GS(sweeps:cr) 46:9.55E-02 - preconditioner: GS(sweeps:cr) 38:9.97E-02 - preconditioner: GS(sweeps:cr) 38:9.76E-02 - preconditioner: GS(sweeps:cr) 48:9.58E-02 CFLp= 5.03E+06 steps= 3 L1(res)= 2.41E-11 1.10E-10 6.01E-10 c.r.: 0.0691 0.0205 0.0541 GCR(projections:sweeps:cr)= 4: 170:6.69E-02 ---------------------------------------------------------------------------------- - preconditioner: GS(sweeps:cr) 44:9.84E-02 - preconditioner: GS(sweeps:cr) 45:9.57E-02 - preconditioner: GS(sweeps:cr) 39:9.61E-02 - preconditioner: GS(sweeps:cr) 39:9.91E-02 CFLp= 7.51E+06 steps= 4 L1(res)= 4.24E-13 2.44E-12 1.97E-11 c.r.: 0.0176 0.0221 0.0327 GCR(projections:sweeps:cr)= 4: 167:2.47E-02 ---------------------------------------------------------------------------------- - preconditioner: GS(sweeps:cr) 43:9.86E-02 - preconditioner: GS(sweeps:cr) 38:9.46E-02 - preconditioner: GS(sweeps:cr) 37:9.73E-02 - preconditioner: GS(sweeps:cr) 33:9.35E-02 CFLp= 9.97E+06 steps= 5 L1(res)= 1.94E-14 9.39E-14 5.31E-13 c.r.: 0.0458 0.0385 0.0270 GCR(projections:sweeps:cr)= 4: 151:3.26E-02 ---------------------------------------------------------------------------------- - preconditioner: GS(sweeps:cr) 44:9.63E-02 - preconditioner: GS(sweeps:cr) 37:9.39E-02 - preconditioner: GS(sweeps:cr) 36:9.67E-02 - preconditioner: GS(sweeps:cr) 39:9.97E-02 CFLp= 1.24E+07 steps= 6 L1(res)= 5.46E-16 2.88E-15 2.37E-14 c.r.: 0.0280 0.0307 0.0447 GCR(projections:sweeps:cr)= 4: 156:4.33E-02 ---------------------------------------------------------------------------------- steps= 6 L1(res)= 5.46E-16 2.88E-15 2.37E-14 Converged End of Jacobian-Free Newton-Krylov Solver -------------------------------------------------------------------- Finished a steady solver... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Computing error norms... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> ------------------------------------------------------------------ Error norms computed over the interior nodes L1, L2, Linf ERRORS for u: 8.287E-05 1.050E-04 3.215E-04 ERRORS for p: 3.314E-03 4.306E-03 1.675E-02 ERRORS for q: 3.368E-03 4.331E-03 1.753E-02 Max error location Max error(u) = 3.215E-04 (x,y) = 4.90617E-01 8.76008E-01 Max error(p) = 1.675E-02 (x,y) = 5.42856E-01 1.58350E-02 Max error(q) = 1.753E-02 (x,y) = 2.91939E-01 7.90833E-01 ------------------------------------------------------------------ Error norms computed over boundary nodes: - Boundary number = 1 - Number of nodes = 201 L1, L2, Linf ERRORS for u: 0.000E+00 0.000E+00 0.000E+00 ERRORS for p: 0.000E+00 0.000E+00 0.000E+00 ERRORS for q: 0.000E+00 0.000E+00 0.000E+00 Max error location Max error(u) = 0.000E+00 (x,y) = 0.00000E+00 1.00000E+00 Max error(p) = 0.000E+00 (x,y) = 0.00000E+00 1.00000E+00 Max error(q) = 0.000E+00 (x,y) = 0.00000E+00 1.00000E+00 ------------------------------------------------------------------ Error norms computed over boundary nodes: - Boundary number = 2 - Number of nodes = 21 L1, L2, Linf ERRORS for u: 0.000E+00 0.000E+00 0.000E+00 ERRORS for p: 0.000E+00 0.000E+00 0.000E+00 ERRORS for q: 2.264E-02 2.539E-02 4.186E-02 Max error location Max error(u) = 0.000E+00 (x,y) = 0.00000E+00 0.00000E+00 Max error(p) = 0.000E+00 (x,y) = 0.00000E+00 0.00000E+00 Max error(q) = 4.186E-02 (x,y) = 6.50000E-01 0.00000E+00 ------------------------------------------------------------------ Error norms computed over boundary nodes: - Boundary number = 3 - Number of nodes = 201 L1, L2, Linf ERRORS for u: 0.000E+00 0.000E+00 0.000E+00 ERRORS for p: 0.000E+00 0.000E+00 0.000E+00 ERRORS for q: 0.000E+00 0.000E+00 0.000E+00 Max error location Max error(u) = 0.000E+00 (x,y) = 1.00000E+00 0.00000E+00 Max error(p) = 0.000E+00 (x,y) = 1.00000E+00 0.00000E+00 Max error(q) = 0.000E+00 (x,y) = 1.00000E+00 0.00000E+00 ------------------------------------------------------------------ Error norms computed over boundary nodes: - Boundary number = 4 - Number of nodes = 21 L1, L2, Linf ERRORS for u: 0.000E+00 0.000E+00 0.000E+00 ERRORS for p: 0.000E+00 0.000E+00 0.000E+00 ERRORS for q: 0.000E+00 0.000E+00 0.000E+00 Max error location Max error(u) = 0.000E+00 (x,y) = 1.00000E+00 1.00000E+00 Max error(p) = 0.000E+00 (x,y) = 1.00000E+00 1.00000E+00 Max error(q) = 0.000E+00 (x,y) = 1.00000E+00 1.00000E+00 ------------------------------------------------------------------- >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> >>> Writing Tecplot files... >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> Writing steady_nc_tecplot.dat ... Writing steady_nc_tecplot_b.dat ... Finished the EDU2D-AdvDiff solver... Bye!