View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0004236 | OpenFOAM | Bug | public | 2025-05-05 16:03 | 2025-05-06 09:24 |
| Reporter | alainislas | Assigned To | henry | ||
| Priority | normal | Severity | crash | Reproducibility | always |
| Status | closed | Resolution | no change required | ||
| Platform | GNU/Linux | OS | Other | OS Version | (please specify) |
| Product Version | 12 | ||||
| Summary | 0004236: Preconditioner gives SIGFPE when using turbulence models in multiphase context | ||||
| Description | When kEpsilon, kOmegaSST, kOmegaSSTSato is used and the alpha fields are initialized via utility such as setField to value 0 or 1 the preconditioner gives floating point exception during usage of preconditioner class. This happens during first iteration. In certain turbulence models such as kEpsilon, kOmegaSST, etc when used in the multiphase context throw an error in the preconditioner. This is due to alpha being zero. Proposal: line 555 in ~/multiphaseEuler/phaseSystem/phaseSystem/phaseSystemSolve.C // phase.maxMin(0, 1); phase.maxMin(SMALL, 1 - SMALL); This error has been reported some years ago in OpenFOAM v1906 | ||||
| Steps To Reproduce | Run tutorial bubbleColumnLaminar with kEpsilon model for example | ||||
| Tags | No tags attached. | ||||
|
|
in fvSolution for example, use "(k.*|epsilon.*)" { solver PBiCGStab; preconditioner DILU; tolerance 1e-12; relTol 0; maxIter 20; minIter 1; } |
|
|
The kEpsilon, kOmegaSST and kOmegaSSTSato models are for turbulence in continuous phases, they do not provide support for the turbulence of dispersed phases, particularly not in the limit of the phase-fraction ->0. If you want to use these models for cases with phase inversion you will need to add terms to handle the low phase-fraction region, possibly using specific fvModels or by deriving special multi-phase versions and this will also resolve the problem you are having with the numerical problems associated with the matrix -> 0 where the phase-fraction -> 0. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2025-05-05 16:03 | alainislas | New Issue | |
| 2025-05-05 16:07 | alainislas | Note Added: 0013565 | |
| 2025-05-06 09:24 | henry | Assigned To | => henry |
| 2025-05-06 09:24 | henry | Status | new => closed |
| 2025-05-06 09:24 | henry | Resolution | open => no change required |
| 2025-05-06 09:24 | henry | Note Added: 0013567 |
