Volume output fields: COORDINATES, SOLUTION, PRIMITIVE
------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Two dimensional problem.
11153 grid points.
11387 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).
1013 triangles.
10374 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 1013 TRIANGLE volume elements.
There has been a re-orientation of 10374 QUADRILATERAL volume elements.
There has been a re-orientation of 545 LINE surface elements.
Identifying edges and vertices.
Setting the control volume structure.
Area of the computational grid: 999.877.
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.1197| 90|
| CV Face Area Aspect Ratio| 1.00061| 32.974|
| CV Sub-Volume Ratio| 1.00002| 6.32859|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 11153| 1/1.00| 10|
| 1| 2707| 1/4.12| 7.38992|
| 2| 555| 1/4.88| 5.01918|
+-------------------------------------------+
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
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 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.7527e-02| 0.070615| 2.738611| 1.108464| 5.599347| 0.058184|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 8.7293e-02| -0.454267| 2.190680| 1.555679| 5.064740| 0.060630|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 8.5950e-02| -0.565341| 2.070468| 1.485403| 4.949698| 0.060586|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 8.4781e-02| -0.726390| 1.920921| 1.383781| 4.787685| 0.072845|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 8.4518e-02| -0.883308| 1.823260| 1.266042| 4.630445| 0.126180|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 8.5095e-02| -1.042499| 1.755083| 1.151079| 4.468758| 0.301945|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 8.5206e-02| -1.231529| 1.700482| 1.105648| 4.268760| 0.849538|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 8.5016e-02| -1.327203| 1.639352| 1.129519| 4.117576| 2.769795|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 8.5026e-02| -1.267223| 1.565550| 1.141885| 4.144438| 46.482609|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 8.5160e-02| -1.220125| 1.493905| 1.121919| 4.191590| -5.798185|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -1.22013| < -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.083415 Cd = -0.014386
Cmz = 0.080305
Cx = -0.014386 Cy = 0.083415
Pressure Contribution
Cl_p = 0.083415 Cd_p = -0.014386
Cmz_p = 0.080305
Cx_p = -0.014386 Cy_p = 0.083415
Viscous Contribution
Cl_v = 0.0 Cd_v = 0.0
Cmz_v = 0.0
Cx_v = 0.0 Cy_v = 0.0
real 0m3.401s
user 0m3.364s
sys 0m0.977s
+ 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/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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_Surf_0.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: AFLR4_Mesh_Surf_1.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: AFLR4_Mesh_Surf_2.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: AFLR4_Mesh_Surf_3.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
Number of surface nodes - 6330
Number of surface elements - 12620
Getting volume mesh
AFLR2 : ---------------------------------------
AFLR2 : AFLR2 LIBRARY
AFLR2 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR2 : TRIA/QUAD GRID GENERATOR
AFLR2 : Version Number 9.14.3
AFLR2 : Version Date 06/09/22 @ 05:16PM
AFLR2 : Compile OS Darwin 19.6.0 x86_64
AFLR2 : Compile Date 07/02/22 @ 03:34PM
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 10.22.23
AFLR4 : Version Date 07/14/22 @ 02:23PM
AFLR4 : Compile OS Darwin 19.6.0 x86_64
AFLR4 : Compile Date 07/18/22 @ 06:42PM
AFLR4 : Copyright 1994-2021, D.L. Marcum
AFLR4 : ---------------------------------------
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.22
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
OVERALL : CPU Time = 0.372 seconds
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.22
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
OVERALL : CPU Time = 0.295 seconds
Volume mesh:
Number of nodes = 39165
Number of elements = 225907
Number of triangles = 12620
Number of quadrilatarals = 0
Number of tetrahedrals = 213287
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.3.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.3.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-2022, 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.
39165 grid points.
213287 volume elements.
4 surface markers.
6726 boundary elements in index 0 (Marker = BC_1).
5418 boundary elements in index 1 (Marker = BC_2).
72 boundary elements in index 2 (Marker = BC_3).
404 boundary elements in index 3 (Marker = BC_4).
213287 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.08442e+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: 4867.4. Mean K: 61.4052. Standard deviation K: 267.141.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 17.0386| 84.7494|
| CV Face Area Aspect Ratio| 1.16458| 28384.3|
| CV Sub-Volume Ratio| 1.03674| 352199|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 39165| 1/1.00| 10|
| 1| 6423| 1/6.1| 8.21053|
| 2| 994| 1/6.46| 6.61219|
+-------------------------------------------+
Finding max control volume width.
Semi-span length = 7.90569 m.
Wetted area = 118.039 m^2.
Area projection in the x-plane = 11.1556 m^2, y-plane = 4.97851 m^2, z-plane = 55.3964 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
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 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.1876e+00| -0.066767| 2.295343| 2.035828| 5.424129| 1.969084|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 1.2125e+00| -0.008267| 2.465034| 2.199175| 5.490089| 2.673815|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 1.2169e+00| 0.028570| 2.552634| 2.439862| 5.520834| 4.098747|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 1.1977e+00| -0.022034| 2.577674| 2.567271| 5.454930| 6.057893|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 1.1864e+00| -0.070810| 2.576919| 2.622447| 5.394307| 7.772831|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 1.1736e+00| -0.087069| 2.577170| 2.646484| 5.390691| 8.984026|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 1.2000e+00| -0.105594| 2.563264| 2.662266| 5.385968| 9.816240|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 1.1930e+00| -0.148979| 2.521215| 2.679520| 5.341336| 10.274841|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 1.1852e+00| -0.205735| 2.453495| 2.672519| 5.270764| 10.337581|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 1.1794e+00| -0.255776| 2.398191| 2.634909| 5.213805| 10.124386|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -0.255776| < -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.793676 Cd = 0.078393
Cmx = -0.000359 Cmy = -0.203012 Cmz = 0.000248
Cx = 0.064529 Cy = 0.000299 Cz = 0.794924
Pressure Contribution
Cl_p = 0.793676 Cd_p = 0.078393
Cmx_p = -0.000359 Cmy_p = -0.203012 Cmz_p = 0.000248
Cx_p = 0.064529 Cy_p = 0.000299 Cz_p = 0.794924
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 0m36.871s
user 0m35.528s
sys 0m1.082s
+ 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/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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.3.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.3.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-2022, 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.
6320 grid points.
12346 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).
12346 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.3051|
| CV Face Area Aspect Ratio| 1.0235| 15.7138|
| CV Sub-Volume Ratio| 1| 7.93643|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 6320| 1/1.00| 10|
| 1| 1790| 1/3.53| 7.98287|
| 2| 504| 1/3.55| 6.35388|
+-------------------------------------------+
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
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 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| 7.7673e-02| -1.163171| 1.073770| 1.341444| 4.312165| 0.062120|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 7.3567e-02| -1.561828| 0.808636| 0.967487| 3.922716| 0.139775|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 7.1549e-02| -1.601698| 0.943793| 1.009036| 3.876382| 0.373708|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 7.0572e-02| -1.659374| 1.008221| 1.015635| 3.799753| 0.952935|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 7.0986e-02| -1.663585| 1.011627| 0.956487| 3.788246| 2.344781|
| 0| 0| 5| 0.0000e+00| 0.0000e+00| 6.9806e-02| -1.633451| 0.991734| 0.872954| 3.829708| 6.259549|
| 0| 0| 6| 0.0000e+00| 0.0000e+00| 6.9528e-02| -1.626398| 0.990123| 0.847075| 3.843511| 26.349556|
| 0| 0| 7| 0.0000e+00| 0.0000e+00| 6.9182e-02| -1.643439| 1.009125| 0.895796| 3.825758| -56.104876|
| 0| 0| 8| 0.0000e+00| 0.0000e+00| 6.8702e-02| -1.659565| 1.020774| 0.939260| 3.804303| -24.431167|
| 0| 0| 9| 0.0000e+00| 0.0000e+00| 6.8472e-02| -1.676757| 1.012633| 0.938320| 3.783310| -21.456835|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 10) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| -1.67676| < -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.805231 Cd = -0.037528
Cmz = 0.414714
Cx = -0.037528 Cy = 0.805231
Pressure Contribution
Cl_p = 0.805231 Cd_p = -0.037528
Cmz_p = 0.414714
Cx_p = -0.037528 Cy_p = 0.805231
Viscous Contribution
Cl_v = 0.0 Cd_v = 0.0
Cmz_v = 0.0
Cx_v = 0.0 Cy_v = 0.0
real 0m2.760s
user 0m2.774s
sys 0m0.822s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_Delaundo_PyTest.py passed (as expected)
=================================================
pointwise: /Applications/Pointwise/PointwiseV18.5R1/pointwise
=================================================
+ echo 'su2_and_pointwise_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWIN64/ESP/DARWIN64/examples/regressionTest/pyCAPSlog.txt
su2_and_pointwise_PyTest.py test;
+ python -u su2_and_pointwise_PyTest.py -outLevel=0
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWIN64/ESP/DARWIN64/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
Writing global Glyph inputs...
Writing egads file '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWIN64/ESP/DARWIN64/examples/pyCAPS/SU2PointwiseAnalysisTest/Scratch/pointwise/caps.egads'....
Number of command line arguments = 2
Sat Nov 12 2022 11:32:47AM EST
Pointwise V18.5
Loading user defaults from capsUserDefaults.glf
GeomToMesh: Defaults
Connector level
InitDim = 10
MaxDim = 15
MinDim = 3
TurnAngle = 0.000000
Deviation = 0.000000
SplitAngle = 0.000000
JoinCons = 0
ProxGrowthRate = 1.300000
SourceSpacing = 0
TurnAngleHard = 70.000000
EqualSpacing = 0
Domain level
Algorithm = Delaunay
FullLayers = 0
MaxLayers = 0
GrowthRate = 1.300000
IsoType = Triangle
TRexType = Triangle
TRexARLimit = 200.000000
TRexAngleBC = 0.000000
Decay = 0.500000
MinEdge = 0.000000
MaxEdge = 0.000000
Adapt = 0
WallSpacing = 0.000000
StrDomConvertARTrigger = 0.000000
Block level
Algorithm = Delaunay
VoxelLayers = 3
boundaryDecay = 0.500000
collisionBuffer = 0.500000
maxSkewAngle = 180.000000
TRexSkewDelay = 0
edgeMaxGrowthRate = 1.800000
fullLayers = 0
maxLayers = 0
growthRate = 1.300000
TRexType = TetPyramid
volInitialize = 1
General
SkipMeshing = 1
CAESolver = UGRID
outerBoxScale = 0
sourceBoxLengthScale = 0.000000
sourceBoxDirection = 1.000000 0.000000 0.000000
sourceBoxAngle = 0.000000
sourceGrowthFactor = 10.000000
sourcePCDFile =
ModelSize = 0
writeGMA = 2.0
assembleTolMult = 2.0
modelOrientIntMeshVolume = 1
Elevate On Export
degree = Q1
costThreshold = 0.8
maxIncAngle = 175.0
relax = 0.05
smoothingPasses = 1000
WCNWeight = 0.5
WCNMode = Calculate
writeVTU = false
CADFile = caps.egads
Determining model size from file.
Info: Model Size adjusted to 100
Computing maximum DB edge tolerance.
Quilt list has 19 entries.
Quilt quilt-1: edge tolerance = 5e-07
Quilt quilt-2: edge tolerance = 6.19984e-05
Quilt quilt-3: edge tolerance = 6.19984e-05
Quilt quilt-4: edge tolerance = 6.19984e-05
Quilt quilt-5: edge tolerance = 6.19984e-05
Quilt quilt-6: edge tolerance = 6.19984e-05
Quilt quilt-7: edge tolerance = 6.19984e-05
Quilt quilt-8: edge tolerance = 6.19984e-05
Quilt quilt-9: edge tolerance = 6.19984e-05
Quilt quilt-10: edge tolerance = 8.32902e-05
Quilt quilt-11: edge tolerance = 8.32902e-05
Quilt quilt-12: edge tolerance = 8.32902e-05
Quilt quilt-13: edge tolerance = 8.32902e-05
Quilt quilt-14: edge tolerance = 8.32902e-05
Quilt quilt-15: edge tolerance = 8.32902e-05
Quilt quilt-16: edge tolerance = 8.32902e-05
Quilt quilt-17: edge tolerance = 8.32902e-05
Quilt quilt-18: edge tolerance = 6.19601e-06
Quilt quilt-19: edge tolerance = 6.19601e-06
Maximum DB edge tolerance = 8.32902e-05
Performing model assembly pass on geometry entities.
domParams(MinEdge) defined as 1.665804e-04 due to model assembly tolerance
Looking for quilt names in geometry.
Quilt list has 19 entries.
Quilt list has 19 entries.
Original connector list has 36 entries.
Minimum connector length = 0.00316783
Tolerance = 8.32902e-05
Performing merge on connectors with equal dimensions and common end points.
After regular merge, connector list has 36 entries.
Updating conData
Min Edge Tol: 5.000e-07 on con-1
Max Edge Tol: 8.329e-05 on con-27
Adjusting connector dimensions based on average spacings.
Checking dimension for connector 1/36, name = con-1
Checking dimension for connector 2/36, name = con-2
Changing dimension for connector 2/36 from 10 to 3
Checking dimension for connector 3/36, name = con-3
Checking dimension for connector 4/36, name = con-4
Changing dimension for connector 4/36 from 10 to 3
Checking dimension for connector 5/36, name = con-5
Changing dimension for connector 5/36 from 10 to 4
Checking dimension for connector 6/36, name = con-6
Checking dimension for connector 7/36, name = con-7
Changing dimension for connector 7/36 from 10 to 3
Checking dimension for connector 8/36, name = con-8
Checking dimension for connector 9/36, name = con-9
Changing dimension for connector 9/36 from 10 to 3
Checking dimension for connector 10/36, name = con-10
Changing dimension for connector 10/36 from 10 to 3
Checking dimension for connector 11/36, name = con-11
Changing dimension for connector 11/36 from 10 to 3
Checking dimension for connector 12/36, name = con-12
Checking dimension for connector 13/36, name = con-13
Changing dimension for connector 13/36 from 10 to 4
Checking dimension for connector 14/36, name = con-14
Checking dimension for connector 15/36, name = con-16
Changing dimension for connector 15/36 from 10 to 3
Checking dimension for connector 16/36, name = con-17
Checking dimension for connector 17/36, name = con-20
Changing dimension for connector 17/36 from 10 to 3
Checking dimension for connector 18/36, name = con-21
Checking dimension for connector 19/36, name = con-24
Changing dimension for connector 19/36 from 10 to 3
Checking dimension for connector 20/36, name = con-25
Changing dimension for connector 20/36 from 10 to 4
Checking dimension for connector 21/36, name = con-26
Checking dimension for connector 22/36, name = con-27
Changing dimension for connector 22/36 from 10 to 3
Checking dimension for connector 23/36, name = con-28
Checking dimension for connector 24/36, name = con-29
Changing dimension for connector 24/36 from 10 to 3
Checking dimension for connector 25/36, name = con-30
Changing dimension for connector 25/36 from 10 to 3
Checking dimension for connector 26/36, name = con-31
Changing dimension for connector 26/36 from 10 to 3
Checking dimension for connector 27/36, name = con-32
Checking dimension for connector 28/36, name = con-33
Changing dimension for connector 28/36 from 10 to 4
Checking dimension for connector 29/36, name = con-34
Checking dimension for connector 30/36, name = con-36
Changing dimension for connector 30/36 from 10 to 3
Checking dimension for connector 31/36, name = con-37
Checking dimension for connector 32/36, name = con-40
Changing dimension for connector 32/36 from 10 to 3
Checking dimension for connector 33/36, name = con-41
Checking dimension for connector 34/36, name = con-44
Changing dimension for connector 34/36 from 10 to 3
Checking dimension for connector 35/36, name = con-45
Checking dimension for connector 36/36, name = con-46
Processing spacing values for connector 1/36, name = con-1
Processing spacing values for connector 2/36, name = con-2
Processing spacing values for connector 3/36, name = con-3
Processing spacing values for connector 4/36, name = con-4
Processing spacing values for connector 5/36, name = con-5
Processing spacing values for connector 6/36, name = con-6
Processing spacing values for connector 7/36, name = con-7
Processing spacing values for connector 8/36, name = con-8
Processing spacing values for connector 9/36, name = con-9
Processing spacing values for connector 10/36, name = con-10
Processing spacing values for connector 11/36, name = con-11
Processing spacing values for connector 12/36, name = con-12
Processing spacing values for connector 13/36, name = con-13
Processing spacing values for connector 14/36, name = con-14
Processing spacing values for connector 15/36, name = con-16
Processing spacing values for connector 16/36, name = con-17
Processing spacing values for connector 17/36, name = con-20
Processing spacing values for connector 18/36, name = con-21
Processing spacing values for connector 19/36, name = con-24
Processing spacing values for connector 20/36, name = con-25
Processing spacing values for connector 21/36, name = con-26
Processing spacing values for connector 22/36, name = con-27
Processing spacing values for connector 23/36, name = con-28
Processing spacing values for connector 24/36, name = con-29
Processing spacing values for connector 25/36, name = con-30
Processing spacing values for connector 26/36, name = con-31
Processing spacing values for connector 27/36, name = con-32
Processing spacing values for connector 28/36, name = con-33
Processing spacing values for connector 29/36, name = con-34
Processing spacing values for connector 30/36, name = con-36
Processing spacing values for connector 31/36, name = con-37
Processing spacing values for connector 32/36, name = con-40
Processing spacing values for connector 33/36, name = con-41
Processing spacing values for connector 34/36, name = con-44
Processing spacing values for connector 35/36, name = con-45
Processing spacing values for connector 36/36, name = con-46
Number of unique endpoints = 24
Checking for baffle intersections.
Performing Local Proximity Refinement
Adjusting connector dimensions based on proximity to other connectors.
Proximity test for connector 1/36, con-1.
Connector dimension changed from 10 to 15, avg. spacing changed from 1.57135 to 1.01015
Proximity test for connector 2/36, con-2.
Connector dimension changed from 3 to 9, avg. spacing changed from 2.5 to 0.625
Proximity test for connector 3/36, con-3.
Connector dimension changed from 10 to 15, avg. spacing changed from 1.57135 to 1.01015
Proximity test for connector 4/36, con-4.
Connector dimension changed from 3 to 9, avg. spacing changed from 2.5 to 0.625
Proximity test for connector 5/36, con-5.
Connector dimension changed from 4 to 15, avg. spacing changed from 1.33869 to 0.286862
Proximity test for connector 6/36, con-6.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.821627 to 0.528189
Proximity test for connector 7/36, con-7.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.985229 to 0.140747
Proximity test for connector 8/36, con-8.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.907322 to 0.583278
Proximity test for connector 9/36, con-9.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.985232 to 0.140747
Proximity test for connector 10/36, con-10.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.985232 to 0.140747
Proximity test for connector 11/36, con-11.
Proximity test for connector 12/36, con-12.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.907322 to 0.583278
Proximity test for connector 13/36, con-13.
Connector dimension changed from 4 to 15, avg. spacing changed from 1.26705 to 0.271511
Proximity test for connector 14/36, con-14.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.821631 to 0.528191
Proximity test for connector 15/36, con-16.
Proximity test for connector 16/36, con-17.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.821627 to 0.528189
Proximity test for connector 17/36, con-20.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.985229 to 0.140747
Proximity test for connector 18/36, con-21.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.821631 to 0.528191
Proximity test for connector 19/36, con-24.
Proximity test for connector 20/36, con-25.
Connector dimension changed from 4 to 15, avg. spacing changed from 0.535476 to 0.114745
Proximity test for connector 21/36, con-26.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.328651 to 0.211275
Proximity test for connector 22/36, con-27.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.394092 to 0.0562988
Proximity test for connector 23/36, con-28.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.362929 to 0.233311
Proximity test for connector 24/36, con-29.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.394096 to 0.0562994
Proximity test for connector 25/36, con-30.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.394096 to 0.0562994
Proximity test for connector 26/36, con-31.
Proximity test for connector 27/36, con-32.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.362929 to 0.233311
Proximity test for connector 28/36, con-33.
Connector dimension changed from 4 to 15, avg. spacing changed from 0.506821 to 0.108604
Proximity test for connector 29/36, con-34.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.328652 to 0.211277
Proximity test for connector 30/36, con-36.
Proximity test for connector 31/36, con-37.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.328651 to 0.211275
Proximity test for connector 32/36, con-40.
Connector dimension changed from 3 to 15, avg. spacing changed from 0.394092 to 0.0562988
Proximity test for connector 33/36, con-41.
Connector dimension changed from 10 to 15, avg. spacing changed from 0.328652 to 0.211277
Proximity test for connector 34/36, con-44.
Proximity test for connector 35/36, con-45.
Connector dimension changed from 10 to 15, avg. spacing changed from 27.9253 to 17.952
Proximity test for connector 36/36, con-46.
Connector dimension changed from 10 to 15, avg. spacing changed from 27.9253 to 17.952
Adjusting node spacing from geometry
Adjusting connector dimensions based on end spacings.
Connector con-1 current dim = 15
Connector con-2 current dim = 9
Connector con-3 current dim = 15
Connector con-4 current dim = 9
Connector con-5 current dim = 15
Connector con-6 current dim = 15
Connector con-7 current dim = 15
Connector con-8 current dim = 15
Connector con-9 current dim = 15
Connector con-10 current dim = 15
Connector con-11 current dim = 3
Connector con-12 current dim = 15
Connector con-13 current dim = 15
Connector con-14 current dim = 15
Connector con-16 current dim = 3
Connector con-17 current dim = 15
Connector con-20 current dim = 15
Connector con-21 current dim = 15
Connector con-24 current dim = 3
Connector con-25 current dim = 15
Connector con-26 current dim = 15
Connector con-27 current dim = 15
Connector con-28 current dim = 15
Connector con-29 current dim = 15
Connector con-30 current dim = 15
Connector con-31 current dim = 3
Connector con-32 current dim = 15
Connector con-33 current dim = 15
Connector con-34 current dim = 15
Connector con-36 current dim = 3
Connector con-37 current dim = 15
Connector con-40 current dim = 15
Connector con-41 current dim = 15
Connector con-44 current dim = 3
Connector con-45 current dim = 15
Connector con-46 current dim = 15
Domain list has 19 entries.
Looking for periodic domain attributes from geometry.
Tolerance = 8.32902e-05
Number of periodic target domains = 0
Unstructured domain 1 of 19:
boundary name = Wake.
Unstructured domain 2 of 19:
boundary name = Wing1.
Creating boundary name Wing1
Unstructured domain 3 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 4 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 5 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 6 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 7 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 8 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 9 of 19:
boundary name = Wing1.
Wing1 already in boundary name list.
Unstructured domain 10 of 19:
boundary name = Wing2.
Creating boundary name Wing2
Unstructured domain 11 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 12 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 13 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 14 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 15 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 16 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 17 of 19:
boundary name = Wing2.
Wing2 already in boundary name list.
Unstructured domain 18 of 19:
boundary name = Farfield.
Creating boundary name Farfield
Unstructured domain 19 of 19:
boundary name = Farfield.
Farfield already in boundary name list.
Checking for blunt flag attribute on geometry faces.
Number of baffles = 1
Adjusting connector dimensions and distributions for wake treatment.
Number of exit connectors = 0
Number of wake edge connectors = 4
Number of TE connectors = 0
Number of exit nodes = 0
Number of TE nodes = 0
New minimum spacing for Wake treatment = 1.0e20
WARNING: T-Rex wall boundary con-1 not found!
WARNING: T-Rex wall boundary con-2 not found!
WARNING: T-Rex wall boundary con-3 not found!
WARNING: T-Rex wall boundary con-4 not found!
Adding Edge connector con-1 to redimension list.
Adding Edge connector con-2 to redimension list.
Adding Edge connector con-3 to redimension list.
Adding Edge connector con-4 to redimension list.
Checking connector con-1
Checking connector con-2
Checking connector con-3
Checking connector con-4
Performing initialization pass on all unstructured domains.
Unstructured domain list for assembly has 19 entries.
Baffle domain dom-1
Number of baffles = 1
Structured domain list for assembly has 0 entries.
Periodic target domain list for assembly has 0 entries.
Creating unstructured block from 18 domains.
Block = auto-blk
Setting up unstructured block.
Retrieving block attributes from geometry.
Initializing unstructured block.
CAE Solver set to UGRID
GeomToMesh finished!
Sat Nov 12 2022 11:32:51AM EST
Run Time: 00 Hours 00 Minutes 04 Seconds
Pointwise save filename = caps.GeomToMesh.pw
Sat Nov 12 2022 11:32:51AM EST
Exporting in UGRID format to file basename = caps.GeomToMesh
Mapped 104 / 104 curves to EGADS originals
GMA contains 36 edge groups
GMA contains 19 face groups
Wrote geometry-mesh associativity to caps.GeomToMesh.gma
GeomToMesh I/O finished!
Sat Nov 12 2022 11:32:54AM EST
Export Time: 00 Hours 00 Minutes 03 Seconds
Volume mesh:
Number of nodes = 16843
Number of elements = 98599
Number of triangles = 2986
Number of quadrilatarals = 0
Number of tetrahedrals = 98599
Number of pyramids = 0
Number of prisms = 0
Number of hexahedrals = 0
Body = 1
Number of nodes = 548
Number of elements = 1092
Number of tris = 1092
Number of quad = 0
Body = 2
Number of nodes = 556
Number of elements = 1108
Number of tris = 1108
Number of quad = 0
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
Body = 3
Number of nodes = 147
Number of elements = 248
Number of tris = 248
Number of quad = 0
Body = 4
Number of nodes = 395
Number of elements = 786
Number of tris = 786
Number of quad = 0
No project name ("Proj_Name") provided - A volume mesh will not be written out
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.3.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.3.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-2022, 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 40 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 14.1421 m.
Surface(s) where the force coefficients are evaluated and
their reference origin for moment computation:
- BC_2 (0, 0, 0).
- BC_3 (0, 0, 0) m.
Surface(s) plotted in the output file: BC_2, BC_3.
Input mesh file name: ../pointwise/caps.GeomToMesh.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|
| 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_Pointwise.dat.
Surface file name: surface_flow_pyCAPS_SU2_Pointwise.
Volume file name: flow_pyCAPS_SU2_Pointwise.
Restart file name: restart_flow_pyCAPS_SU2_Pointwise.dat.
------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
| Marker Type| Marker Name|
+-----------------------------------------------------------------------+
| Euler wall| BC_2|
| | BC_3|
+-----------------------------------------------------------------------+
| Far-field| BC_1|
+-----------------------------------------------------------------------+
-------------------- 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.
16843 grid points.
98599 volume elements.
3 surface markers.
786 boundary elements in index 0 (Marker = BC_1).
1092 boundary elements in index 1 (Marker = BC_2).
1108 boundary elements in index 2 (Marker = BC_3).
98599 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.11302e+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: 6570.86. Mean K: 103.512. Standard deviation K: 506.754.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 42.667| 82.4507|
| CV Face Area Aspect Ratio| 1.20723| 93.3104|
| CV Sub-Volume Ratio| 1.00035| 2062.86|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 16843| 1/1.00| 10|
| 1| 2644| 1/6.37| 8.09168|
| 2| 387| 1/6.83| 6.39656|
+-------------------------------------------+
Finding max control volume width.
Semi-span length = 7.07107 m.
Wetted area = 90.089 m^2.
Area projection in the x-plane = 8.32094 m^2, y-plane = 3.6838 m^2, z-plane = 42.3969 m^2.
Max. coordinate in the x-direction = 12.3619 m, y-direction = 7.07107 m, z-direction = 0.510646 m.
Min. coordinate in the x-direction = 0 m, y-direction = -7.07107 m, z-direction = -0.368435 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker BC_2 is NOT a single plane.
Boundary marker BC_3 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
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 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| 6.2889e-01| 0.431318| 2.847144| 2.407979| 5.922214| 2.162720|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 5.5537e-01| 0.156042| 2.634244| 2.370361| 5.638293| 2.663044|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 5.3154e-01| 0.114372| 2.631966| 2.489176| 5.586483| 3.627433|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 5.2028e-01| 0.053009| 2.612274| 2.549130| 5.521984| 4.729599|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 5.1353e-01| 0.014424| 2.611922| 2.564028| 5.489762| 5.481903|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 5) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| 0.0144238| < -14| No|
+-----------------------------------------------------------------------+
-------------------------------------------------------------------------
+-----------------------------------------------------------------------+
| File Writing Summary | Filename |
+-----------------------------------------------------------------------+
|SU2 binary restart |restart_flow_pyCAPS_SU2_Pointwise.dat|
+-----------------------------------------------------------------------+
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Pointwise.dat).
|CSV file |surface_flow_pyCAPS_SU2_Pointwise.csv|
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Pointwise.dat).
|Tecplot binary |flow_pyCAPS_SU2_Pointwise.szplt |
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Pointwise.dat).
|Tecplot binary surface |surface_flow_pyCAPS_SU2_Pointwise.szplt|
Writing the forces breakdown file (forces_breakdown_pyCAPS_SU2_Pointwise.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.660206 Cd = 0.120434
Cmx = 0.00478 Cmy = -0.197025 Cmz = 0.000155
Cx = 0.108893 Cy = 0.000977 Cz = 0.662207
real 0m15.072s
user 0m13.567s
sys 0m2.113s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_and_pointwise_PyTest.py passed (as expected)
=================================================
=================================================
+ echo 'su2_and_Tetgen_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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 = 267
Number of elements = 530
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 530
Number of quadrilateral elements = 0
----------------------------
Total number of nodes = 3139
Total number of elements = 6266
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 - 3139, elements - 6266
Getting volume mesh
Generating volume mesh using TetGen.....
Tetgen input string = pYq1.500/0.000T1.00e-16A
Delaunizing vertices...
Delaunay seconds: 0.048695
Creating surface mesh ...
Surface mesh seconds: 0.012742
Recovering boundaries...
Boundary recovery seconds: 0.122433
Removing exterior tetrahedra ...
Spreading region attributes.
Exterior tets removal seconds: 0.013819
Suppressing Steiner points ...
Steiner suppression seconds: 1.5e-05
Recovering Delaunayness...
Delaunay recovery seconds: 0.015093
Refining mesh...
4189 insertions, added 2138 points, 127097 tetrahedra in queue.
1394 insertions, added 409 points, 130640 tetrahedra in queue.
1859 insertions, added 449 points, 127916 tetrahedra in queue.
2478 insertions, added 473 points, 111923 tetrahedra in queue.
3303 insertions, added 484 points, 65315 tetrahedra in queue.
4403 insertions, added 433 points, 571 tetrahedra in queue.
Refinement seconds: 1.25713
Smoothing vertices...
Mesh smoothing seconds: 0.84849
Improving mesh...
Mesh improvement seconds: 0.039707
Jettisoning redundant points.
Writing nodes.
Writing elements.
Writing faces.
Writing edges.
Output seconds: 0.004877
Total running seconds: 2.36346
Statistics:
Input points: 3139
Input facets: 6266
Input segments: 9399
Input holes: 2
Input regions: 0
Mesh points: 8039
Mesh tetrahedra: 39668
Mesh faces: 82470
Mesh faces on exterior boundary: 6268
Mesh faces on input facets: 6268
Mesh edges on input segments: 9400
Steiner points on input segments: 1
Steiner points inside domain: 4899
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.3.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.3.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-2022, 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.
8039 grid points.
39668 volume elements.
3 surface markers.
4468 boundary elements in index 0 (Marker = BC_1).
1270 boundary elements in index 1 (Marker = BC_2).
530 boundary elements in index 2 (Marker = BC_3).
39668 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.08934e+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.561. Standard deviation K: 37.5336.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 1.83144| 84.8022|
| CV Face Area Aspect Ratio| 1.55015| 7101.39|
| CV Sub-Volume Ratio| 1.05639| 198266|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 8039| 1/1.00| 10|
| 1| 1338| 1/6.01| 8.25105|
| 2| 214| 1/6.25| 6.7182|
+-------------------------------------------+
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.58502 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
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 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.0699e-01| 0.619384| 3.051344| 2.094590| 6.110279| 1.812161|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 2.0584e-01| 0.785542| 3.264371| 2.974872| 6.285549| 2.181906|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 2.0392e-01| 0.724139| 3.334285| 3.249736| 6.229534| 3.226648|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 2.0709e-01| 0.621976| 3.305627| 3.334988| 6.113861| 4.648342|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 2.2073e-01| 0.625390| 3.207240| 3.303991| 6.098630| 5.785533|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 5) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| 0.62539| < -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.759441 Cd = 0.131266
Cmx = 0.001646 Cmy = -0.215157 Cmz = 0.000487
Cx = 0.117991 Cy = 0.00137 Cz = 0.761616
real 3m32.285s
user 3m31.414s
sys 0m1.079s
+ 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/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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 = 3808
Number of elements = 7616
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 7616
Number of quadrilateral elements = 0
Body 2 (of 2)
Number of nodes = 5654
Number of elements = 11304
Number of node elements = 0
Number of line elements = 0
Number of triangle elements = 11304
Number of quadrilateral elements = 0
----------------------------
Total number of nodes = 9462
Total number of elements = 18920
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 - 9462, elements - 18920
Getting volume mesh
Generating volume mesh using TetGen.....
Tetgen input string = pYq1.500/0.000T1.00e-16A
Delaunizing vertices...
Delaunay seconds: 0.146408
Creating surface mesh ...
Surface mesh seconds: 0.03721
Recovering boundaries...
Warning: 6 segments are not recovered.
Boundary recovery seconds: 0.317661
Removing exterior tetrahedra ...
Spreading region attributes.
Exterior tets removal seconds: 0.028852
Suppressing Steiner points ...
Steiner suppression seconds: 1.5e-05
Recovering Delaunayness...
Delaunay recovery seconds: 0.039834
Refining mesh...
12618 insertions, added 5266 points, 45 tetrahedra in queue.
Refinement seconds: 1.26823
Smoothing vertices...
Mesh smoothing seconds: 1.14649
Improving mesh...
Mesh improvement seconds: 0.096275
Jettisoning redundant points.
Writing nodes.
Writing elements.
Writing faces.
Writing edges.
Output seconds: 0.019246
Total running seconds: 3.11064
Statistics:
Input points: 9462
Input facets: 18920
Input segments: 28380
Input holes: 1
Input regions: 0
Mesh points: 14871
Mesh tetrahedra: 61451
Mesh faces: 132365
Mesh faces on exterior boundary: 18926
Mesh faces on input facets: 18926
Mesh edges on input segments: 28383
Steiner points on input segments: 3
Steiner points inside domain: 5406
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 - x43A
Done getting CFD boundary conditions
Warning: The su2 cfg file will be overwritten!
Write SU2 configuration file for version "BlackBird (7.3.1) "
Running SU2......
-------------------------------------------------------------------------
| ___ _ _ ___ |
| / __| | | |_ ) Release 7.3.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-2022, 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.
14871 grid points.
61451 volume elements.
2 surface markers.
7622 boundary elements in index 0 (Marker = BC_1).
11304 boundary elements in index 1 (Marker = BC_2).
61451 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.67744e+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.525068. Standard deviation K: 2.57981.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
| Mesh Quality Metric| Minimum| Maximum|
+--------------------------------------------------------------+
| Orthogonality Angle (deg.)| 2.06744| 83.7887|
| CV Face Area Aspect Ratio| 1.43723| 2173.76|
| CV Sub-Volume Ratio| 1| 173135|
+--------------------------------------------------------------+
Setting the multigrid structure.
+-------------------------------------------+
| MG Level| CVs|Aggl. Rate| CFL|
+-------------------------------------------+
| 0| 14871| 1/1.00| 10|
| 1| 2668| 1/5.57| 8.46005|
| 2| 485| 1/5.5| 7.18869|
+-------------------------------------------+
Finding max control volume width.
Semi-span length = 30 m.
Wetted area = 13597.1 m^2.
Area projection in the x-plane = 761.677 m^2, y-plane = 2062.72 m^2, z-plane = 5463.61 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
+----------------------------------------------------------------------------------------------------------------------------------------------+
| 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.9618e-01| 2.731069| 5.665107| 4.728786| 8.730545| -0.436821|
| 0| 0| 1| 0.0000e+00| 0.0000e+00| 3.8640e-01| 3.111151| 6.144144| 5.696527| 9.110995| -0.142019|
| 0| 0| 2| 0.0000e+00| 0.0000e+00| 3.8241e-01| 3.211606| 6.272645| 5.788253| 9.178383| -0.024212|
| 0| 0| 3| 0.0000e+00| 0.0000e+00| 3.7760e-01| 3.212276| 6.263360| 5.685266| 9.173095| 0.000170|
| 0| 0| 4| 0.0000e+00| 0.0000e+00| 3.7525e-01| 3.198810| 6.249774| 5.549315| 9.170668| -0.026654|
----------------------------- Solver Exit -------------------------------
Maximum number of iterations reached (ITER = 5) before convergence.
+-----------------------------------------------------------------------+
| Convergence Field | Value | Criterion | Converged |
+-----------------------------------------------------------------------+
| rms[Rho]| 3.19881| < -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 135 non-physical points in the solution.
Warning: 506 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 = -3.127208 Cd = 117.324459
Cmx = -34.067503 Cmy = 9269.22877 Cmz = 116.064843
Cx = 117.327335 Cy = 1.889135 Cz = 3.017366
real 0m23.243s
user 0m37.822s
sys 0m1.082s
+ status=0
+ set +x
=================================================
CAPS verification (using pyCAPS) case su2_X43a_PyTest.py passed (as expected)
=================================================
./execute_PyTestRegression.sh: line 403: ulimit: stack size: cannot modify limit: Operation not permitted
flowCart: /Users/jenkins/util/cart3d/cart3d_v1.5.7_OSX64_ICC--20.10.28/bin/OSX64_ICC//flowCart
=================================================
+ echo 'cart3d_PyTest.py test;'
+ tee -a /Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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 0m22.347s
user 0m37.276s
sys 0m0.854s
+ 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/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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
(11:38:09) Running aero.csh on 2 core(s) in M0.5901A2B0_DP1
(11:38:13) 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
(11:38:14) Running aero.csh on 2 core(s) in M0.5901A2B0_DP1
(11:38:17) 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
(11:38:18) Linearizing cut-cells
(11:38:18) ... done
o Number of simultaneous xsensit jobs: 1
(11:38:18) Computing gradients
(11:38:22) ... 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
(11:38:23) Running aero.csh on 2 core(s) in M0.5901A1.98427B0_DP1
(11:38:25) 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
(11:38:26) Running aero.csh on 2 core(s) in M0.5901A1.98427B0_DP1
(11:38:30) 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
(11:38:30) Linearizing cut-cells
(11:38:30) ... done
o Number of simultaneous xsensit jobs: 1
(11:38:30) Computing gradients
(11:38:34) ... 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
(11:38:35) Running aero.csh on 2 core(s) in M0.5901A-0.0109918B0_DP1
(11:38:37) aero.csh done
o case M0.5901A-0.0109918B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [6.05653687e-06] CL 0.00246101
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:38:39) Running aero.csh on 2 core(s) in M0.5901A-0.0109918B0_DP1
(11:38:42) 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
(11:38:42) Linearizing cut-cells
(11:38:42) ... done
o Number of simultaneous xsensit jobs: 1
(11:38:42) Computing gradients
(11:38:46) ... done
o case M0.5901A-0.0109918B0_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
(11:38:47) Running aero.csh on 2 core(s) in M0.5901A-0.0514929B0_DP1
(11:38:49) aero.csh done
o case M0.5901A-0.0514929B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [2.59817201e-10] CL 1.612e-05
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:38:50) Running aero.csh on 2 core(s) in M0.5901A-0.0514929B0_DP1
(11:38: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
(11:38:54) Linearizing cut-cells
(11:38:54) ... done
o Number of simultaneous xsensit jobs: 1
(11:38:54) Computing gradients
(11:38:58) ... done
o case M0.5901A-0.0514929B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[2.03612486e-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
(11:38:59) Running aero.csh on 2 core(s) in M0.5901A-0.0517603B0_DP1
(11:39:01) aero.csh done
o case M0.5901A-0.0517603B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [9.26120217e-15] CL 1e-07
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:39:02) Running aero.csh on 2 core(s) in M0.5901A-0.0517603B0_DP1
(11:39:05) 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
(11:39:06) Linearizing cut-cells
(11:39:06) ... done
o Number of simultaneous xsensit jobs: 1
(11:39:06) Computing gradients
(11:39:10) ... done
o case M0.5901A-0.0517603B0_DP1, averaging data over 1 cycle(s)
>>> CL2_alpha [[1.21394333e-08]]
========================================
At iterate 4 f= 9.26120D-15 |proj g|= 1.21394D-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.261D-15
F = 9.2612021713000013E-015
CONVERGENCE: NORM_OF_PROJECTED_GRADIENT_<=_PGTOL
Optimization Complete
-----------------------------------
Optimized value: [-0.05176027]
real 1m52.467s
user 0m52.352s
sys 0m9.295s
+ 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/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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
(11:39:16) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:39:34) 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
(11:39:37) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:40:03) 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
(11:40:03) Linearizing cut-cells
(11:40:05) ... done
o Number of simultaneous xsensit jobs: 1
(11:40:05) Computing gradients
(11:40:09) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[0.00209615]]
========================================
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.09615D-03
----------------------------------------
>>> twist [3.99790385]
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
(11:40:14) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:40:31) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.00411103] CL 0.0641104
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:40:34) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:41:00) 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
(11:41:01) Linearizing cut-cells
(11:41:03) ... done
o Number of simultaneous xsensit jobs: 1
(11:41:03) Computing gradients
(11:41:07) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[0.00209539]]
========================================
At iterate 1 f= 4.11103D-03 |proj g|= 2.09539D-03
----------------------------------------
>>> twist [-1.83549221]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.627738
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
(11:41:11) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:41:28) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [0.00079649] CL -0.0282221
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:41:31) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:41:58) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 37721
Estimated cutter memory: 0.003772 GB / DV
Estimated xsensit memory: 0.071670 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(11:41:58) Linearizing cut-cells
(11:42:00) ... done
o Number of simultaneous xsensit jobs: 1
(11:42:00) Computing gradients
(11:42:04) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[-0.00094963]]
========================================
At iterate 2 f= 7.96486D-04 |proj g|= 9.49631D-04
----------------------------------------
>>> twist [-0.01627177]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.625631
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
(11:42:08) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:42:25) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [1.46441729e-06] CL 0.00122349
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:42:29) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:42:55) 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
(11:42:55) Linearizing cut-cells
(11:42:57) ... done
o Number of simultaneous xsensit jobs: 1
(11:42:57) Computing gradients
(11:43:01) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[4.03789358e-05]]
========================================
At iterate 3 f= 1.46442D-06 |proj g|= 4.03789D-05
----------------------------------------
>>> twist [-0.09047123]
Mapping capsGroup attributes ................
Number of unique capsGroup attributes = 1
Name = Wing, index = 1
Body size = 3.625707
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
(11:43:06) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:43:24) aero.csh done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2 [7.13761415e-07] CL 0.00085978
----------------------------------------
========================================
o Evaluating Cart3D functionals and/or gradients
ATTENTION: using static geometry ../inputs/Components.i.tri
(11:43:27) Running aero.csh on 2 core(s) in M0.5901A0B0_DP1
(11:43:53) aero.csh done
Max memory allowed: 32 GB
Flow hexes: 37754
Estimated cutter memory: 0.003775 GB / DV
Estimated xsensit memory: 0.071733 GB / DV
Available threads: 2
# DVs: 1
o Number of simultaneous cutter jobs: 1
(11:43:53) Linearizing cut-cells
(11:43:55) ... done
o Number of simultaneous xsensit jobs: 1
(11:43:55) Computing gradients
(11:43:59) ... done
o case M0.5901A0B0_DP1, averaging data over 1 cycle(s)
>>> CL2_twist [[2.79531675e-05]]
========================================
At iterate 4 f= 7.13761D-07 |proj g|= 2.79532D-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.795D-05 7.138D-07
F = 7.1376141541999997E-007
CONVERGENCE: REL_REDUCTION_OF_F_<=_FACTR*EPSMCH
Optimization Complete
-----------------------------------
Optimized values: [-0.09047123]
real 4m49.641s
user 7m9.456s
sys 0m11.148s
+ 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/DARWIN64/ESP/DARWIN64/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/DARWIN64/ESP/DARWIN64/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: pyCAPS_AFLR4_AFLR3_Surf_0.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: pyCAPS_AFLR4_AFLR3_Surf_1.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: pyCAPS_AFLR4_AFLR3_Surf_2.dat ....
Finished writing TECPLOT file
Writing TECPLOT file: pyCAPS_AFLR4_AFLR3_Surf_3.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
Number of surface nodes - 7964
Number of surface elements - 15886
Getting volume mesh
AFLR2 : ---------------------------------------
AFLR2 : AFLR2 LIBRARY
AFLR2 : ADVANCING-FRONT/LOCAL-RECONNECTION
AFLR2 : TRIA/QUAD GRID GENERATOR
AFLR2 : Version Number 9.14.3
AFLR2 : Version Date 06/09/22 @ 05:16PM
AFLR2 : Compile OS Darwin 19.6.0 x86_64
AFLR2 : Compile Date 07/02/22 @ 03:34PM
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 10.22.23
AFLR4 : Version Date 07/14/22 @ 02:23PM
AFLR4 : Compile OS Darwin 19.6.0 x86_64
AFLR4 : Compile Date 07/18/22 @ 06:42PM
AFLR4 : Copyright 1994-2021, D.L. Marcum
AFLR4 : ---------------------------------------
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.22
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
OVERALL : CPU Time = 0.455 seconds
EGADS : ---------------------------------------
EGADS : Engineering Geometry Aircraft Design System
EGADS : Version 1.22
EGADS : OpenCASCADE Version 7.6.0
EGADS : ---------------------------------------
OVERALL : CPU Time = 0.311 seconds
Volume mesh:
Number of nodes = 57557
Number of elements = 333420
Number of triangles = 15886
Number of quadrilatarals = 0
Number of tetrahedrals = 317534
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
Warning: No boundary condition specified for capsGroup Wake!
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 sensitivity file for body 4, File - Rubberize/model.tec.4.sd1
Writing rubber.data
Writing surface file for 1 (of 4) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 4) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 4) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 4) 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 4) on to new body
Projecting tessellation 2 (of 4) on to new body
Projecting tessellation 3 (of 4) on to new body
Projecting tessellation 4 (of 4) 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 sensitivity file for body 4, File - Rubberize/model.tec.4.sd1
Writing surface file for 1 (of 4) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 4) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 4) in body 1 File - Flow/fun3dAFLRTest_body1.dat
Writing surface file for 1 (of 4) 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 0m11.841s
user 0m10.888s
sys 0m0.817s
+ 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/DARWIN64')
[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/DARWIN64')
[Clang] -> found 0 issues (skipped 0 duplicates)
[Clang] Parsing console log (workspace: '/Users/jenkins/workspace/ESP_UndefinedCaps/7.6/CAPS/CFD/DARWIN64')
[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/DARWIN64')
[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] Reference build recorder is not configured
[Static Analysis] Obtaining reference build from same job (7.6,CAPS/CFD,DARWIN64,macys)
[Static Analysis] Using reference build 'ESP_UndefinedCaps/CASREV=7.6,ESPTEST=CAPS/CFD,ESP_ARCH=DARWIN64,buildnode=macys #171' 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=DARWIN64,buildnode=macys #172'.
[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/d4a6861f/images/svgs/logo.svg){kind=logo}
