volume mesh generation in pHYP

[bhargavisrinivas666]

I'm trying to generate volume mesh for wing+body using pyHYP .I'm able to run pyhyp script without errors,but the grid is only generating for some distance instead of given march distance. this is the script

from pyhyp import pyHyp

options = {
    # ---------------------------
    #   General options
    # ---------------------------
    "inputFile": "wb_1.6ar_0.67s.cgns",
    "fileType": "CGNS",
    "unattachedEdgesAreSymmetry": False,
    "outerFaceBC": "farfield",
    "autoConnect": True,
    "BC": {1: {"iHigh": "zSymm"},5: {"jLow": "zSymm"},6: {"jHigh": "zSymm"},7: {"jHigh": "zSymm"},8: {"jLow": "zSymm"},9: {"jLow": "zSymm"},10: {"jLow": "zSymm"},11: {"iHigh": "zSymm"},12: {"jHigh": "zSymm"},14: {"iHigh": "zSymm"}, 15: {"iHigh": "zSymm"}, 17: {"iLow": "zSymm","iHigh": "zSymm"},25: {"jLow": "zSymm"},26: {"jLow": "zSymm"},27: {"jHigh": "zSymm"}},
    
    "families": "wall",

    # ---------------------------
    #   Grid Parameters
    # ---------------------------
    "N":129,
    "s0": 1e-6,
    "marchDist" : 22.5,
    "nConstantStart":100,
    
   
    # ---------------------------
    #   Pseudo Grid Parameters
    # ---------------------------
    "ps0": -1.0,
    "pGridRatio":-1.0,
    "cMax":1.0,
   
    # ---------------------------
    #   Smoothing parameters
    # ---------------------------
    "epsE": 3.0,
    "epsI": 6.0,
    "theta": 3.0,
    "volCoef": 0.15,
    "volBlend": 0.0005,
    "volSmoothIter": 100,
    "volSmoothSchedule": [[0, 10], [0.4, 50], [1.0, 100]],
    "slExp":0.073,
  
    # ---------------------------
    #   Solution Parameters
    # ---------------------------
    "kspRelTol": 1e-10,
    "kspMaxIts": 1500,
    "kspSubspaceSize": 120,
    "writeMetrics":True,
}

hyp = pyHyp(options=options)
hyp.run()
hyp.writeCGNS("wb_vol.cgns")

what changes can avoid this.Thanks in advance

[sseraj]

Using "nConstantStart":100 means that the first 100 layers will have a height of 1e-6. This leaves only 28 layers to go from 1e-4 to your marchDist of 22.5. You should try reducing nConstantStart to less than 10 or increase N if you really need 100 constant offwall layers.

[bhargavisrinivas666]

Thank you for the response .If I'm reducing nConstantStart,I'm getting negative volume cells and then segmentation fault error.

[gawng]

I recommend visualizing the mesh to see where the cells are going negative. Then change the options dictionary accordingly.

[bhargavisrinivas666]

When I'm getting segmentation fault error ,I'm unable to visualize it.

[eirikurj]

For future reference, instead of an image, please post the full terminal output as text.
Also please check and make sure that all surface normals are pointing out.

[bhargavisrinivas666]

Yes,I'm sure that all surface normals are pointing out

This is the terminal output for reference ,I want to know why I'm getting NaN ,what it means?

#--------------------#
 Total Nodes:   52982 
 Unique Nodes:  50816 
 Total Faces:   50462 
#--------------------#
 Normal orientation check ...
 Normals are consistent!
 Determining topology ...
 Topology complete.
