[davidicl]

Hello everyone,

I'm trying to simulate a wavemaker with the interFoam solver (air phase + water phase) considering a laminar flow.
I have implemented my own input kinematics in the waves2foam package based on the uniform velocity profile of a piston-type wavemaker.
The velocity BC at the wavemaker is a no-slip.
My mesh is a 3m long wave flume, with a water depth of 0.6m.
The resolution is 5cm, except around the free surface where I had to refine the mesh to capture the free surface position accurately (resolution of 2mm).
I also had to refine the mesh at the wavemaker to account for the local high velocity gradients due to the no-slip boundary condition an the paddle (resolution of 5mm).

My simulation always crashes at the same time (t = 0.40s).
I didn't use adaptative time-stepping, and reducing the time step doesn't avoid the crash.
Also, the courant number doesn't go up before it crashes.

My boundary conditions are:

U:

inlet
{
type waveVelocity;
refValue uniform ( 0 0 0 );
refGradient uniform ( 0 0 0 );
valueFraction uniform 1;
value uniform ( 0 0 0 );
}
bottom
{
type fixedValue;
value uniform ( 0 0 0 );
}
outlet
{
type fixedValue;
value uniform ( 0 0 0 );
}
atmosphere
{
type pressureInletOutletVelocity;
value uniform ( 0 0 0 );
}
front
{
type slip;
}
back
{
type slip;
}


p_rgh:

front
{
type zeroGradient;
}
back
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type totalPressure;
rho none;
psi none;
gamma 1;
p0 uniform 0;
value uniform 0;
}


alpha1:

inlet
{
type waveAlpha;
refValue uniform 0;
refGrad uniform 0;
valueFraction uniform 1;
value uniform 0;
}
bottom
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
front
{
type zeroGradient;
}
back
{
type zeroGradient;
}


Here are the last 3 iterations before it crashes, could anyone help ?


MULES: Solving for alpha1
Phase-1 volume fraction = 0.884502 Min(alpha1) = -0.000918702 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884501 Min(alpha1) = -0.000975817 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884501 Min(alpha1) = -0.00102995 Max(alpha1) = 1
Relaxing time: 0.04 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000459341, Final residual = 9.65422e-15, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000712516, Final residual = 2.7008e-14, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.00157154, Final residual = 3.37197e-14, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00131087, Final residual = 8.49745e-13, No Iterations 26
GAMG: Solving for p_rgh, Initial residual = 5.6643e-05, Final residual = 8.0606e-13, No Iterations 19
GAMG: Solving for p_rgh, Initial residual = 5.58024e-06, Final residual = 9.06201e-13, No Iterations 15
GAMG: Solving for p_rgh, Initial residual = 3.8308e-07, Final residual = 3.82669e-13, No Iterations 12
time step continuity errors : sum local = 2.65472e-16, global = -1.89997e-17, cumulative = -1.20753e-15
GAMG: Solving for p_rgh, Initial residual = 2.30848e-05, Final residual = 9.23408e-13, No Iterations 18
GAMG: Solving for p_rgh, Initial residual = 2.16917e-06, Final residual = 6.50879e-13, No Iterations 13
GAMG: Solving for p_rgh, Initial residual = 2.36853e-07, Final residual = 9.31239e-13, No Iterations 11
GAMG: Solving for p_rgh, Initial residual = 1.9956e-08, Final residual = 7.22653e-13, No Iterations 8
time step continuity errors : sum local = 5.01315e-16, global = -7.7352e-17, cumulative = -1.28489e-15
GAMG: Solving for p_rgh, Initial residual = 8.72431e-07, Final residual = 9.6311e-13, No Iterations 12
GAMG: Solving for p_rgh, Initial residual = 5.31826e-08, Final residual = 8.48754e-13, No Iterations 9
GAMG: Solving for p_rgh, Initial residual = 5.95305e-09, Final residual = 7.37499e-13, No Iterations 7
GAMG: Solving for p_rgh, Initial residual = 5.83901e-10, Final residual = 6.19992e-13, No Iterations 4
time step continuity errors : sum local = 4.301e-16, global = -2.95754e-17, cumulative = -1.31446e-15
ExecutionTime = 905.95 s ClockTime = 911 s

