CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Free Surface Ship Flow

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree31Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 24, 2014, 01:42
Default
  #281
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
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.
Sachin m is offline   Reply With Quote

Old   December 30, 2014, 23:43
Default Flow around ship hull
  #282
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
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:rintStack(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)
Sachin m is offline   Reply With Quote

Old   December 31, 2014, 03:26
Default
  #283
Member
 
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15
vince_44 is on a distinguished road
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
vince_44 is offline   Reply With Quote

Old   March 31, 2015, 06:30
Default Trim for the hull (Positioning the geometry in the Domain)
  #284
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
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?
Sachin m is offline   Reply With Quote

Old   March 31, 2015, 06:34
Default
  #285
New Member
 
saman
Join Date: Apr 2012
Posts: 4
Rep Power: 14
samy_20042004 is on a distinguished road
Quote:
Originally Posted by Sachin m View Post
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?
samy_20042004 is offline   Reply With Quote

Old   April 2, 2015, 02:13
Default
  #286
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
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)
Sachin m is offline   Reply With Quote

Old   April 2, 2015, 02:28
Default
  #287
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
And also i would like to know how to get the wetted surface area in paraview
Sachin m is offline   Reply With Quote

Old   April 2, 2015, 02:29
Default
  #288
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
n thanx for the reply saman
Sachin m is offline   Reply With Quote

Old   April 2, 2015, 03:34
Default
  #289
New Member
 
romainRH
Join Date: Jul 2012
Posts: 24
Rep Power: 13
romainRH is on a distinguished road
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 likes this.
romainRH is offline   Reply With Quote

Old   April 2, 2015, 04:32
Default
  #290
New Member
 
saman
Join Date: Apr 2012
Posts: 4
Rep Power: 14
samy_20042004 is on a distinguished road
Quote:
Originally Posted by Sachin m View Post
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)
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
Sachin m likes this.
samy_20042004 is offline   Reply With Quote

Old   April 7, 2015, 09:09
Default
  #291
New Member
 
Join Date: Feb 2015
Posts: 18
Rep Power: 11
jojosaxo is on a distinguished road
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
jojosaxo is offline   Reply With Quote

Old   April 7, 2015, 09:13
Default
  #292
New Member
 
romainRH
Join Date: Jul 2012
Posts: 24
Rep Power: 13
romainRH is on a distinguished road
Hi jojosaxo

Can you be more precise about your problem? When does it occur? Does it start or crash before the first iteration?
romainRH is offline   Reply With Quote

Old   April 7, 2015, 09:33
Default
  #293
New Member
 
Join Date: Feb 2015
Posts: 18
Rep Power: 11
jojosaxo is on a distinguished road
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.
jojosaxo is offline   Reply With Quote

Old   April 7, 2015, 09:35
Default
  #294
New Member
 
James hakes
Join Date: Aug 2013
Posts: 5
Rep Power: 12
jdchakes is on a distinguished road
Try mapfields from your LTS sim and use as initial conditions
jdchakes is offline   Reply With Quote

Old   April 7, 2015, 09:55
Default
  #295
New Member
 
romainRH
Join Date: Jul 2012
Posts: 24
Rep Power: 13
romainRH is on a distinguished road
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 :/
romainRH is offline   Reply With Quote

Old   April 7, 2015, 10:43
Default
  #296
New Member
 
Join Date: Feb 2015
Posts: 18
Rep Power: 11
jojosaxo is on a distinguished road
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...
jojosaxo is offline   Reply With Quote

Old   April 7, 2015, 10:58
Default
  #297
New Member
 
romainRH
Join Date: Jul 2012
Posts: 24
Rep Power: 13
romainRH is on a distinguished road
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?
romainRH is offline   Reply With Quote

Old   April 7, 2015, 11:13
Default
  #298
New Member
 
Join Date: Feb 2015
Posts: 18
Rep Power: 11
jojosaxo is on a distinguished road
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!
jojosaxo is offline   Reply With Quote

Old   April 7, 2015, 11:47
Default
  #299
New Member
 
romainRH
Join Date: Jul 2012
Posts: 24
Rep Power: 13
romainRH is on a distinguished road
I don't have it now i have to look in my old work.
But SuperBee is still available in v2.3.1 ^^
romainRH is offline   Reply With Quote

Old   April 10, 2015, 01:50
Default
  #300
Member
 
Sachin
Join Date: Aug 2014
Location: India
Posts: 84
Rep Power: 11
Sachin m is on a distinguished road
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?
Sachin m is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Free-Surface Ship Flow - Boundary Conditions James Date CFX 1 February 19, 2013 05:42
ship free-surface analysis Andrea Mercuri Siemens 0 September 28, 2004 11:01
Free Surface Flow for Ship sam FLUENT 6 October 24, 2003 05:29
viscous free surface flow past a ship hull lololo Main CFD Forum 0 June 12, 2002 23:02
meshing for surface ship flow boris FLUENT 0 April 24, 2002 20:27


All times are GMT -4. The time now is 14:04.