#--------------------#
Grid Ratio:  1.1251 
#--------------------#
#-------------------------------------------------------------------------------------------------------------------#
# 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.5     2    18      0  0.282  1.01279  0.98644  0.42704  0.277E-14  0.633E-06  0.106E-05  0.0141  0.0000 
      3    1.2     2    18      0  0.300  1.02580  0.98402  0.27494  0.184E-14  0.801E-06  0.241E-05  0.0178  1.7367 
      4    1.5     1    18      0  0.307  1.02896  0.98254  0.24218  0.186E-14  0.901E-06  0.321E-05  0.0200  1.1286 
      5    1.9     2    18      0  0.318  1.04695  0.97349  0.23294  0.235E-14  0.113E-05  0.512E-05  0.0253  1.3145 
      6    2.2     1    18      0  0.322  1.05317  0.97091  0.22859  0.273E-14  0.127E-05  0.625E-05  0.0281  1.1941 
      7    2.4     1    18      0  0.327  1.06609  0.96533  0.22466  0.341E-14  0.142E-05  0.752E-05  0.0315  1.2838 
      8    2.6     1    18      0  0.331  1.08314  0.95858  0.22323  0.423E-14  0.160E-05  0.894E-05  0.0355  1.2710 
      9    2.9     1    18      0  0.335  1.10449  0.95101  0.22432  0.521E-14  0.180E-05  0.105E-04  0.0399  1.2612 
     10    3.1     1    18      0  0.339  1.13143  0.94263  0.22789  0.638E-14  0.203E-05  0.123E-04  0.0448  1.2549 
     11    3.5     2    18      0  0.347  1.21047  0.92344  0.23439  0.777E-14  0.228E-05  0.167E-04  0.0565  1.2541 
     12    3.8     1    19      0  0.350  1.25002  0.91721  0.24607  0.908E-14  0.257E-05  0.189E-04  0.0565  1.1832 
     13    4.0     1    19      0  0.354  1.27928  0.91497  0.25846  0.112E-13  0.289E-05  0.215E-04  0.0634  1.2503 
     14    4.3     1    19      0  0.357  1.35664  0.90399  0.26979  0.137E-13  0.325E-05  0.244E-04  0.0712  1.2563 
     15    4.5     1    19      0  0.360  1.46460  0.89247  0.28320  0.168E-13  0.366E-05  0.276E-04  0.0800  1.2668 
     16    4.8     1    19      0  0.364  1.60803  0.88071  0.29806  0.206E-13  0.411E-05  0.313E-04  0.0898  1.2700 
     17    5.2     2    19      0  0.370  1.92265  0.84861  0.31457  0.253E-13  0.463E-05  0.400E-04  0.1130  1.2644 
     18    5.5     1    19      0  0.374  1.88525  0.83089  0.33402  0.301E-13  0.521E-05  0.447E-04  0.1205  1.1873 
     19    5.7     1    19      0  0.377  1.69282  0.83577  0.35194  0.375E-13  0.586E-05  0.499E-04  0.1466  1.2273 
     20    6.0     1    19      0  0.380  1.63768  0.81812  0.36611  0.458E-13  0.659E-05  0.557E-04  0.1722  1.2286 
     21    6.2     1    19      0  0.383  1.72753  0.80313  0.38015  0.559E-13  0.742E-05  0.623E-04  0.1944  1.2395 
     22    6.5     1    20      0  0.386  1.80817  0.79100  0.39348  0.681E-13  0.834E-05  0.697E-04  0.2225  1.2354 
     23    7.0     2    21      0  0.393  1.85155  0.76576  0.40630  0.826E-13  0.939E-05  0.875E-04  0.3014  1.2297 
     24    7.3     1    20      0  0.396  1.84648  0.75948  0.40456  0.978E-13  0.106E-04  0.969E-04  0.3070  1.1972 
     25    7.6     1    21      0  0.399  1.74095  0.77050  0.38910  0.120E-12  0.119E-04  0.107E-03  0.3562  1.1862 
     26    7.8     1    21      0  0.402  1.76512  0.75959  0.39047  0.144E-12  0.134E-04  0.119E-03  0.4065  1.1925 
     27    8.1     1    22      0  0.405  1.80038  0.74697  0.40632  0.172E-12  0.150E-04  0.133E-03  0.4562  1.2002 
     28    8.4     1    22      0  0.409  1.82579  0.73778  0.41812  0.206E-12  0.169E-04  0.148E-03  0.5065  1.1975 
     29    8.9     2    23      0  0.415  1.87201  0.72086  0.41852  0.244E-12  0.190E-04  0.184E-03  0.6150  1.1964 
     30    9.2     1    22      0  0.418  1.87546  0.72018  0.41904  0.286E-12  0.214E-04  0.203E-03  0.6004  1.1866 
     31    9.5     1    23      0  0.422  1.75480  0.73805  0.41953  0.347E-12  0.241E-04  0.224E-03  0.6567  1.1779 
     32    9.8     1    24      0  0.425  1.77231  0.72813  0.42020  0.411E-12  0.271E-04  0.248E-03  0.7174  1.1910 
     33   10.1     1    24      0  0.428  1.83722  0.72145  0.42100  0.486E-12  0.305E-04  0.275E-03  0.7832  1.1940 
     34   10.4     1    25      0  0.432  1.88163  0.71630  0.42193  0.574E-12  0.343E-04  0.306E-03  0.8549  1.1947 
     35   11.0     2    26      0  0.438  1.77897  0.70537  0.42304  0.677E-12  0.386E-04  0.378E-03  1.0000  1.1959 
     36   11.3     1    26      0  0.442  1.76102  0.70991  0.42433  0.796E-12  0.434E-04  0.417E-03  0.9962  1.1801 
     37   11.6     1    26      0  0.445  1.68920  0.72331  0.42556  0.967E-12  0.447E-04  0.456E-03  1.0000  1.1774 
     38   11.9     1    26      0  0.448  1.60474  0.73463  0.42696  0.118E-11  0.458E-04  0.497E-03  1.0000  1.1774 
     39   12.5     2    25      0  0.453  1.50978  0.75755  0.42855  0.144E-11  0.468E-04  0.580E-03  1.0000  1.1815 
     40   12.8     1    25      0  0.455  1.47435  0.77143  0.43030  0.176E-11  0.470E-04  0.622E-03  1.0000  1.1712 
     41   13.4     2    24      0  0.459  1.43370  0.79217  0.43225  0.217E-11  0.473E-04  0.706E-03  1.0000  1.1837 
     42   14.0     2    24      0  0.463  1.39176  0.81067  0.43440  0.268E-11  0.475E-04  0.790E-03  1.0000  1.1787 
     43   14.8     3    24      0  0.469  1.35170  0.83376  0.43679  0.331E-11  0.478E-04  0.917E-03  1.0000  1.1841 
     44   15.4     2    24      0  0.472  1.32217  0.84657  0.43939  0.410E-11  0.479E-04  0.100E-02  1.0000  1.1786 
     45   16.2     3    23      0  0.476  1.29349  0.86204  0.44198  0.511E-11  0.476E-04  0.113E-02  1.0000  1.1887 
     46   17.3     4    23      0  0.481  1.25817  0.88057  0.39853  0.638E-11  0.471E-04  0.130E-02  1.0000  1.1880 
     47   18.2     3    22      0  0.484  1.23923  0.89047  0.23519  0.797E-11  0.469E-04  0.142E-02  1.0000  1.1982 
     48   19.9     5    22      0  0.489  1.21362  0.90449  0.02483  0.984E-11  0.465E-04  0.163E-02  1.0000  1.2298 
     49   21.3     5    22      1  0.494  1.19397  0.91530 -0.15498  0.118E-10  0.462E-04  0.184E-02  1.0000  1.2313 
