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

Free Surface Ship Flow

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

Like Tree12Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   February 13, 2008, 15:19
Default Hi, I'm new to OpenFOAM and
  #1
New Member
 
Tim Franke
Join Date: Mar 2009
Posts: 7
Rep Power: 7
timfranke is on a distinguished road
Hi,

I'm new to OpenFOAM and I'm interested in
calculating the flow around ships incl. the
free surface.
Is there already something "ready to use"
for that application?
The 2nd step might be considering dynamic trim
and sinkage, so 6DOF capabilities will be needed.
Is something like that also available?

Best regards,
Tim
timfranke is offline   Reply With Quote

Old   February 14, 2008, 06:25
Default Hi Tim I am facing similar
  #2
New Member
 
George Kokotis
Join Date: Mar 2009
Posts: 9
Rep Power: 7
gkokgr is on a distinguished road
Hi Tim

I am facing similar issues (and I'm also at the starting point...). I would like to investigate OpenFOAM's potential to deal with this. I cannot help with solid answers ATM, but here are some thoughts:

1. "...The 2nd step might be considering dynamic trim and sinkage, so 6DOF capabilities will be needed....". Are you going after the seakeeping problem/behavior, the increased resistance (i.e. "increased resistance when traveling in a seaway") or to detect differences in the wake pattern (if they are at all detectable...)?
(Since the first step is actually a "leap"*, going after the seakeeping problem in this fashion seems a bit too much...the traditional method is to employ frequency domain methods and from what I understand they yield decent results).
[*: external viscous problem on an unbounded half-space with a priori unknown boundary - the free surface - with non linear boundary conditions...]

2. I think (I am still going through the documentation...) that interFoam may be a suitable solver. A question I still have is about geometry handling and continuity attainability on mesh surfaces. Is it C2?

3. Choice of specific problem formulation? Assuming the "simplest" situation - trying to establish total resistance and form of wake behind a known hull form moving at constant velocity, different options exist:
> Treat this a priori as a steady-state case - assuming that it exists (actually the "steady-state" may be actually an observed average of a fluctuating solution of the true unsteady case as t->oo) or
> Treat the unsteady/transient problem (undisturbed, flat, free surface, hull starting to move from rest @ t=0 to target velocity) while time gets sufficiently large...

4. Proper Domain decomposition? Since this is an external problem, only a portion (containing the hull) of the half space may be examined (and meshed...). One needs to carefully choose the extent of the treated sub-domain and apply suitable matching conditions at the artificial boundaries to reach a realistic simulation.

Best Regards
GK
gkokgr is offline   Reply With Quote

Old   February 14, 2008, 10:04
Default Hi George, at the first ste
  #3
New Member
 
Tim Franke
Join Date: Mar 2009
Posts: 7
Rep Power: 7
timfranke is on a distinguished road
Hi George,

at the first step I'm interested in predicting
the resistance of a ship. To get that value
right the ship has to have to possibility to
adjust its position depending on the flow forces
that act on the hull.
So 6dof capabilities are only needed for resistance prediction.

Best regards,
Tim
timfranke is offline   Reply With Quote

Old   February 15, 2008, 14:30
Default To my knowledge OpenFOAM doesn
  #4
New Member
 
Eddi
Join Date: Mar 2009
Posts: 1
Rep Power: 0
eddi is on a distinguished road
To my knowledge OpenFOAM doesn't include a surface following mesher (VOF with interface reconstruction) as would be necessary for efficient solution.

Or do you plan to use a phase fraction approach similar to the breaking dam example in the OpenFOAM manual? How would you balance the high/low mesh resolution close to the free surface and to the hull on the one hand side and in the distance on the oher hand side, all this while the free surface moves?

Eddi
eddi is offline   Reply With Quote

Old   February 18, 2008, 01:57
Default Ship flow simulations are defi
  #5
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 196
Blog Entries: 1
Rep Power: 8
egp is on a distinguished road
Ship flow simulations are definatey possible with rasInterFoam.

Here is a simple example for the Wigley Hull at Fr=0.316. The plot shows wave-elevation contours.


Tsasan and james.conger like this.
egp is offline   Reply With Quote

Old   February 18, 2008, 02:07
Default One more figure: gamma contou
  #6
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 196
Blog Entries: 1
Rep Power: 8
egp is on a distinguished road
One more figure: gamma contours at a cross-plane near the bow.


james.conger likes this.
egp is offline   Reply With Quote

Old   February 18, 2008, 05:17
Default Hi Eric, could you please
  #7
Member
 
Francesco Boschetto
Join Date: Mar 2009
Location: Italy
Posts: 56
Rep Power: 7
francesco_b is on a distinguished road
Hi Eric,

could you please tell us which is your test-case setup? BC, mesh, etc... Is your hull fixed in the mesh? If not how can you move it?

I'm not able to do ship simulations and it would be something really interesting, please give me any hints to achieve your results

Thanks in advance

Francesco
francesco_b is offline   Reply With Quote

Old   February 18, 2008, 14:04
Default Francesco, You can download
  #8
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 196
Blog Entries: 1
Rep Power: 8
egp is on a distinguished road
Francesco,

You can download the Wigley Hull example at http://idisk.mac.com/egpaterson-Public. Look for the file, wigley.tar.gz.

In naval architecture terminology, this problem is set up for the steady resistance problem, i.e., no motion, and no incident waves. However, I am currently working on these issues.

Eric
egp is offline   Reply With Quote

Old   February 21, 2008, 10:40
Default Hi Eric, Thank you for the
  #9
Member
 
Francesco Boschetto
Join Date: Mar 2009
Location: Italy
Posts: 56
Rep Power: 7
francesco_b is on a distinguished road
Hi Eric,

Thank you for the example, very interesting! I'm interested in coupling interFoam with 6DOF solver and icoDyMFoam to make the simulation more realistic and interesting, I'll try to work on it, if you have any suggestions, you are welcome

Regards

Francesco
francesco_b is offline   Reply With Quote

Old   February 25, 2008, 07:43
Default Hi Eric, Really nice pics a
  #10
Member
 
Pablo
Join Date: Mar 2009
Posts: 95
Rep Power: 7
pablodecastillo is on a distinguished road
Hi Eric,

Really nice pics and interesting results, it is opening a lot of questions, i miss a pic with drag force time history convergence (do you have any?, or routine to get forces in every step?)
Did you try with diferents froude numbers, transom stern?

I am new in OF (i installed the past week), i would like to join to find the best parameters to free surface about bodies. Let me known what mesher are you using.

Thanks,

Pablo
pablodecastillo is offline   Reply With Quote

Old   February 25, 2008, 09:00
Default hi i want to modify the cod
  #11
Senior Member
 
weihong yao
Join Date: Mar 2009
Posts: 117
Rep Power: 7
ivanyao is on a distinguished road
hi
i want to modify the code,how can i get the source code? 
ivanyao is offline   Reply With Quote

Old   November 12, 2008, 03:18
Default Hi, I am currently doing a
  #12
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 299
Rep Power: 8
openfoam_user is on a distinguished road
Hi,

I am currently doing a simulation around a planing hull (15 meters) with an incoming flow velocity of 10 m/s.

The mesh was done using GridPro. The 'checkMesh' is OK.

I am using the rasnterFoam solver.

I have re-used the settings of E. Patterson (model scale of wigley hull) with a maximum Courant number of 0.8.

After 10-15 the run crashes (bounding k and epsilon).

So I have decieded to run the case setting off the turbulence (RASProperties and turbulenceProperties files). And the calculation runs well, but very slowly.

So, is it possible to resume how to define (formula) the turbulence parameters for such a case :
- k
- epsilon
- nuTilda
- nut

And what about the initialization of the flow domain ?

Thanks for helping me !

Stephane.
openfoam_user is offline   Reply With Quote

Old   March 25, 2009, 09:30
Unhappy does not work
  #13
New Member
 
yann
Join Date: Mar 2009
Posts: 2
Rep Power: 0
yannocean is on a distinguished road
hello,
I have download the wigley test case and when I launch interFoam or rasinterFoam, I have the following message. Thers is someone who can help me?


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4.1 |
| \\ / A nd | Web:
http://www.openfoam.org |
|
\\/ M anipulation | |
\*---------------------------------------------------------------------------*/

Exec : foamToFieldview9 ./ wigley
Date : Mar 25 2009
Time : 14:16:00
Host : hpxo.crihan.fr
PID : 18985
Root : /state/partition1/home/hp08003/yandri01/TEST/
Case : wigley
Nprocs : 1
Create time

Create mesh for time = 0
#0 Foam::error:rintStack(Foam::Ostream&) in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigSegv::sigSegvHandler(int) in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt in "/lib64/tls/libc.so.6"
#3 Foam:olyMesh::initMesh() in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#4 Foam:olyMesh:olyMesh(Foam::IOobject const&) in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libOpenFOAM.so"
#5 Foam::fvMesh::fvMesh(Foam::IOobject const&) in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 main in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/foamToFieldview9"
#7 __libc_start_main in "/lib64/tls/libc.so.6"
#8 Foam::regIOobject::readIfModified() in "/share/apps/soft/OpenFOAM/OpenFOAM-1.4.1/applications/bin/linux64GccDPOpt/foamToFieldview9"
yannocean is offline   Reply With Quote

Old   April 3, 2009, 04:50
Default
  #14
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 299
Rep Power: 8
openfoam_user is on a distinguished road
Hi OpenFOAM users,

I have problems in running free surface computations around ships using grids generated with ICEMCFD Hexa and the rasInterFoam solver . Very soon messages about bounding epsilon and bounding k appear and then the case crashes.

ICEMCFD Hexa is not able to keep orthogonality at the walls (with a kind of elliptic smoothing).

With GridPro meshes I succeed in running this kind of cases. Unfortunatly I do not have GridPro. I think that the mesh of the wigley hull of Paterson was done using GridPro. The wigley hull runs well.

Has someone experienced similar problems with ICEMCFD meshes ?

Regards,

Stephane
openfoam_user is offline   Reply With Quote

Old   May 14, 2009, 06:31
Default
  #15
Member
 
matteo lombardi
Join Date: Apr 2009
Posts: 58
Rep Power: 7
matteoL is on a distinguished road
Hi stephane,
I am currently having the same issues you posted. I produce my meshes with ICEM (structured hexa) and altough they seems ok and runs well with CFX, they diverge very quickly with Openfoam. I tried both with interFoam and rasinterFoam and in both cases the simulations diverges (let's say that interfoam is a bit more forgiving and tend to diverge more slowly, while rasinterfoam explodes a bit faster..)

Postprocessing the results, it can clearly be observed that the critical area where the simulation diverges are the one where the mesh is not orthogonal (constant/poly/sets/nonorthoFaces... it can be plotted via paraview and VTK with setSet commands.. btw you can observe it also in icem plotting the maxortho quality criteria over a value of about 65.)

Theorically Openfoam should be able to take care of the non-orthogonality via the "nNonOrthogonalCorrectors" parameters but I have tried many value (from 0 to 6) and it doesn't seems to change at all..

I have tried with Patterson parameters and with many others, but i only manage to delay a bit the divergence at best..

Any ideas?

I don't think it is a problem related to the VOF approach since in some cases the critical area are well under the freesurface (thus completely immerged in water and thus sould be the same of a single-phase simulation in that area), but it seems very strange to me that thus this is a general problem because if so it would affect almost all cases...

Is the nonorthogonal correction for some reason not suited to hexa meshes?

Or if the checkmesh reports some non-orthogonal cells, are those over the limit of Openfoam and thus we cannot expect to recover them via the correction?

suggestions?

thansk,
Matteo
matteoL is offline   Reply With Quote

Old   May 14, 2009, 07:05
Default
  #16
Senior Member
 
stephane sanchi
Join Date: Mar 2009
Posts: 299
Rep Power: 8
openfoam_user is on a distinguished road
Hi Matteo, hi OpenFOAM readers,

Until now I haven't solve the problem. Cases computed with ICEMCFD Hexa grids (wigley hull, Series 60) and Eric Paterson settings always diverge !

I don't know which grid software is using E. Paterson. Maybe Gridgen, maybe GridPro !

I would be nice to have the opinion of Eric Paterson about our specific problems. Hoping Eric paterson is reading this thread.

Matteo, I know that you are mainly interested in computing simulations with this kind of geometries, but try to use a more rounded boat to begin. Let me know about an alternative geometry.

Regards,
Stephane.
openfoam_user is offline   Reply With Quote

Old   May 15, 2009, 06:35
Default
  #17
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 196
Blog Entries: 1
Rep Power: 8
egp is on a distinguished road
Hi Stephane & Matteo,

Is it possible for you to post your case at http://idisk.mac.com/egpaterson-Public ? If so, I'll try to take a look and see if it is something obvious.

It is probably dangerous to think that the settings from a simple problem like the Wigley hull will be uniformly applicable to all problems. You may need to turn on the limiters for the laplacian and grad calculations, e.g.,

gradSchemes
{
default faceLimited leastSquares 0.5;
}
laplacianSchemes
{
default Gauss linear limited 0.33;
}
snGradSchemes
{
default limited 0.33;
}

and also set divSchemes to upwind until you establish something close the correct pressure and wave field.

Good luck, Eric

P.S. I primarily use Pointwise & Gridgen, from Pointwise Inc. Pointwise has an OpenFOAM parser which makes it particularly nice.
lth and Rafael_Coelho like this.
egp is offline   Reply With Quote

Old   May 18, 2009, 10:07
Default
  #18
Member
 
matteo lombardi
Join Date: Apr 2009
Posts: 58
Rep Power: 7
matteoL is on a distinguished road
Hi Eric,
thank you very much for the post!

I have managed to find a way in Icem to control a bit the orthogonality and now checkmesh doesn't complain anymore.

Unfortunately with my old settings the simulation used to diverge as well.

I have then added the schemes correction suggested by Eric and now finally I manage to run it! I still have to check results and so on, but at least it is not exploding!
thank you very much Eric!


A couple more questions:
1)are the limited schemes a sort of under-relaxation and thus at some point I have to switch them off to obtain the correct result or are they just just numerical correction and thus I can leave them on? If the second case: why not leaving them like that as default? and also, I have seen on the user-guide they are some sort of orthogonal correction: what's the difference between these correction and the the non-orthogonal correction step in the fvSolution?

Thanks again,
Matteo
matteoL is offline   Reply With Quote

Old   July 7, 2009, 11:54
Default Wigley hull and series 60 test case
  #19
Member
 
matteo lombardi
Join Date: Apr 2009
Posts: 58
Rep Power: 7
matteoL is on a distinguished road
Hello,
I am wondering whether anyone has managed to validate rasInterfoam on any test case (i.e. wigley hull or similar) since I have tried to do it with the Series60 case but I am getting results that are not very close at all with the experimental one.

In brief:
-Openfoam predicts a value of pressure drag and viscous drag that is about twice as big as the experimental and CFX one (the CFX one run on the same grid)
-My openfoam case starts as first-order-Upwind and at a certain time switches to second order-SFCD-corrected etc etc..
-I have tried on different grid, a medium one (500 000 cells) and a more refined one (1.5 millions cell). Mesh check is fine and the meshes are produced with ICEM and are hexa-structured
-As turbulence model i am using only Kepsilon since SST was not converging.

Thus, my question is:
Has anyone obtained any good agreement on any test case with Openfoam-rasInterfoam?Prof. Patterson showed before some qualitative very interesting results, are those also quantitevely correct? how do they compare with experimental data?

Finally:
anyone can point me to some web page or similar where i can find some experimental data about the wigley hull (waterline, watersurface elevation, measurememt conditions, sink and trim)? I cannot find them anywhere.

Thanks,
best regards,
Matteo
matteoL is offline   Reply With Quote

Old   November 23, 2009, 14:34
Default wigley simulation with OF1.6
  #20
New Member
 
GRD
Join Date: Jun 2009
Posts: 28
Rep Power: 7
g.redondo is on a distinguished road
Interfoam has been updated in version 1.6 to use p instead of pd. I've been trying to update the wigley case to 1.6 without success. Did anyone tried to run it in 1.6? Any help would be greatly appreciated.

Gonzalo
g.redondo is offline   Reply With Quote

Reply

Thread Tools
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 On
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 CD-adapco 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 18:52.