waveFoam + k-OmegaSSTStable: Wave Height Only 86% of Target at 1 Wavelength — Why?

[LJM-97]

Hello everyone,

How does the improved k-ω SST turbulence model kOmegaSSTStablehere in olaFlow perform with respect to wave attenuation?

I installed it in OpenFOAM v2012 and generated waves using waveFoam with relaxation zones.

I found that the wave height is not much different from the case without the new model: at one wavelength downstream, the wave height is only about 86% of the target input height, and there are no structures in the domain. Is this behavior normal?

I’ll post the specific configuration files (turbulenceProperties, fvSchemes, fvSolution) below. Could someone explain the cause and suggest solutions?

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

simulationType RAS;

RAS
{
RASModel kOmegaSSTStable;//kOmegaSSTBuoyancy;
turbulence on;
printCoeffs on;
}

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

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

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
}

divSchemes
{
default Gauss linear;

div(rhoPhi,U)		Gauss SFCD;
div(phi,alpha)		Gauss vanLeer;
div(phirb,alpha)	Gauss interfaceCompression;

div(rhoPhi,k)		Gauss upwind;
div(rhoPhi,omega)	Gauss upwind;
div(phi,k)			Gauss upwind;
div(phi,omega)		Gauss upwind;
div((phi*interpolate(rho)),k)		Gauss upwind;
div((phi*interpolate(rho)),omega)	Gauss upwind;

div((muEff*dev(T(grad(U)))))		Gauss linear;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;

}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

wallDist
{
method meshWave;
}
fluxRequired
{
default no;
p_rgh;
pcorr;
alpha.water;
}

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

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

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

    MULESCorr       yes;
    nLimiterIter    10;

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

"cellDisplacement.*"
{
    solver          PCG;
    preconditioner  DIC;
    tolerance       1e-06;
    relTol          0;
    maxIter         100;
}

"pcorr.*"
{
    solver          PCG;
    preconditioner  DIC;
    tolerance       1e-10;
    relTol          0;
}

p_rgh
{
    solver          PCG;
    preconditioner  DIC;
    tolerance       1e-07;
    relTol          0.01;
}

p_rghFinal
{
    solver          PCG;
    preconditioner  DIC;
    tolerance       1e-07;
    relTol          0;
}

"(k|omega)"
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance       1e-07;
    relTol          0;
}

"(k|omega)Final"
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance       1e-08;
    relTol          0;
}

U
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance       1e-06;
    relTol          0;
}

UFinal PBiCG
{
    preconditioner   DILU;
    tolerance        1e-06;
    relTol           0;
};

}

PIMPLE
{
momentumPredictor no;
nOuterCorrectors 2;
nCorrectors 3;
nNonOrthogonalCorrectors 1;

ddtCorr                    yes;
correctPhi                 yes;

moveMeshOuterCorrectors yes;

}

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

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