#-------------------------------------------------------------------------------------------------------------------#
# Grid | CPU  | Sub | KSP | nAvg |  Sl  | Sensor | Sensor | Min     | Min     |  deltaS  | March    | cMax  | Ratio |
# Lvl  | Time | Its | Its |      |      | Max    | Min    | Quality | Volume  |          | Distance |       | kMax  |
#-------------------------------------------------------------------------------------------------------------------#
     50   22.7     5    22      1  0.498  1.17788  0.92366 -0.21937  0.141E-10  0.460E-04  0.204E-02  1.0000  1.2224 
     51   24.3     6    22      0  0.502  1.16076  0.93226 -0.26442  0.168E-10  0.458E-04  0.229E-02  1.0000  1.1998 
     52   26.2     7    21      0  0.506  1.14624  0.93954 -0.12789  0.200E-10  0.458E-04  0.257E-02  1.0000  1.1999 
     53   28.5     8    21      0  0.511  1.13363  0.94635  0.18639  0.236E-10  0.458E-04  0.290E-02  1.0000  1.2031 
     54   30.8     9    21      0  0.515  1.12386  0.95250  0.19510  0.277E-10  0.458E-04  0.326E-02  1.0000  1.2080 
     55   33.3    10    21      0  0.520  1.11859  0.95778  0.27834  0.325E-10  0.460E-04  0.367E-02  1.0000  1.2368 
     56   36.2    11    21      0  0.524  1.12205  0.96253  0.41333  0.381E-10  0.462E-04  0.412E-02  1.0000  1.2431 
     57   39.5    13    21      0  0.529  1.11063  0.96654  0.43227  0.446E-10  0.468E-04  0.466E-02  1.0000  1.2476 
     58   43.1    14    21      0  0.534  1.10177  0.97013  0.43165  0.521E-10  0.478E-04  0.525E-02  1.0000  1.2407 
     59   46.9    15    21      0  0.539  1.10016  0.97286  0.43070  0.607E-10  0.490E-04  0.589E-02  1.0000  1.2390 
     60   51.6    17    21      0  0.544  1.09894  0.97563  0.42968  0.707E-10  0.504E-04  0.665E-02  1.0000  1.2530 
     61   56.3    18    21      0  0.548  1.09649  0.97776  0.42861  0.826E-10  0.520E-04  0.746E-02  1.0000  1.2701 
     62   61.5    20    21      0  0.553  1.09434  0.97986  0.13110  0.967E-10  0.539E-04  0.841E-02  1.0000  1.2892 
     63   67.1    21    21      0  0.558  1.09270  0.98115 -0.29189  0.113E-09  0.562E-04  0.943E-02  1.0000  1.3020 
     64   73.2    23    21      0  0.563  1.09107  0.98219 -0.14737  0.133E-09  0.591E-04  0.106E-01  1.0000  1.3079 
     65   79.8    25    22      0  0.568  1.08922  0.96995  0.07976  0.156E-09  0.623E-04  0.120E-01  1.0000  1.3170 
     66   86.8    26    22      0  0.573  1.08778  0.98423  0.06832  0.181E-09  0.658E-04  0.134E-01  1.0000  1.3221 
     67   94.3    28    22      0  0.578  1.08173  0.98549  0.05708  0.210E-09  0.672E-04  0.151E-01  1.0000  1.3272 
     68  103.0    32    22      0  0.583  1.07018  0.98721  0.06060  0.244E-09  0.664E-04  0.170E-01  1.0000  1.3218 
     69  112.5    36    21      0  0.588  1.06486  0.98872  0.05004  0.286E-09  0.659E-04  0.191E-01  1.0000  1.3291 
     70  123.8    41    21      0  0.594  1.05371  0.99006 -0.54981  0.336E-09  0.634E-04  0.215E-01  1.0000  1.3330 
     71  136.7    49    21      0  0.599  1.04525  0.98872 -0.43406  0.318E-09  0.598E-04  0.242E-01  1.0000  1.3214 
     72  152.3    59    21      0  0.604  1.03748  0.99244 -0.26974  0.452E-09  0.568E-04  0.272E-01  1.0000  1.3255 
     73  170.2    69    21      0  0.610  1.03063  0.99359 -0.26769  0.513E-09  0.546E-04  0.306E-01  1.0000  1.3146 
     74  191.0    80    21      2  0.615  1.02607  0.99394 -0.48963  0.514E-09  0.521E-04  0.344E-01  1.0000  1.3153 
     75  215.9    95    21      4  0.620  1.02053  0.99426 -0.78555  0.381E-09  0.506E-04  0.388E-01  1.0000  1.3054 
     76  244.6   109    21      5  0.626  1.01794  0.99514 -1.00000  0.858E-10  0.493E-04  0.436E-01  1.0000  1.3010 
     77  277.1   125    21      7  0.632  1.02472  0.99505 -1.00000 -0.305E-09  0.484E-04  0.490E-01  1.0000  1.2951 
     78  314.8   143    22     10  0.637  1.02881  0.99480 -1.00000 -0.906E-09  0.483E-04  0.552E-01  1.0000  1.3625 
     79  358.5   160    23     13  0.643  1.02895  0.99539 -1.00000 -0.196E-08  0.489E-04  0.621E-01  1.0000  1.3957 
     80  413.5   194    24     26  0.648  1.02572  0.99650 -1.00000 -0.464E-08  0.422E-04  0.698E-01  1.0000  1.4385 
     81  494.3   247    31     20  0.654  1.02873  0.93885 -1.00000 -0.303E-07  0.462E-04  0.786E-01  1.0000  1.4697 
     82  568.5   232    23     11  0.660  1.03435  0.96935 -1.00000 -0.205E-07  0.467E-04  0.884E-01  1.0000  1.5696 
     83  647.3   268    22     15  0.666  1.01049  0.97439 -1.00000 -0.106E-06  0.547E-04  0.995E-01  1.0000  7.6808 
     84  710.3   227    22     14  0.672  1.02478  0.94155 -1.00000 -0.103E-05  0.873E-04  0.112E+00  1.0000  2.3729 
     85  750.2   135    27     12  0.678  1.02008  0.97084 -1.00000 -0.456E-06  0.117E-03  0.126E+00  1.0000  8.1356 
     86  909.0    93     0  50804  0.684  0.10000  0.10000      NaN        NaN  0.136E-02  0.143E+00     NaN     Inf 
     87  909.6     8     0  50804  0.690  0.10000  0.10000      NaN        NaN  0.349E-02  0.160E+00     NaN  0.0000 
     88  910.0     5     0  50804  0.697  0.10000  0.10000      NaN        NaN  0.629E-02  0.182E+00     NaN  0.0000 
     89  910.2     3     0  50804  0.703  0.10000  0.10000      NaN        NaN  0.895E-02  0.203E+00     NaN  0.0000 
     90  910.5     3     0  50804  0.710  0.10000  0.10000      NaN        NaN  0.127E-01  0.234E+00     NaN  0.0000 
     91  910.7     2     0  50804  0.716  0.10000  0.10000      NaN        NaN  0.161E-01  0.261E+00     NaN  0.0000 
     92  910.8     2     0  50804  0.722  0.10000  0.10000      NaN        NaN  0.204E-01  0.295E+00     NaN  0.0000 
     93  911.0     2     0  50804  0.730  0.10000  0.10000      NaN        NaN  0.259E-01  0.338E+00     NaN  0.0000 
     94  911.1     1     0  50804  0.734  0.10000  0.10000      NaN        NaN  0.291E-01  0.364E+00     NaN  0.0000 
     95  911.3     2     0  50804  0.743  0.10000  0.10000      NaN        NaN  0.368E-01  0.426E+00     NaN  0.0000 
     96  911.4     1     0  50804  0.747  0.10000  0.10000      NaN        NaN  0.414E-01  0.463E+00     NaN  0.0000 
     97  911.6     2     0  50804  0.757  0.10000  0.10000      NaN        NaN  0.524E-01  0.551E+00     NaN  0.0000 
     98  911.7     1     0  50804  0.762  0.10000  0.10000      NaN        NaN  0.590E-01  0.603E+00     NaN  0.0000 
     99  911.9     1     0  50804  0.768  0.10000  0.10000      NaN        NaN  0.664E-01  0.662E+00     NaN  0.0000 
