View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0003561 | OpenFOAM | Bug | public | 2020-09-29 12:29 | 2020-09-29 15:52 |
Reporter | rrossi | Assigned To | henry | ||
Priority | high | Severity | minor | Reproducibility | always |
Status | closed | Resolution | no change required | ||
Platform | Linux 64 bit | OS | Centos | OS Version | 7 |
Summary | 0003561: multiphaseEulerFoam non-zero velocities whre alpha is zero | ||||
Description | I'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. | ||||
Tags | No tags attached. | ||||
|
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? |
|
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. |
|
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. |
|
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). |
|
It is not possible to define the velocity of something that isn't there, zero or otherwise. |
|
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? |
|
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. |
|
User support request. |
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 |