search

MicroNanoFlows / OpenFOAM-2.4.0-MNF

OpenFOAM 2.4.0 plus the MicroNanoFlow Group Codes
http://www.micronanoflows.ac.uk/
Other
49 stars 36 forks source link

ab initio potentials based collision submodel #6

Open zhulianhua opened 5 years ago

zhulianhua commented 5 years ago

Hi, I've implemented the ab initio potentials based collision submodel AbInitio. It would be interesting if it can be integrated into the OpenFOAM-MNF. #3

craigwhite87 commented 5 years ago

I have also been working with Prof. Sharipov to introduce the ab initio collision model. Your implementation looks slightly more polished than mine currently. Do you have results comparing the heat flux and shear stress to Prof. Sharipov's Couette and Fourier flow cases? I have achieved agreement to within fractions of a percent of each of his, with the last step to achieving this being the interpolation between sigmaT values. We should also check with Prof. Sharipov that he is happy for it to be included in release of the code along with the look up tables.

zhulianhua commented 5 years ago

Hi Criag, Glad to hear you are already working on that. I've made validation of the He4 Couette flow at delta = 1 at various equilibrium temperature from 1K to 3000K, with one of his paper in 2018. The validation results are given as iPython Notebook. (https://github.com/zhulianhua/AbInitio/blob/master/demo/validate_AbInitio.ipynb). For that case, I didn't use the interpolation of sigmaT, but I made this hidden option (by default disabled) available already in the code (https://github.com/zhulianhua/AbInitio/commit/705c2fb1601eb2d4a55f1fe97de65a22efc1cd54). The relative difference in the non-dimensional shear stress is less than 0.051% with respect to the result in his paper.

zhulianhua commented 5 years ago

One case for the validtion is given here (https://github.com/zhulianhua/AbInitio/tree/master/demo/Couette_He4_delta_1_10K)