Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 5
Name = 2DSlice, index = 1
Name = TunnelWall, index = 2
Name = OutFlow, index = 3
Name = InFlow, index = 4
Name = Airfoil, index = 5
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 4
Name = TunnelWall, index = 1
Name = OutFlow, index = 2
Name = InFlow, index = 3
Name = Airfoil, index = 4
Getting mesh sizing parameters
Mesh sizing name - 2DSlice
Mesh sizing name - Airfoil
Mesh sizing name - InFlow
Mesh sizing name - OutFlow
Mesh sizing name - TunnelWall
Done getting mesh sizing parameters
Getting 2D mesh for body 1 (of 1)
AFLR2 : ---------------------------------------
AFLR2 : AFLR2 LIBRARY
AFLR2 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR2 : TRIA/QUAD GRID GENERATOR
AFLR2 : Version Number 9.14.6
AFLR2 : Version Date 11/20/22 @ 09:21PM
AFLR2 : Compile OS Darwin 22.3.0 arm64
AFLR2 : Compile Date 03/08/23 @ 05:04PM
AFLR2 : Copyright 1994-2021, D.L. Marcum
AFLR2 : ---------------------------------------
AFLR2 IC : INPUT DATA CHECK
AFLR2 IC : Nodes, Faces = 224 0
AFLR2 : CPU Time = 0.000 seconds
AFLR2 EG : EDGE GRID RE-GENERATION
AFLR2 EG : Nodes, Edges = 224 224
AFLR2 : CPU Time = 0.000 seconds
AFLR2 IG : INITIAL GRID GENERATION
AFLR2 IG : Nodes, Faces = 0 2
AFLR2 IG : Nodes, Faces = 46 94
AFLR2 IG : Nodes, Faces = 92 186
AFLR2 IG : Nodes, Faces = 138 278
AFLR2 IG : Nodes, Faces = 184 370
AFLR2 IG : Nodes, Faces = 196 394
AFLR2 IG : Nodes, Faces = 224 450
AFLR2 IG : Nodes, Faces = 228 450
AFLR2 IG : Nodes, Faces = 224 224
AFLR2 : CPU Time = 0.000 seconds
AFLR2 GG : FIELD GRID GENERATION
AFLR2 GG : Nodes, Faces = 224 224
AFLR2 GG : Nodes, Faces = 440 656
AFLR2 GG : Nodes, Faces = 649 1074
AFLR2 GG : Nodes, Faces = 853 1482
AFLR2 GG : Nodes, Faces = 1050 1876
AFLR2 GG : Nodes, Faces = 1238 2252
AFLR2 GG : Nodes, Faces = 1421 2618
AFLR2 GG : Nodes, Faces = 1594 2964
AFLR2 GG : Nodes, Faces = 1760 3296
AFLR2 GG : Nodes, Faces = 1915 3606
AFLR2 GG : Nodes, Faces = 2063 3902
AFLR2 GG : Nodes, Faces = 2202 4180
AFLR2 GG : Nodes, Faces = 2314 4404
AFLR2 GG : Nodes, Faces = 2423 4622
AFLR2 GG : Nodes, Faces = 2528 4832
AFLR2 GG : Nodes, Faces = 2634 5044
AFLR2 GG : Nodes, Faces = 2732 5240
AFLR2 GG : Nodes, Faces = 2831 5438
AFLR2 GG : Nodes, Faces = 2927 5630
AFLR2 GG : Nodes, Faces = 3018 5812
AFLR2 GG : Nodes, Faces = 3109 5994
AFLR2 GG : Nodes, Faces = 3202 6180
AFLR2 GG : Nodes, Faces = 3295 6366
AFLR2 GG : Nodes, Faces = 3386 6548
AFLR2 GG : Nodes, Faces = 3464 6704
AFLR2 GG : Nodes, Faces = 3535 6846
AFLR2 GG : Nodes, Faces = 3594 6964
AFLR2 GG : Nodes, Faces = 3639 7054
AFLR2 GG : Nodes, Faces = 3671 7118
AFLR2 GG : Nodes, Faces = 3691 7158
AFLR2 GG : Nodes, Faces = 3705 7186
AFLR2 GG : Nodes, Faces = 3712 7200
AFLR2 GG : Nodes, Faces = 3715 7206
AFLR2 : CPU Time = 0.004 seconds
AFLR2 QI : QUALITY IMPROVEMENT
AFLR2 QI : Nodes, Faces = 3715 7206
AFLR2 : CPU Time = 0.000 seconds
AFLR2 : DONE
Writing AFLR3 file ....
Finished writing AFLR3 file
Number of nodes = 3715
Number of elements = 7430
Number of tris = 7206
Number of quad = 0
Writing SU2 file ....
SU2 expects 2D meshes be in the x-y plane... attempting to rotate mesh!
Swapping z and y coordinates!
Finished writing SU2 file
Getting CFD boundary conditions
Boundary condition name - Airfoil
Boundary condition name - InFlow
Boundary condition name - OutFlow
Boundary condition name - TunnelWall
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.5.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.5.1 "Blackbird" |
| \__ \ |_| |/ / |
| |___/\___//___| Suite (Computational Fluid Dynamics Code) |
| |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io |
| |
| The SU2 Project is maintained by the SU2 Foundation |
| (http://su2foundation.org) |
-------------------------------------------------------------------------
| Copyright 2012-2023, SU2 Contributors |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 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 |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Parsing config file for zone 0
----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 0.4.
Angle of attack (AoA): 0 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 1 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 1 m.
Reference origin for moment evaluation is (0, 0, 0).
Surface(s) where the force coefficients are evaluated: BC_5.
Surface(s) plotted in the output file: BC_5, BC_2.
Input mesh file name: ../aflr2/aflr2.su2
--------------- Space Numerical Integration ( Zone 0 ) ------------------
Roe (with entropy fix = 0.001) solver for the flow inviscid terms.
Standard Roe without low-dissipation function.
Second order integration in space, with slope limiter.
Venkatakrishnan slope-limiting method, with constant: 0.05.
The reference element size is: 1.
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.
--------------- Time Numerical Integration ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-06.
Max number of linear iterations: 10.
V Multigrid Cycle, with 2 multigrid levels.
Damping factor for the residual restriction: 0.75.
Damping factor for the correction prolongation: 0.75.
No CFL adaptation.
+-------------------------------------------+
| MG Level| Presmooth|PostSmooth|CorrectSmo|
+-------------------------------------------+
| 0| 1| 0| 0|
| 1| 2| 0| 0|
| 2| 3| 0| 0|
+-------------------------------------------+
Courant-Friedrichs-Lewy number: 10
------------------ Convergence Criteria ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 10.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-14.
Cauchy series min. value: 1e-10.
Number of Cauchy elements: 100.
Begin windowed time average at iteration 0.
-------------------- Output Information ( Zone 0 ) ----------------------
File writing frequency:
+------------------------------------+
| File| Frequency|
+------------------------------------+
| RESTART| 250|
| SURFACE_CSV| 250|
| TECPLOT| 250|
| SURFACE_TECPLOT| 250|
+------------------------------------+
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is CSV (.csv).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown_pyCAPS_SU2_aflr2.dat.
Surface file name: surface_flow_pyCAPS_SU2_aflr2.
Volume file name: flow_pyCAPS_SU2_aflr2.
Restart file name: restart_flow_pyCAPS_SU2_aflr2.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_5|
| | BC_2|
+-----------------------------------------------------------------------+
| Inlet boundary| BC_4|
+-----------------------------------------------------------------------+
| Outlet boundary| BC_3|
+-----------------------------------------------------------------------+
-------------------- Output Preprocessing ( Zone 0 ) --------------------
WARNING: SURFACE_PRESSURE_DROP can only be computed for at least 2 surfaces (outlet, inlet, ...)
Screen output fields: TIME_ITER, OUTER_ITER, INNER_ITER, CUR_TIME, TIME_STEP, WALL_TIME, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY, EFFICIENCY
History output group(s): ITER, TIME_DOMAIN, WALL_TIME, RMS_RES, AERO_COEFF
Convergence field(s): RMS_DENSITY
Warning: No (valid) fields chosen for time convergence monitoring. Time convergence monitoring inactive.
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Two dimensional problem.
3715 grid points.
7206 volume elements.
4 surface markers.
98 boundary elements in index 0 (Marker = BC_2).
24 boundary elements in index 1 (Marker = BC_3).
24 boundary elements in index 2 (Marker = BC_4).
78 boundary elements in index 3 (Marker = BC_5).
7206 triangles.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
There has been a re-orientation of 7206 TRIANGLE volume elements.
There has been a re-orientation of 73 LINE surface elements.
Identifying edges and vertices.
Setting the control volume structure.
Area of the computational grid: 999.917.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 177.06. Mean K: 1.97231. Standard deviation K: 13.0434.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 63.4343| 90|
| CV Face Area Aspect Ratio| 1| 5.23543|
| CV Sub-Volume Ratio| 1| 3.11804|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 3715| 1/1.00| 10|
| 1| 1048| 1/3.54| 7.96696|
| 2| 298| 1/3.52| 6.37252|
+-------------------------------------------+
Finding max control volume width.
Wetted area = 2.06136 m.
Area projection in the x-plane = 0.178383 m, y-plane = 0.994925 m.
Max. coordinate in the x-direction = 1 m, y-direction = 0.138389 m.
Min. coordinate in the x-direction = -0.000250564 m, y-direction = -0.0151999 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_5 is NOT a single straight.
Boundary marker BC_2 is NOT a single straight.
Computing wall distances.
-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.
-- Models:
+------------------------------------------------------------------------------+
| Viscosity Model| Conductivity Model| Fluid Model|
+------------------------------------------------------------------------------+
| -| -| STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Gas Constant| 287.058| 1| N.m/kg.K| 287.058|
| Spec. Heat Ratio| -| -| -| 1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Static Pressure| 101325| 1| Pa| 101325|
| Density| 1.22498| 1| kg/m^3| 1.22498|
| Temperature| 288.15| 1| K| 288.15|
| Total Energy| 216054| 1| m^2/s^2| 216054|
| Velocity-X| 136.119| 1| m/s| 136.119|
| Velocity-Y| 0| 1| m/s| 0|
| Velocity Magnitude| 136.119| 1| m/s| 136.119|
+------------------------------------------------------------------------------+
| Mach Number| -| -| -| 0.4|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.
Initialize Jacobian structure (Euler). MG level: 1.
Initialize Jacobian structure (Euler). MG level: 2.
------------------- Numerics Preprocessing ( Zone 0 ) -------------------
----------------- Integration Preprocessing ( Zone 0 ) ------------------
------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.
------------------------------ Begin Solver -----------------------------
Simulation Run using the Single-zone Driver
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
+----------------------------------------------------------------------------------------------------------------------------------------------+
| Time_Iter| Outer_Iter| Inner_Iter| Cur_Time| Time_Step| Time(sec)| rms[Rho]| rms[RhoU]| rms[RhoV]| rms[RhoE]| CEff|
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 0| 0| 0| 0.0000e+00| 0.0000e+00| 2.5309e-02| -0.385358| 2.453598| 1.625650| 5.369499| 0.360138|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 2.9025e-02| -0.512634| 2.087781| 1.374671| 4.973554| 1.183036|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 2.9912e-02| -0.792582| 1.733529| 1.142729| 4.657622| 3.896320|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 2.9686e-02| -0.880738| 1.486539| 1.023158| 4.334759| 14.592265|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 3.0788e-02| -0.987528| 1.351478| 0.932491| 4.191092| 111.166022|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 3.0358e-02| -1.104038| 1.188715| 0.874764| 3.951825| -483.924340|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 3.1048e-02| -1.178622| 1.124587| 0.844310| 3.820426| 127.331200|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 3.1387e-02| -1.243517| 1.097463| 0.824869| 3.768858| 40.162273|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 3.1249e-02| -1.287932| 1.075843| 0.810639| 3.733440| 25.154653|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 3.1643e-02| -1.300689| 1.046363| 0.775886| 3.684245| 21.097801|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -1.30069| < -14| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_pyCAPS_SU2_aflr2.dat |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_aflr2.dat).
|CSV file |surface_flow_pyCAPS_SU2_aflr2.csv |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_aflr2.dat).
|Tecplot binary |flow_pyCAPS_SU2_aflr2.szplt |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_aflr2.dat).
|Tecplot binary surface |surface_flow_pyCAPS_SU2_aflr2.szplt|
Error in tecFileWriterClose: Not all node map values for zone 1 have been written.
Currently need 352 more values.0: Error finishing Tecplot file output.
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_aflr2.dat).
+-----------------------------------------------------------------------+
------------------------- Solver Postprocessing -------------------------
Deleted CNumerics container.
Deleted CIntegration container.
Deleted CSolver container.
Deleted CIteration container.
Deleted CInterface container.
Deleted CGeometry container.
Deleted CFreeFormDefBox class.
Deleted CSurfaceMovement class.
Deleted CVolumetricMovement class.
Deleted CConfig container.
Deleted nInst container.
Deleted COutput class.
-------------------------------------------------------------------------
------------------------- Exit Success (SU2_CFD) ------------------------
Total Force - Pressure + Viscous
Cl = 0.650457 Cd = 0.030831
Cmz = 0.373226
Cx = 0.030831 Cy = 0.650457
Pressure Contribution
Cl_p = 0.650457 Cd_p = 0.030831
Cmz_p = 0.373226
Cx_p = 0.030831 Cy_p = 0.650457
Viscous Contribution
Cl_v = 0.0 Cd_v = 0.0
Cmz_v = 0.0
Cx_v = 0.0 Cy_v = 0.0
real 0m2.969s
user 0m1.790s
sys 0m0.562s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_AFLR2_NodeDist_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'su2_and_AFLR2_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
su2_and_AFLR2_PyTest.py test;
+ python -u su2_and_AFLR2_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Saving geometry
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 5
Name = 2DSlice, index = 1
Name = TunnelWall, index = 2
Name = OutFlow, index = 3
Name = InFlow, index = 4
Name = Airfoil, index = 5
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 4
Name = TunnelWall, index = 1
Name = OutFlow, index = 2
Name = InFlow, index = 3
Name = Airfoil, index = 4
Getting mesh sizing parameters
Mesh sizing name - 2DSlice
Mesh sizing name - Airfoil
Mesh sizing name - InFlow
Mesh sizing name - OutFlow
Mesh sizing name - TunnelWall
Done getting mesh sizing parameters
Getting 2D mesh for body 1 (of 1)
AFLR2 : ---------------------------------------
AFLR2 : AFLR2 LIBRARY
AFLR2 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR2 : TRIA/QUAD GRID GENERATOR
AFLR2 : Version Number 9.14.6
AFLR2 : Version Date 11/20/22 @ 09:21PM
AFLR2 : Compile OS Darwin 22.3.0 arm64
AFLR2 : Compile Date 03/08/23 @ 05:04PM
AFLR2 : Copyright 1994-2021, D.L. Marcum
AFLR2 : ---------------------------------------
AFLR2 IC : INPUT DATA CHECK
AFLR2 IC : Nodes, Faces = 545 0
AFLR2 : CPU Time = 0.000 seconds
AFLR2 EG : EDGE GRID RE-GENERATION
AFLR2 EG : Nodes, Edges = 545 545
AFLR2 : CPU Time = 0.000 seconds
AFLR2 IG : INITIAL GRID GENERATION
AFLR2 IG : Nodes, Faces = 0 2
AFLR2 IG : Nodes, Faces = 110 222
AFLR2 IG : Nodes, Faces = 220 442
AFLR2 IG : Nodes, Faces = 330 662
AFLR2 IG : Nodes, Faces = 440 882
AFLR2 IG : Nodes, Faces = 517 1036
AFLR2 IG : Nodes, Faces = 545 1092
AFLR2 IG : Nodes, Faces = 549 1092
AFLR2 IG : Nodes, Faces = 545 545
AFLR2 : CPU Time = 0.001 seconds
AFLR2 GG : FIELD GRID GENERATION
AFLR2 GGp: Nodes, Faces = 545 545
AFLR2 GGp: Nodes, Faces = 1083 1621
AFLR2 GGp: Nodes, Faces = 1616 2687
AFLR2 GGp: Nodes, Faces = 2144 3743
AFLR2 GGp: Nodes, Faces = 2661 4777
AFLR2 GGp: Nodes, Faces = 3171 5797
AFLR2 GGp: Nodes, Faces = 3672 6799
AFLR2 GGp: Nodes, Faces = 4158 7771
AFLR2 GGp: Nodes, Faces = 4653 8761
AFLR2 GGp: Nodes, Faces = 5142 9739
AFLR2 GGp: Nodes, Faces = 5608 10671
AFLR2 GGp: Nodes, Faces = 6055 11565
AFLR2 GGp: Nodes, Faces = 6482 12419
AFLR2 GGp: Nodes, Faces = 6880 13215
AFLR2 GGp: Nodes, Faces = 7257 13969
AFLR2 GGp: Nodes, Faces = 7614 14683
AFLR2 GGp: Nodes, Faces = 7947 15349
AFLR2 GGp: Nodes, Faces = 8276 16007
AFLR2 GGp: Nodes, Faces = 8575 16605
AFLR2 GGp: Nodes, Faces = 8869 17193
AFLR2 GGp: Nodes, Faces = 9140 17735
AFLR2 GGp: Nodes, Faces = 9401 18257
AFLR2 GGp: Nodes, Faces = 9654 18763
AFLR2 GGp: Nodes, Faces = 9887 19229
AFLR2 GGp: Nodes, Faces = 10110 19675
AFLR2 GGp: Nodes, Faces = 10315 20085
AFLR2 GGp: Nodes, Faces = 10497 20449
AFLR2 GGp: Nodes, Faces = 10656 20767
AFLR2 GGp: Nodes, Faces = 10791 21037
AFLR2 GGp: Nodes, Faces = 10919 21293
AFLR2 GGp: Nodes, Faces = 11022 21499
AFLR2 GGp: Nodes, Faces = 11092 21639
AFLR2 GGp: Nodes, Faces = 11145 21745
AFLR2 GGp: Nodes, Faces = 11156 21767
AFLR2 : CPU Time = 0.029 seconds
AFLR2 QI : QUALITY IMPROVEMENT
AFLR2 QI : Nodes, Faces = 11156 21767
AFLR2 QI : Nodes, Faces = 11155 819
AFLR2 QI : Nodes, Faces = 11155 10473
AFLR2 QI : Nodes, Faces = 11155 21765
AFLR2 QI : Nodes, Faces = 11154 811
AFLR2 QI : Nodes, Faces = 11154 10476
AFLR2 QI : Nodes, Faces = 11154 21763
AFLR2 QI : Nodes, Faces = 11154 805
AFLR2 QI : Nodes, Faces = 11154 10479
AFLR2 QI : Nodes, Faces = 11154 21763
AFLR2 QI : Nodes, Faces = 11154 911
AFLR2 QI : Nodes, Faces = 11154 10426
AFLR2 QI : Nodes, Faces = 11154 21763
AFLR2 QI : Nodes, Faces = 11153 915
AFLR2 QI : Nodes, Faces = 11153 10423
AFLR2 QI : Nodes, Faces = 11153 917
AFLR2 QI : Nodes, Faces = 11153 10423
AFLR2 : CPU Time = 0.015 seconds
AFLR2 : DONE
EGADS Internal: Face 1, Side 11077 11153 complete [but 18248] (EG_makeConnect)!
EGADS Internal: Face 1, Side 10985 11077 complete [but 18248] (EG_makeConnect)!
EGADS Internal: Face 1, Side 10985 11152 complete [but 18249] (EG_makeConnect)!
EGADS Internal: Face 1, Side 11152 11153 complete [but 19022] (EG_makeConnect)!
Writing AFLR3 file ....
Finished writing AFLR3 file
Number of nodes = 11153
Number of elements = 11885
Number of tris = 917
Number of quad = 10423
Writing TECPLOT file: pyCAPS_su2_aflr2.dat ....
Finished writing TECPLOT file
EGADS Internal: Face 1, Side 11077 11153 complete [but 18248] (EG_makeConnect)!
EGADS Internal: Face 1, Side 10985 11077 complete [but 18248] (EG_makeConnect)!
EGADS Internal: Face 1, Side 10985 11152 complete [but 18249] (EG_makeConnect)!
EGADS Internal: Face 1, Side 11152 11153 complete [but 19022] (EG_makeConnect)!
EGADS Internal: Face 1, Side 11077 11153 complete [but 18248] (EG_makeConnect)!
EGADS Internal: Face 1, Side 10985 11077 complete [but 18248] (EG_makeConnect)!
EGADS Internal: Face 1, Side 10985 11152 complete [but 18249] (EG_makeConnect)!
EGADS Internal: Face 1, Side 11152 11153 complete [but 19022] (EG_makeConnect)!
Writing SU2 file ....
SU2 expects 2D meshes be in the x-y plane... attempting to rotate mesh!
Swapping z and y coordinates!
Finished writing SU2 file
Getting CFD boundary conditions
Boundary condition name - Airfoil
Boundary condition name - InFlow
Boundary condition name - OutFlow
Boundary condition name - TunnelWall
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.5.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.5.1 "Blackbird" |
| \__ \ |_| |/ / |
| |___/\___//___| Suite (Computational Fluid Dynamics Code) |
| |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io |
| |
| The SU2 Project is maintained by the SU2 Foundation |
| (http://su2foundation.org) |
-------------------------------------------------------------------------
| Copyright 2012-2023, SU2 Contributors |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 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 |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Parsing config file for zone 0
----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 0.4.
Angle of attack (AoA): 0 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 1 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 1 m.
Reference origin for moment evaluation is (0, 0, 0).
Surface(s) where the force coefficients are evaluated: BC_5.
Surface(s) plotted in the output file: BC_5, BC_2.
Input mesh file name: ../aflr2/aflr2.su2
--------------- Space Numerical Integration ( Zone 0 ) ------------------
Roe (with entropy fix = 0.001) solver for the flow inviscid terms.
Standard Roe without low-dissipation function.
Second order integration in space, with slope limiter.
Venkatakrishnan slope-limiting method, with constant: 0.05.
The reference element size is: 1.
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.
--------------- Time Numerical Integration ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-06.
Max number of linear iterations: 10.
V Multigrid Cycle, with 2 multigrid levels.
Damping factor for the residual restriction: 0.75.
Damping factor for the correction prolongation: 0.75.
No CFL adaptation.
+-------------------------------------------+
| MG Level| Presmooth|PostSmooth|CorrectSmo|
+-------------------------------------------+
| 0| 1| 0| 0|
| 1| 2| 0| 0|
| 2| 3| 0| 0|
+-------------------------------------------+
Courant-Friedrichs-Lewy number: 10
------------------ Convergence Criteria ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 10.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-14.
Cauchy series min. value: 1e-10.
Number of Cauchy elements: 100.
Begin windowed time average at iteration 0.
-------------------- Output Information ( Zone 0 ) ----------------------
File writing frequency:
+------------------------------------+
| File| Frequency|
+------------------------------------+
| RESTART| 250|
| SURFACE_CSV| 250|
| TECPLOT| 250|
| SURFACE_TECPLOT| 250|
+------------------------------------+
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is CSV (.csv).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown_pyCAPS_su2_aflr2.dat.
Surface file name: surface_flow_pyCAPS_su2_aflr2.
Volume file name: flow_pyCAPS_su2_aflr2.
Restart file name: restart_flow_pyCAPS_su2_aflr2.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_5|
| | BC_2|
+-----------------------------------------------------------------------+
| Inlet boundary| BC_4|
+-----------------------------------------------------------------------+
| Outlet boundary| BC_3|
+-----------------------------------------------------------------------+
-------------------- Output Preprocessing ( Zone 0 ) --------------------
WARNING: SURFACE_PRESSURE_DROP can only be computed for at least 2 surfaces (outlet, inlet, ...)
Screen output fields: TIME_ITER, OUTER_ITER, INNER_ITER, CUR_TIME, TIME_STEP, WALL_TIME, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY, EFFICIENCY
History output group(s): ITER, TIME_DOMAIN, WALL_TIME, RMS_RES, AERO_COEFF
Convergence field(s): RMS_DENSITY
Warning: No (valid) fields chosen for time convergence monitoring. Time convergence monitoring inactive.
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Two dimensional problem.
11153 grid points.
11340 volume elements.
4 surface markers.
98 boundary elements in index 0 (Marker = BC_2).
24 boundary elements in index 1 (Marker = BC_3).
24 boundary elements in index 2 (Marker = BC_4).
399 boundary elements in index 3 (Marker = BC_5).
917 triangles.
10423 quadrilaterals.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
There has been a re-orientation of 917 TRIANGLE volume elements.
There has been a re-orientation of 10423 QUADRILATERAL volume elements.
There has been a re-orientation of 73 LINE surface elements.
Identifying edges and vertices.
Setting the control volume structure.
Area of the computational grid: 999.878.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 385.665. Mean K: 2.15881. Standard deviation K: 17.7215.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 49.3184| 90|
| CV Face Area Aspect Ratio| 1.00035| 33.4982|
| CV Sub-Volume Ratio| 1.00001| 7.89914|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 11153| 1/1.00| 10|
| 1| 2702| 1/4.13| 7.38309|
| 2| 549| 1/4.92| 4.99198|
+-------------------------------------------+
Finding max control volume width.
Wetted area = 2.06265 m.
Area projection in the x-plane = 0.211192 m, y-plane = 0.997011 m.
Max. coordinate in the x-direction = 1 m, y-direction = 0.102037 m.
Min. coordinate in the x-direction = 0.000323229 m, y-direction = -0.108492 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_5 is NOT a single straight.
Boundary marker BC_2 is NOT a single straight.
Computing wall distances.
-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.
-- Models:
+------------------------------------------------------------------------------+
| Viscosity Model| Conductivity Model| Fluid Model|
+------------------------------------------------------------------------------+
| -| -| STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Gas Constant| 287.058| 1| N.m/kg.K| 287.058|
| Spec. Heat Ratio| -| -| -| 1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Static Pressure| 101325| 1| Pa| 101325|
| Density| 1.22498| 1| kg/m^3| 1.22498|
| Temperature| 288.15| 1| K| 288.15|
| Total Energy| 216054| 1| m^2/s^2| 216054|
| Velocity-X| 136.119| 1| m/s| 136.119|
| Velocity-Y| 0| 1| m/s| 0|
| Velocity Magnitude| 136.119| 1| m/s| 136.119|
+------------------------------------------------------------------------------+
| Mach Number| -| -| -| 0.4|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.
Initialize Jacobian structure (Euler). MG level: 1.
Initialize Jacobian structure (Euler). MG level: 2.
------------------- Numerics Preprocessing ( Zone 0 ) -------------------
----------------- Integration Preprocessing ( Zone 0 ) ------------------
------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.
------------------------------ Begin Solver -----------------------------
Simulation Run using the Single-zone Driver
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
+----------------------------------------------------------------------------------------------------------------------------------------------+
| Time_Iter| Outer_Iter| Inner_Iter| Cur_Time| Time_Step| Time(sec)| rms[Rho]| rms[RhoU]| rms[RhoV]| rms[RhoE]| CEff|
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 0| 0| 0| 0.0000e+00| 0.0000e+00| 8.9282e-02| 0.071067| 2.739763| 1.616077| 5.599701| 0.055267|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 8.7200e-02| -0.447576| 2.193837| 1.571349| 5.070382| 0.046028|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 8.8400e-02| -0.561159| 2.074597| 1.490746| 4.953698| 0.033910|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 8.7675e-02| -0.723326| 1.921155| 1.390756| 4.790278| 0.049695|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 8.6405e-02| -0.876572| 1.825653| 1.275093| 4.637161| 0.133976|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 8.6665e-02| -1.029177| 1.759196| 1.156078| 4.482282| 0.368542|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 8.6312e-02| -1.215334| 1.704909| 1.114966| 4.286060| 0.974759|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 8.6022e-02| -1.329970| 1.644937| 1.145030| 4.116982| 2.716332|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 8.5969e-02| -1.270672| 1.570739| 1.155400| 4.139480| 16.569216|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 8.5801e-02| -1.221243| 1.499646| 1.130319| 4.189559| -6.277653|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -1.22124| < -14| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_pyCAPS_su2_aflr2.dat |
Writing the forces breakdown file (forces_breakdown_pyCAPS_su2_aflr2.dat).
|CSV file |surface_flow_pyCAPS_su2_aflr2.csv |
Writing the forces breakdown file (forces_breakdown_pyCAPS_su2_aflr2.dat).
|Tecplot binary |flow_pyCAPS_su2_aflr2.szplt |
Writing the forces breakdown file (forces_breakdown_pyCAPS_su2_aflr2.dat).
|Tecplot binary surface |surface_flow_pyCAPS_su2_aflr2.szplt|
Error in tecFileWriterClose: Not all node map values for zone 1 have been written.
Currently need 994 more values.0: Error finishing Tecplot file output.
Writing the forces breakdown file (forces_breakdown_pyCAPS_su2_aflr2.dat).
+-----------------------------------------------------------------------+
------------------------- Solver Postprocessing -------------------------
Deleted CNumerics container.
Deleted CIntegration container.
Deleted CSolver container.
Deleted CIteration container.
Deleted CInterface container.
Deleted CGeometry container.
Deleted CFreeFormDefBox class.
Deleted CSurfaceMovement class.
Deleted CVolumetricMovement class.
Deleted CConfig container.
Deleted nInst container.
Deleted COutput class.
-------------------------------------------------------------------------
------------------------- Exit Success (SU2_CFD) ------------------------
Total Force - Pressure + Viscous
Cl = 0.072253 Cd = -0.01151
Cmz = 0.070975
Cx = -0.01151 Cy = 0.072253
Pressure Contribution
Cl_p = 0.072253 Cd_p = -0.01151
Cmz_p = 0.070975
Cx_p = -0.01151 Cy_p = 0.072253
Viscous Contribution
Cl_v = 0.0 Cd_v = 0.0
Cmz_v = 0.0
Cx_v = 0.0 Cy_v = 0.0
real 0m3.049s
user 0m2.914s
sys 0m0.622s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_AFLR2_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'su2_and_AFLR4_AFLR3_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
su2_and_AFLR4_AFLR3_PyTest.py test;
+ python -u su2_and_AFLR4_AFLR3_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 4
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Wake, index = 3
Name = Farfield, index = 4
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 4
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Wake, index = 3
Name = Farfield, index = 4
AFLR4 : Proximity BG Surface Grid Generation Skipped
AFLR4 : No Modifications Required
Warning: No capsGroup/capsIgnore attribute found on edge 1 of face 1, unable to assign a boundary index value
Warning: No capsGroup/capsIgnore attribute found on edge 2 of face 1, unable to assign a boundary index value
Warning: No capsGroup/capsIgnore attribute found on edge 3 of face 1, unable to assign a boundary index value
Warning: No capsGroup/capsIgnore attribute found on edge 4 of face 1, unable to assign a boundary index value
Writing TECPLOT file: AFLR4_Mesh.dat ....
Finished writing TECPLOT file
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 4
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Wake, index = 3
Name = Farfield, index = 4
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 4
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Wake, index = 3
Name = Farfield, index = 4
Getting volume mesh
AFLR2 : ---------------------------------------
AFLR2 : AFLR2 LIBRARY
AFLR2 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR2 : TRIA/QUAD GRID GENERATOR
AFLR2 : Version Number 9.14.6
AFLR2 : Version Date 11/20/22 @ 09:21PM
AFLR2 : Compile OS Darwin 22.3.0 arm64
AFLR2 : Compile Date 03/08/23 @ 05:04PM
AFLR2 : Copyright 1994-2021, D.L. Marcum
AFLR2 : ---------------------------------------
AFLR4 : ---------------------------------------
AFLR4 : AFLR4 LIBRARY
AFLR4 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR4 : TRIA/QUAD SURFACE GRID GENERATOR
AFLR4 : Version Number 11.4.5
AFLR4 : Version Date 04/20/23 @ 05:37PM
AFLR4 : Compile OS Darwin 22.4.0 arm64
AFLR4 : Compile Date 04/20/23 @ 05:37PM
AFLR4 : Copyright 1994-2021, D.L. Marcum
AFLR4 : ---------------------------------------
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.23
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.23
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
Volume mesh:
Number of nodes = 39166
Number of elements = 225915
Number of triangles = 12446
Number of quadrilatarals = 0
Number of tetrahedrals = 213469
Number of pyramids = 0
Number of prisms = 0
Number of hexahedrals = 0
Writing TECPLOT file: AFLR3_Mesh.dat ....
Finished writing TECPLOT file
Writing SU2 file ....
Finished writing SU2 file
Getting CFD boundary conditions
Boundary condition name - Farfield
Boundary condition name - Wake
Boundary condition name - Wing1
Boundary condition name - Wing2
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.5.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.5.1 "Blackbird" |
| \__ \ |_| |/ / |
| |___/\___//___| Suite (Computational Fluid Dynamics Code) |
| |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io |
| |
| The SU2 Project is maintained by the SU2 Foundation |
| (http://su2foundation.org) |
-------------------------------------------------------------------------
| Copyright 2012-2023, SU2 Contributors |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 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 |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Parsing config file for zone 0
----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 0.5901.
Angle of attack (AoA): 1 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 50 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 15.8114 m.
Surface(s) where the force coefficients are evaluated and
their reference origin for moment computation:
- BC_1 (0, 0, 0).
- BC_2 (0, 0, 0) m.
Surface(s) plotted in the output file: BC_1, BC_2.
Input mesh file name: ../aflr3/aflr3_0.su2
--------------- Space Numerical Integration ( Zone 0 ) ------------------
Roe (with entropy fix = 0.001) solver for the flow inviscid terms.
Standard Roe without low-dissipation function.
Second order integration in space, with slope limiter.
Venkatakrishnan slope-limiting method, with constant: 0.05.
The reference element size is: 1.
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.
--------------- Time Numerical Integration ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-06.
Max number of linear iterations: 10.
V Multigrid Cycle, with 2 multigrid levels.
Damping factor for the residual restriction: 0.75.
Damping factor for the correction prolongation: 0.75.
No CFL adaptation.
+-------------------------------------------+
| MG Level| Presmooth|PostSmooth|CorrectSmo|
+-------------------------------------------+
| 0| 1| 0| 0|
| 1| 2| 0| 0|
| 2| 3| 0| 0|
+-------------------------------------------+
Courant-Friedrichs-Lewy number: 10
------------------ Convergence Criteria ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 10.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-8.
Cauchy series min. value: 1e-10.
Number of Cauchy elements: 100.
Begin windowed time average at iteration 0.
-------------------- Output Information ( Zone 0 ) ----------------------
File writing frequency:
+------------------------------------+
| File| Frequency|
+------------------------------------+
| RESTART| 250|
| SURFACE_CSV| 250|
| PARAVIEW| 250|
| SURFACE_PARAVIEW| 250|
+------------------------------------+
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is CSV (.csv).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown_su2AFLRTest.dat.
Surface file name: surface_flow_su2AFLRTest.
Volume file name: flow_su2AFLRTest.
Restart file name: restart_flow_su2AFLRTest.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_1|
| | BC_2|
+-----------------------------------------------------------------------+
| Far-field| BC_4|
+-----------------------------------------------------------------------+
| Internal boundary| BC_3|
+-----------------------------------------------------------------------+
-------------------- Output Preprocessing ( Zone 0 ) --------------------
WARNING: SURFACE_PRESSURE_DROP can only be computed for at least 2 surfaces (outlet, inlet, ...)
Screen output fields: TIME_ITER, OUTER_ITER, INNER_ITER, CUR_TIME, TIME_STEP, WALL_TIME, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY, EFFICIENCY
History output group(s): ITER, TIME_DOMAIN, WALL_TIME, RMS_RES, AERO_COEFF
Convergence field(s): RMS_DENSITY
Warning: No (valid) fields chosen for time convergence monitoring. Time convergence monitoring inactive.
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Three dimensional problem.
39166 grid points.
213469 volume elements.
4 surface markers.
6594 boundary elements in index 0 (Marker = BC_1).
5358 boundary elements in index 1 (Marker = BC_2).
74 boundary elements in index 2 (Marker = BC_3).
420 boundary elements in index 3 (Marker = BC_4).
213469 tetrahedra.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
All volume elements are correctly orientend.
All surface elements are correctly orientend.
Identifying edges and vertices.
Setting the control volume structure.
Volume of the computational grid: 2.08667e+06.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 4873.1. Mean K: 61.8682. Standard deviation K: 266.942.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 21.3732| 85.3104|
| CV Face Area Aspect Ratio| 1.08044| 9403.53|
| CV Sub-Volume Ratio| 1.0189| 157083|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 39166| 1/1.00| 10|
| 1| 6505| 1/6.02| 8.24525|
| 2| 1027| 1/6.33| 6.68451|
+-------------------------------------------+
Finding max control volume width.
Semi-span length = 7.90569 m.
Wetted area = 118.038 m^2.
Area projection in the x-plane = 11.1419 m^2, y-plane = 4.96714 m^2, z-plane = 55.4182 m^2.
Max. coordinate in the x-direction = 12.6407 m, y-direction = 7.90569 m, z-direction = 0.582636 m.
Min. coordinate in the x-direction = 0 m, y-direction = -7.90569 m, z-direction = -0.41192 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_1 is NOT a single plane.
Boundary marker BC_2 is NOT a single plane.
Computing wall distances.
-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.
-- Models:
+------------------------------------------------------------------------------+
| Viscosity Model| Conductivity Model| Fluid Model|
+------------------------------------------------------------------------------+
| -| -| STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Gas Constant| 287.058| 1| N.m/kg.K| 287.058|
| Spec. Heat Ratio| -| -| -| 1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Static Pressure| 101325| 1| Pa| 101325|
| Density| 1.22498| 1| kg/m^3| 1.22498|
| Temperature| 288.15| 1| K| 288.15|
| Total Energy| 226952| 1| m^2/s^2| 226952|
| Velocity-X| 200.779| 1| m/s| 200.779|
| Velocity-Y| 0| 1| m/s| 0|
| Velocity-Z| 3.50461| 1| m/s| 3.50461|
| Velocity Magnitude| 200.809| 1| m/s| 200.809|
+------------------------------------------------------------------------------+
| Mach Number| -| -| -| 0.5901|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.
Initialize Jacobian structure (Euler). MG level: 1.
Initialize Jacobian structure (Euler). MG level: 2.
------------------- Numerics Preprocessing ( Zone 0 ) -------------------
----------------- Integration Preprocessing ( Zone 0 ) ------------------
------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.
------------------------------ Begin Solver -----------------------------
Simulation Run using the Single-zone Driver
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
+----------------------------------------------------------------------------------------------------------------------------------------------+
| Time_Iter| Outer_Iter| Inner_Iter| Cur_Time| Time_Step| Time(sec)| rms[Rho]| rms[RhoU]| rms[RhoV]| rms[RhoE]| CEff|
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 0| 0| 0| 0.0000e+00| 0.0000e+00| 1.0676e+00| -0.063778| 2.298484| 2.034108| 5.427118| 1.983988|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 1.0480e+00| -0.017700| 2.469822| 2.196237| 5.484157| 2.692721|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 1.0116e+00| 0.035482| 2.574097| 2.432587| 5.535685| 4.130422|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 1.0071e+00| -0.008022| 2.598875| 2.568018| 5.481619| 6.089410|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 1.0087e+00| -0.062467| 2.590538| 2.625188| 5.414457| 7.802347|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 1.0100e+00| -0.085318| 2.586172| 2.642507| 5.398217| 9.008384|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 1.0114e+00| -0.100720| 2.573638| 2.649855| 5.396450| 9.881394|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 1.0053e+00| -0.123403| 2.529904| 2.652341| 5.371566| 10.357805|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 9.9716e-01| -0.170356| 2.460750| 2.642108| 5.313757| 10.422913|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 9.9100e-01| -0.245372| 2.388233| 2.605782| 5.230218| 10.195974|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -0.245372| < -8| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_su2AFLRTest.dat |
Writing the forces breakdown file (forces_breakdown_su2AFLRTest.dat).
|CSV file |surface_flow_su2AFLRTest.csv |
Writing the forces breakdown file (forces_breakdown_su2AFLRTest.dat).
|Paraview |flow_su2AFLRTest.vtu |
Writing the forces breakdown file (forces_breakdown_su2AFLRTest.dat).
|Paraview surface |surface_flow_su2AFLRTest.vtu |
Writing the forces breakdown file (forces_breakdown_su2AFLRTest.dat).
+-----------------------------------------------------------------------+
------------------------- Solver Postprocessing -------------------------
Deleted CNumerics container.
Deleted CIntegration container.
Deleted CSolver container.
Deleted CIteration container.
Deleted CInterface container.
Deleted CGeometry container.
Deleted CFreeFormDefBox class.
Deleted CSurfaceMovement class.
Deleted CVolumetricMovement class.
Deleted CConfig container.
Deleted nInst container.
Deleted COutput class.
-------------------------------------------------------------------------
------------------------- Exit Success (SU2_CFD) ------------------------
Total Force - Pressure + Viscous
Cl = 0.792328 Cd = 0.07771
Cmx = 0.000776 Cmy = -0.20156 Cmz = 0.00016
Cx = 0.06387 Cy = 0.000359 Cz = 0.793563
Pressure Contribution
Cl_p = 0.792328 Cd_p = 0.07771
Cmx_p = 0.000776 Cmy_p = -0.20156 Cmz_p = 0.00016
Cx_p = 0.06387 Cy_p = 0.000359 Cz_p = 0.793563
Viscous Contribution
Cl_v = 0.0 Cd_v = 0.0
Cmx_v = 0.0 Cmy_v = 0.0 Cmz_v = 0.0
Cx_v = 0.0 Cy_v = 0.0 Cz_v = 0.0
real 0m31.179s
user 0m29.209s
sys 0m0.770s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_AFLR4_AFLR3_PyTest.py passed (as expected)
=================================================
delaundo: /Users/jenkins/util/delaundo/delaundo/src/delaundo
=================================================
+ echo 'su2_and_Delaundo_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
su2_and_Delaundo_PyTest.py test;
+ python -u su2_and_Delaundo_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 5
Name = 2DSlice, index = 1
Name = TunnelWall, index = 2
Name = OutFlow, index = 3
Name = InFlow, index = 4
Name = Airfoil, index = 5
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 5
Name = 2DSlice, index = 1
Name = TunnelWall, index = 2
Name = OutFlow, index = 3
Name = InFlow, index = 4
Name = Airfoil, index = 5
Getting mesh sizing parameters
Mesh sizing name - 2DSlice
Mesh sizing name - Airfoil
Mesh sizing name - InFlow
Mesh sizing name - OutFlow
Mesh sizing name - TunnelWall
Done getting mesh sizing parameters
Getting edge discretization for body 1
Delaundo expects 2D meshes be in the x-y plane... attempting to rotate mesh through node swapping!
Swapping z and y coordinates!
Writing delaundo control file - delaundo.ctr
Writing out *.pts file
Reading delaundo mesh file - delaundoMesh.mesh
Swapping y and z coordinates!
Writing TECPLOT file: delaundoMesh.dat ....
Finished writing TECPLOT file
Writing SU2 file ....
SU2 expects 2D meshes be in the x-y plane... attempting to rotate mesh!
Swapping z and y coordinates!
Finished writing SU2 file
Getting CFD boundary conditions
Boundary condition name - Airfoil
Boundary condition name - InFlow
Boundary condition name - OutFlow
Boundary condition name - TunnelWall
Warning: No boundary condition specified for capsGroup 2DSlice!
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.5.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.5.1 "Blackbird" |
| \__ \ |_| |/ / |
| |___/\___//___| Suite (Computational Fluid Dynamics Code) |
| |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io |
| |
| The SU2 Project is maintained by the SU2 Foundation |
| (http://su2foundation.org) |
-------------------------------------------------------------------------
| Copyright 2012-2023, SU2 Contributors |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 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 |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Parsing config file for zone 0
----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 0.4.
Angle of attack (AoA): 0 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 1 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 1 m.
Reference origin for moment evaluation is (0, 0, 0).
Surface(s) where the force coefficients are evaluated: BC_5.
Surface(s) plotted in the output file: BC_5, BC_2.
Input mesh file name: ../delaundo/delaundoMesh.su2
--------------- Space Numerical Integration ( Zone 0 ) ------------------
Roe (with entropy fix = 0.001) solver for the flow inviscid terms.
Standard Roe without low-dissipation function.
Second order integration in space, with slope limiter.
Venkatakrishnan slope-limiting method, with constant: 0.05.
The reference element size is: 1.
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.
--------------- Time Numerical Integration ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-06.
Max number of linear iterations: 10.
V Multigrid Cycle, with 2 multigrid levels.
Damping factor for the residual restriction: 0.75.
Damping factor for the correction prolongation: 0.75.
No CFL adaptation.
+-------------------------------------------+
| MG Level| Presmooth|PostSmooth|CorrectSmo|
+-------------------------------------------+
| 0| 1| 0| 0|
| 1| 2| 0| 0|
| 2| 3| 0| 0|
+-------------------------------------------+
Courant-Friedrichs-Lewy number: 10
------------------ Convergence Criteria ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 10.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-14.
Cauchy series min. value: 1e-10.
Number of Cauchy elements: 100.
Begin windowed time average at iteration 0.
-------------------- Output Information ( Zone 0 ) ----------------------
File writing frequency:
+------------------------------------+
| File| Frequency|
+------------------------------------+
| RESTART| 250|
| SURFACE_CSV| 250|
| TECPLOT| 250|
| SURFACE_TECPLOT| 250|
+------------------------------------+
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is CSV (.csv).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown_delaundoMesh.dat.
Surface file name: surface_flow_delaundoMesh.
Volume file name: flow_delaundoMesh.
Restart file name: restart_flow_delaundoMesh.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_5|
| | BC_2|
+-----------------------------------------------------------------------+
| Inlet boundary| BC_4|
+-----------------------------------------------------------------------+
| Outlet boundary| BC_3|
+-----------------------------------------------------------------------+
-------------------- Output Preprocessing ( Zone 0 ) --------------------
WARNING: SURFACE_PRESSURE_DROP can only be computed for at least 2 surfaces (outlet, inlet, ...)
Screen output fields: TIME_ITER, OUTER_ITER, INNER_ITER, CUR_TIME, TIME_STEP, WALL_TIME, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY, EFFICIENCY
History output group(s): ITER, TIME_DOMAIN, WALL_TIME, RMS_RES, AERO_COEFF
Convergence field(s): RMS_DENSITY
Warning: No (valid) fields chosen for time convergence monitoring. Time convergence monitoring inactive.
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Two dimensional problem.
6317 grid points.
12340 volume elements.
5 surface markers.
24 boundary elements in index 0 (Marker = BC_2).
24 boundary elements in index 1 (Marker = BC_3).
24 boundary elements in index 2 (Marker = BC_2).
24 boundary elements in index 3 (Marker = BC_4).
198 boundary elements in index 4 (Marker = BC_5).
12340 triangles.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
All volume elements are correctly orientend.
There has been a re-orientation of 294 LINE surface elements.
Identifying edges and vertices.
Setting the control volume structure.
Area of the computational grid: 999.917.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 1792.95. Mean K: 10.8142. Standard deviation K: 104.33.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 45.3524| 89.2595|
| CV Face Area Aspect Ratio| 1.03245| 15.7138|
| CV Sub-Volume Ratio| 1| 7.93643|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 6317| 1/1.00| 10|
| 1| 1792| 1/3.53| 7.98923|
| 2| 509| 1/3.52| 6.38684|
+-------------------------------------------+
Finding max control volume width.
Wetted area = 2.07206 m.
Area projection in the x-plane = 0.179342 m, y-plane = 1.00053 m.
Max. coordinate in the x-direction = 1 m, y-direction = 0.138492 m.
Min. coordinate in the x-direction = -0.00104578 m, y-direction = -0.0151891 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_5 is NOT a single straight.
Boundary marker BC_2 is a single straight.
Computing wall distances.
-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.
-- Models:
+------------------------------------------------------------------------------+
| Viscosity Model| Conductivity Model| Fluid Model|
+------------------------------------------------------------------------------+
| -| -| STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Gas Constant| 287.058| 1| N.m/kg.K| 287.058|
| Spec. Heat Ratio| -| -| -| 1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Static Pressure| 101325| 1| Pa| 101325|
| Density| 1.22498| 1| kg/m^3| 1.22498|
| Temperature| 288.15| 1| K| 288.15|
| Total Energy| 216054| 1| m^2/s^2| 216054|
| Velocity-X| 136.119| 1| m/s| 136.119|
| Velocity-Y| 0| 1| m/s| 0|
| Velocity Magnitude| 136.119| 1| m/s| 136.119|
+------------------------------------------------------------------------------+
| Mach Number| -| -| -| 0.4|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.
Initialize Jacobian structure (Euler). MG level: 1.
Initialize Jacobian structure (Euler). MG level: 2.
------------------- Numerics Preprocessing ( Zone 0 ) -------------------
----------------- Integration Preprocessing ( Zone 0 ) ------------------
------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.
------------------------------ Begin Solver -----------------------------
Simulation Run using the Single-zone Driver
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
+----------------------------------------------------------------------------------------------------------------------------------------------+
| Time_Iter| Outer_Iter| Inner_Iter| Cur_Time| Time_Step| Time(sec)| rms[Rho]| rms[RhoU]| rms[RhoV]| rms[RhoE]| CEff|
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 0| 0| 0| 0.0000e+00| 0.0000e+00| 5.4710e-02| -1.163068| 1.073873| 1.341547| 4.312268| 0.062082|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 4.7828e-02| -1.558471| 0.813117| 0.969484| 3.926501| 0.143132|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 5.4079e-02| -1.599269| 0.950037| 1.012154| 3.879885| 0.386300|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 5.2406e-02| -1.661284| 1.005970| 1.023067| 3.797319| 0.981260|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 5.2430e-02| -1.658421| 1.009907| 0.965279| 3.792658| 2.399755|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 5.1711e-02| -1.630602| 0.996801| 0.859556| 3.833204| 6.390579|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 5.2281e-02| -1.628009| 1.002770| 0.822847| 3.846297| 27.872536|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 5.5859e-02| -1.646113| 1.024424| 0.882502| 3.830202| -48.293527|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 5.5849e-02| -1.672152| 1.033657| 0.939200| 3.797228| -22.113788|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 5.6181e-02| -1.695131| 1.018255| 0.944150| 3.765919| -19.281368|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -1.69513| < -14| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_delaundoMesh.dat |
Writing the forces breakdown file (forces_breakdown_delaundoMesh.dat).
|CSV file |surface_flow_delaundoMesh.csv |
Writing the forces breakdown file (forces_breakdown_delaundoMesh.dat).
|Tecplot binary |flow_delaundoMesh.szplt |
Writing the forces breakdown file (forces_breakdown_delaundoMesh.dat).
|Tecplot binary surface |surface_flow_delaundoMesh.szplt |
Error in tecFileWriterClose: Not all node map values for zone 1 have been written.
Currently need 492 more values.0: Error finishing Tecplot file output.
Writing the forces breakdown file (forces_breakdown_delaundoMesh.dat).
+-----------------------------------------------------------------------+
------------------------- Solver Postprocessing -------------------------
Deleted CNumerics container.
Deleted CIntegration container.
Deleted CSolver container.
Deleted CIteration container.
Deleted CInterface container.
Deleted CGeometry container.
Deleted CFreeFormDefBox class.
Deleted CSurfaceMovement class.
Deleted CVolumetricMovement class.
Deleted CConfig container.
Deleted nInst container.
Deleted COutput class.
-------------------------------------------------------------------------
------------------------- Exit Success (SU2_CFD) ------------------------
Total Force - Pressure + Viscous
Cl = 0.805707 Cd = -0.041787
Cmz = 0.408515
Cx = -0.041787 Cy = 0.805707
Pressure Contribution
Cl_p = 0.805707 Cd_p = -0.041787
Cmz_p = 0.408515
Cx_p = -0.041787 Cy_p = 0.805707
Viscous Contribution
Cl_v = 0.0 Cd_v = 0.0
Cmz_v = 0.0
Cx_v = 0.0 Cy_v = 0.0
real 0m1.817s
user 0m1.816s
sys 0m0.452s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_Delaundo_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'su2_and_Tetgen_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
su2_and_Tetgen_PyTest.py test;
+ python -u su2_and_Tetgen_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
Getting surface mesh for body 1 (of 3)
Getting surface mesh for body 2 (of 3)
Getting surface mesh for body 3 (of 3)
Body 1 (of 3)
Number of nodes = 2236
Number of elements = 4468
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 4468
Number of quadrilateral elements = 0
Body 2 (of 3)
Number of nodes = 636
Number of elements = 1268
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 1268
Number of quadrilateral elements = 0
Body 3 (of 3)
Number of nodes = 258
Number of elements = 512
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 512
Number of quadrilateral elements = 0
----------------------------
Total number of nodes = 3130
Total number of elements = 6248
Writing TECPLOT file: egadsTessMesh_Surf_0.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: egadsTessMesh_Surf_1.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: egadsTessMesh_Surf_2.dat ....
Finished writing TECPLOT file
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
Number of surface: nodes - 3130, elements - 6248
Getting volume mesh
Generating volume mesh using TetGen.....
Tetgen input string = pYq1.500/0.000T1.00e-16A
Delaunizing vertices...
Delaunay seconds: 0.048278
Creating surface mesh ...
Surface mesh seconds: 0.012797
Recovering boundaries...
Boundary recovery seconds: 0.141893
Removing exterior tetrahedra ...
Spreading region attributes.
Exterior tets removal seconds: 0.011412
Suppressing Steiner points ...
Steiner suppression seconds: 2.6e-05
Recovering Delaunayness...
Delaunay recovery seconds: 0.021118
Refining mesh...
4180 insertions, added 2061 points, 130359 tetrahedra in queue.
1391 insertions, added 425 points, 134059 tetrahedra in queue.
1855 insertions, added 465 points, 131734 tetrahedra in queue.
2472 insertions, added 531 points, 115562 tetrahedra in queue.
3296 insertions, added 526 points, 67935 tetrahedra in queue.
4393 insertions, added 433 points, 372 tetrahedra in queue.
5856 insertions, added 571 points, 794 tetrahedra in queue.
Refinement seconds: 1.47694
Smoothing vertices...
Mesh smoothing seconds: 0.92792
Improving mesh...
Mesh improvement seconds: 0.045388
Jettisoning redundant points.
Writing nodes.
Writing elements.
Writing faces.
Writing edges.
Output seconds: 0.002561
Total running seconds: 2.68909
Statistics:
Input points: 3130
Input facets: 6248
Input segments: 9372
Input holes: 2
Input regions: 0
Mesh points: 8154
Mesh tetrahedra: 40351
Mesh faces: 83827
Mesh faces on exterior boundary: 6250
Mesh faces on input facets: 6250
Mesh edges on input segments: 9373
Steiner points on input segments: 1
Steiner points inside domain: 5023
Done meshing using TetGen!
Volume mesher did not preserve surface elements - data transfer will NOT be possible.
Writing SU2 file ....
Finished writing SU2 file
Getting CFD boundary conditions
Boundary condition name - Farfield
Boundary condition name - Wing1
Boundary condition name - Wing2
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.5.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.5.1 "Blackbird" |
| \__ \ |_| |/ / |
| |___/\___//___| Suite (Computational Fluid Dynamics Code) |
| |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io |
| |
| The SU2 Project is maintained by the SU2 Foundation |
| (http://su2foundation.org) |
-------------------------------------------------------------------------
| Copyright 2012-2023, SU2 Contributors |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 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 |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Parsing config file for zone 0
----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 0.5901.
Angle of attack (AoA): 1 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 50 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 15.8114 m.
Surface(s) where the force coefficients are evaluated and
their reference origin for moment computation:
- BC_1 (0, 0, 0).
- BC_2 (0, 0, 0) m.
Surface(s) plotted in the output file: BC_1, BC_2.
Input mesh file name: ../myMesh/tetgen_0.su2
--------------- Space Numerical Integration ( Zone 0 ) ------------------
Roe (with entropy fix = 0.001) solver for the flow inviscid terms.
Standard Roe without low-dissipation function.
Second order integration in space, with slope limiter.
Venkatakrishnan slope-limiting method, with constant: 0.05.
The reference element size is: 1.
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.
--------------- Time Numerical Integration ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-06.
Max number of linear iterations: 10.
V Multigrid Cycle, with 2 multigrid levels.
Damping factor for the residual restriction: 0.75.
Damping factor for the correction prolongation: 0.75.
No CFL adaptation.
+-------------------------------------------+
| MG Level| Presmooth|PostSmooth|CorrectSmo|
+-------------------------------------------+
| 0| 1| 0| 0|
| 1| 2| 0| 0|
| 2| 3| 0| 0|
+-------------------------------------------+
Courant-Friedrichs-Lewy number: 10
------------------ Convergence Criteria ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 5.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-14.
Cauchy series min. value: 1e-10.
Number of Cauchy elements: 100.
Begin windowed time average at iteration 0.
-------------------- Output Information ( Zone 0 ) ----------------------
File writing frequency:
+------------------------------------+
| File| Frequency|
+------------------------------------+
| RESTART| 250|
| SURFACE_CSV| 250|
| PARAVIEW| 250|
| SURFACE_PARAVIEW| 250|
+------------------------------------+
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is CSV (.csv).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown_pyCAPS_SU2_Tetgen.dat.
Surface file name: surface_flow_pyCAPS_SU2_Tetgen.
Volume file name: flow_pyCAPS_SU2_Tetgen.
Restart file name: restart_flow_pyCAPS_SU2_Tetgen.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_1|
| | BC_2|
+-----------------------------------------------------------------------+
| Far-field| BC_3|
+-----------------------------------------------------------------------+
-------------------- Output Preprocessing ( Zone 0 ) --------------------
WARNING: SURFACE_PRESSURE_DROP can only be computed for at least 2 surfaces (outlet, inlet, ...)
Screen output fields: TIME_ITER, OUTER_ITER, INNER_ITER, CUR_TIME, TIME_STEP, WALL_TIME, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY, EFFICIENCY
History output group(s): ITER, TIME_DOMAIN, WALL_TIME, RMS_RES, AERO_COEFF
Convergence field(s): RMS_DENSITY
Warning: No (valid) fields chosen for time convergence monitoring. Time convergence monitoring inactive.
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Three dimensional problem.
8154 grid points.
40351 volume elements.
3 surface markers.
4468 boundary elements in index 0 (Marker = BC_1).
1270 boundary elements in index 1 (Marker = BC_2).
512 boundary elements in index 2 (Marker = BC_3).
40351 tetrahedra.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
All volume elements are correctly orientend.
All surface elements are correctly orientend.
Identifying edges and vertices.
Setting the control volume structure.
Volume of the computational grid: 2.08741e+06.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 414.134. Mean K: 30.6158. Standard deviation K: 37.5195.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 1.93303| 83.8522|
| CV Face Area Aspect Ratio| 1.66492| 8220.52|
| CV Sub-Volume Ratio| 1.0945| 249103|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 8154| 1/1.00| 10|
| 1| 1365| 1/5.97| 8.26695|
| 2| 222| 1/6.15| 6.76877|
+-------------------------------------------+
Finding max control volume width.
Semi-span length = 7.90569 m.
Wetted area = 120.046 m^2.
Area projection in the x-plane = 10.7914 m^2, y-plane = 4.5843 m^2, z-plane = 57.287 m^2.
Max. coordinate in the x-direction = 12.6407 m, y-direction = 7.90569 m, z-direction = 0.557773 m.
Min. coordinate in the x-direction = -0.00409054 m, y-direction = -7.90569 m, z-direction = -0.41192 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_1 is NOT a single plane.
Boundary marker BC_2 is NOT a single plane.
Computing wall distances.
-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.
-- Models:
+------------------------------------------------------------------------------+
| Viscosity Model| Conductivity Model| Fluid Model|
+------------------------------------------------------------------------------+
| -| -| STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Gas Constant| 287.058| 1| N.m/kg.K| 287.058|
| Spec. Heat Ratio| -| -| -| 1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Static Pressure| 101325| 1| Pa| 101325|
| Density| 1.22498| 1| kg/m^3| 1.22498|
| Temperature| 288.15| 1| K| 288.15|
| Total Energy| 226952| 1| m^2/s^2| 226952|
| Velocity-X| 200.779| 1| m/s| 200.779|
| Velocity-Y| 0| 1| m/s| 0|
| Velocity-Z| 3.50461| 1| m/s| 3.50461|
| Velocity Magnitude| 200.809| 1| m/s| 200.809|
+------------------------------------------------------------------------------+
| Mach Number| -| -| -| 0.5901|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.
Initialize Jacobian structure (Euler). MG level: 1.
Initialize Jacobian structure (Euler). MG level: 2.
------------------- Numerics Preprocessing ( Zone 0 ) -------------------
----------------- Integration Preprocessing ( Zone 0 ) ------------------
------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.
------------------------------ Begin Solver -----------------------------
Simulation Run using the Single-zone Driver
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
+----------------------------------------------------------------------------------------------------------------------------------------------+
| Time_Iter| Outer_Iter| Inner_Iter| Cur_Time| Time_Step| Time(sec)| rms[Rho]| rms[RhoU]| rms[RhoV]| rms[RhoE]| CEff|
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 0| 0| 0| 0.0000e+00| 0.0000e+00| 1.6174e-01| 0.616284| 3.048260| 2.091357| 6.107180| 1.815389|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 1.5905e-01| 0.741001| 3.258088| 2.952168| 6.243294| 2.156707|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 1.6392e-01| 0.650783| 3.322878| 3.201040| 6.147243| 3.157692|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 1.6417e-01| 0.569960| 3.293032| 3.250810| 6.036360| 4.529265|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 1.6448e-01| 0.533226| 3.191940| 3.195181| 5.994483| 5.686317|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 5) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| 0.533226| < -14| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_pyCAPS_SU2_Tetgen.dat |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Tetgen.dat).
|CSV file |surface_flow_pyCAPS_SU2_Tetgen.csv |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Tetgen.dat).
|Paraview |flow_pyCAPS_SU2_Tetgen.vtu |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Tetgen.dat).
|Paraview surface |surface_flow_pyCAPS_SU2_Tetgen.vtu |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Tetgen.dat).
+-----------------------------------------------------------------------+
------------------------- Solver Postprocessing -------------------------
Deleted CNumerics container.
Deleted CIntegration container.
Deleted CSolver container.
Deleted CIteration container.
Deleted CInterface container.
Deleted CGeometry container.
Deleted CFreeFormDefBox class.
Deleted CSurfaceMovement class.
Deleted CVolumetricMovement class.
Deleted CConfig container.
Deleted nInst container.
Deleted COutput class.
-------------------------------------------------------------------------
------------------------- Exit Success (SU2_CFD) ------------------------
Total Force - Pressure + Viscous
Cl = 0.757437 Cd = 0.133203
Cmx = -0.003663 Cmy = -0.214942 Cmz = 0.000732
Cx = 0.119964 Cy = 0.000385 Cz = 0.759646
real 0m7.970s
user 0m7.197s
sys 0m0.420s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_Tetgen_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'su2_X43a_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
su2_X43a_PyTest.py test;
+ python -u su2_X43a_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 2
Name = x43A, index = 1
Name = Farfield, index = 2
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 0
Getting surface mesh for body 1 (of 2)
Getting surface mesh for body 2 (of 2)
Body 1 (of 2)
Number of nodes = 3801
Number of elements = 7602
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 7602
Number of quadrilateral elements = 0
Body 2 (of 2)
Number of nodes = 5678
Number of elements = 11352
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 11352
Number of quadrilateral elements = 0
----------------------------
Total number of nodes = 9479
Total number of elements = 18954
Writing TECPLOT file: egadsTessMesh_Surf_0.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: egadsTessMesh_Surf_1.dat ....
Finished writing TECPLOT file
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 2
Name = x43A, index = 1
Name = Farfield, index = 2
Number of surface: nodes - 9479, elements - 18954
Getting volume mesh
Generating volume mesh using TetGen.....
Tetgen input string = pYq1.500/0.000T1.00e-16A
Delaunizing vertices...
Delaunay seconds: 0.12771
Creating surface mesh ...
Surface mesh seconds: 0.025213
Recovering boundaries...
Boundary recovery seconds: 0.224894
Removing exterior tetrahedra ...
Spreading region attributes.
Exterior tets removal seconds: 0.020409
Recovering Delaunayness...
Delaunay recovery seconds: 0.052709
Refining mesh...
12635 insertions, added 5200 points, 13361 tetrahedra in queue.
Refinement seconds: 1.31939
Smoothing vertices...
Mesh smoothing seconds: 1.11519
Improving mesh...
Mesh improvement seconds: 0.06313
Writing nodes.
Writing elements.
Writing faces.
Writing edges.
Output seconds: 0.016564
Total running seconds: 2.96624
Statistics:
Input points: 9479
Input facets: 18954
Input segments: 28431
Input holes: 1
Input regions: 0
Mesh points: 14860
Mesh tetrahedra: 61188
Mesh faces: 131853
Mesh faces on exterior boundary: 18954
Mesh faces on input facets: 18954
Mesh edges on input segments: 28431
Steiner points inside domain: 5381
Done meshing using TetGen!
Writing SU2 file ....
Finished writing SU2 file
Getting CFD boundary conditions
Boundary condition name - Farfield
Boundary condition name - x43A
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.5.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.5.1 "Blackbird" |
| \__ \ |_| |/ / |
| |___/\___//___| Suite (Computational Fluid Dynamics Code) |
| |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io |
| |
| The SU2 Project is maintained by the SU2 Foundation |
| (http://su2foundation.org) |
-------------------------------------------------------------------------
| Copyright 2012-2023, SU2 Contributors |
| |
| SU2 is free software; you can redistribute it and/or |
| modify it under the terms of the GNU Lesser General Public |
| License as published by the Free Software Foundation; either |
| version 2.1 of the License, or (at your option) any later version. |
| |
| SU2 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 |
| Lesser General Public License for more details. |
| |
| You should have received a copy of the GNU Lesser General Public |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>. |
-------------------------------------------------------------------------
Parsing config file for zone 0
----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 3.5.
Angle of attack (AoA): 3 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 1 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 1 m.
Reference origin for moment evaluation is (0, 0, 0).
Surface(s) where the force coefficients are evaluated: BC_1.
Surface(s) plotted in the output file: BC_1.
Input mesh file name: ../myMesh/tetgen_0.su2
--------------- Space Numerical Integration ( Zone 0 ) ------------------
Roe (with entropy fix = 0.001) solver for the flow inviscid terms.
Standard Roe without low-dissipation function.
Second order integration in space, with slope limiter.
Venkatakrishnan slope-limiting method, with constant: 0.05.
The reference element size is: 1.
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.
--------------- Time Numerical Integration ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-06.
Max number of linear iterations: 10.
V Multigrid Cycle, with 2 multigrid levels.
Damping factor for the residual restriction: 0.75.
Damping factor for the correction prolongation: 0.75.
No CFL adaptation.
+-------------------------------------------+
| MG Level| Presmooth|PostSmooth|CorrectSmo|
+-------------------------------------------+
| 0| 1| 0| 0|
| 1| 2| 0| 0|
| 2| 3| 0| 0|
+-------------------------------------------+
Courant-Friedrichs-Lewy number: 10
------------------ Convergence Criteria ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 5.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-14.
Cauchy series min. value: 1e-10.
Number of Cauchy elements: 100.
Begin windowed time average at iteration 0.
-------------------- Output Information ( Zone 0 ) ----------------------
File writing frequency:
+------------------------------------+
| File| Frequency|
+------------------------------------+
| RESTART| 250|
| SURFACE_CSV| 250|
| PARAVIEW| 250|
| SURFACE_PARAVIEW| 250|
+------------------------------------+
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is CSV (.csv).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown_x43a_Test.dat.
Surface file name: surface_flow_x43a_Test.
Volume file name: flow_x43a_Test.
Restart file name: restart_flow_x43a_Test.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_1|
+-----------------------------------------------------------------------+
| Far-field| BC_2|
+-----------------------------------------------------------------------+
-------------------- Output Preprocessing ( Zone 0 ) --------------------
WARNING: SURFACE_PRESSURE_DROP can only be computed for at least 2 surfaces (outlet, inlet, ...)
Screen output fields: TIME_ITER, OUTER_ITER, INNER_ITER, CUR_TIME, TIME_STEP, WALL_TIME, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY, EFFICIENCY
History output group(s): ITER, TIME_DOMAIN, WALL_TIME, RMS_RES, AERO_COEFF
Convergence field(s): RMS_DENSITY
Warning: No (valid) fields chosen for time convergence monitoring. Time convergence monitoring inactive.
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Three dimensional problem.
14860 grid points.
61188 volume elements.
2 surface markers.
7602 boundary elements in index 0 (Marker = BC_1).
11352 boundary elements in index 1 (Marker = BC_2).
61188 tetrahedra.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
All volume elements are correctly orientend.
All surface elements are correctly orientend.
Identifying edges and vertices.
Setting the control volume structure.
Volume of the computational grid: 2.67745e+08.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 95.6599. Mean K: 0.525142. Standard deviation K: 2.58752.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 4.9354| 84.3083|
| CV Face Area Aspect Ratio| 1.26054| 6343.09|
| CV Sub-Volume Ratio| 1| 800252|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 14860| 1/1.00| 10|
| 1| 2681| 1/5.54| 8.47586|
| 2| 494| 1/5.43| 7.23467|
+-------------------------------------------+
Finding max control volume width.
Semi-span length = 30 m.
Wetted area = 13592.9 m^2.
Area projection in the x-plane = 761.521 m^2, y-plane = 2060.65 m^2, z-plane = 5463.74 m^2.
Max. coordinate in the x-direction = 150 m, y-direction = 30 m, z-direction = 19 m.
Min. coordinate in the x-direction = 0 m, y-direction = -30 m, z-direction = -10.3 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_1 is NOT a single plane.
Computing wall distances.
-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.
-- Models:
+------------------------------------------------------------------------------+
| Viscosity Model| Conductivity Model| Fluid Model|
+------------------------------------------------------------------------------+
| -| -| STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Gas Constant| 287.058| 1| N.m/kg.K| 287.058|
| Spec. Heat Ratio| -| -| -| 1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
| Name| Dim. value| Ref. value| Unit|Non-dim. value|
+------------------------------------------------------------------------------+
| Static Pressure| 101325| 1| Pa| 101325|
| Density| 1.22498| 1| kg/m^3| 1.22498|
| Temperature| 288.15| 1| K| 288.15|
| Total Energy| 916077| 1| m^2/s^2| 916077|
| Velocity-X| 1189.41| 1| m/s| 1189.41|
| Velocity-Y| 0| 1| m/s| 0|
| Velocity-Z| 62.3342| 1| m/s| 62.3342|
| Velocity Magnitude| 1191.04| 1| m/s| 1191.04|
+------------------------------------------------------------------------------+
| Mach Number| -| -| -| 3.5|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.
Initialize Jacobian structure (Euler). MG level: 1.
Initialize Jacobian structure (Euler). MG level: 2.
------------------- Numerics Preprocessing ( Zone 0 ) -------------------
----------------- Integration Preprocessing ( Zone 0 ) ------------------
------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.
------------------------------ Begin Solver -----------------------------
Simulation Run using the Single-zone Driver
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
WARNING: SU2 was not compiled for an AVX-capable architecture. Performance could be better,
see https://su2code.github.io/docs_v7/Build-SU2-Linux-MacOS/#compiler-optimizations
+----------------------------------------------------------------------------------------------------------------------------------------------+
| Time_Iter| Outer_Iter| Inner_Iter| Cur_Time| Time_Step| Time(sec)| rms[Rho]| rms[RhoU]| rms[RhoV]| rms[RhoE]| CEff|
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 0| 0| 0| 0.0000e+00| 0.0000e+00| 3.3797e-01| 2.731357| 5.665667| 4.731819| 8.730832| -0.417340|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 3.1920e-01| 3.146139| 6.181490| 5.754536| 9.163855| -0.162226|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 3.1893e-01| 3.242298| 6.314030| 5.766862| 9.216747| -0.072314|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 3.1808e-01| 3.234919| 6.292726| 5.661514| 9.202116| -0.058811|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 3.2129e-01| 3.185506| 6.236716| 5.511324| 9.155240| -0.078638|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 5) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| 3.18551| < -14| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_x43a_Test.dat |
Writing the forces breakdown file (forces_breakdown_x43a_Test.dat).
|CSV file |surface_flow_x43a_Test.csv |
Writing the forces breakdown file (forces_breakdown_x43a_Test.dat).
|Paraview |flow_x43a_Test.vtu |
Writing the forces breakdown file (forces_breakdown_x43a_Test.dat).
|Paraview surface |surface_flow_x43a_Test.vtu |
Writing the forces breakdown file (forces_breakdown_x43a_Test.dat).
+-----------------------------------------------------------------------+
Warning: there are 120 non-physical points in the solution.
Warning: 57 reconstructed states for upwinding are non-physical.
------------------------- Solver Postprocessing -------------------------
Deleted CNumerics container.
Deleted CIntegration container.
Deleted CSolver container.
Deleted CIteration container.
Deleted CInterface container.
Deleted CGeometry container.
Deleted CFreeFormDefBox class.
Deleted CSurfaceMovement class.
Deleted CVolumetricMovement class.
Deleted CConfig container.
Deleted nInst container.
Deleted COutput class.
-------------------------------------------------------------------------
------------------------- Exit Success (SU2_CFD) ------------------------
Total Force - Pressure + Viscous
Cl = -9.346996 Cd = 118.86092
Cmx = -1.671905 Cmy = 9708.921181 Cmz = 254.31444
Cx = 119.187209 Cy = 2.633602 Cz = -3.113487
real 0m17.397s
user 0m21.346s
sys 0m0.665s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_X43a_PyTest.py passed (as expected)
=================================================
./execute_PyTestRegression.sh: line 404: ulimit: stack size: cannot modify limit: Operation not permitted
flowCart: /Users/jenkins/util/cart3d/cart3d_v1.5.9_OSX64_ICC--22.11.16/bin/OSX64_ICC//flowCart
=================================================
+ echo 'cart3d_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
cart3d_PyTest.py test;
+ python -u cart3d_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Loading AIM
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.796888
Tessellating body 1 with MaxEdge = 0.094922 Sag = 0.003797 Angle = 15.000000
Executing: autoInputs -r 30.000000 -nDiv 6 -maxR 9 > autoInputs.out
Executing: ./aero.csh > aero.out
C_A 0.050108722
C_Y -0.0015626888
C_N 0.0065057682
C_D 0.050305245
C_S -0.0015626888
C_L 0.0047530359
C_l 0.00381162
C_m -0.00346146
C_n -0.0401149
C_M_x 0.00240931
C_M_y -0.00346146
C_M_z -0.0402235
real 0m43.801s
user 0m35.795s
sys 0m1.708s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case cart3d_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'cart3d_OpenMDAO_3_alpha_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
cart3d_OpenMDAO_3_alpha_PyTest.py test;
+ python -u cart3d_OpenMDAO_3_alpha_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
----------------------------------------
>>> alpha [2.]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Getting CFD design variables.......
Number of design variables - 1
Design Variable name - alpha
Done getting CFD design variables
Getting CFD functional.......
Number of design variables - 1
Objective name - CL2
Done getting CFD functional
Body size = 3.625615
Tessellating body 1 with MaxEdge = 0.090640 Sag = 0.003626 Angle = 15.000000
Executing: autoInputs -r 30.000000 -nDiv 5 -maxR 7 -mesh2d > autoInputs.out
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:23:15) Running aero.csh on 2 core(s) in M0.5901A2B0_DP1
(14:23:22) aero.csh done
o case M0.5901A2B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.01631562] CL 0.12773309
----------------------------------------
----------------------------------------
>>> alpha [2.]
>>> CL2 [0.01631562] CL 0.12773309
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:23:24) Running aero.csh on 2 core(s) in M0.5901A2B0_DP1
(14:23:32) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 2188
Estimated cutter memory: 0.000219 GB / DV
Estimated xsensit memory: 0.004157 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:23:32) Linearizing cut-cells
(14:23:32) ... done
o Number of simultaneous xsensit jobs: 1
(14:23:32) Computing gradients
(14:23:36) ... done
o case M0.5901A2B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[0.01572627]]
========================================
RUNNING THE L-BFGS-B CODE
* * *
Machine precision = 2.220D-16
N = 1 M = 10
At X0 0 variables are exactly at the bounds
At iterate 0 f= 1.63156D-02 |proj g|= 1.57263D-02
----------------------------------------
>>> alpha [1.98427373]
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:23:38) Running aero.csh on 2 core(s) in M0.5901A1.98427B0_DP1
(14:23:51) aero.csh done
o case M0.5901A1.98427B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.01606724] CL 0.12675617
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:23:53) Running aero.csh on 2 core(s) in M0.5901A1.98427B0_DP1
(14:24:01) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 2188
Estimated cutter memory: 0.000219 GB / DV
Estimated xsensit memory: 0.004157 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:24:02) Linearizing cut-cells
(14:24:02) ... done
o Number of simultaneous xsensit jobs: 1
(14:24:02) Computing gradients
(14:24:06) ... done
o case M0.5901A1.98427B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[0.01560329]]
========================================
At iterate 1 f= 1.60672D-02 |proj g|= 1.56033D-02
----------------------------------------
>>> alpha [-0.01099175]
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:24:07) Running aero.csh on 2 core(s) in M0.5901A-0.0109917B0_DP1
(14:24:14) aero.csh done
o case M0.5901A-0.0109917B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [6.05653711e-06] CL 0.00246101
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:24:16) Running aero.csh on 2 core(s) in M0.5901A-0.0109917B0_DP1
(14:24:23) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 2188
Estimated cutter memory: 0.000219 GB / DV
Estimated xsensit memory: 0.004157 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:24:23) Linearizing cut-cells
(14:24:23) ... done
o Number of simultaneous xsensit jobs: 1
(14:24:23) Computing gradients
(14:24:28) ... done
o case M0.5901A-0.0109917B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[0.00031042]]
========================================
At iterate 2 f= 6.05654D-06 |proj g|= 3.10424D-04
----------------------------------------
>>> alpha [-0.05149286]
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:24:29) Running aero.csh on 2 core(s) in M0.5901A-0.0514929B0_DP1
(14:24:35) aero.csh done
o case M0.5901A-0.0514929B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [2.59817062e-10] CL 1.612e-05
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:24:36) Running aero.csh on 2 core(s) in M0.5901A-0.0514929B0_DP1
(14:24:53) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 2188
Estimated cutter memory: 0.000219 GB / DV
Estimated xsensit memory: 0.004157 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:24:54) Linearizing cut-cells
(14:24:54) ... done
o Number of simultaneous xsensit jobs: 1
(14:24:54) Computing gradients
(14:24:58) ... done
o case M0.5901A-0.0514929B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[2.03612433e-06]]
========================================
At iterate 3 f= 2.59817D-10 |proj g|= 2.03612D-06
----------------------------------------
>>> alpha [-0.05176027]
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:24:59) Running aero.csh on 2 core(s) in M0.5901A-0.0517603B0_DP1
(14:25:23) aero.csh done
o case M0.5901A-0.0517603B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [9.25766089e-15] CL 1e-07
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:25:24) Running aero.csh on 2 core(s) in M0.5901A-0.0517603B0_DP1
(14:25:48) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 2188
Estimated cutter memory: 0.000219 GB / DV
Estimated xsensit memory: 0.004157 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:25:48) Linearizing cut-cells
(14:25:48) ... done
o Number of simultaneous xsensit jobs: 1
(14:25:49) Computing gradients
(14:25:53) ... done
o case M0.5901A-0.0517603B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[1.21372408e-08]]
========================================
At iterate 4 f= 9.25766D-15 |proj g|= 1.21372D-08
* * *
Tit = total number of iterations
Tnf = total number of function evaluations
Tnint = total number of segments explored during Cauchy searches
Skip = number of BFGS updates skipped
Nact = number of active bounds at final generalized Cauchy point
Projg = norm of the final projected gradient
F = final function value
* * *
N Tit Tnf Tnint Skip Nact Projg F
1 4 5 4 0 0 1.214D-08 9.258D-15
F = 9.2576608860000000E-015
CONVERGENCE: NORM_OF_PROJECTED_GRADIENT_<=_PGTOL
Optimization Complete
-----------------------------------
Optimized value: [-0.05176027]
real 2m45.615s
user 1m10.384s
sys 0m14.876s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case cart3d_OpenMDAO_3_alpha_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'cart3d_OpenMDAO_3_twist_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
cart3d_OpenMDAO_3_twist_PyTest.py test;
+ python -u cart3d_OpenMDAO_3_twist_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
----------------------------------------
>>> twist [4.]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Getting CFD design variables.......
Number of design variables - 1
Design Variable name - twist
Done getting CFD design variables
Getting CFD functional.......
Number of design variables - 1
Objective name - CL2
Done getting CFD functional
Body size = 3.630958
Tessellating body 1 with MaxEdge = 0.018155 Sag = 0.000908 Angle = 10.000000
Executing: autoInputs -r 5.123450 -nDiv 5 -maxR 7 -symmY -halfBody > autoInputs.out
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:25:59) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:26:30) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.00411365] CL 0.0641354
----------------------------------------
----------------------------------------
>>> twist [4.]
>>> CL2 [0.00411365] CL 0.0641354
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:26:33) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:27:16) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 37528
Estimated cutter memory: 0.003753 GB / DV
Estimated xsensit memory: 0.071303 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:27:16) Linearizing cut-cells
(14:27:18) ... done
o Number of simultaneous xsensit jobs: 1
(14:27:18) Computing gradients
(14:27:22) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[0.00210028]]
========================================
RUNNING THE L-BFGS-B CODE
* * *
Machine precision = 2.220D-16
N = 1 M = 10
At X0 0 variables are exactly at the bounds
At iterate 0 f= 4.11365D-03 |proj g|= 2.10028D-03
----------------------------------------
>>> twist [3.99789972]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.630955
Tessellating body 1 with MaxEdge = 0.018155 Sag = 0.000908 Angle = 10.000000
Executing: autoInputs -r 5.123450 -nDiv 5 -maxR 7 -symmY -halfBody > autoInputs.out
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:27:27) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:28:12) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.00411054] CL 0.0641174
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:28:14) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:29:01) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 37527
Estimated cutter memory: 0.003753 GB / DV
Estimated xsensit memory: 0.071301 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:29:01) Linearizing cut-cells
(14:29:03) ... done
o Number of simultaneous xsensit jobs: 1
(14:29:03) Computing gradients
(14:29:07) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[0.00209952]]
========================================
At iterate 1 f= 4.11054D-03 |proj g|= 2.09952D-03
----------------------------------------
>>> twist [-1.82264925]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.627721
Tessellating body 1 with MaxEdge = 0.018139 Sag = 0.000907 Angle = 10.000000
Executing: autoInputs -r 5.123450 -nDiv 5 -maxR 7 -symmY -halfBody > autoInputs.out
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:29:12) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:29:40) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.00078736] CL -0.0280601
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:29:42) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:30:14) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 37740
Estimated cutter memory: 0.003774 GB / DV
Estimated xsensit memory: 0.071706 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:30:15) Linearizing cut-cells
(14:30:17) ... done
o Number of simultaneous xsensit jobs: 1
(14:30:17) Computing gradients
(14:30:21) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[-0.00094377]]
========================================
At iterate 2 f= 7.87362D-04 |proj g|= 9.43772D-04
----------------------------------------
>>> twist [-0.01760921]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.625633
Tessellating body 1 with MaxEdge = 0.018128 Sag = 0.000906 Angle = 10.000000
Executing: autoInputs -r 5.123450 -nDiv 5 -maxR 7 -symmY -halfBody > autoInputs.out
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:30:25) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:30:53) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [1.54163358e-06] CL 0.00123669
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:30:55) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:31:37) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 38636
Estimated cutter memory: 0.003864 GB / DV
Estimated xsensit memory: 0.073408 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:31:38) Linearizing cut-cells
(14:31:40) ... done
o Number of simultaneous xsensit jobs: 1
(14:31:40) Computing gradients
(14:31:44) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[4.14686939e-05]]
========================================
At iterate 3 f= 1.54163D-06 |proj g|= 4.14687D-05
----------------------------------------
>>> twist [-0.09358318]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.625710
Tessellating body 1 with MaxEdge = 0.018129 Sag = 0.000906 Angle = 10.000000
Executing: autoInputs -r 5.123450 -nDiv 5 -maxR 7 -symmY -halfBody > autoInputs.out
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:31:49) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:32:27) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [6.09330132e-07] CL 0.00078708
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(14:32:29) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(14:33:42) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 37755
Estimated cutter memory: 0.003776 GB / DV
Estimated xsensit memory: 0.071735 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(14:33:43) Linearizing cut-cells
(14:33:45) ... done
o Number of simultaneous xsensit jobs: 1
(14:33:45) Computing gradients
(14:33:49) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[2.58505396e-05]]
========================================
At iterate 4 f= 6.09330D-07 |proj g|= 2.58505D-05
* * *
Tit = total number of iterations
Tnf = total number of function evaluations
Tnint = total number of segments explored during Cauchy searches
Skip = number of BFGS updates skipped
Nact = number of active bounds at final generalized Cauchy point
Projg = norm of the final projected gradient
F = final function value
* * *
N Tit Tnf Tnint Skip Nact Projg F
1 4 5 4 0 0 2.585D-05 6.093D-07
F = 6.0933013220999997E-007
CONVERGENCE: REL_REDUCTION_OF_F_<=_FACTR*EPSMCH
Optimization Complete
-----------------------------------
Optimized values: [-0.09358318]
real 7m56.146s
user 6m41.294s
sys 0m21.916s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case cart3d_OpenMDAO_3_twist_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'fun3d_Morph_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
fun3d_Morph_PyTest.py test;
+ python -u fun3d_Morph_PyTest.py -outLevel=0 -noAnalysis
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1/ESP/DARWINM1/examples/regressionTest/pyCAPSlog.txt
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
AFLR4 : Proximity BG Surface Grid Generation Skipped
AFLR4 : No Modifications Required
Writing TECPLOT file: pyCAPS_AFLR4_AFLR3.dat ....
Finished writing TECPLOT file
Mapping capsMesh attributes ................
Number of unique capsMesh attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 3
Name = Wing1, index = 1
Name = Wing2, index = 2
Name = Farfield, index = 3
Getting volume mesh
AFLR2 : ---------------------------------------
AFLR2 : AFLR2 LIBRARY
AFLR2 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR2 : TRIA/QUAD GRID GENERATOR
AFLR2 : Version Number 9.14.6
AFLR2 : Version Date 11/20/22 @ 09:21PM
AFLR2 : Compile OS Darwin 22.3.0 arm64
AFLR2 : Compile Date 03/08/23 @ 05:04PM
AFLR2 : Copyright 1994-2021, D.L. Marcum
AFLR2 : ---------------------------------------
AFLR4 : ---------------------------------------
AFLR4 : AFLR4 LIBRARY
AFLR4 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR4 : TRIA/QUAD SURFACE GRID GENERATOR
AFLR4 : Version Number 11.4.5
AFLR4 : Version Date 04/20/23 @ 05:37PM
AFLR4 : Compile OS Darwin 22.4.0 arm64
AFLR4 : Compile Date 04/20/23 @ 05:37PM
AFLR4 : Copyright 1994-2021, D.L. Marcum
AFLR4 : ---------------------------------------
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.23
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.23
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
Volume mesh:
Number of nodes = 83060
Number of elements = 484549
Number of triangles = 15680
Number of quadrilatarals = 0
Number of tetrahedrals = 468869
Number of pyramids = 0
Number of prisms = 0
Number of hexahedrals = 0
Writing TECPLOT file: pyCAPS_AFLR4_AFLR3_VolMesh.dat ....
Finished writing TECPLOT file
Getting CFD boundary conditions
Boundary condition name - Farfield
Boundary condition name - Wing1
Boundary condition name - Wing2
Done getting CFD boundary conditions
Getting CFD design variables.......
Number of design variables - 2
Design Variable name - Alpha
Design Variable name - twist
Done getting CFD design variables
Getting CFD functional.......
Number of design variables - 2
Objective name - Cd2
Objective name - Const
Done getting CFD functional
Python library was linked, but will not be used!
Creating FUN3D directory tree
DesignVariable = Alpha
DesignVariable = twist
Writing sensitivity file for body 1, File - Rubberize/model.tec.1.sd1
Writing sensitivity file for body 2, File - Rubberize/model.tec.2.sd1
Writing sensitivity file for body 3, File - Rubberize/model.tec.3.sd1
Writing rubber.data
Writing surface file for 1 (of 3) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 3) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 3) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing MAPBC file ....
Finished writing MAPBC file
Warning: The fun3d.nml file will be overwritten!
Writing fun3d.nml
Info: No recognized data transfer names found.
Reading rubber.data
Projecting tessellation 1 (of 3) on to new body
Projecting tessellation 2 (of 3) on to new body
Projecting tessellation 3 (of 3) on to new body
Python library was linked, but will not be used!
Creating FUN3D directory tree
DesignVariable = Alpha
DesignVariable = twist
Writing sensitivity file for body 1, File - Rubberize/model.tec.1.sd1
Writing sensitivity file for body 2, File - Rubberize/model.tec.2.sd1
Writing sensitivity file for body 3, File - Rubberize/model.tec.3.sd1
Writing surface file for 1 (of 3) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 3) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 3) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing MAPBC file ....
Finished writing MAPBC file
Warning: The fun3d.nml file will be overwritten!
Writing fun3d.nml
Info: No recognized data transfer names found.
Reading rubber.data
real 0m13.697s
user 0m11.788s
sys 0m0.560s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case fun3d_Morph_PyTest.py passed (as expected)
=================================================
*************************************************
*************************************************
=================================================
Did not run examples for:
Fun3D
=================================================
All tests pass!
=================================================
=================================================
[GNU C Compiler (gcc)] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1')
[GNU C Compiler (gcc)] Successfully parsed console log
[GNU C Compiler (gcc)] -> found 0 issues (skipped 0 duplicates)
[GNU C Compiler (gcc)] Skipping post processing
[GNU C Compiler (gcc)] No filter has been set, publishing all 0 issues
[GNU C Compiler (gcc)] Repository miner is not configured, skipping repository mining
[Clang] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1')
[Clang] -> found 0 issues (skipped 0 duplicates)
[Clang] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1')
[Clang] -> found 0 issues (skipped 0 duplicates)
[Clang] Successfully parsed console log
[Clang] -> found 0 issues (skipped 0 duplicates)
[Clang] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1')
[Clang] Skipping post processing
[Clang] No filter has been set, publishing all 0 issues
[Clang] Repository miner is not configured, skipping repository mining
[Static Analysis] Successfully parsed console log
[Static Analysis] -> found 0 issues (skipped 0 duplicates)
[Static Analysis] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1')
[Static Analysis] Successfully parsed console log
[Static Analysis] -> found 0 issues (skipped 0 duplicates)
[Static Analysis] Skipping post processing
[Static Analysis] No filter has been set, publishing all 0 issues
[Static Analysis] Repository miner is not configured, skipping repository mining
[Static Analysis] Successfully parsed console log
[Static Analysis] -> found 0 issues (skipped 0 duplicates)
[Static Analysis] Successfully parsed console log
[Static Analysis] -> found 0 issues (skipped 0 duplicates)
[Static Analysis] Successfully parsed console log
[Static Analysis] -> found 0 issues (skipped 0 duplicates)
[Static Analysis] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWINM1')
[Static Analysis] Skipping post processing
[Static Analysis] No filter has been set, publishing all 0 issues
[Static Analysis] Repository miner is not configured, skipping repository mining
[Static Analysis] Reference build recorder is not configured
[Static Analysis] Obtaining reference build from same job (7.6,CAPS/CFD,DARWINM1,viggen)
[Static Analysis] Using reference build 'ESP_UndefinedCaps/CASREV=7.6,ESPTEST=CAPS/CFD,ESP_ARCH=DARWINM1,buildnode=viggen #179' to compute new, fixed, and outstanding issues
[Static Analysis] Issues delta (vs. reference build): outstanding: 0, new: 0, fixed: 0
[Static Analysis] Evaluating quality gates
[Static Analysis] -> PASSED - Total (any severity): 0 - Quality Gate: 1
[Static Analysis] -> All quality gates have been passed
[Static Analysis] Health report is disabled - skipping
[Static Analysis] Created analysis result for 0 issues (found 0 new issues, fixed 0 issues)
[Static Analysis] Attaching ResultAction with ID 'analysis' to build 'ESP_UndefinedCaps/CASREV=7.6,ESPTEST=CAPS/CFD,ESP_ARCH=DARWINM1,buildnode=viggen #180'.
[Checks API] No suitable checks publisher found.
[PostBuildScript] - [INFO] Executing post build scripts.
[PostBuildScript] - [INFO] build step #0 should only be executed on MATRIX
[WS-CLEANUP] Deleting project workspace...
[WS-CLEANUP] done
Finished: SUCCESS
![[Jenkins]](/csi/static/e031249a/images/svgs/logo.svg){kind=logo}
