Segmentation fault
Gittwell opened this issue · comments
Hello. I'm trying to use pyhyp for a complex wing configuration. What could be causing this error: Segmentation fault? Quality mesh more the 0.7. What could be wrong? In what direction should I think? I use docker.
ERROR 7: VSPAERO Viewer Not Found.
Expected here: /home/mdolabuser/packages/OpenVSP/python/openvsp/openvsp/vspviewer
#--------------------#
Total Nodes: 36557
Unique Nodes: 33713
Total Faces: 33536
#--------------------#
Normal orientation check ...
Normals are consistent!
Determining topology ...
Topology complete.
#--------------------#
Grid Ratio: 1.7306
#--------------------#
#-------------------------------------------------------------------------------------------------------------------#
# Grid | CPU | Sub | KSP | nAvg | Sl | Sensor | Sensor | Min | Min | deltaS | March | cMax | Ratio |
# Lvl | Time | Its | Its | | | Max | Min | Quality | Volume | | Distance | | kMax |
#-------------------------------------------------------------------------------------------------------------------#
2 0.3 2 14 61 0.019 1.00001 0.99989 0.03577 0.198E-16 0.100E-08 0.137E-08 0.0000 0.0000
3 0.7 2 14 61 0.023 1.00002 0.99903 0.03552 0.368E-15 0.173E-08 0.410E-08 0.0000 24.5351
4 0.8 1 14 61 0.024 1.00003 0.99840 0.03542 0.119E-14 0.299E-08 0.583E-08 0.0000 3.2697
5 1.1 2 14 61 0.027 1.00005 0.99340 0.03539 0.386E-14 0.518E-08 0.140E-07 0.0000 3.2786
6 1.4 2 14 61 0.030 1.00008 0.98651 0.03536 0.813E-14 0.897E-08 0.282E-07 0.0000 2.1136
7 1.6 1 14 61 0.032 1.00012 0.97773 0.03532 0.171E-13 0.155E-07 0.371E-07 0.0000 2.1076
8 1.9 2 14 61 0.035 1.00020 0.94985 0.03529 0.390E-13 0.269E-07 0.795E-07 0.0001 2.2893
9 2.2 2 14 61 0.039 1.00035 0.90997 0.03525 0.752E-13 0.465E-07 0.153E-06 0.0002 1.9290
10 2.4 1 13 61 0.041 1.00069 0.85486 0.03521 0.146E-12 0.805E-07 0.199E-06 0.0002 1.9409
11 2.7 2 14 61 0.045 1.00202 0.73178 0.03516 0.298E-12 0.139E-06 0.419E-06 0.0007 2.0663
12 3.0 2 14 61 0.050 1.00543 0.57029 0.03511 0.560E-12 0.241E-06 0.799E-06 0.0012 1.8683
13 3.2 1 14 61 0.052 1.01237 0.40303 0.03505 0.104E-11 0.417E-06 0.104E-05 0.0012 1.8839
14 3.5 2 15 61 0.058 1.03184 0.21303 0.03499 0.179E-11 0.722E-06 0.218E-05 0.0035 1.9798
15 3.8 2 17 61 0.064 1.05095 0.11586 0.03493 0.317E-11 0.125E-05 0.415E-05 0.0061 1.8327
16 4.0 1 17 61 0.067 1.06841 0.10000 0.03485 0.537E-11 0.216E-05 0.540E-05 0.0061 1.8399
17 4.4 2 30 61 0.075 1.08002 0.10000 -0.01915 0.913E-11 0.374E-05 0.113E-04 0.0180 1.9021
18 4.9 2 40 61 0.082 1.22213 0.10000 -0.02140 0.185E-10 0.647E-05 0.215E-04 0.0307 1.7985
19 5.2 1 31 61 0.085 5.08805 0.10000 0.03460 0.359E-10 0.112E-04 0.280E-04 0.0304 1.9837
20 6.3 2 92 62 0.095 ******** 0.10000 -1.00000 -0.711E-04 0.194E-04 0.586E-04 0.0977 *******
21 6.6 2 22 110 0.105 ******** 0.10000 -1.00000 -0.105E+01 0.335E-04 0.112E-03 0.1902 *******
22 15.3 1 1500 164 0.109 ******** 0.10000 -1.00000 -0.756E+18 0.581E-04 0.145E-03 0.1370 *******
23 15.5 2 0 33709 0.118 0.10000 0.10000 NaN NaN 0.100E-03 0.239E-03 NaN Inf
24 15.6 2 0 33709 0.132 0.10000 0.10000 NaN NaN 0.174E-03 0.513E-03 NaN 0.0000
25 15.7 2 0 33709 0.146 0.10000 0.10000 NaN NaN 0.301E-03 0.988E-03 NaN 0.0000
26 15.8 1 0 33709 0.152 0.10000 0.10000 NaN NaN 0.521E-03 0.129E-02 NaN 0.0000
27 15.9 2 0 33709 0.170 0.10000 0.10000 NaN NaN 0.901E-03 0.271E-02 NaN 0.0000
28 16.0 2 0 33709 0.187 0.10000 0.10000 NaN NaN 0.156E-02 0.517E-02 NaN 0.0000
29 16.1 1 0 33709 0.194 0.10000 0.10000 NaN NaN 0.270E-02 0.673E-02 NaN 0.0000
30 16.2 2 0 33709 0.217 0.10000 0.10000 NaN NaN 0.467E-02 0.141E-01 NaN 0.0000
31 16.4 2 0 33709 0.239 0.10000 0.10000 NaN NaN 0.808E-02 0.269E-01 NaN 0.0000
32 16.4 1 0 33709 0.249 0.10000 0.10000 NaN NaN 0.140E-01 0.349E-01 NaN 0.0000
33 16.6 2 0 33709 0.278 0.10000 0.10000 NaN NaN 0.242E-01 0.731E-01 NaN 0.0000
34 16.7 2 0 33709 0.306 0.10000 0.10000 NaN NaN 0.419E-01 0.139E+00 NaN 0.0000
35 16.7 1 0 33709 0.319 0.10000 0.10000 NaN NaN 0.725E-01 0.181E+00 NaN 0.0000
36 16.9 2 0 33709 0.356 0.10000 0.10000 NaN NaN 0.125E+00 0.379E+00 NaN 0.0000
37 17.0 2 0 33709 0.392 0.10000 0.10000 NaN NaN 0.217E+00 0.722E+00 NaN 0.0000
38 17.0 1 0 33709 0.408 0.10000 0.10000 NaN NaN 0.376E+00 0.939E+00 NaN 0.0000
39 17.2 2 0 33709 0.456 0.10000 0.10000 NaN NaN 0.650E+00 0.197E+01 NaN 0.0000
40 17.3 2 0 33709 0.502 0.10000 0.10000 NaN NaN 0.113E+01 0.374E+01 NaN 0.0000
41 17.3 1 0 33709 0.522 0.10000 0.10000 NaN NaN 0.195E+01 0.487E+01 NaN 0.0000
42 17.5 2 0 33709 0.583 0.10000 0.10000 NaN NaN 0.337E+01 0.102E+02 NaN 0.0000
43 17.6 2 0 33709 0.642 0.10000 0.10000 NaN NaN 0.583E+01 0.194E+02 NaN 0.0000
44 17.6 1 0 33709 0.668 0.10000 0.10000 NaN NaN 0.101E+02 0.252E+02 NaN 0.0000
45 17.8 2 0 33709 0.746 0.10000 0.10000 NaN NaN 0.175E+02 0.528E+02 NaN 0.0000
46 17.9 2 0 33709 0.822 0.10000 0.10000 NaN NaN 0.302E+02 0.100E+03 NaN 0.0000
47 18.0 1 0 33709 0.855 0.10000 0.10000 NaN NaN 0.523E+02 0.131E+03 NaN 0.0000
48 18.1 2 0 33709 0.955 0.10000 0.10000 NaN NaN 0.906E+02 0.274E+03 NaN 0.0000
49 18.2 2 0 33709 1.052 0.10000 0.10000 NaN NaN 0.157E+03 0.521E+03 NaN 0.0000
Segmentation fault
Please post your pyhyp extrusion options (it should be the python dictionary). We do not have enough information to help you with this. Your output text says your mesh volumes are negative
My pyhyp options.
options = {
# ---------------------------
# General options
# ---------------------------
"inputFile": "wing_p2_2.cgns",
"fileType": "CGNS",
"unattachedEdgesAreSymmetry": True,
"outerFaceBC": "farfield",
"autoConnect": True,
"BC": {},
"families": "wall",
# ---------------------------
# Grid Parameters
# ---------------------------
"N": 49,
"s0": 1e-9,
"marchDist": 23.2 * 16,
# ---------------------------
# Pseudo Grid Parameters
# ---------------------------
"ps0": -1.0,
"pGridRatio": -1.0,
"cMax": 1.0,
# ---------------------------
# Smoothing parameters
# ---------------------------
"epsE": 1.0,
"epsI": 2.0,
"theta": 3.0,
"volCoef": 0.2,
"volBlend": 0.0005,
"volSmoothIter": 20,
# ---------------------------
# Solution Parameters
# ---------------------------
"kspRelTol": 1e-10,
"kspMaxIts": 1500,
"kspSubspaceSize": 50,
}
Great, now I recommend you read the documentation for the smoothing parameters options as that will probably solve your problem. Let us know how it goes.
I would also recommend checking the normal direction of your surface in ICEM. They should be pointing out towards the farfield.