Zentari |
April 4, 2017 09:59 |
Resistance prediction in shallow water using OpenFoam
Dear OpenFOAM Community and forum users,
I would be very thankfull for any advice. I am currently dealing with ship resistance in shallow waters, using a starCCM+ mesh. The domain is pretty simple, (Inlet / Outlet / Bottom / Side1 & 2 / Hull / Atmosphere) meshed around the KCS HULL Model, approximately 6x10^6 Cells. I use interFoam ver. 4.10, using the DTC Interfoam tutorial. No Mesh morphing, the ship is fixed. After a few timesteps, the simulation collapses with following message:
Code:
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 7.89558e-07, Final residual = 3.7645e-10, No Iterations 1
Phase-1 volume fraction = 0.284503 Min(alpha.water) = -0.150908 Max(alpha.water) = 17.5372
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.284503 Min(alpha.water) = -0.150907 Max(alpha.water) = 17.5372
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00458485, No Iterations 5
time step continuity errors : sum local = 7.14223e+63, global = 3.62167e+62, cumulative = 3.62167e+62
GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 9.96527e-08, No Iterations 516
time step continuity errors : sum local = 1.72037e+90, global = 3.60294e+87, cumulative = 3.60294e+87
[17] [18] #0 #0 Foam::error::printStack(Foam::Ostream&)Foam::error::printStack(Foam::Ostream&) addr2line failed
[17] #1 Foam::sigFpe::sigHandler(int) addr2line failed
[18] #1 Foam::sigFpe::sigHandler(int) addr2line failed
[17] #2 ? addr2line failed
[18] #2 ? addr2line failed
[17] #3 void Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&) addr2line failed
[18] #3 void Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&) addr2line failed
[17] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) addr2line failed
[18] #4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::magSqr<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> > const&) addr2line failed
[17] #5 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > >, Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() addr2line failed
[18] #5 Foam::kOmegaSST<Foam::eddyViscosity<Foam::RASModel<Foam::IncompressibleTurbulenceModel<Foam::transportModel> > >, Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::correct() addr2line failed
[17] #6 ? addr2line failed
[18] #6 ? addr2line failed
[17] #7 __libc_start_main addr2line failed
[18] #7 __libc_start_main addr2line failed
[17] #8 ? addr2line failed
[18] #8 ? addr2line failed
addr2line failed
Do you have any Ideas?
Kind regards.
|