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   April 13, 2011, 16:41
Default
  #101
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 17
pablodecastillo is on a distinguished road
I ask about Quality mesh because it is very important in multiphase, in fact, he wrote "after 2 days of running in an 8 core pc, i'm stil at 0.15secs just 0.15 seconds, sounds bad quality mesh.

In my experience run with laminar helps to discard problems in turbulent model.
pablodecastillo is offline   Reply With Quote

Old   April 15, 2011, 04:44
Default
  #102
New Member
 
Ippokratis
Join Date: Nov 2010
Location: Athens, Greece
Posts: 13
Rep Power: 15
chripp is on a distinguished road
Hello everyone,

Thanks for your replies! My mesh quality is really smooth now, i've already faced that problem before. I have a question though, as you can see in the attached images, the mesh grading on my case is different from the wigley's. On my mesh is dense only round the ship, while wigleys mesh denses gradually. Is this difference is responsible for the small timestep issue i'm dealing with or it's completely irrelevant?

Is it easy to explain why and how i will test the laminar flow? I'm not really familiar with this.

Right now i'm running 2 cases. One with the controlDict i've posted before (default wigley case) and one with the changes Niels suggested. The first is becoming unstable now (i get some negative prices) but generally until that it was right. The other one is really stable but the drag force is much larger than expected. Of course it's early to say because both of them are still below 1 second.

The timestep on the first case is small [it's running on an 8-core pc for 4 days now and it's still at 0.45 sec]. The other one is running on an 4-core pc and the timestep is quite small, bigger than the first of course, but stable. [about a day to reach 0,035 sec]

Thanks again for your time!
Best regards
Ippokratis

mymesh.jpg wigleymesh.jpg

Edit: I've finally tested the laminar flow, it's going quite better than before [little bigger timestep and the forces are stable]. Although the drag force is bigger than expected but it's early to tell. Another question though, in the courant number section of the user guide, it mentions that the dx is defined from some cell number to the distance. Is this about the cells in the x axis to the distance of the domain i've defined in the blockMeshDict or something else?

Last edited by chripp; April 17, 2011 at 19:20.
chripp is offline   Reply With Quote

Old   April 18, 2011, 14:12
Default
  #103
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 17
pablodecastillo is on a distinguished road
Hi Ralph,

In my opinion in forcesCalc.h, the force must be calculated like const volScalarField& p_ = mesh_.lookupObject<volScalarField>("p");, and not using p_rgh, if not you never find the balance.

I hope that it is helping.
pablodecastillo is offline   Reply With Quote

