View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0001850 | OpenFOAM | Bug | public | 2015-09-08 07:56 | 2015-09-09 19:46 |
| Reporter | hannes | Assigned To | henry | ||
| Priority | normal | Severity | major | Reproducibility | random |
| Status | resolved | Resolution | fixed | ||
| Platform | GNU/Linux | OS | Ubuntu | OS Version | 14.04 |
| Summary | 0001850: dynOneEqEddy LES model causes NaN in first time step | ||||
| Description | In the first time step after starting a LES simulation with the dynOneEqEddy model, the residuals in U equation are displayed as NaN. The issue holds for OpenFOAM-2.4.x. | ||||
| Steps To Reproduce | The issue occurs spuriously, depending on the prescribed initial conditions. | ||||
| Additional Information | The problem could be traced into the constructor of the dynOneEqEddy class. It was caused by negative values in the field KK with very small magnitude. While these were already clipped in the "correct"-method (line 153), this clipping procedure was not done in the constructor (where the same expression was repeated). The following modification fixes the issue: @@ -133,7 +133,8 @@ dynOneEqEddy::dynOneEqEddy { bound(k_, kMin_); - const volScalarField KK(0.5*(filter_(magSqr(U)) - magSqr(filter_(U)))); + volScalarField KK(0.5*(filter_(magSqr(U)) - magSqr(filter_(U)))); + KK.max(dimensionedScalar("small", KK.dimensions(), SMALL)); updateSubGridScaleFields(symm(fvc::grad(U)), KK); printCoeffs(); | ||||
| Tags | No tags attached. | ||||
