[JFM]

Good day All

I am attempting to create a LTSInterFoam model as the interFoam models I currently have developed are taking an inordinate amount of time to finish (two months so far on 8/16 nodes with approximately 500,000 cells). Anyways the error output I am getting is:

Time = 1

Flow time scale min/max = 0.0223363, 1
Smoothed flow time scale min/max = 0.0223363, 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191047  Min(alpha1) = -93.8154  Max(alpha1) = 26.8419
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191234  Min(alpha1) = -2410  Max(alpha1) = 8705.97
MULES: Solving for alpha1
Phase-1 volume fraction = 0.191417  Min(alpha1) = -808793  Max(alpha1) = 223717
[0] #0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #2   in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #4  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[0] #6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[0] #7  Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[0] #8  
[0]  in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[0] #9  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[0] #10  
[0]  in "/opt/openfoam220/platforms/linux64GccDPOpt/bin/LTSInterFoam"
[Ubuntu:29885] *** Process received signal ***
[Ubuntu:29885] Signal: Floating point exception (8)
[Ubuntu:29885] Signal code:  (-6)
[Ubuntu:29885] Failing at address: 0x3e8000074bd
[Ubuntu:29885] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fc8032b04a0]
[Ubuntu:29885] [ 1] /lib/x86_64-linux-gnu/libc.so.6(gsignal+0x35) [0x7fc8032b0425]
[Ubuntu:29885] [ 2] /lib/x86_64-linux-gnu/libc.so.6(+0x364a0) [0x7fc8032b04a0]
[Ubuntu:29885] [ 3] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xd2) [0x7fc8043953c2]
[Ubuntu:29885] [ 4] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_RNS1_IS8_EESD_h+0xf1d) [0x7fc804397cfd]
[Ubuntu:29885] [ 5] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x4d3) [0x7fc804398c43]
[Ubuntu:29885] [ 6] /opt/openfoam220/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x137) [0x7fc805e86e17]
[Ubuntu:29885] [ 7] LTSInterFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x152) [0x4562a2]
[Ubuntu:29885] [ 8] LTSInterFoam() [0x42c93a]
[Ubuntu:29885] [ 9] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xed) [0x7fc80329b76d]
[Ubuntu:29885] [10] LTSInterFoam() [0x43145d]
[Ubuntu:29885] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 29885 on node Ubuntu exited on signal 8 (Floating point exception).
--------------------------------------------------------------------------

I have been using the following fvSolutions file:

solvers
{
	pcorr
	{
		solver				PCG;

		preconditioner
		{
			preconditioner  		GAMG;

			smoother				DICGaussSeidel;
			agglomerator			faceAreaPair;
			mergeLevels				1;
			nCellsInCoarsestLevel	10;
			cacheAgglomeration		true;

			tolerance				1e-6;
			relTol					0;
		};

		tolerance		1e-6;
		relTol			0;
	};

	p_rgh
	{
        solver          			GAMG;

        smoother        			GaussSeidel;
        agglomerator   				faceAreaPair;
        mergeLevels     			1;
        nCellsInCoarsestLevel 		10;
        cacheAgglomeration 			true;

        tolerance       			1e-7;
        relTol          			0.01;
    };

    p_rghFinal
    {
        $p_rgh;
        tolerance       1e-7;
        relTol          0;
    }

	U
	{
		solver			smoothSolver;	//PBiCG;
		smoother		GaussSeidel;	//preconditioner	DILU;
		tolerance		1e-06;
		relTol			0;
		nSweeps			1;
	}

	UFinal
	{
		solver			smoothSolver;	//$U;
		smoother		GaussSeidel;
		tolerance		1e-6;
		relTol			0;
		nSweeps			1;
	}

/*	"U.*"
	{
		solver		smoothSolver;
		
		smoother	GaussSeidel;
		nSweeps		1;
		
		tolerance	1e-8;
		relTol		0;
	};
*/

}

PIMPLE
{
    momentumPredictor	no;	// yes;

    nCorrectors			2;	// 1;
    nNonOrthogonalCorrectors 0;

    nAlphaCorr      1;
    nAlphaSubCycles		3; //1;
    cAlpha          1;
    correctPhi		no;

    maxCo				0.75;
    maxAlphaCo			0.75;
    nAlphaSweepIter		1;

    rDeltaTSmoothingCoeff	0.1;
    rDeltaTDampingCoeff		1;
    maxDeltaT				1;
}

relaxationFactors
{
    fields
    {
    }
    equations
    {
    	"U.*"				1;
    }
}


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

And fvSchemes file:

ddtSchemes
{
    default         localEuler rDeltaT;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    div(rho*phi,U)  Gauss linearUpwind grad(U);
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
//	div(phi,k)      Gauss upwind;
//	div(phi,omega)  Gauss upwind;
    div((muEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p_rgh;
    pcorr;
    alpha1;
}


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

I would appreciate if someone could have a look at these and let me know what I am missing. I suspect this is an issue with alpha?

Regards
JFM