[Emanueledes]

Hi guys,

I'm a beginner in the use of openFoam so I'm having some problems. I'm trying to simulate the slamming of an hull on the water by using interDyMFoam for my graduation project. I have placed the hull at a certain height from the water and let it move at the velocity of 1 m/s in the negative z direction to evaluate pressure on the hull due to the impact on the water. Now, when I run interDyMFoam, i got this error after some timesteps:

Courant Number mean: 0.00212778 max: 6.57496
Interface Courant Number mean: 1.95423e-05 max: 0.0120632
deltaT = 0.000951527
Time = 0.00565446

PIMPLE: iteration 1
DICPCG: Solving for cellMotionUz, Initial residual = 0.00875879, Final residual = 9.55665e-09, No Iterations 44
Execution time for mesh.update() = 9.36 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.0598646, No Iterations 4
time step continuity errors : sum local = 8.80194e-07, global = 3.47383e-08, cumulative = -2.26463e-08
smoothSolver: Solving for alpha.water, Initial residual = 1.0246e-05, Final residual = 4.76097e-11, No Iterations 1
Phase-1 volume fraction = 0.412872 Min(alpha.water) = -1.62871e-88 Max(alpha.water) = 1
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.412872 Min(alpha.water) = -3.47749e-16 Max(alpha.water) = 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::GAMGSolver::scale(Foam::Field&, Foam::Field&, Foam::lduMatrix const&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList&, Foam::PtrList&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const at ??:?
#6 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) at ??:?
#7 Foam::fvMatrix::solve(Foam::dictionary const&) in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/interDyMFoam"
#8 ? in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/interDyMFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10 ? in "/opt/openfoam4/platforms/linux64GccDPInt32Opt/bin/interDyMFoam"

I tried to run moveDynamicMesh and I noticed that some bad cells appeared (non-orthogonalities, negative volumes ecc..) after some timesteps, but the initial mesh is perfect. Maybe it can be a problem of the solver but I've read that velocityLaplacian is the best for these cases. Can someone help me to solve this problem? I would be really grateful. I attach my controlDict and dynamicMeshDict files. Thank you all.

Emanuele