[rishabhk28]

I'm trying to simulate fluid flow over a hull to calculate it's resistace, the issue is i'm getting unusual high value of resistance on the hull patch.

I'm using interFoam on openFOAM2212. I have around 6 mil cells with very poor generatoin of layers, I tried a lot to improve the layer additions using snappyhexmesh but to no avail.

I need suggestions, why my force values on the half patch are very high ?


controldict
================================================
application interFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 4000;
deltaT 1;
writeControl timeStep;
writeInterval 100;
purgeWrite 0;
writeFormat binary;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
functions
{
forces
{
type forces;
libs (forces);
patches (hull);
log on;
writeControl timeStep;
writeInterval 1;
CofR (0 0 0);
}

}

FVschemes
==================================================
ddtSchemes
{
default localEuler;
}

gradSchemes
{
default Gauss linear;
limitedGrad cellLimited Gauss linear 1;
}

divSchemes
{
div(rhoPhi,U) Gauss linearUpwind grad(U);
div(phi,alpha) Gauss vanLeer;
div(phirb,alpha) Gauss linear;
div(phi,k) Gauss linearUpwind limitedGrad;
div(phi,omega) Gauss linearUpwind limitedGrad;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}


FVSolution
================================================== =====

solvers
{
"alpha.water.*"
{
nAlphaCorr 2;
nAlphaSubCycles 1;
cAlpha 1;
icAlpha 0;

MULESCorr yes;
nLimiterIter 10;
alphaApplyPrevCorr yes;

solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-8;
relTol 0;
minIter 1;
}

"pcorr.*"
{
solver PCG;
preconditioner
{
preconditioner GAMG;
smoother GaussSeidel;
tolerance 1e-5;
relTol 0;
};
tolerance 1e-5;
relTol 0;
};

p_rgh
{
solver GAMG;
smoother DIC;
tolerance 1e-7;
relTol 0.01;
};

p_rghFinal
{
$p_rgh;
relTol 0;
}

"(U|k|omega).*"
{
solver smoothSolver;
smoother symGaussSeidel;
nSweeps 1;
tolerance 1e-7;
relTol 0.1;
minIter 1;
};
}

PIMPLE
{
momentumPredictor no;

nOuterCorrectors 1;
nCorrectors 2;
nNonOrthogonalCorrectors 1;

maxCo 10;
maxAlphaCo 5;

rDeltaTSmoothingCoeff 0.05;
rDeltaTDampingCoeff 0.5;
nAlphaSpreadIter 0;
nAlphaSweepIter 0;
maxDeltaT 1;
}

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

cache
{
grad(U);
}