#-------------------------------------------------------------------------------------------------------------------#
# Grid | CPU  | Sub | KSP | nAvg |  Sl  | Sensor | Sensor | Min     | Min     |  deltaS  | March    | cMax  | Ratio |
# Lvl  | Time | Its | Its |      |      | Max    | Min    | Quality | Volume  |          | Distance |       | kMax  |
#-------------------------------------------------------------------------------------------------------------------#
    100  912.0     2     0  50804  0.779  0.10000  0.10000      NaN        NaN  0.820E-01  0.803E+00     NaN  0.0000 
    101  912.1     1     0  50804  0.785  0.10000  0.10000      NaN        NaN  0.923E-01  0.885E+00     NaN  0.0000 
    102  912.3     1     0  50804  0.790  0.10000  0.10000      NaN        NaN  0.104E+00  0.978E+00     NaN  0.0000 
    103  912.4     1     0  50804  0.796  0.10000  0.10000      NaN        NaN  0.117E+00  0.108E+01     NaN  0.0000 
    104  912.5     1     0  50804  0.803  0.10000  0.10000      NaN        NaN  0.131E+00  0.120E+01     NaN  0.0000 
    105  912.6     1     0  50804  0.809  0.10000  0.10000      NaN        NaN  0.148E+00  0.133E+01     NaN  0.0000 
    106  912.8     2     0  50804  0.822  0.10000  0.10000      NaN        NaN  0.166E+00  0.164E+01     NaN  0.0000 
    107  912.9     1     0  50804  0.828  0.10000  0.10000      NaN        NaN  0.187E+00  0.181E+01     NaN  0.0000 
    108  913.0     1     0  50804  0.834  0.10000  0.10000      NaN        NaN  0.211E+00  0.200E+01     NaN  0.0000 
    109  913.1     1     0  50804  0.840  0.10000  0.10000      NaN        NaN  0.237E+00  0.221E+01     NaN  0.0000 
    110  913.2     1     0  50804  0.847  0.10000  0.10000      NaN        NaN  0.266E+00  0.244E+01     NaN  0.0000 
    111  913.3     1     0  50804  0.853  0.10000  0.10000      NaN        NaN  0.300E+00  0.271E+01     NaN  0.0000 
    112  913.5     2     0  50804  0.867  0.10000  0.10000      NaN        NaN  0.337E+00  0.335E+01     NaN  0.0000 
    113  913.7     1     0  50804  0.873  0.10000  0.10000      NaN        NaN  0.380E+00  0.369E+01     NaN  0.0000 
    114  913.8     1     0  50804  0.880  0.10000  0.10000      NaN        NaN  0.427E+00  0.407E+01     NaN  0.0000 
    115  913.9     1     0  50804  0.886  0.10000  0.10000      NaN        NaN  0.480E+00  0.449E+01     NaN  0.0000 
    116  914.0     1     0  50804  0.893  0.10000  0.10000      NaN        NaN  0.540E+00  0.497E+01     NaN  0.0000 
    117  914.1     1     0  50804  0.900  0.10000  0.10000      NaN        NaN  0.608E+00  0.551E+01     NaN  0.0000 
    118  914.3     2     0  50804  0.914  0.10000  0.10000      NaN        NaN  0.684E+00  0.681E+01     NaN  0.0000 
    119  914.4     1     0  50804  0.921  0.10000  0.10000      NaN        NaN  0.770E+00  0.749E+01     NaN  0.0000 
    120  914.5     1     0  50804  0.928  0.10000  0.10000      NaN        NaN  0.866E+00  0.826E+01     NaN  0.0000 
    121  914.6     1     0  50804  0.935  0.10000  0.10000      NaN        NaN  0.974E+00  0.913E+01     NaN  0.0000 
    122  914.8     1     0  50804  0.942  0.10000  0.10000      NaN        NaN  0.110E+01  0.101E+02     NaN  0.0000 
    123  914.9     1     0  50804  0.949  0.10000  0.10000      NaN        NaN  0.123E+01  0.112E+02     NaN  0.0000 
    124  915.1     2     0  50804  0.964  0.10000  0.10000      NaN        NaN  0.139E+01  0.138E+02     NaN  0.0000 
    125  915.2     1     0  50804  0.971  0.10000  0.10000      NaN        NaN  0.156E+01  0.152E+02     NaN  0.0000 
    126  915.3     1     0  50804  0.978  0.10000  0.10000      NaN        NaN  0.176E+01  0.168E+02     NaN  0.0000 
    127  915.4     1     0  50804  0.986  0.10000  0.10000      NaN        NaN  0.198E+01  0.185E+02     NaN  0.0000 
    128  915.6     1     0  50804  0.993  0.10000  0.10000      NaN        NaN  0.222E+01  0.205E+02     NaN  0.0000 
    129  915.7     1     0  50804  1.001  0.10000  0.10000      NaN        NaN  0.250E+01  0.227E+02     NaN  0.0000 

