command: ./cs_solver --mpi --app-name FLUID *************************************************************** (R) Code_Saturne Version 7.0.2 Copyright (C) 1998-2021 EDF S.A., France revision d029e16a2 build Thu Nov 25 22:29:31 2021 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: Thu Dec 2 11:26:48 2021 System: Linux 5.10.60.1-microsoft-standard-WSL2 (Ubuntu 20.04 LTS) Machine: CIL-CND137BL3Y Processor: model name : 11th Gen Intel(R) Core(TM) i9-11950H @ 2.60GHz Memory: 25390 MB User: yikun () Directory: /mnt/c/Users/y.wang/Documents/Projects/20211126_code_saturne_demo/02_Thermal_Coupling/Coupling/4_2D_DISKS/RESU_COUPLING/20211202-1126/FLUID MPI ranks: 2 (appnum attribute: 1) MPI_COMM_WORLD size: 4 OpenMP threads: 4 Processors/node: 16 Compilers used for build: C compiler: gcc (Ubuntu 9.3.0-17ubuntu1~20.04) 9.3.0 C++ compiler: g++ (Ubuntu 9.3.0-17ubuntu1~20.04) 9.3.0 Fortran compiler: GNU Fortran (Ubuntu 9.3.0-17ubuntu1~20.04) 9.3.0 MPI version 3.1 (Open MPI 4.0.3) I/O read method: collective MPI-IO (explicit offsets) I/O write method: collective MPI-IO (explicit offsets) I/O rank step: 1 External libraries for partitioning: ParMETIS 4.0.3 PT-SCOTCH 6.1.0 =============================================================== CALCULATION PREPARATION ======================= =========================================================== No error detected during the data verification cs_user_parameters.f90 and others). =========================================================== CALCULATION PARAMETERS SUMMARY ============================== ----------------------------------------------------------- ** DIMENSIONS ---------- --- Physics NVAR = 8 (Nb variables ) NSCAL = 1 (Nb scalars ) NSCAUS = 0 (Nb user scalars ) NSCAPP = 1 (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 k 1 0.1000E+01 epsilon 1 0.1000E+01 TempC 1 0.1000E+01 Hydrostatic Pres 0 0.1000E+01 ---------------------------- ISTAT = 0 ou 1 (1 for unsteady ) CDTVAR > 0 (time step multiplier ) --- Order of base time stepping scheme ISCHTP = 1 (1: order 1; 2: order 2 ) ----------------------------------------------------------- ** STOKES ------ -- Continuous phase: ISTMPF = 1 (time scheme for flow (0: explicit (THETFL = 0 ) (1: std scheme (Saturne 1.0 ) (2: 2nd-order (THETFL = 0.5 ) 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.00000E+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 ) ----------------------------------------------------------- ** SCALARS ------- ITBRRB = 0 (T or H reconstruction at bdy) ----------------------------------------------------- Variable Number ITURT ----------------------------------------------------- TempC 1 0 ------------------------------------------------------------- ------------------------------------ Variable Number RVARFL ------------------------------------ TempC 1 0.8000E+00 ------------------------------------------- ------------------------------------------------------- Variable Number ICLVFL SCAMIN SCAMAX ------------------------------------------------------- TempC 1 -1 -0.2731E+03 0.1000E+13 ------------------------------------------------------- ------------------------------------------------------------- For each scalar, the number indicates it's rank in the list of all scalars. User scalars are placed first, from 1 to NSCAUS. Specific physics scalars are placed at the end, from NSCAUS+1 to NSCAPP+NSCAUS=NSCAL. is_temerature = 0 or 1 (use Cp or not ) VISLS0 = >0 (Reference viscosity ) SIGMAS = >0 (Schmidt ) RVARFL = >0 (Rf, cf variance dissipation ) ICLVFL = 0, 1 or 2 (Variance clipping mode ) SCAMIN = (Min authorized value ) SCAMAX = (Max authorized value ) For variances, SCAMIN is ignored and SCAMAX is used only if ICLVFL = 2 ------------------------------------------------------ Scalar THETSS THETVS ------------------------------------------------------ 1 0.0000E+00 0.0000E+00 ------------------------------------------------------ THETSS = (theta for source terms ) ((1+theta).new-theta.old ) THETVS = (theta for scalar diffusivity ((1+theta).new-theta.old ) ----------------------------------------------------------- ** 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 = 100 (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 = 1 (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) ----------------------------------------------------------- Applications accessible through MPI: ------------------------------------ 1; type: "SYRTHES 4" case name: "SOLID" lead rank: 0; n_ranks: 2 (group synchronized) 2; type: "Code_Saturne 7.0.2" (this instance) case name: "FLUID" lead rank: 2; n_ranks: 2 (group synchronized) SYRTHES coupling 0: initializing MPI communication ... [ok] Local ranks = [2..3], distant ranks = [0..1]. SYRTHES coupling: coupling id: 0 version: "SYRTHES 4" local name: "SOLID" distant application name: "SOLID" MPI application id: 0 MPI root rank: 0 number of MPI ranks: 2 Postprocessing output writers: ------------------------------ -1: name: results directory: postprocessing format: EnSight Gold options: separate_meshes time dependency: fixed mesh output: every 25 time steps, at end -5: name: directory: monitoring format: time_plot options: CSV time dependency: fixed mesh output: every time step -6: name: directory: profiles format: plot options: time dependency: fixed mesh output: at end -7: name: histograms directory: histograms format: histogram options: txt time dependency: fixed mesh output: at end Reading file: mesh_input.csm Finished reading: mesh_input.csm No "partition_input/domain_number_2" file available; ---------------------------------------------------------- Partitioning 3866 cells to 2 domains on 2 ranks (ParMETIS_V3_PartKway). Total number of faces on parallel boundaries: 41 wall-clock time: 0.001289 s Number of cells per domain (histogramm): [ 1931 ; 1932 [ = 1 [ 1932 ; 1933 [ = 0 [ 1933 ; 1934 [ = 0 [ 1934 ; 1935 ] = 1 Writing file: partition_output/domain_number_2 Finished writing: partition_output/domain_number_2 Partitioning finished (0.0217 s) ---------------------------------------------------------- Halo construction with standard neighborhood ============================================ Face interfaces creation Vertex interfaces creation Halo creation Halo definition Local halo definition Distant halo creation Updating the faces -> cells connectivity Histogram of the number of cells per rank: minimum value = 1931 maximum value = 1935 1 : [ 1931 ; 1932 [ = 1 2 : [ 1932 ; 1933 [ = 0 3 : [ 1933 ; 1934 [ = 0 4 : [ 1934 ; 1935 ] = 1 ---------------------------------------------------------- Histogram of the number of standard + halo cells per rank: minimum value = 1972 maximum value = 1976 1 : [ 1972 ; 1973 [ = 1 2 : [ 1973 ; 1974 [ = 0 3 : [ 1974 ; 1975 [ = 0 4 : [ 1975 ; 1976 ] = 1 ---------------------------------------------------------- Histogram of the number of ghost cells per rank: minimum value = 41 maximum value = 41 1 : [ 41 ; 41 ] = 2 ---------------------------------------------------------- Histogram of the number of interior faces per rank: minimum value = 3105 maximum value = 3119 1 : [ 3105 ; 3107 [ = 1 2 : [ 3107 ; 3110 [ = 0 3 : [ 3110 ; 3113 [ = 0 4 : [ 3113 ; 3116 [ = 0 5 : [ 3116 ; 3119 ] = 1 ---------------------------------------------------------- Histogram of the number of boundary faces per rank: minimum value = 4050 maximum value = 4066 1 : [ 4050 ; 4053 [ = 1 2 : [ 4053 ; 4056 [ = 0 3 : [ 4056 ; 4059 [ = 0 4 : [ 4059 ; 4062 [ = 0 5 : [ 4062 ; 4066 ] = 1 ---------------------------------------------------------- Histogram of the number of interior faces per cell: minimum value = 3 maximum value = 4 1 : [ 3 ; 4 ] = 3866 ---------------------------------------------------------- Histogram of the number of neighboring domains per rank: minimum value = 1 maximum value = 1 1 : [ 1 ; 1 ] = 2 ---------------------------------------------------------- Global definition of the number of elements (cells, vertices, faces...) Synchronizing cell families 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 Numbering for interior faces: type: threads minimum maximum mean number of threads: 4 4 4 number of exclusive groups: 3 3 3 number of elements in group 0: 2169 2257 2213 number of elements in group 1: 717 744 730 number of elements in group 2: 131 206 168 estimated thread imbalance: 0.555 0.670 0.612 Numbering for boundary faces: type: threads minimum maximum mean number of threads: 4 4 4 number of exclusive groups: 1 1 1 number of elements in group 0: 4050 4066 4058 estimated thread imbalance: 0.001 0.001 0.001 ---------------------------------------------------------- Histogram of volume mesh matrix bandwidth per rank: minimum value = 1865 maximum value = 1880 1 : [ 1865 ; 1868 [ = 1 2 : [ 1868 ; 1871 [ = 0 3 : [ 1871 ; 1874 [ = 0 4 : [ 1874 ; 1877 [ = 0 5 : [ 1877 ; 1880 ] = 1 Histogram of volume mesh matrix profile/lines per rank: minimum value = 511 maximum value = 536 1 : [ 511 ; 516 [ = 1 2 : [ 516 ; 521 [ = 0 3 : [ 521 ; 526 [ = 0 4 : [ 526 ; 531 [ = 0 5 : [ 531 ; 536 ] = 1 Mesh coordinates: minimum and maximum X : -3.0000001e-01 3.0000001e-01 Y : -1.7000000e+00 -1.1000000e+00 Z : 0.0000000e+00 9.9999998e-03 Mesh Number of cells: 3866 Number of interior faces: 6183 Number of boundary faces: 8116 Number of vertices: 5398 Groups: "1" cells: 3866 boundary faces: 384 "2" boundary faces: 3866 "3" boundary faces: 3866 --- Information on the volumes Minimum control volume = 1.4004559e-07 Maximum control volume = 1.0120791e-06 Total volume for the domain = 1.8852644e-03 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 (0.000817 s) --- Information on volume zones Volume zone "cells" id = 0 Number of cells = 3866 Volume = 0.0018853 Center of gravity = (-3.1372e-07, -1.4, 0.005) Surface = -1 (not computed) Volume zone "all_cells" id = 1 Number of cells = 3866 Volume = 0.0018853 Center of gravity = (-3.1372e-07, -1.4, 0.005) Surface = -1 (not computed) --- Information on boundary zones Boundary zone "boundary_faces" id = 0 Number of faces = 8116 Surface = 0.41475 Center of gravity = (-3.8205e-09, -1.4, 0.005) Perimeter = -1 (not computed) Boundary zone "BC_1" id = 1 Number of faces = 384 Surface = 0.037692 Center of gravity = (-7.6184e-10, -1.4, 0.005) Perimeter = -1 (not computed) Boundary zone "BC_2" id = 2 Number of faces = 7732 Surface = 0.37705 Center of gravity = (-4.1262e-09, -1.4, 0.005) Perimeter = -1 (not computed) ALMAX = 0.12354E+00 (Characteristic length ) ALMAX is the cubic root of the domain volume. ALMAX is the length used to initialize the turbulence. -------------------------------------------------------------------------------- ** VARIABLES INITIALIZATION ------------------------ ----------------------------------------- Property Min. value Max. value ----------------------------------------- density 0.1855E+06 0.1855E+06 molecular_viscos 0.1830E-04 0.1830E-04 turbulent_viscos 0.5612E+03 0.5612E+03 ----------------------------------------- --- Diffusivity: ----------------------------------------------- Scalar Number Min. value Max. value ----------------------------------------------- TempC 1 0.2495E-01 0.2495E-01 ----------------------------------------------- ** INFORMATION ON BOUNDARY FACES TYPE ---------------------------------- ------------------------------------------------------------------------- Boundary type Code Nb faces ------------------------------------------------------------------------- Inlet 2 0 Smooth wall 5 384 Rough wall 6 0 Symmetry 4 7732 Free outlet 3 0 Free inlet 14 0 Convective inlet 16 0 Free surface 15 0 Undefined 1 0 ------------------------------------------------------------------------- ** Clippings for computed fields ----------------------------- field initial min initial max clips to min clips to max - ------- -------------- -------------- ------------ ------------ a k a epsilon a TempC ** 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 k 0.0006 0.0006 0.0006 0.0006 v epsilon 1.0707e-05 1.0707e-05 1.0707e-05 1.0707e-05 v TempC 20 20 20 20 v Density 1.8546e+05 1.8546e+05 1.8546e+05 1.8546e+05 v TurbVisc 561.2 561.2 561.2 561.2 v CourantNb 0 0 0 0 v FourierNb 0 0 0 0 v total_pressure 1.0134e+05 1.0134e+05 1.0134e+05 1.0134e+05 v volume_forces[X] 0 0 0 0 v volume_forces[Y] 0 0 0 0 v volume_forces[Z] 0 0 0 0 v ǁvolume_forcesǁ 0 0 0 0 ** Field values on boundary_faces ------------------------------ field minimum maximum set mean spatial mean - -------------------- -------------- -------------- -------------- -------------- v Boundary temperature 20 20 20 20 v Boundary Density 1.8546e+05 1.8546e+05 1.8546e+05 1.8546e+05 v Yplus 0 0 0 0 =============================================================== MAIN CALCULATION ================ =============================================================== ===============================================================