64th IUVSTA workshop - arrata and comments
Dear participants of the 64th IUVSTA workshop,
Channel - 2D planarLet me discuss the channel cases first. The molecular species for all benchmark problems is helium. The gas temperature is 300K.
Prof. Sharipov used the hard sphere model while I used the variable hard sphere model (VHS) as a model for the helium molecules. The hard sphere model is the most simple model that is not valid for real gases. The VHS model, developed by Dr. Bird, is more sophisticated as the collision cross section depends on the relative velocity. This diference in the model parameters exlains the deviation in the calculated mas flow rate for benchmark cases 1-15.
The VHS model was also used for the benchmark cases 16-20. After checking the configuration file, I realized that I did not define the geometry correctly. The ratio between channel width and channel diameter was 0.5 instead of 1. The obtained results for a ratio of 1 and molecules modeled as hard spheres is shown in the table below. The computation time is given for the most demanding cases (16-20). The remaning cases can be solved in less time, e.g. all cases with delta=100 can be solved in less time than case 20.
First of all I want to make clear that PI-DSMC produces results that are in very good agreement with the results obtained by Prof. Sharipov.
The results of Prof. Sharipov in turn agree very well with the results presented by Prof. Titarev who used numerical methods to solve the kinetic equation.