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

Help solving LTSInterFoam-wigley hull

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

Like Tree2Likes
  • 1 Post By Dr.Faustus
  • 1 Post By jhoepken

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   July 16, 2012, 09:36
Default Help solving LTSInterFoam-wigley hull
  #1
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
Hi to everybody...
as i've wrote in another post i'm a rookie with open foam... so i've problems with very first steps...
now i'm trying to lounce a wigley-hull project
(i've followed this tutorial http://www.hydroniumion.de/general/s...mesh-tutorial/)

but it appears this error:
Code:
--> FOAM FATAL IO ERROR: 
cannot find file

file: /media/04095B560627AED8/tutorial_foam/wigley2/0/alpha1 at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 73.

FOAM exiting
probably there's some mistakes in alpha1.org but i've checked it and found nothing strange except the version of openfoam (i've installed the last 2.1.1).
Any idea?
Thanks
Attached Files
File Type: txt alpha1.org.txt (1.3 KB, 23 views)
Dr.Faustus is offline   Reply With Quote

Old   July 16, 2012, 13:28
Default
  #2
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
Looks like it can't find your alpha1 file. Did you copy alpha1.org to alpha1 before running setFields?
mgdenno is offline   Reply With Quote

Old   July 17, 2012, 05:40
Default
  #3
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
No sorry in fact now it's works... just erasing the .org extension.
thanks
Maanieh.h likes this.
Dr.Faustus is offline   Reply With Quote

Old   July 17, 2012, 05:53
Default
  #4
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 159
Rep Power: 16
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
Just a tip for the alpha1: Keep 0/alpha1.org and just copy it to 0/alpha1. Afterwards you may run setFields without any issues. You do not have to delete the nonuniform lists from 0/alpha1 if you change your mesh, since you have an untouched copy (0/alpha1.org) that you can just recopy again.
Dr.Faustus likes this.
jhoepken is offline   Reply With Quote

Old   July 17, 2012, 12:25
Default
  #5
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
Could you have a look to my alpha1 ? i've made some changes... but i'm not so conviced by the outlet properties... that in the default file i set as zeroGradient ()

i've setted the sides and the bottom as walls (the domain is quite large)
Attached Files
File Type: txt alpha1.txt (1.3 KB, 23 views)
Dr.Faustus is offline   Reply With Quote

Old   July 18, 2012, 08:08
Default
  #6
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
There's something i've missed in my project as you can see in the image

something it's going better from last running but it seems to be still some mistake...
Now it's go to convergence (i've stopped near 10^-3)
but it seems there's some errors in inlet or U vector... seem that water came in from one side...
Attached Images
File Type: jpg Screenshot-ParaView 3.12.0 64-bit-1.jpg (28.4 KB, 90 views)
Attached Files
File Type: txt U.txt (1.4 KB, 19 views)
File Type: txt omega.txt (1.5 KB, 10 views)
File Type: txt k.txt (1.4 KB, 7 views)
File Type: txt g.txt (936 Bytes, 7 views)

Last edited by Dr.Faustus; July 18, 2012 at 10:30.
Dr.Faustus is offline   Reply With Quote

Old   July 19, 2012, 07:56
Default
  #7
Member
 
Jan
Join Date: Dec 2009
Location: Berlin
Posts: 50
Rep Power: 19
SirWombat is on a distinguished road
Send a message via Skype™ to SirWombat
Where's the bottom in those files ?
__________________
~~~_/)~~~
SirWombat is offline   Reply With Quote

Old   July 19, 2012, 09:51
Default
  #8
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
The mesh is oriented as in the reality: the first vertical face is the inlet (the opposite is the outlet) sides are parallel at the axe of the hull... The top orizontal plan is the atmosfear and the opposite is the bottom....

I've checked from blockmesh and the boundaries seems ok...
But from results it seems that's an opening in the bottom and the inleti is setted on side...
There's a way to visualize the boundaries?
Dr.Faustus is offline   Reply With Quote

Old   July 19, 2012, 11:00
Default
  #9
Member
 
Jan
Join Date: Dec 2009
Location: Berlin
Posts: 50
Rep Power: 19
SirWombat is on a distinguished road
Send a message via Skype™ to SirWombat
ok, i had a look at the wigley-tutorial, where the bottom is defined together with the sides as a "symmetryPlane" in the "blockMeshDict". if you used the same setup that should be ok.

you did not post the boundary setup for "p_rgh". would be interesting to see what that looks like. the rest seems to be ok

The only thing i noticed is, that you set the velocity to stream from outlet to inlet: U_inlet=fixedValue;(-0.4 0 0). Be sure to have "p_rgh" defined accordingly!

Normally p would be "bouyantPressure" where the water comes in. And "zeroGradient" where the water leaves. in your case you will probably need the following setup ( just the opposite of what is defined in the wigley-hull tutorial):

p_rgh:
Code:
...

boundaryField
{
    outlet
    {
        type            buoyantPressure;
        value           uniform 0;
    }

    inlet
    {
        type            zeroGradient;
    }

...
Hope thats the clue here!

Greetings,
Jan
__________________
~~~_/)~~~
SirWombat is offline   Reply With Quote

Old   July 19, 2012, 11:52
Default
  #10
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
Many thanks
In my p_rgh is exactly the opposite (it cames from tutorial)


anyway, sorry but i've to clarify some aspect:
the velocity vector is negative accordling with the direction of x of the reference system..
Ive tried to set it positively but comes out a floating point error...
maybe the velocity is to be setted positively and the floating point error is dued to another error

Last edited by Dr.Faustus; July 23, 2012 at 10:08.
Dr.Faustus is offline   Reply With Quote

Old   July 23, 2012, 11:30
Default
  #11
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
Ok, i've solved the great part of the troubles...
now i can see something...

Now i've to improve the mesh.. that seems vary bad....
The last thing of solving i haven't understood is how i can predict the free survace level
Attached Images
File Type: jpg Screenshot-ParaView 3.12.0 64-bit.jpg (23.6 KB, 78 views)

Last edited by Dr.Faustus; July 25, 2012 at 07:08.
Dr.Faustus is offline   Reply With Quote

Old   July 25, 2012, 16:00
Default
  #12
Senior Member
 
sail's Avatar
 
Vieri Abolaffio
Join Date: Jul 2010
Location: Always on the move.
Posts: 308
Rep Power: 16
sail is on a distinguished road
Ciao Ale.

to see the free surface level in paraview use the contour of VOF = 0.5. that will give you the free surface.

to see the height, on the same filter apply the calculator using z as the value.

let me know if you need further help.
__________________
http://www.leadingedge.it/
Naval architecture and CFD consultancy
sail is offline   Reply With Quote

Old   July 26, 2012, 02:52
Default
  #13
New Member
 
Alessandro Pani
Join Date: Jul 2012
Posts: 21
Rep Power: 13
Dr.Faustus is on a distinguished road
Ciao bello com'è?

I'm going to do some runs and let you know if it works
however i've to predict the height of free surface before running...
the model has no dof so i've to modify the mesh for having a fixed dept
maybe it's better to understand before how this solver calculates the free surface
Dr.Faustus 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 07:20
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 05:03
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 06:24
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 12:16
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


All times are GMT -4. The time now is 20:37.