Courant Number mean: 0.00122961 max: 0.218001
Interface Courant Number mean: 8.42295e-05 max: 0.113623
Time = 0.4055

MULES: Solving for alpha1
Phase-1 volume fraction = 0.884501 Min(alpha1) = -0.00108301 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.8845 Min(alpha1) = -0.00113319 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.8845 Min(alpha1) = -0.00118055 Max(alpha1) = 1
Relaxing time: 0.04 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000458296, Final residual = 1.93453e-14, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000708513, Final residual = 3.09945e-14, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.000680929, Final residual = 2.04159e-13, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.00145775, Final residual = 7.57413e-13, No Iterations 105
GAMG: Solving for p_rgh, Initial residual = 0.000112549, Final residual = 7.91032e-13, No Iterations 86
GAMG: Solving for p_rgh, Initial residual = 1.2736e-05, Final residual = 8.64368e-13, No Iterations 87
GAMG: Solving for p_rgh, Initial residual = 1.70598e-06, Final residual = 8.5032e-13, No Iterations 73
time step continuity errors : sum local = 5.88191e-16, global = -9.4146e-17, cumulative = -1.40861e-15
GAMG: Solving for p_rgh, Initial residual = 4.04471e-05, Final residual = 7.85174e-13, No Iterations 86
GAMG: Solving for p_rgh, Initial residual = 6.2229e-06, Final residual = 8.98764e-13, No Iterations 82
GAMG: Solving for p_rgh, Initial residual = 7.91941e-07, Final residual = 9.64218e-13, No Iterations 70
GAMG: Solving for p_rgh, Initial residual = 1.15067e-07, Final residual = 8.55931e-13, No Iterations 60
time step continuity errors : sum local = 5.92062e-16, global = 8.55776e-17, cumulative = -1.32303e-15
GAMG: Solving for p_rgh, Initial residual = 6.83839e-07, Final residual = 8.45255e-13, No Iterations 70
GAMG: Solving for p_rgh, Initial residual = 7.13413e-08, Final residual = 9.04201e-13, No Iterations 52
GAMG: Solving for p_rgh, Initial residual = 3.75992e-09, Final residual = 8.70047e-13, No Iterations 23
GAMG: Solving for p_rgh, Initial residual = 5.88567e-10, Final residual = 8.6016e-13, No Iterations 29
time step continuity errors : sum local = 5.94978e-16, global = -4.85994e-17, cumulative = -1.37163e-15
ExecutionTime = 908.83 s ClockTime = 914 s

Courant Number mean: 0.00122992 max: 0.21805
Interface Courant Number mean: 8.40171e-05 max: 0.113785
Time = 0.406

MULES: Solving for alpha1
Phase-1 volume fraction = 0.884499 Min(alpha1) = -0.00123552 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884499 Min(alpha1) = -0.00128781 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.884499 Min(alpha1) = -0.00133746 Max(alpha1) = 1
Relaxing time: 0.04 s
DILUPBiCG: Solving for Ux, Initial residual = 0.000459398, Final residual = 1.06474e-14, No Iterations 4
DILUPBiCG: Solving for Uy, Initial residual = 0.000707475, Final residual = 2.65673e-14, No Iterations 4
DILUPBiCG: Solving for Uz, Initial residual = 0.00125732, Final residual = 1.15904e-13, No Iterations 4
[3] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #1 Foam::sigFpe::sigHandler(int) in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #2 in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #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 "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::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 "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
[3] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) in "/home/david/OpenFOAM/openfoam221/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
[3] #7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) in "/home/david/OpenFOAM/david-2.2.1/platforms/linux64GccDPOpt/bin/waveFoam"
[3] #8
[3] in "/home/david/OpenFOAM/david-2.2.1/platforms/linux64GccDPOpt/bin/waveFoam"
[3] #9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
[3] #10
[3] in "/home/david/OpenFOAM/david-2.2.1/platforms/linux64GccDPOpt/bin/waveFoam"

Thanks !
David