View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0003997 | OpenFOAM | Bug | public | 2023-07-20 00:33 | 2023-07-20 10:52 |
| Reporter | alejandro.lopez | Assigned To | henry | ||
| Priority | high | Severity | major | Reproducibility | always |
| Status | closed | Resolution | no change required | ||
| Platform | GNU/Linux | OS | Ubuntu | OS Version | 22.04 LTS |
| Product Version | 10 | ||||
| Fixed in Version | 11 | ||||
| Summary | 0003997: fvModels momentumSource not updating in transient simulations with denseParticleFoam | ||||
| Description | When using momentumSource in fvModels and running a steady state simulation with simpleFoam, it updates fine and the momentum is kept constant. However, when that same momentumSource is set in the denseParticleFoam solver, although available, it does not update. It seems the fvModels source is not considered un UcEqn.H. The initial velocity mapped from the steady state simulation is lost as time progresses and tends to zero. | ||||
| Steps To Reproduce | Set up a periodic pipe simulation with denseParticleFoam and set an initial velocity for the field. Include a momentumSource through fvModels in an attempt to keep that velocity and you will notice the velocity will tend to zero. | ||||
| Additional Information | I have modified the code in UcEqn.H to look like this, recompiled a new solver and it now seems to work. fvVectorMatrix UcEqn ( fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc) - fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) + continuousPhaseTurbulence->divDevTau(Uc) == (1.0/rhoc)*cloudSU + fvModels.source(alphac, Uc) ); The addition is + fvModels.source(alphac, Uc) I am not sure if I am missing something else. It would help a lot if you could update or clarify. The case I have is slightly big to be uploades here but happy to share it with you if needed. | ||||
| Tags | No tags attached. | ||||
|
|
Try the incompressibleDenseParticleFluid solver module in OpenFOAM-11 and OpenFOAM-dev which supersedes denseParticleFoam with many improvements including support for fvModels. |
|
|
Thanks for your quick reply, Henry, I will give it a try asap.. For the time being however, I have openfoam 10 painfully installed in a HPC so it would help a lot if I could at least confirm that the inclussion of that term is correct and that I missed nothing else. Regarding version 11 it seems the source structure has changed quite a bit so I expect the cases from the previous version wont run without making some changes. Is that correct? |
|
|
I do not have a suitable test-case to test your proposed change and you did not provide one. I don't know what changes might be needed to your case to get it to run in OpenFOAM-11, either way you should make them as the new algorithm for the dense particle drag on the continuous phase is significantly better. |
|
|
Already fixed in OpenFOAM-11 and OpenFOAM-dev. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2023-07-20 00:33 | alejandro.lopez | New Issue | |
| 2023-07-20 08:27 | henry | Note Added: 0013078 | |
| 2023-07-20 10:31 | alejandro.lopez | Note Added: 0013079 | |
| 2023-07-20 10:52 | henry | Note Added: 0013080 | |
| 2023-07-20 10:52 | henry | Assigned To | => henry |
| 2023-07-20 10:52 | henry | Status | new => closed |
| 2023-07-20 10:52 | henry | Resolution | open => no change required |
| 2023-07-20 10:52 | henry | Fixed in Version | => 11 |
| 2023-07-20 10:52 | henry | Note Added: 0013081 |
