[kenrichaiz]

Dear Niels,
Really appreciate for being responsive. I've successfully compiled the solver. Btw, I'm using OpenFOAM 3.0.1. However, i've encounter this error when I try to use the waveFoam solver in the provided waveFlume tutorial, also when running the same solver in damBreak. Really hope you can help me out on this because it is going to be extremely helpful in my research. Thanks mate.

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.0.1-119cac7e8750
Exec : waveFoam
Date : May 30 2016
Time : 14:47:24
Host : "OECU"
PID : 28150
Case : /home/oecu/OpenFOAM/oecu-3.0.1/waves2Foam/tutorials/waveFoam/waveFlume
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

// using new solver syntax:
pcorr
{
solver GAMG;
tolerance 1e-07;
relTol 0;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

// using new solver syntax:
p_rgh
{
solver GAMG;
tolerance 1e-07;
relTol 0;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

// using new solver syntax:
p_rghFinal
{
solver GAMG;
tolerance 1e-08;
relTol 0;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

// using new solver syntax:
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-09;
relTol 0;
}

// using new solver syntax:
UFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-09;
relTol 0;
}

// using new solver syntax:
gamma
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
}


PIMPLE: Operating solver in PISO mode


Reading g

Reading waveProperties

Reading waveProperties
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present

No finite volume options present

GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.0011976
Time = 0.0011976

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0
Phase-1 volume fraction = 0.666667 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.666667 Min(alpha.water) = 0 Max(alpha.water) = 1
Relaxing time: 0.01 s


--> FOAM FATAL IO ERROR:
keyword div(((rho*nuEff)*dev2(T(grad(U))))) is undefined in dictionary "/home/oecu/OpenFOAM/oecu-3.0.1/waves2Foam/tutorials/waveFoam/waveFlume/system/fvSchemes.divSchemes"

file: /home/oecu/OpenFOAM/oecu-3.0.1/waves2Foam/tutorials/waveFoam/waveFlume/system/fvSchemes.divSchemes from line 33 to line 38.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 442.

FOAM exiting

[kenrichaiz]

After adding div(((rho*nuEff)*dev2(T(grad(U))))) to /home/oecu/OpenFOAM/oecu-3.0.1/waves2Foam/tutorials/waveFoam/waveFlume/system/fvSchemes.divSchemes, I also have another error as below:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 3.0.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 3.0.1-119cac7e8750
Exec : waveFoam
Date : May 30 2016
Time : 15:13:22
Host : "OECU"
PID : 28309
Case : /home/oecu/OpenFOAM/oecu-3.0.1/waves2Foam/tutorials/waveFoam/waveFlume
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

// using new solver syntax:
pcorr
{
solver GAMG;
tolerance 1e-07;
relTol 0;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

// using new solver syntax:
p_rgh
{
solver GAMG;
tolerance 1e-07;
relTol 0;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

// using new solver syntax:
p_rghFinal
{
solver GAMG;
tolerance 1e-08;
relTol 0;
smoother DIC;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

// using new solver syntax:
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-09;
relTol 0;
}

// using new solver syntax:
UFinal
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-09;
relTol 0;
}

// using new solver syntax:
gamma
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0;
}


PIMPLE: Operating solver in PISO mode


Reading g

Reading waveProperties

Reading waveProperties
Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
No MRF models present

No finite volume options present

GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG: Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Courant Number mean: 0 max: 0

Starting time loop

Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.0011976
Time = 0.0011976

PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0, Final residual = 0, No Iterations 0
Phase-1 volume fraction = 0.666667 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.666667 Min(alpha.water) = 0 Max(alpha.water) = 1
Relaxing time: 0.02 s


--> FOAM FATAL IO ERROR:
unknown div scheme

Valid div schemes are :
1(Gauss)

file: /home/oecu/OpenFOAM/oecu-3.0.1/waves2Foam/tutorials/waveFoam/waveFlume/system/fvSchemes.divSchemes.div(((rho*nuEff)*dev2(T(grad( U))))) at line 0.

From function divScheme<Type>::New(const fvMesh&, Istream&)
in file /home/openfoam/OpenFOAM/OpenFOAM-3.0.1/src/finiteVolume/lnInclude/divScheme.C at line 82.

FOAM exiting


FYI this is my fvSchemes:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
// default cellLimited fourth 1;
default Gauss linear;
grad(U) Gauss linear;
grad(alpha1) Gauss linear;
}

divSchemes
{
// default Gauss linear;
div(rhoPhi,U) Gauss limitedLinearV 1;
// div(phi,gamma) Gauss vanLeer;
// div(rho*phi,U) Gauss MUSCL;
div((muEff*dev(T(grad(U))))) Gauss linear;
div(phi,alpha) Gauss MUSCL;
div(phirb,alpha) Gauss interfaceCompression;
// div(rhoPhi,U) Gauss limitedLinearV 1;
div(((rho*nuEff)*dev2(T(grad(U)))));
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p_rgh;
pcorr;
alpha.water;
}

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

[kenrichaiz]

oh i didnt add Gauss linear to it. Sorry and thank you very much. Hahaha