Old   April 19, 2011, 06:56
Default
  #104
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 18
egp is on a distinguished road
Who is coming to the 6th OpenFOAM Workshop (http://www.openfoamworkshop.org)? We have just posted the Technical Program for the workshop.

Dr. S.E. Kim will be giving a talk titled, ""Taming OpenFOAM for Ship Hydrodynamics Applications."

Dr. Kevin Maki will be teaching a class in the Training Sessions titled, "Ship Resistance and Propulsion Simulations with OpenFOAM."

Also, the Ship Hydrodynamics Special Interest Group will be meeting on Thursday 16 June as part of the Workshop.

Last edited by egp; April 20, 2011 at 07:19.
egp is offline   Reply With Quote

Old   May 6, 2011, 12:58
Default
  #105
Member
 
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15
vince_44 is on a distinguished road
Hi

I'm a new FOAM user. I'm testing a KRISO hull and when I run InterFoam, I have the following message:
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 void Foam::MULES::limiter<Foam::geometricOneField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 void Foam::MULES::explicitSolve<Foam::geometricOneField , Foam::zeroField, Foam::zeroField>(Foam::geometricOneField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/opt/openfoam171/lib/linux64GccDPOpt/libfiniteVolume.so"
#6
in "/opt/openfoam171/applications/bin/linux64GccDPOpt/interFoam"
#7 __libc_start_main in "/lib/libc.so.6"
#8
in "/opt/openfoam171/applications/bin/linux64GccDPOpt/interFoam"
Exception en point flottant


Any idea about the problem? It's a meshing problem or rather a problem with the solver?

Vincent
vince_44 is offline   Reply With Quote

Old   May 6, 2011, 17:31
Default
  #106
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16
Ralph M is on a distinguished road
Hello Vincent,

In general pressure can cause interfoam to blow up for cases which involve ships... I'm not sure if that's the cause with your problem (search for "floating point" or something like SIGFPE)
A solution for this is to relaxate the pressure which is done with the adopted interFoam solver called shipFoam. This was developed for OF1.6; with a group of people we're trying to get it to work for OF1.7. Go to the shipHydromechanics group for more info (and maybe join the club? )

Cheers,

Ralph
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html
Ralph M is offline   Reply With Quote

Old   May 7, 2011, 11:48
Default
  #107
Member
 
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15
vince_44 is on a distinguished road
Thanks Ralph for your answer. I download OF 1.6 but when I'm trying to unpack the file, there is an error. So for the moment, I stay with OF1.7, where and how I can install ShipFoam for OF1.7?

Cheers

Vincent
vince_44 is offline   Reply With Quote

Old   May 7, 2011, 14:39
Default
  #108
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16
Ralph M is on a distinguished road
You can find the download link in the ship hydromechanics group (follow the link in my signature). Be aware that the tool isn't working as supposed to be; there are still some bugs in it.

Regards,

Ralph
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html
Ralph M is offline   Reply With Quote

Old   May 16, 2011, 03:30
Default
  #109
Member
 
norman1981's Avatar
 
Norman Del Puppo
Join Date: Mar 2009
Location: Hinwil, CH
Posts: 57
Rep Power: 17
norman1981 is on a distinguished road
Quote:
Originally Posted by egp View Post
Who is coming to the 6th OpenFOAM Workshop (http://www.openfoamworkshop.org)? We have just posted the Technical Program for the workshop.

Dr. S.E. Kim will be giving a talk titled, ""Taming OpenFOAM for Ship Hydrodynamics Applications."

Dr. Kevin Maki will be teaching a class in the Training Sessions titled, "Ship Resistance and Propulsion Simulations with OpenFOAM."

Also, the Ship Hydrodynamics Special Interest Group will be meeting on Thursday 16 June as part of the Workshop.
Dear Eric,

I am happy to say that I'll come to the 6th OpenFOAM Workshop . I read the Technical Program and it looks really promising!

Regards

Norman
norman1981 is offline   Reply With Quote

Old   May 16, 2011, 09:19
Default
  #110
Member
 
Albert Tong
Join Date: Dec 2010
Location: Perth, WA, Australia
Posts: 76
Blog Entries: 1
Rep Power: 15
tfuwa is on a distinguished road
Hi Foamers,

I follow this Thread and swak4Foam at openfoamwiki (http://openfoamwiki.net/index.php/Contrib/swak4Foam) to re-simulate the groovyWaveTank case kindly provided by Eric. While the case run smoothly after changing gamma to alpha1 and pd to p_rgh (reason could be found in this thread), the surface elevation increases gradually to the limit of the tank. I believe the outflow B.C. is not very good but am not sure. Anybody has opinion about this issue? Can you please give me any suggestions to improve the simulation?

Cheers,
Albert
Attached Images
File Type: jpg 1.jpg (8.5 KB, 98 views)

Last edited by tfuwa; May 16, 2011 at 23:36.
tfuwa is offline   Reply With Quote

Old   May 19, 2011, 13:59
Default
  #111
Senior Member
 
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 16
Ralph M is on a distinguished road
Hi Albert,

InterFoam causes sometimes the same issue with a domain which fills itself or has a leak in it. Probably the BC between the versions are not very compatible. What kind of conditions for p, alpha and U did you use for your case?
__________________
CFD for marine applications? Go to http://www.marinecfd.com/ and join the OF Ship Hydromechanics Group: http://www.cfd-online.com/Forums/gro...mechanics.html
Ralph M is offline   Reply With Quote

Old   May 19, 2011, 21:46
Default
  #112
Member
 
Albert Tong
Join Date: Dec 2010
Location: Perth, WA, Australia
Posts: 76
Blog Entries: 1
Rep Power: 15
tfuwa is on a distinguished road
Hi Ralph,

Thanks for your reply. I thought interFoam is quite mature, and I am surprised when you point out it has a problem like this. Anyway, I hope the increase of water surface in my case caused by initial variables or B.C. settings, rather than by code.

Please find attached the whole case (including p_rgh, alpha and U, of course), which can be run directly by "blockMesh, setFields, and interFoam" under OF17x, provided groovyBC has already been installed .

Kind regards,
Albert
Attached Files
File Type: gz groovyWaveTank171.tar.gz (5.8 KB, 106 views)
tfuwa is offline   Reply With Quote

Old   June 16, 2011, 10:10
Default
  #113
New Member
 
Jianxi Yao
Join Date: Apr 2011
Posts: 17
Rep Power: 15
jianxiyao is on a distinguished road
hello, which tool you use to generate the pictures?
jianxiyao is offline   Reply With Quote

Old   August 1, 2011, 05:04
Default
  #114
Member
 
Ben Vernieres
Join Date: Jul 2011
Location: Valencia, Spain
Posts: 42
Rep Power: 15
bouclette is on a distinguished road
Send a message via Skype™ to bouclette
Hi All,

I'm currently trying to run a simple yacht hull at Fn 0.35 using the Wigley setup in the OF 2.0 tutorials and I'm getting the usual parasitic waves in the refined region of the domain.

The mesh is of OK quality (nothing amazing though) with max skew at 2.8 and max orthogonality at 51.6 and average at 12.1. Mesh size about 2.2M cells.

When running without the refinement, I don't get those waves except really close to the hull and I guess this is due to diffusion because of coarse mesh and not mesh quality.

I've been trying different settings playing with Co, nAlphaSubCycles and nAlphaCorr without much success. I also tried different interpolation scheemes without much success neither.

Even with the Wigley case I'm getting those waves when refining the mesh supplied (increasing the box 2 refinement level to 3 instead of 2). I assume the original case is coarse enough so that those waves are not visible but it is not treating them.

Has anyone managed to get some nice results without having to make a beautiful structured mesh? We are currently using FineMarine and with the unstructured mesh from Hexpress having local refinements we are getting real nice results both in terms of pretty pictures and drag data. I wish I could do the same with OF!

Attached are 2 shots of what I'm getting after 3000 iterations with LTSInterFoam.

Regards,

Ben
Attached Images
File Type: jpg Pic1.jpg (55.6 KB, 155 views)
File Type: jpg Pic2.jpg (74.2 KB, 142 views)

Last edited by bouclette; August 1, 2011 at 05:04. Reason: spelling...
bouclette is offline   Reply With Quote

Old   August 1, 2011, 16:18
Default
  #115
Senior Member
 
Dave
Join Date: Jul 2010
Posts: 100
Rep Power: 16
daveatstyacht is on a distinguished road
bouclette,
Looking at your pictures you appear to have a refinement zone that is at the free surface that only covers a portion of your domain. I have noticed this kind of behavior tends to occur in a number of cfd programs due to the change in mesh density. As an experiment it may be worth trying a run with an equal refinement across the whole domain free surface and seeing if that eliminates the problem (ie extending the refinement to the inlet and side of the domain). If you think it is due to reflection you might try extending your domain further aft and applying a coarser mesh there to intentionally dampen out waves.
Regards,
Dave
daveatstyacht is offline   Reply With Quote

Old   August 1, 2011, 17:22
Default
  #116
Member
 
Ben Vernieres
Join Date: Jul 2011
Location: Valencia, Spain
Posts: 42
Rep Power: 15
bouclette is on a distinguished road
Send a message via Skype™ to bouclette
Hi Dave,

Thanks for your reply.

I have tried with no local refinement around the hull and coarser mesh (237k cells) and it does get rid of those waves away from the boat but I still have some dodgy waves close to the hull, especially at bow wave. I will attach a picture later.

Could it be possible to have reflection even with symmetryPlane? my side walls are symmetryPlane.

I guess my issue is partly due to the mesh since there are some settings that do work very well with a structured mesh and do not work when using an unstructured mesh, even if it is of OK quality, for example the series 60 container ship ( https://documents.epfl.ch/users/l/lo...i-OF2.0.tar.gz ) settings will not run at all with a sHM mesh but work perfectly with a structured one (my residuals and forces blow up after 600iterations, problem coming from a random cell that is not far away from a perfect cube).

I don't know if there is a setting that allows to compensate for the unstructured mesh intrinsic defects... Would setting nNonOrthogonalCorrectors (in in fvSolution) to 1 or 2 help ? I will try that next.

Regards,

Ben
bouclette is offline   Reply With Quote

Old   August 1, 2011, 18:17
Default
  #117
Senior Member
 
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 17
pablodecastillo is on a distinguished road
Did you try to solve with interFoam?, i think that it is going to work better or at least you can get more info obseving the courant numbers residuals etc.....
pablodecastillo is offline   Reply With Quote

Old   August 1, 2011, 19:52
Default
  #118
Senior Member
 
Dave
Join Date: Jul 2010
Posts: 100
Rep Power: 16
daveatstyacht is on a distinguished road
Quote:
Originally Posted by bouclette View Post
Could it be possible to have reflection even with symmetryPlane? my side walls are symmetryPlane.
I would avoid using a symmetry plane as your far wall BC since if you think about it physically that is saying you have another boat mirrored on the other side (of course if your wall is sufficiently far enough away it should be fine since the wave pattern is constrained).

The non-ortho corrector will improve the calculation of the fluxes across cell boundaries of non-orthogonal cells which is not a bad thing if you are interested in improving accuracy. I don't think it will help though in addressing this issue.

I have not used the LTS version of interFoam yet so I can't comment on the settings to use with regards to it, but I second the thought of trying a run with interFoam to see if it is specific to your setup with regards to the LTS variant or not. I have only ever seen the artificial waves form at the edges of refinement zones but not nearly as severe as those in your pictures. A higher refinement across the free surface region may help alleviate the issue by at least reducing the amplitude of the waves.

Regards,
Dave

Last edited by daveatstyacht; August 1, 2011 at 19:54. Reason: spelling
daveatstyacht is offline   Reply With Quote

Old   August 1, 2011, 20:02
Default
  #119
Member
 
Ben Vernieres
Join Date: Jul 2011
Location: Valencia, Spain
Posts: 42
Rep Power: 15
bouclette is on a distinguished road
Send a message via Skype™ to bouclette
Hi Pablo, Dave,

Attached is a picture of the coarse mesh that is refined on the whole domain at level 2 in sHM.

Solver was interFoam and the time is 6s so waves just start reaching the outlet and the run isn't converged yet.

As you can see, there are those dodgy waves along the hull and in the wake. I have set the nNonOrthogonalCorrectors to 1 and it didn't seem to have improved anything.

I will run the case longer over night and see what happens but on the last iterations, Co mean was at 0.0051 and Co max at 0.35 for a timestep of 0.0045s so it still seems to be going ok.

Will post an update maņana.

Regards,

Ben
Attached Images
File Type: jpg Pic3.jpg (48.6 KB, 121 views)
bouclette is offline   Reply With Quote

Old   August 2, 2011, 04:23
Default
  #120
Member
 
vincent
Join Date: Apr 2011
Posts: 45
Rep Power: 15
vince_44 is on a distinguished road
Hello

I use LTSInterFoam too. For a fishing vessel, I have resistant results closed to commercial CFD code and analytical results. But the free surface is not really nice (I join some picture).
On the another case (multihull), the free surface is more realistic and results close to another CFD code.
I think the problem is in the meshing. Any idea?

best regard

Vincent
Attached Images
File Type: jpg chalutier_v=6nds_sl_01.jpg (76.7 KB, 197 views)
File Type: jpg MES_02.jpg (40.5 KB, 187 views)
vince_44 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 09:11.