View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0004017 | OpenFOAM | Bug | public | 2023-09-28 15:51 | 2023-09-30 23:44 |
| Reporter | MatthiasK | Assigned To | henry | ||
| Priority | normal | Severity | crash | Reproducibility | always |
| Status | closed | Resolution | suspended | ||
| Platform | Linux | OS | Ubuntu | OS Version | 22.04 |
| Product Version | 10 | ||||
| Summary | 0004017: multiphaseEulerFoam not working with dynamic mesh and CrankNicolson ddtScheme | ||||
| Description | I ran the OpenFOAM-10 tutorial case: OpenFOAM-10/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVessel2D In system/fvSchemes I changed the ddtScheme from "Euler" to "CrankNicolson 0". The two schemes should be equivalent. With "Euler" the simulation runs without error. With "CrankNicolson 0" the following error occurs: --> FOAM FATAL ERROR: LHS and RHS of - have different dimensions dimensions : [1 2 -2 0 0 0 0] - [1 2 -3 0 0 0 0] From function Foam::dimensionSet Foam::operator-(const Foam::dimensionSet&, const Foam::dimensionSet&) in file dimensionSet/dimensionSet.C at line 527. FOAM aborting The tutorial case: OpenFOAM-10/tutorials/multiphase/multiphaseEulerFoam/laminar/injection works with both ddtSchemes. The error seems to occur only when using a dynamic mesh. | ||||
| Steps To Reproduce | Open the tutorial case "OpenFOAM-10/tutorials/multiphase/multiphaseEulerFoam/laminar/mixerVessel2D" In system/fvSchemes change the ddtScheme from "Euler" to "CrankNicolson 0". | ||||
| Tags | No tags attached. | ||||
|
|
The CrankNicolson scheme is not currently supported in multiphaseEuler. It is complicated to add support for CrankNicolson with MULES and sub-cycling, see twoPhaseSolver/alphaPredictor.C for details of how it was done for 2-phase VoF, and it would be even more complicated for multiphaseEuler. If you would be interested in sponsoring this development please contact CFD Direct to discuss funding. |