Even,I'm unable to visualize the output file

[sseraj]

You can use the nTruncate option to stop the extrusion before the NaNs occur but after the negative volumes occur (for example, nTruncate: 83). This will allow you to visualize the negative volumes.

[bhargavisrinivas666]

yes,I'm able to visualize,Now how can I eliminate negative volumes and NaN

…

On Thu, 18 Apr, 2024, 03:22 Sabet Seraj, ***@***.***> wrote: You can use the nTruncate option to stop the extrusion before the NaNs occur but after the negative volumes occur (for example, nTruncate: 83). This will allow you to visualize the negative volumes. — Reply to this email directly, view it on GitHub <#85 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/BCIUJYVUSWJUMI2WNFVVV23Y53VIPAVCNFSM6AAAAABGJ46TKOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANRSGQ3TQOJZGY> . You are receiving this because you authored the thread.Message ID: ***@***.***>

[gawng]

Please attach a picture of the negative volumes

[bhargavisrinivas666]

Usually I check negative volumes in pointwise,But I'm unable to open the output file in pointwise ,I'm attaching paraview results

And I updated few options in pyhyp script and I added both script and terminal output below

rst Imports

from pyhyp import pyHyp

rst general

options = {
# ---------------------------
# General options
# ---------------------------
"inputFile": "wb.cgns",
"fileType": "CGNS",
"unattachedEdgesAreSymmetry": False,
"outerFaceBC": "farfield",
"autoConnect": True,
"BC": {1: {"iHigh": "zSymm"},5: {"jLow": "zSymm"},6: {"jHigh": "zSymm"},7: {"jHigh": "zSymm"},8: {"jLow": "zSymm"},9: {"jLow": "zSymm"},10: {"jLow": "zSymm"},11: {"iHigh": "zSymm"},12: {"jHigh": "zSymm"},14: {"iHigh": "zSymm"}, 15: {"iHigh": "zSymm"}, 17: {"iLow": "zSymm","iHigh": "zSymm"},25: {"jLow": "zSymm"},26: {"jLow": "zSymm"},27: {"jHigh": "zSymm"}},

"families": "wall",
# rst grid
# ---------------------------
#   Grid Parameters
# ---------------------------
"N":129,
"s0": 1e-6,
"marchDist" : 22.5,
"nConstantStart":5,
#"ntruncate":104,
# rst pseudo
# ---------------------------
#   Pseudo Grid Parameters
# ---------------------------
"ps0": -1.0,
"pGridRatio":-1.0,
"cMax":1.0,
# rst smoothing
# ---------------------------
#   Smoothing parameters
# ---------------------------
"epsE": 2.0,
"epsI": 4.0,
"theta": 3.0,
"volCoef": 0.4,
"volBlend": 0.001,
"volSmoothIter": 150,
"volSmoothSchedule": [[0, 150], [0.4, 250],[1.0, 350]],
"slExp":0.073,
"cornerAngle":31.0,
# rst solution
# ---------------------------
#   Solution Parameters
# ---------------------------
"kspRelTol": 1e-10,
"kspMaxIts": 1500,
"kspSubspaceSize": 120,
"writeMetrics":True,

}

rst run pyHyp

hyp = pyHyp(options=options)
hyp.run()
hyp.writeCGNS("wb_vol.cgns")

