I am working on a ship hull using LTSInterFoam. My simulations went well in the beginning until after around 15000 timesteps, the force values boomed up. It was around 9000 N and all of a sudden it shooted up to 4e08 +.
Can anyone tell me what could be the possible reason. The hull is that of a barge.

Flow around ship hull
 

Iam using Open Foam 2.3.1.

Iam trying to simulate a flow around ship using LTSInterFoam. I encountered with the below error. can anyone help me out with this?




Time = 0.001

Flow time scale min/max = 4.879e-94, 1.174e-20
Smoothed flow time scale min/max = 4.879e-94, 3.15e-91
Damped flow time scale min/max = 4.879e-94, 3.15e-91
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 0.001018, Final residual = 5.476e-07, No Iterations 2
Phase-1 volume fraction = 0.6114 Min(alpha1) = -1.91e-08 Max(alpha1) = 2.414
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.6114 Min(alpha1) = -3.396e-05 Max(alpha1) = 2.365
smoothSolver: Solving for alpha.water, Initial residual = 0.001111, Final residual = 5.757e-07, No Iterations 4
Phase-1 volume fraction = 0.6114 Min(alpha1) = -1.565e-08 Max(alpha1) = 3.426
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.6215 Min(alpha1) = -0.418 Max(alpha1) = 10.33
smoothSolver: Solving for alpha.water, Initial residual = 0.162, Final residual = 0.0001505, No Iterations 5
Phase-1 volume fraction = 0.617 Min(alpha1) = -0.0001251 Max(alpha1) = 4.383
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.6113 Min(alpha1) = -12.64 Max(alpha1) = 4.254
smoothSolver: Solving for alpha.water, Initial residual = 0.2362, Final residual = 0.0002031, No Iterations 5
Phase-1 volume fraction = 0.6144 Min(alpha1) = -0.2591 Max(alpha1) = 5.291
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.615 Min(alpha1) = -7.501 Max(alpha1) = 25.08
smoothSolver: Solving for alpha.water, Initial residual = 0.2808, Final residual = 0.0002388, No Iterations 5
Phase-1 volume fraction = 0.618 Min(alpha1) = -0.1152 Max(alpha1) = 6.153
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.6383 Min(alpha1) = -7.343 Max(alpha1) = 57.66
smoothSolver: Solving for omega, Initial residual = 2.776e-05, Final residual = 1.459e-08, No Iterations 6
bounding omega, min: -1.657e+05 max: 2.675e+07 average: 9843
smoothSolver: Solving for k, Initial residual = 0.0004899, Final residual = 1.194e-08, No Iterations 9
bounding k, min: -0.007236 max: 5.074 average: 0.00445
smoothSolver: Solving for Ux, Initial residual = 0.5367, Final residual = 2.726e-18, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.6717, Final residual = 4.852e-20, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.7993, Final residual = 9.98e-34, No Iterations 3
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<d ouble>&, Foam::Field<double> const&) const at ??:?
#6 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::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#9 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10
at ??:?
#11 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#12
at ??:?
Floating point exception (core dumped)

Hi Sachin

Could you send more details about your case (fvScheme and FvSolution dict)... It's difficult to find the problem withtout this.
Have you check your mesh? May be it's due to bad mesh.

Best regards

Vince

Trim for the hull (Positioning the geometry in the Domain)
 

Has anyone tried simulation of a vessel with a trim.??
I would like to simulate a ship with a trim of 7 degrees. Is there anyone who can help me out with this?

yes, I already considered trim on my simulations. what is exactly your problem?

is it possible to give a trim to the hull in open foam itself.
That is, if i have a geometry can i simulate it with different trim using that same geometry by positioning it at an angle using openFoam. Is there a way to do that in OpenFoam.(Other than creating diffferent geometries in a cad software)

And also i would like to know how to get the wetted surface area in paraview

n thanx for the reply saman

Hi Sachin M,
Yes you can move the boat with openfoam. You can take a look at interDYnFoam or dynamicFvMesh

For the wetted surface you can try to use a Quality filter to compute surfaces then use calculator to multiply by the volume fraction

Sachin m

here is the command I use for rotating STL file.
surfaceTransformPoints 1.stl 2.stl -yawPitchRoll "(0 -5 0)"
To see the free surface around your hull you should to use Quality filter as romainRH said, put the alpha1 eqaul to 0.5.

Regards,
Saman

Hi Everyone,
I'm currently trying to make a simulation working for a flow around a motorboat (with free surface and so one) using inter foam.
The reason I use inter foam instead of LTSInterFoam is because I want to go for dynamic simulation of the boat with interDyMFoam later on.
For now, I've been able to run a proper simulation of my boat with LTSInterFoam but I'm unable to do the same with interFoam.

I've asked questions in some other threads, got some useful answer but for now nothing is working.
Does anyone have a bit of experience and would agree to share his source file (BC fv Schemes, fv solution controlDict) with me so as I have a working example.
Obviously I don't need geometry since I understand it may be confidential...

If you people have any tricks or tips I'll be more than happy to hear your advice since I'm a wee bit desperate...

Regards
Jordan

Hi jojosaxo

Can you be more precise about your problem? When does it occur? Does it start or crash before the first iteration?

The simulation start normally, it only crashes after a few hundred of iteration. Generally, pressure is increasing. I also often see time step decreasing rf several order that means that velocity is exploding as well.

Try mapfields from your LTS sim and use as initial conditions

Also check courant number. I think the most relevant scheme is for the freeSurface transport equation. I use SuperBee for alpha1 (of.v2.2). i think it is not the same name now :/

I did not get your answer about the "relevant scheme"...

I've already tried to initialize with LTSinterFoam but if I carry on the computation with interFoam it blows up...

I have found that in fvscheme, the parameter that most influence the behaviour of the simulation, is the numerical scheme used for solving the transport equation of volume fraction.

Look for divSchemes and div (phi, alpha) in fvSchemes

Is it clear now?

Yes I got it!
Unfortunately the scheme you're talking about is not available anymore...
Would you be ok to share your file?
Thank's a lot!

I don't have it now i have to look in my old work.
But SuperBee is still available in v2.3.1 ^^

Thank you romainRH and saman for your valuable replies. I would try out and let you know.

In between is there anyone who has tried incorporating wave conditions in to the ship simulation using LTSInterFoam?