View Issue Details
|ID||Project||Category||View Status||Date Submitted||Last Update|
|0002785||OpenFOAM||[All Projects] Bug||public||2017-12-09 18:47||2017-12-13 12:44|
|Fixed in Version||5.x|
|Summary||0002785: compressibleInterFoam issue with transonic option and mass conservation|
|Description||compressibleInterFoam doesn't work with the "transonic" option on. It gives an error message as:|
incompatible dimensions for operation
[p_rgh[1 -3 -1 0 0 0 0] ] + [p_rgh[0 0 -1 0 0 0 0] ]
Also, I have used the solver to model rebound and collapse of a single bubble and I believe that it still doesn't conserve mass accurately (like previous OF releases).
|Tags||No tags attached.|
> compressibleInterFoam doesn't work with the "transonic" option on
The transonic option is out of date and not tested. I will update it now and run the tutorials with it and push into OpenFOAM-dev.
> Also, I have used the solver to model rebound and collapse of a single bubble and I believe that it still doesn't conserve mass accurately
Mass should be conserved if the pressure/phase-fraction loop converges which is not guaranteed. Have you formulated any methods which would guarantee formal mass conservation even if the equation coupling does not converge?
This is now fixed in OpenFOAM-dev: commit 862fa9e10674d63383a03b2499f53d015443c97c
Please test and if it works for your cases I will make the same changes to OpenFOAM-5.x.
Hi Henry, thank you very much. However, I have just tested the updated version with the transonic option (the only difference was in the pressure equation?) and I got this error message:
cannot scale a matrix containing a faceFluxCorrection
Regarding the mass conservation issue, I have noticed that the solver is working much better in OF5 in comparison to the previous version but it needs a very small Co number for challenging problems like bubble collapse. Also, I have been trying to evaluate the solver against the Rayleigh–Plesset equation (for a bubble rebound problem) but it hasn’t been able to predict the maximum rebound radius of the bubble with an acceptable accuracy. I will run the simulation on a larger domain to see if there is any problem with the pressure boundary condition.
||From the error message it appears that you have not recompiled everything affected by the changes in the commit.|
||Yes, I have recompiled everything and it is now working with the transonic option (I think the issue was with fvMatrices).|
BubbleReb.tar.gz (4,682,545 bytes)
Resolved in OpenFOAM-dev by commit 862fa9e10674d63383a03b2499f53d015443c97c
Resolved in OpenFOAM-5.x by commit 73d878406f6bfb877fbff9206aa449b9cb753f58
|2017-12-09 18:47||arahamz||New Issue|
|2017-12-09 20:07||henry||Note Added: 0009142|
|2017-12-09 21:05||henry||Note Added: 0009143|
|2017-12-13 10:10||arahamz||Note Added: 0009150|
|2017-12-13 10:21||henry||Note Added: 0009151|
|2017-12-13 12:21||arahamz||Note Added: 0009152|
|2017-12-13 12:39||arahamz||File Added: BubbleReb.tar.gz|
|2017-12-13 12:44||henry||Assigned To||=> henry|
|2017-12-13 12:44||henry||Status||new => resolved|
|2017-12-13 12:44||henry||Resolution||open => fixed|
|2017-12-13 12:44||henry||Fixed in Version||=> 5.x|
|2017-12-13 12:44||henry||Note Added: 0009154|