terminal output
#--------------------#
Total Nodes: 52982
Unique Nodes: 50816
Total Faces: 50462
#--------------------#
Normal orientation check ...
Normals are consistent!
Determining topology ...
Topology complete.
#--------------------#
Grid Ratio: 1.1274
#--------------------#
#-------------------------------------------------------------------------------------------------------------------#

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.7 2 19 0 0.292 1.08310 0.95515 0.41396 0.961E-14 0.636E-06 0.106E-05 0.0141 0.0000
3 1.5 2 19 0 0.310 1.18245 0.93475 0.25903 0.708E-14 0.808E-06 0.242E-05 0.0179 2.1727
4 1.8 1 19 0 0.317 1.22069 0.92891 0.22838 0.801E-14 0.911E-06 0.322E-05 0.0201 1.2763
5 2.2 1 19 0 0.322 1.32381 0.90461 0.22271 0.105E-13 0.100E-05 0.413E-05 0.0227 1.3748
6 2.7 1 19 0 0.327 1.49372 0.87845 0.22815 0.129E-13 0.100E-05 0.513E-05 0.0249 1.3206
7 3.1 1 19 0 0.332 1.74487 0.85507 0.24295 0.176E-13 0.113E-05 0.613E-05 0.0249 1.4396
8 3.7 2 19 0 0.340 2.48272 0.80293 0.25955 0.235E-13 0.127E-05 0.853E-05 0.0315 1.4592
9 4.1 1 19 0 0.343 2.43769 0.77989 0.28559 0.288E-13 0.143E-05 0.980E-05 0.0315 1.2402
10 4.4 1 20 0 0.347 1.94669 0.76497 0.30780 0.385E-13 0.162E-05 0.112E-04 0.0386 1.3395
11 4.7 1 21 0 0.350 1.97694 0.74539 0.32723 0.498E-13 0.182E-05 0.129E-04 0.0485 1.3281
12 5.2 1 21 0 0.354 2.00568 0.73228 0.34658 0.640E-13 0.205E-05 0.147E-04 0.0612 1.3229
13 5.6 1 22 0 0.357 1.98555 0.73468 0.30528 0.814E-13 0.231E-05 0.167E-04 0.0724 1.3041
14 6.7 2 23 0 0.364 1.92879 0.69525 0.30508 0.102E-12 0.261E-05 0.217E-04 0.0927 1.2910
15 7.1 1 23 0 0.367 1.86900 0.69608 0.33669 0.122E-12 0.294E-05 0.243E-04 0.0904 1.2248
16 7.4 1 23 0 0.370 1.73498 0.70777 0.37252 0.154E-12 0.332E-05 0.272E-04 0.0990 1.2304
17 7.8 1 24 0 0.373 1.73735 0.69517 0.40188 0.189E-12 0.374E-05 0.305E-04 0.1080 1.2344
18 8.1 1 24 0 0.376 1.74378 0.67969 0.40173 0.230E-12 0.422E-05 0.343E-04 0.1172 1.2383
19 8.5 1 25 0 0.379 1.75163 0.67062 0.40160 0.278E-12 0.475E-05 0.385E-04 0.1265 1.2371
20 9.2 2 26 0 0.386 1.75274 0.65333 0.40153 0.333E-12 0.536E-05 0.486E-04 0.1456 1.2402
21 9.6 1 26 0 0.389 1.73488 0.66095 0.40168 0.388E-12 0.604E-05 0.539E-04 0.1386 1.2033
22 10.0 1 27 0 0.392 1.64389 0.67982 0.40166 0.476E-12 0.681E-05 0.600E-04 0.1489 1.2319
23 10.4 1 27 0 0.395 1.63475 0.67134 0.40177 0.567E-12 0.768E-05 0.668E-04 0.1600 1.2397
24 10.8 1 28 0 0.398 1.69859 0.66658 0.40194 0.673E-12 0.866E-05 0.745E-04 0.1737 1.2406
25 11.2 1 29 0 0.401 1.72603 0.65927 0.40217 0.797E-12 0.976E-05 0.831E-04 0.1936 1.2427
26 12.0 2 30 0 0.408 1.73111 0.65333 0.40252 0.937E-12 0.110E-04 0.104E-03 0.2393 1.2469
27 12.4 1 29 0 0.411 1.69351 0.65780 0.40315 0.108E-11 0.124E-04 0.115E-03 0.2353 1.2146
28 12.7 1 30 0 0.414 1.63197 0.68303 0.40364 0.132E-11 0.140E-04 0.127E-03 0.2610 1.2321
29 13.1 1 31 0 0.417 1.67010 0.67780 0.40437 0.155E-11 0.158E-04 0.141E-03 0.2889 1.2444
30 13.6 1 32 0 0.420 1.67397 0.67301 0.40526 0.182E-11 0.178E-04 0.157E-03 0.3188 1.2450
31 14.0 1 33 0 0.424 1.69387 0.67169 0.40637 0.214E-11 0.200E-04 0.175E-03 0.3510 1.2442
32 14.9 2 34 0 0.430 1.72031 0.66703 0.40773 0.249E-11 0.226E-04 0.217E-03 0.4211 1.2483
33 15.3 1 33 0 0.434 1.72406 0.67349 0.40948 0.288E-11 0.255E-04 0.240E-03 0.4070 1.2238
34 15.8 1 34 0 0.437 1.68431 0.69660 0.41105 0.347E-11 0.287E-04 0.266E-03 0.4437 1.2300
35 16.2 1 35 0 0.440 1.68539 0.69230 0.41317 0.406E-11 0.324E-04 0.294E-03 0.4818 1.2428
36 16.7 1 36 0 0.443 1.69013 0.69100 0.41559 0.475E-11 0.365E-04 0.327E-03 0.5211 1.2427
37 17.3 1 37 0 0.447 1.74621 0.68928 0.41832 0.556E-11 0.411E-04 0.363E-03 0.5613 1.2419
38 18.3 2 39 0 0.454 1.71623 0.68743 0.42143 0.648E-11 0.464E-04 0.451E-03 0.6443 1.2451
39 18.8 1 37 0 0.457 1.76188 0.69129 0.42491 0.750E-11 0.523E-04 0.497E-03 0.6126 1.2258
40 19.3 1 38 0 0.461 1.70171 0.71277 0.42706 0.905E-11 0.590E-04 0.549E-03 0.6764 1.2251
41 19.9 1 39 0 0.464 1.66669 0.71083 0.42977 0.106E-10 0.665E-04 0.608E-03 0.7460 1.2373
42 20.4 1 40 0 0.468 1.71652 0.71016 0.43264 0.124E-10 0.749E-04 0.675E-03 0.8224 1.2376
43 21.0 1 41 0 0.471 1.75790 0.70718 0.43567 0.145E-10 0.845E-04 0.750E-03 0.9058 1.2365
44 22.0 2 40 0 0.478 1.73283 0.70760 0.43885 0.169E-10 0.952E-04 0.921E-03 1.0000 1.2347
45 22.5 1 39 0 0.482 1.70796 0.72457 0.44141 0.202E-10 0.101E-03 0.101E-02 1.0000 1.2140
46 23.0 1 38 0 0.485 1.65168 0.73880 0.44377 0.246E-10 0.104E-03 0.110E-02 1.0000 1.2144
47 24.0 2 37 0 0.490 1.59562 0.76331 0.44610 0.300E-10 0.110E-03 0.129E-02 1.0000 1.2169
48 24.5 1 36 0 0.493 1.58592 0.77428 0.44843 0.361E-10 0.113E-03 0.139E-02 1.0000 1.2258
49 25.6 2 35 0 0.498 1.49516 0.79124 0.44800 0.440E-10 0.120E-03 0.160E-02 1.0000 1.2244
#-------------------------------------------------------------------------------------------------------------------#

