search

OpenFOAM / OpenFOAM-10

OpenFOAM Foundation repository for OpenFOAM version 10
Other
74 stars 41 forks source link

cavitationModels in compressibleInterFoam don't conserve mass #7

Open kincaidkc opened 1 year ago

kincaidkc commented 1 year ago

Hello all,

I've been playing with the cavitation models in compressibleInterFoam within a closed domain as a precursor to a more complex project that will be based off of the same mass transfer framework, i.e. source terms for the alpha and p_rgh equations. As a first test for the stock solver I created a small 10cm cube filled partially with liquid water, and partially with water vapor, with the initial pressure set below the cavitation pressure to trigger mass transfer. Using a functionObject I'm tracking the total mass within the domain. At the beginning of the simulation, there is about 0.5 kg total within the closed volume, but after 5s only 0.38 kg remains. I've tried a wide variety of solution parameters to increase accuracy but can't get better than 93% mass conservation (final mass / initial mass) after the system reaches steady state.

I started a thread on CFD Online detailing the steps I took to try to solve the problem, and attached my baseline case, which can be found here: https://www.cfd-online.com/Forums/openfoam-solving/247487-significant-continuity-errors-compressibleinterfoam-cavitation-models.html

I assume the cavitation models were validated before release, but I can't get them to work correctly and haven't received any feedback from the forums. Any help would be appreciated.

Thanks.

kincaidkc commented 1 year ago

Is there any documentation about the implementation of the source terms in the cavitation model? I can't find any references in the usual places (doxygen, header files, Chalmers student papers, etc). Thanks.