[woodie]

Good day all,

I have a problem with running a simulation from results that were mapped onto the mesh. It seems that the omega field stops the solver from starting the time loop. I am using OpenFoam-v1912, i am simulating a rotor blade with SSTLM turbulence model and SRFSimpleFoam.

Initially, i have obtained a solution on a coarse mesh and mapped it onto the finer one using mapFieldsPar. This seems to be working fine, all the fields are created and ParaView shows no irregularaties.

my procedure:

reconstruct initial solution with reconstructPar -latestTime
map the fields to empty folder "1" with mapFieldsPar ../coarse/ -consistent -sourceTime -latestTime
decompose all mapped fields in folder "1" with decomposePar -fields -latestTime

However, when i start the simulation, it stops before the timeloop even starts; producing this log-file (which does not give any information why it stops):

Exec   : SRFSimpleFoam -parallel
Case   : baseCaseMedium
nProcs : 8
Hosts  :
(
    (b230 8)
)
Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 1


SIMPLE: convergence criteria
    field p      tolerance 1e-05
    field Urel   tolerance 1e-06
    field k      tolerance 1e-06
    field omega  tolerance 1e-06
    field ReThetat       tolerance 1e-06
    field gammaInt       tolerance 1e-06

Reading field p

Reading field Urel

Reading/calculating face flux field phi

AMI: Creating addressing and weights between 59136 source faces and 55552 target faces
AMI: Patch source sum(weights) min:0.999998897 max:1.00000001 average:0.999999999
AMI: Patch target sum(weights) min:0.999986073 max:1.00000009 average:0.999999999
Creating SRF model

Selecting SRFModel rpm
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kOmegaSSTLM
Selecting patchDistMethod meshWave
bounding k, min: 0 max: 16703.522 average: 7.94369825
bounding omega, min: 0 max: 907782106 average: 304329.529
kOmegaSSTLMCoeffs
{
   ...
   ...
}

No finite volume options present

I have boiled the problem down to the omega field; when i replace the omega field in the "1" folder with the initial file from the "0" folder and then do the decomposition, the simulation starts like it should. I believe, fvSchemes and fvSolution are also correct, i will attach them in the next post.

Edit: the same behaviour occurs when running "in serial" as well.

Does anyone have an idea why this is happening?

I hope my question is understandable, if any more information is needed please contact me.

kind regards
Christoph

[woodie]

fvSchemes:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1912                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         cellLimited Gauss linear 1;
    grad(Urel)      cellLimited leastSquares 1;
}

divSchemes
{
    default         none;
    div(phi,Urel)   bounded Gauss linearUpwindV grad(Urel);
    div(phi,k)      bounded Gauss upwind;
    div(phi,omega)  bounded Gauss upwind;
	div(phi,ReThetat) bounded Gauss upwind;
	div(phi,gammaInt) bounded Gauss upwind;

    div((nuEff*dev2(T(grad(Urel))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited 0.333;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited 0.333;
}

wallDist
{
    method meshWave;
}


// ************************************************************************* //

fvSolution:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1912                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
    solver          GAMG;
    smoother		GaussSeidel;
    tolerance       1e-7;		//change to 10-12;
    relTol          0.0;			//change to 0.0 later
	minIter			3;
	nPreSweeps		0;
	nPostSweeps		2;
	cacheAgglomeration	on;
	agglomerator		faceAreaPair;
	nCellsInCoarsestLevel	1000;
	mergeLevels		1;
	maxIter			1500;
    }

    Urel
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-8;
        relTol          0;
	minIter			2;
    }
	
	"(k|ReThetat|gammaInt)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-7;
        relTol          0.00;
		minIter			2;
    }
	
    omega
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-7;
        relTol          0.00;
		minIter			2;
    }
}


SIMPLE
{
    nNonOrthogonalCorrectors 3;
    consistent          yes;
	
	residualControl
	{
		p	1e-5;
		Urel	1e-6;
		k	1e-6;
		omega 	1e-6;
		ReThetat 1e-6;
		gammaInt 1e-6;
	}
}

relaxationFactors
{

	fields
	{
		p	0.8;
	}
    equations
    {
        ".*"            0.75;
    }
}


// ************************************************************************* //

omega BC file, using this instead of the mapped data works just fine:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1912                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 -1 0 0 0 0];


internalField   uniform 14.3325;

boundaryField
{
    BLADE
    {
        blended         1;
        beta1           0.075;
        type            omegaWallFunction;
        value           uniform 14.3325;
    }
    INLET
    {
        type            fixedValue;
        value           uniform 14.3325;
    }
    HUB
    {
        type            zeroGradient;
    }
    OUTLET
    {
        type            inletOutlet;
        inletValue      uniform 14.3325;
        value           uniform 14.3325;
    }
    SIDE_LEFT
    {
        type            cyclicAMI;
        value           uniform 14.3325;
    }
    SIDE_RIGHT
    {
        type            cyclicAMI;
        value           uniform 14.3325;
    }
    TOP
    {
        type            fixedValue;
        value           uniform 14.3325;
    }
}


// ************************************************************************* //