Grid | CPU | Sub | KSP | nAvg | Sl | Sensor | Sensor | Min | Min | deltaS | March | cMax | Ratio |

Lvl | Time | Its | Its | | | Max | Min | Quality | Volume | | Distance | | kMax |

#-------------------------------------------------------------------------------------------------------------------#
50 26.5 2 35 0 0.503 1.51348 0.80561 0.44725 0.534E-10 0.127E-03 0.182E-02 1.0000 1.2250
51 27.4 2 34 0 0.507 1.46600 0.81781 0.44665 0.649E-10 0.134E-03 0.206E-02 1.0000 1.2305
52 28.3 2 33 0 0.511 1.44352 0.82912 0.44612 0.789E-10 0.141E-03 0.231E-02 1.0000 1.2377
53 29.3 2 33 0 0.515 1.44091 0.83739 0.44561 0.961E-10 0.149E-03 0.257E-02 1.0000 1.2491
54 30.2 2 32 0 0.519 1.38611 0.84575 0.44510 0.117E-09 0.157E-03 0.284E-02 1.0000 1.2537
55 31.6 3 31 0 0.525 1.39562 0.85624 0.44461 0.143E-09 0.169E-03 0.328E-02 1.0000 1.2581
56 32.4 2 31 0 0.528 1.35884 0.86185 0.44433 0.175E-09 0.178E-03 0.359E-02 1.0000 1.2636
57 33.8 3 30 0 0.533 1.35691 0.87055 0.44382 0.215E-09 0.191E-03 0.409E-02 1.0000 1.2719
58 35.1 3 30 0 0.538 1.33033 0.87880 0.44324 0.265E-09 0.199E-03 0.461E-02 1.0000 1.2701
59 36.4 3 29 0 0.542 1.31322 0.88801 0.44253 0.327E-09 0.207E-03 0.515E-02 1.0000 1.2725
60 38.2 4 28 0 0.548 1.28292 0.89633 0.44181 0.406E-09 0.217E-03 0.591E-02 1.0000 1.2765
61 39.5 3 28 0 0.552 1.27299 0.90258 0.44120 0.504E-09 0.224E-03 0.650E-02 1.0000 1.2767
62 41.1 4 27 0 0.556 1.24355 0.90864 0.44048 0.629E-09 0.233E-03 0.732E-02 1.0000 1.2764
63 43.0 5 27 0 0.562 1.24151 0.91597 0.43983 0.786E-09 0.245E-03 0.838E-02 1.0000 1.2805
64 44.9 5 27 0 0.567 1.21971 0.92109 0.43918 0.974E-09 0.255E-03 0.949E-02 1.0000 1.2854
65 47.0 5 26 0 0.572 1.19992 0.92633 0.43854 0.117E-08 0.265E-03 0.107E-01 1.0000 1.2998
66 49.0 5 26 0 0.576 1.19714 0.93021 0.43786 0.142E-08 0.274E-03 0.118E-01 1.0000 1.3075
67 51.7 6 26 0 0.581 1.17718 0.93512 0.43717 0.172E-08 0.283E-03 0.133E-01 1.0000 1.3175
68 54.2 7 26 0 0.587 1.15220 0.94075 0.43652 0.208E-08 0.292E-03 0.151E-01 1.0000 1.3232
69 56.8 7 25 0 0.592 1.12584 0.94413 0.43577 0.252E-08 0.297E-03 0.170E-01 1.0000 1.3210
70 59.6 8 25 0 0.597 1.11277 0.94372 0.43493 0.305E-08 0.300E-03 0.191E-01 1.0000 1.3194
71 62.8 9 25 0 0.602 1.14140 0.94086 0.43413 0.371E-08 0.303E-03 0.215E-01 1.0000 1.3121
72 66.7 11 24 0 0.608 1.12666 0.94017 0.43336 0.451E-08 0.303E-03 0.244E-01 1.0000 1.3376
73 70.5 11 24 0 0.613 1.10820 0.94890 0.43260 0.551E-08 0.296E-03 0.274E-01 1.0000 1.3250
74 75.4 14 25 0 0.618 1.08193 0.95863 0.43170 0.669E-08 0.270E-03 0.309E-01 1.0000 1.3029
75 82.0 18 25 0 0.623 1.05759 0.97091 0.39965 0.789E-08 0.219E-03 0.348E-01 1.0000 1.2823
76 91.8 26 25 0 0.629 1.04073 0.97986 0.21860 0.888E-08 0.179E-03 0.392E-01 1.0000 1.2821
77 103.8 31 25 0 0.634 1.03701 0.98259 0.17172 0.971E-08 0.184E-03 0.442E-01 1.0000 1.3007
78 117.2 34 26 0 0.640 1.04096 0.98514 0.16711 0.114E-07 0.187E-03 0.498E-01 1.0000 1.3379
79 132.3 39 25 0 0.646 1.03777 0.98793 0.17319 0.148E-07 0.179E-03 0.561E-01 1.0000 1.3598
80 151.5 44 26 0 0.651 1.03409 0.98865 0.09943 0.166E-07 0.189E-03 0.633E-01 1.0000 1.4019
81 171.8 47 26 0 0.657 1.03105 0.98981 0.09261 0.211E-07 0.191E-03 0.713E-01 1.0000 1.3637
82 193.1 55 25 0 0.663 1.02359 0.99153 0.19563 0.261E-07 0.183E-03 0.804E-01 1.0000 1.3502
83 217.8 64 25 0 0.669 1.02071 0.99273 0.19109 0.338E-07 0.179E-03 0.907E-01 1.0000 1.3408
84 243.7 72 25 0 0.674 1.01663 0.99361 0.18224 0.419E-07 0.181E-03 0.102E+00 1.0000 1.3290
85 273.2 80 24 1 0.680 1.01514 0.99414 0.18072 0.512E-07 0.184E-03 0.115E+00 1.0000 1.3203
86 305.6 89 24 2 0.686 1.01266 0.99462 0.08419 0.594E-07 0.188E-03 0.130E+00 1.0000 1.3126
87 339.6 97 24 2 0.692 1.01179 0.99343 -0.06065 0.647E-07 0.198E-03 0.146E+00 1.0000 1.3051
88 375.7 103 23 2 0.698 1.01046 0.99336 -0.31999 0.696E-07 0.211E-03 0.165E+00 1.0000 1.3005
89 425.2 151 20 0 0.705 1.01679 0.99669 -1.00000 -0.221E-06 0.121E-03 0.186E+00 1.0000 1.2941
90 494.7 227 20 1 0.711 1.00468 0.99566 -1.00000 -0.995E-06 0.123E-03 0.210E+00 1.0000 1.3325
91 559.2 209 20 9 0.717 1.00666 0.96400 -1.00000 -0.165E-05 0.187E-03 0.236E+00 1.0000 1.3350
92 594.9 105 28 3 0.723 1.01145 0.95113 -1.00000 -0.108E-05 0.332E-03 0.266E+00 1.0000 1.3348
93 629.2 96 25 3 0.730 1.05030 0.93498 -1.00000 -0.788E-06 0.586E-03 0.300E+00 1.0000 2.3338
94 657.7 74 29 2 0.736 1.02139 0.97669 -1.00000 -0.701E-06 0.586E-03 0.339E+00 1.0000 1.7638
95 690.5 83 27 4 0.743 1.02468 0.97351 -1.00000 -0.594E-06 0.594E-03 0.382E+00 1.0000 1.6117
96 726.5 90 25 5 0.749 1.03622 0.96696 -1.00000 -0.150E-05 0.619E-03 0.430E+00 1.0000 1.5060
97 766.8 98 27 3 0.756 1.06831 0.98239 -1.00000 -0.712E-05 0.642E-03 0.485E+00 1.0000 2.5406
98 808.6 109 28 5 0.762 1.05518 0.92894 -1.00000 -0.442E-05 0.625E-03 0.547E+00 1.0000 2.0337
99 880.8 163 29 27 0.769 1.20683 0.80634 -1.00000 -0.616E-03 0.464E-03 0.617E+00 1.0000 8.3679
#-------------------------------------------------------------------------------------------------------------------#

