command: ./cs_solver --trace --logp *************************************************************** (R) Code_Saturne Version 6.3.0 Copyright (C) 1998-2020 EDF S.A., France revision 2410f48a4 build 2021年04月08日 星期四 22时08分52秒 MPI version 3.1 (Open MPI 2.1.1) The Code_Saturne CFD tool is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. The Code_Saturne CFD tool is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. *************************************************************** Local case configuration: Date: 2021年04月13日 星期二 21时11分56秒 System: Linux 4.15.0-132-generic (Ubuntu 18.04.5 LTS) Machine: hp-HP-Z6-G4-Workstation Processor: model name : Intel(R) Xeon(R) Silver 4110 CPU @ 2.10GHz Memory: 64122 MB User: wu (wu) Directory: /home/wu/Desktop/槽道湍流/LES_channelflow/test2/CASE1/RESU/20210413-2111 OpenMP threads: 1 Processors/node: 32 Compilers used for build: C compiler: gcc (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0 C++ compiler: g++ (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0 Fortran compiler: GNU Fortran (Ubuntu 7.5.0-3ubuntu1~18.04) 7.5.0 I/O read method: standard input and output, serial access I/O write method: standard input and output, serial access =============================================================== CALCULATION PREPARATION ======================= =========================================================== @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ @ @@ WARNING : WHEN READING INPUT DATA @ ========= @ CHOICE OF TIME-SCHEME @ @ TIME-SCHEME FOR VELOCITY IS SECOND ORDER @ (THETAV = 0.50E+00) @ SOME TERMS ARE HOWEVER FIRST ORDER IN TIME WITH @ THE FOLLOWING SETTINGS: @ @ parameters ISTMPF ISNO2T ISTO2T IROEXT IVIEXT ICPEXT @ Values entered 2 1 0 0 0 0 @ @ computation will go on. @ @ Check the input data given through the User Interface @ or in cs_user_parameters.f90. @ @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ @ No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 6 (Nb variables ) NSCAL = 0 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 0 (Nb specific physics scalars ) ----------------------------------------------------------- ** HOMOGENEOUS MIXTURE MODEL VoF ---------------------------------------- IVOFMT = 0 ( 0 : disabled ) ( > 0: enabled ) ----------------------------------------------------------- ** TIME STEPPING ------------- --- Per-variable properties ------------------------------------ Variable ISTAT CDTVAR ------------------------------------ Velocity 1 0.1000E+01 Pressure 0 0.1000E+01 Wall distance 0 0.1000E+01 Wall Y+ 0 0.1000E+01 ---------------------------- ISTAT = 0 ou 1 (1 for unsteady ) CDTVAR > 0 (time step multiplier ) --- Order of base time stepping scheme ISCHTP = 2 (1: order 1; 2: order 2 ) ----------------------------------------------------------- ** STOKES ------ -- Continuous phase: ISTMPF = 2 (time scheme for flow (0: explicit (THETFL = 0 ) (1: std scheme (Saturne 1.0 ) (2: 2nd-order (THETFL = 0.5 ) THETFL = 0.50000E+00 (theta for mass flow ) THETVI = 0.00000E+00 (theta for total viscosity ((1+theta).new-theta.old THETCP = -0.99900E+03 (specific heat theta-scheme ((1+theta).new-theta.old THETSN = 0.50000E+00 (Nav-Stokes S.T. theta scheme) ((1+theta).new-theta.old THETST = 0.00000E+00 (Turbulence S.T. theta-scheme) ((1+theta).new-theta.old EPSUP = 0.10000E-04 (Velocity/pressure coupling stop test ) ----------------------------------------------------------- ** WALL DISTANCE COMPUTATION ------------------------- ICDPAR = -1 ( 1: std, reread if restart (-1: std, recomputed if restrt ( 2: old, reread if restart (-2: old, recomputed if restrt ----------------------------------------------------------- ** CALCULATION MANAGEMENT ---------------------- --- Restarted calculation ISUITE = 0 (1: restarted calculation ) ILEAUX = 1 (1: read restart/auxiliary ) IECAUX = 1 (1: write checkpoint/auxiliary) --- Calculation time The numbering of time steps and the measure of simulated physical time are absolute values, and not values relative to the current calculation. NTPABS = 0 (Initial time step) NTMABS = 10000 (Final time step required) ----------------------------------------------------------- ** INPUT-OUTPUT ------------ --- Restart file NTSUIT = 0 (Checkpoint frequency ) --- Probe history files NTHIST = 1 (Output frequency ) FRHIST = -.10000E+01 (Output frequency (s) ) -- -- -- --- run_solver.log files NTLIST = 100 (Output frequency ) -- -- -- --- Additional post-processing variables (ipstdv) ipstfo = 1 (Force exerted by the fluid on the boundary) ipstyp = 1 (y+ at boundary) ipsttp = 0 (T+ at boundary) ipstft = 1 (Thermal flux at boundary) ipstnu = 0 (Dimensionless thermal flux at boundary) ----------------------------------------------------------- ** ALE METHOD (MOVING MESH) ----------- IALE = 0 (1: activated,2: CDO ) NALINF = 0 (Fluid initialization iterations) IFLXMW = 0 (ALE mass flux computation 0: thanks to vertices 1: thanks to mesh velocity) ----------------------------------------------------------- Postprocessing output writers: ------------------------------ -1: name: results directory: postprocessing format: EnSight Gold options: separate_meshes time dependency: fixed mesh output: at calculation end -5: name: directory: monitoring format: time_plot options: time dependency: fixed mesh output: every 1 time steps -6: name: directory: profiles format: plot options: time dependency: fixed mesh output: at calculation end -7: name: histograms directory: histograms format: histogram options: txt time dependency: fixed mesh output: at calculation end Reading file: mesh_input.csm Finished reading: mesh_input.csm ------------------------------------------------------- Joining number 1: Periodicity type: translation Transformation matrix: 1 0 0 0 0 1 0 0 0 0 1 3 Selection criteria: "Front or Back" Parameters for the joining operation: Shortest incident edge fraction: 0.10000 Maximum angle between joined face planes: 25.00000 Advanced joining parameters: Verbosity level: 1 Visualization level: 1 Deepest level reachable in tree building: 30 Max boxes by leaf: 25 Max ratio of linked boxes / init. boxes: 5.00000 Max ratio of boxes for distribution: 2.00000 Merge step tolerance multiplier: 1.00000 Pre-merge factor: 0.05000 Tolerance computation mode: 1 Intersection computation mode: 1 Max. number of equiv. breaks: 500 Max. number of subfaces by face: 200 Before joining Number of cells: 1080000 Number of interior faces: 3208200 Number of boundary faces: 63600 Number of vertices: 1112111 Global number of boundary faces selected for joining: 24000 Element selection successfully done. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 224159 Vertex-Vertex intersections: 224159 Other intersections: 0 Joining operation is conforming. Edge intersections and vertex creation done. Pre-merge for 12221 global element couples. No need to merge vertices. Edge removed for 0 faces (global). Join mesh cleaning done. Edge removed for 0 faces (global). Join mesh cleaning done. Merge vertices and mesh update done. Global number of faces after splitting: 36000 Global configuration after the joining operation: Global number of border faces to add: 0 Mesh cleaning done for degenerate faces. Global number of cleaned interior faces: 0 Global number of cleaned border faces: 0 After joining Number of cells: 1080000 Number of interior faces: 3232200 Number of boundary faces: 39600 Number of vertices: 1112111 Joining 1 completed (4.58 s). ------------------------------------------------------- Joining number 2: Periodicity type: translation Transformation matrix: 1 0 0 6 0 1 0 0 0 0 1 0 Selection criteria: "Inlet or Outlet" Parameters for the joining operation: Shortest incident edge fraction: 0.10000 Maximum angle between joined face planes: 25.00000 Advanced joining parameters: Verbosity level: 1 Visualization level: 1 Deepest level reachable in tree building: 30 Max boxes by leaf: 25 Max ratio of linked boxes / init. boxes: 5.00000 Max ratio of boxes for distribution: 2.00000 Merge step tolerance multiplier: 1.00000 Pre-merge factor: 0.05000 Tolerance computation mode: 1 Intersection computation mode: 1 Max. number of equiv. breaks: 500 Max. number of subfaces by face: 200 Before joining Number of cells: 1080000 Number of interior faces: 3232200 Number of boundary faces: 39600 Number of vertices: 1112111 Global number of boundary faces selected for joining: 18000 Element selection successfully done. Determination of possible face intersections: bounding-box tree layout: 3D Global number of intersections detected: 162821 Vertex-Vertex intersections: 162821 Other intersections: 0 Joining operation is conforming. Edge intersections and vertex creation done. Pre-merge for 9191 global element couples. No need to merge vertices. Edge removed for 0 faces (global). Join mesh cleaning done. Edge removed for 0 faces (global). Join mesh cleaning done. Merge vertices and mesh update done. Global number of faces after splitting: 27000 Global configuration after the joining operation: Global number of border faces to add: 0 Mesh cleaning done for degenerate faces. Global number of cleaned interior faces: 0 Global number of cleaned border faces: 0 After joining Number of cells: 1080000 Number of interior faces: 3250200 Number of boundary faces: 21600 Number of vertices: 1112111 Joining 2 completed (3.99 s). All joining operations successfully finished: Wall clock time: 8.57 ---------------------------------------------------------- Composing periodicities Halo construction with standard neighborhood ============================================ Face interfaces creation Definition of periodic vertices Vertex interfaces creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Number of cells: 1080000 ---------------------------------------------------------- Number of cells + halo cells: 1122000 ---------------------------------------------------------- Number of ghost cells: 42000 ---------------------------------------------------------- Number of interior faces: 3250200 ---------------------------------------------------------- Number of boundary faces: 21600 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 5 maximum value = 6 1 : [ 5 ; 6 ] = 1080000 ---------------------------------------------------------- Synchronizing cell families Writing file: mesh_output.csm Finished writing: mesh_output.csm Renumbering mesh: renumbering for cells: pre-numbering: no renumbering cells adjacent to ghost cells last: no numbering: no renumbering renumbering for interior faces: cell adjacency pre-ordering: lowest id first faces adjacent to ghost cells last: no numbering: multipass renumbering for boundary faces: numbering: no shared cell across threads renumbering for vertices: numbering: no renumbering Numbering for cells: type: default number of halo-adjacent elements: 1080000 Numbering for interior faces: type: default number of halo-adjacent elements: 3250200 Numbering for boundary faces: type: default number of halo-adjacent elements: 21600 ---------------------------------------------------------- Matrix bandwidth for volume mesh : 1113000 Matrix profile/lines for volume mesh : 31996 Mesh coordinates: minimum and maximum X : -8.8817842e-16 6.0000000e+00 Y : 0.0000000e+00 2.0000000e+00 Z : -4.4408921e-16 3.0000000e+00 Mesh Number of cells: 1080000 Number of interior faces: 3250200 Number of boundary faces: 21600 Number of vertices: 1112111 Periodic faces (which are also interior faces): Periodicity 1: 12000 face couples Periodicity 2: 9000 face couples Groups: "Back" interior faces: 12000 "Down" boundary faces: 10800 "Front" interior faces: 12000 "Inlet" interior faces: 9000 "Outlet" interior faces: 9000 "Up" boundary faces: 10800 --- Information on the volumes Minimum control volume = 4.1666667e-06 Maximum control volume = 6.0087473e-05 Total volume for the domain = 3.6000000e+01 Criterion 1: Orthogonality: Number of bad cells detected: 0 --> 0 % Criterion 2: Offset: Number of bad cells detected: 0 --> 0 % Criterion 3: Least-Squares Gradient Quality: Number of bad cells detected: 0 --> 0 % Criterion 4: Cells Volume Ratio: Number of bad cells detected: 0 --> 0 % Criterion 5: Guilt by Association: Number of bad cells detected: 0 --> 0 % Computing geometric quantities (3.87 s) --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 1080000 Volume = 36 Center of gravity = (3, 1, 1.5) Surface = -1 (not computed) Volume zone "all_cells" id = 1 Number of cells = 1080000 Volume = 36 Center of gravity = (3, 1, 1.5) Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 21600 Surface = 36 Center of gravity = (3, 1, 1.5) Perimeter = -1 (not computed) Boundary zone "Wall" id = 1 Number of faces = 21600 Surface = 36 Center of gravity = (3, 1, 1.5) Perimeter = -1 (not computed) Boundary zone "Inlet" id = 2 Number of faces = 0 Surface = 0 Center of gravity = (0, 0, 0) Perimeter = -1 (not computed) Boundary zone "Outlet" id = 3 Number of faces = 0 Surface = 0 Center of gravity = (0, 0, 0) Perimeter = -1 (not computed) Boundary zone "BackFront" id = 4 Number of faces = 0 Surface = 0 Center of gravity = (0, 0, 0) Perimeter = -1 (not computed) ALMAX = 0.33019E+01 (Characteristic length ) ALMAX is the cubic root of the domain volume. Definition of the LES inflow for zone "Inlet" Number of structures: 50 -------------------------------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.1000E+01 0.1000E+01 molecular_viscos 0.2000E-04 0.2000E-04 turbulent_viscos 0.0000E+00 0.0000E+00 ----------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 0 Smooth wall 5 21600 Rough wall 6 0 Symmetry 4 0 Free outlet 3 0 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- ** Field values on cells --------------------- field minimum maximum set mean spatial mean - -------------- -------------- -------------- -------------- -------------- v Velocity[X] 0 0 0 0 v Velocity[Y] 0 0 0 0 v Velocity[Z] 0 0 0 0 v ǁVelocityǁ 0 0 0 0 v Pressure 0 0 0 0 v TurbVisc 0 0 0 0 v CourantNb 0 0 0 0 v FourierNb 0 0 0 0 v total_pressure 1.0132e+05 1.0132e+05 1.0132e+05 1.0132e+05 v Wall distance 0 0 0 0 v Wall Y+ 0 0 0 0 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - ----- -------------- -------------- -------------- -------------- v Yplus 0 0 0 0 =============================================================== MAIN CALCULATION ================ =============================================================== =============================================================== INSTANT 0.500000000E-03 TIME STEP NUMBER 1 ============================================================= SIGSEGV signal (forbidden memory area access) intercepted! Call stack: 1: 0x7f1cf9905040 <+0x3f040> (libc.so.6) 2: 0x5612af3d225e (cs_solver) 3: 0x7f1cfa48dacc (libsaturne-6.3.so) 4: 0x7f1cf9df9f09 (libsaturne-6.3.so) 5: 0x7f1cfa0852ba (libsaturne-6.3.so) 6: 0x7f1cf9d81cab (libsaturne-6.3.so) 7: 0x7f1cfb751838 <+0x4838> (libcs_solver-6.3.so) 8: 0x7f1cfb751b26 (libcs_solver-6.3.so) 9: 0x7f1cf98e7bf7 <__libc_start_main+0xe7> (libc.so.6) 10: 0x5612af3d1c1a <_start+0x2a> (cs_solver) End of stack