View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0000156 | OpenFOAM | Bug | public | 2011-03-04 19:42 | 2011-03-04 22:50 |
Reporter | Assigned To | henry | |||
Priority | normal | Severity | minor | Reproducibility | always |
Status | closed | Resolution | fixed | ||
Platform | Linux | OS | Ubuntu | OS Version | 10.10 |
Summary | 0000156: High viscosity gives wrong result in multiphase problem | ||||
Description | The InterFoam solver seems to be giving the wrong answer when the viscosity is large. To demonstrate this I generated a test case called waterDrop which I uploaded as a zip file. It is essentially identical to the damBreak tutorial except that the liquid is now in a rectangular blob suspended above the ground by about 0.218 meters. At t=0 it is dropped and falls down. The expected fall time can be obtained from the usual physics formula: t = sqrt(2*h/g), which in this case is about 0.21 seconds. In the simulation I varied the viscosity of the water, not of the air. It works as expected for low viscosity, but surprisingly the fall time is found to increase when the water viscosity is increased to a large value. As uploaded, the viscosity of the "water" has been increased from 1e-6 to 1e0 and the fall time has increased to about 0.78 seconds. When further increased to 1e2 the water actually floats upwards and disappears out the atmosphere patch! | ||||
Steps To Reproduce | The viscosity nu for phase1 is set in the file constant/transportProperties. To run and generate log file, enter the commands: blockMesh, checkMesh, setFields, interFoam | tee log1. | ||||
Additional Information | I have tried adjusting most of the parameters in the fvSolution file with little success. Setting momentumPredictor to yes helps somewhat as does increasing nCorrectors but drastically increases the run time and doesn't completely fix the problem. Decreasing Courant number in the controlDict file also has some beneficial effect. I'm hoping that someone will know how to solve this problem, or else can confirm my suspicion that this is a symptom of some kind of bug in the software. Thanks - Paul | ||||
Tags | interFoam, multiphase, viscosity | ||||
2011-03-04 19:42
|
|
|
Increasing the viscosity ratio increases the stiffness of the problem and in order to converge the velocity in high viscosity region which is accelerating will indeed require an increase in the number of PISO iterations: 10 to 20 are required for a liquid viscosity of 0.1. With a viscosity of 1e2 you may require hundreds of PISO iterations or it may not even converge and you may have to change to PIMPLE and apply under-relaxation, either way it would be EXTREMELY expensive and it is not at all clear that VoF would be appropriate anyway. |
Date Modified | Username | Field | Change |
---|---|---|---|
2011-03-04 19:42 |
|
New Issue | |
2011-03-04 19:42 |
|
File Added: waterDrop.zip | |
2011-03-04 19:47 |
|
Tag Attached: multiphase | |
2011-03-04 19:47 |
|
Tag Attached: viscosity | |
2011-03-04 19:47 |
|
Tag Attached: interFoam | |
2011-03-04 22:50 | henry | Note Added: 0000276 | |
2011-03-04 22:50 | henry | Status | new => closed |
2011-03-04 22:50 | henry | Assigned To | => henry |
2011-03-04 22:50 | henry | Resolution | open => fixed |