View Issue Details

IDProjectCategoryView StatusLast Update
0003561OpenFOAMBugpublic2020-09-29 15:52
Reporterrrossi Assigned Tohenry  
PriorityhighSeverityminorReproducibilityalways
Status closedResolutionno change required 
PlatformLinux 64 bitOSCentosOS Version7
Summary0003561: multiphaseEulerFoam non-zero velocities whre alpha is zero
DescriptionI'm testing the new multiphaseEulerFoam solver and I'm seeing extreme non zero-velocities for the nth-phase in regions where the corresponding volume fraction is zero, causing the time step size of the simulation to decrease considerably due to the Courant number restriction.

Beside my own simulation, I've seen the same behavior in the solver tutorials, for example the bubbleColumn one.

Is this the correct behavior for the solver? Shouldn't the fluid/flow properties of the corresponding phase be zero if the volume fraction is zero?

I guess there might be some stability/numerical needs behind this, but would be great to have some feedback.
TagsNo tags attached.

Activities

henry

2020-09-29 13:18

manager   ~0011546

I don't see a problem with the bubbleColumn, could you be more specific about which bubbleColumn case and at what time you see what level of extreme non zero-velocities?

rrossi

2020-09-29 14:19

reporter   ~0011547

Sorry, I didn't mean there are extreme velocities in the bubbleColumn test case (laminar one) but only non zero-velocities in the region occupied by air where the water fraction is zero (see picture attached).

The extreme values are found in my test case where, of course, they might be due to issues with the setup but I'm wondering if this is the correct behavior of the solver.
pic.png (199,783 bytes)   
pic.png (199,783 bytes)   

henry

2020-09-29 14:49

manager   ~0011548

Why should the water velocity be zero where the water fraction is zero? Formally it is unknown but in the simulations it will approach the velocity of the continuous phase as the phase-fraction -> 0 due to the drag.

rrossi

2020-09-29 14:55

reporter   ~0011549

Multiphase flows are not my main field of expertise, but I was guessing there can't be, physically, any velocity associated with one phase which is not present at that particular location (i.e. the phase/volume fraction is zero).

henry

2020-09-29 15:02

manager   ~0011550

It is not possible to define the velocity of something that isn't there, zero or otherwise.

rrossi

2020-09-29 15:44

reporter   ~0011552

Ok, but wouldn't be zero-padding in cells where the volume fraction is zero within a given threshold be a better option to avoid velocity values higher than physical ones (as in the air region of the picture above where maximum water velocity is about 0.5 where there is no water while the actual maximum value is about 0.3 where the volume fraction is not zero) and relax the Co constraint or would this lead to numerical/stability issues?

henry

2020-09-29 15:52

manager   ~0011553

If you need the phase velocity to be exactly zero when the phase-fraction is zero you will need to add a damping model to do this. e.g. via an fvOption.

henry

2020-09-29 15:52

manager   ~0011554

User support request.

Issue History

Date Modified Username Field Change
2020-09-29 12:29 rrossi New Issue
2020-09-29 13:18 henry Note Added: 0011546
2020-09-29 14:19 rrossi File Added: pic.png
2020-09-29 14:19 rrossi Note Added: 0011547
2020-09-29 14:49 henry Note Added: 0011548
2020-09-29 14:55 rrossi Note Added: 0011549
2020-09-29 15:02 henry Note Added: 0011550
2020-09-29 15:44 rrossi Note Added: 0011552
2020-09-29 15:52 henry Note Added: 0011553
2020-09-29 15:52 henry Assigned To => henry
2020-09-29 15:52 henry Status new => closed
2020-09-29 15:52 henry Resolution open => no change required
2020-09-29 15:52 henry Note Added: 0011554