Grid | CPU | Sub | KSP | nAvg | Sl | Sensor | Sensor | Min | Min | deltaS | March | cMax | Ratio |

Lvl | Time | Its | Its | | | Max | Min | Quality | Volume | | Distance | | kMax |

#-------------------------------------------------------------------------------------------------------------------#
100 965.6 201 28 9 0.776 1.03060 0.97029 -1.00000 -0.731E-04 0.446E-03 0.695E+00 1.0000 2.7459
101 1059.9 229 27 8 0.783 1.01665 0.97755 -1.00000 -0.561E-04 0.435E-03 0.784E+00 1.0000 3.4234
102 1157.6 252 27 17 0.790 1.01137 0.98835 -1.00000 -0.104E-03 0.465E-03 0.884E+00 1.0000 2.0325
103 1263.9 262 26 17 0.796 1.01741 0.97877 -1.00000 -0.127E-03 0.522E-03 0.996E+00 1.0000 4.0848
104 1356.5 254 25 6 0.803 1.09338 0.95825 -1.00000 -0.961E-04 0.545E-03 0.112E+01 1.0000 2.1616

[bhargavisrinivas666]

Can anyone help me with this?Why I'm getting Negative Volumes and NaN's

…

On Fri, 19 Apr, 2024, 02:02 Galen Ng, ***@***.***> wrote: Send a picture of the negative volumes — Reply to this email directly, view it on GitHub <#85 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/BCIUJYQKIE7UQTC2AOV5U7DY6AUWNAVCNFSM6AAAAABGJ46TKOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANRVGI3DOMJYGE> . You are receiving this because you authored the thread.Message ID: ***@***.***>