mdolab / pyhyp

pyHyp generates volume meshes from surface meshes using hyperbolic marching.
Other
45 stars 38 forks source link

volume mesh generation in pHYP #85

Closed bhargavisrinivas666 closed 4 months ago

bhargavisrinivas666 commented 5 months ago

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 commented 5 months ago

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 commented 5 months ago

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

gawng commented 5 months ago

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

bhargavisrinivas666 commented 5 months ago

When I'm getting segmentation fault error ,I'm unable to visualize it. Screenshot 2024-04-17 114706

eirikurj commented 5 months ago

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 commented 5 months ago

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 commented 5 months ago

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 commented 5 months ago

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 https://github.com/mdolab/pyhyp/issues/85#issuecomment-2062478996, or unsubscribe https://github.com/notifications/unsubscribe-auth/BCIUJYVUSWJUMI2WNFVVV23Y53VIPAVCNFSM6AAAAABGJ46TKOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANRSGQ3TQOJZGY . You are receiving this because you authored the thread.Message ID: @.***>

gawng commented 5 months ago

Please attach a picture of the negative volumes

bhargavisrinivas666 commented 5 months ago

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

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 commented 5 months ago

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 https://github.com/mdolab/pyhyp/issues/85#issuecomment-2065267181, or unsubscribe https://github.com/notifications/unsubscribe-auth/BCIUJYQKIE7UQTC2AOV5U7DY6AUWNAVCNFSM6AAAAABGJ46TKOVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANRVGI3DOMJYGE . You are receiving this because you authored the thread.Message ID: @.***>