View Issue Details

IDProjectCategoryView StatusLast Update
0002824OpenFOAMBugpublic2018-02-07 15:39
Reporterarohner Assigned Tohenry  
Status resolvedResolutionfixed 
PlatformUnixOSWindows Sub (Ubuntu)OS Version16.04.3 LTS
Fixed in Versiondev 
Summary0002824: violation of mass conservation in compressibleInterDyMFoam
DescriptionComparing a simple case, where the system is closed with just a piston compressing the gas inside, shows wrong behavior in pressure and density. They both do not rise as expected due to mass loss in the system. With rhoPimpleDyMFoam I do get the correct behavior.
Steps To ReproduceRunning the attached case, which has been adapted from a previous issue (
on the same topic shows the mentioned behavior. The case has been slightly changed, so that the whole system is closed, hence no inlet/outlet, just a piston compressing the air inside the volume.
Pressure and density do not rise as expected due to mass loss in the system.
TagsNo tags attached.



2018-02-05 08:33

reporter (119,222 bytes)


2018-02-05 08:52

manager   ~0009252

Have you tested OpenFOAM-dev?


2018-02-05 09:04

reporter   ~0009253

Hi Henry,

Thanks for the prompt answer. No, I didn't check it with the -dev version. Have there been any recent updates regarding this solver?


2018-02-05 09:13

manager   ~0009254

Yes there have been changes which might affect this behavior. In anycase all tests and changes will be made to OpenFOAM-dev so it is necessary to first test this version.

> shows wrong behavior in pressure and density.

Could you provide more details? What is wrong about the behavior in pressure and density?


2018-02-05 13:04

reporter   ~0009255

When I evaluate the System after runtime I should get a density of around 1.2kg/m3, from which the solution is still far away with around 1.7kg/m3. This is due to mass loss observed by integrating the density. I wrote this issue here because I couldn't find any problems regarding convergence. Pressure seems to converge fine.


2018-02-05 13:37

manager   ~0009256

I am surprised that mass loss is causing an increase in density, could you clarify?


2018-02-05 13:53

reporter   ~0009257

Nono. To clarify. I observe that density as well as pressure does not increase as expected. I mean the setup can be hand calculated using isotropic compression. Reason for this behavior is that mass conservations is violated by losing mass during the process of compressing. I do just not find a reason - at least not by checking the setup and all the solver settings - why there is this mass loss.


2018-02-05 13:54

reporter   ~0009258

Oh. the number is 1.17kg/m3.


2018-02-05 14:05

manager   ~0009259

Note that compressibleInterDyMFoam operates a corrected incompressible algorithms in order to handle the large phase density differences in two-phase sytems. This means that mass conservation is only obtained by converging the PIMPLE loop sufficiently, proved it does converge. Try tightening tolerances and running more PIMPLE iterations to see if this reduces the convervation error.

rhoPimpleDyMFoam runs a mass conservative compressible algorithm which is appropriate for single-phase, and as such converves mass simply by converging the pressure solution irrespective of convergence of the PIMPLE loop.


2018-02-05 14:56

reporter   ~0009260

Thanks for the insights. This makes things clearer. Unfortunately I had no success with tightening tolerances nor highering PIMPLE interations.


2018-02-05 15:29

manager   ~0009261

I have isolated the problem and will investigate further.


2018-02-07 15:39

manager   ~0009275

Resolved by commit 852726c25fbe6b4bd4e947bec7ff32bc92bb103e

Issue History

Date Modified Username Field Change
2018-02-05 08:33 arohner New Issue
2018-02-05 08:33 arohner File Added:
2018-02-05 08:52 henry Note Added: 0009252
2018-02-05 09:04 arohner Note Added: 0009253
2018-02-05 09:13 henry Note Added: 0009254
2018-02-05 13:04 arohner Note Added: 0009255
2018-02-05 13:37 henry Note Added: 0009256
2018-02-05 13:53 arohner Note Added: 0009257
2018-02-05 13:54 arohner Note Added: 0009258
2018-02-05 14:05 henry Note Added: 0009259
2018-02-05 14:56 arohner Note Added: 0009260
2018-02-05 15:29 henry Note Added: 0009261
2018-02-07 15:39 henry Assigned To => henry
2018-02-07 15:39 henry Status new => resolved
2018-02-07 15:39 henry Resolution open => fixed
2018-02-07 15:39 henry Fixed in Version => dev
2018-02-07 15:39 henry Note Added: 0009275