 |
 |
 |
|
|
|
[Sponsors] | ||||
April 11, 2010, 20:42
|
|
#121 |
|
Member
Alex
Join Date: Apr 2010
Posts: 32
Rep Power: 16  |
Hi Juliuslein
Thanks for your answer  Code:
alex@iskandhar:~/OpenFOAM/alex-1.6/Breeder_1.6/shipHydrodynamicIG/tutorials/shipFoam/roll$ ls
0 Allrun VTK constant log.blockMesh log.foamToVTK log.setFields log.shipFoam log.snappyHexMesh system
alex@iskandhar:~/OpenFOAM/alex-1.6/Breeder_1.6/shipHydrodynamicIG/tutorials/shipFoam/roll$ cat constant/polyMesh/boundary
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
2
(
Walls
{
type patch;
nFaces 11725;
startFace 269150;
}
Top
{
type patch;
nFaces 1225;
startFace 280875;
}
)
// ************************************************************************* //
alex@iskandhar:~/OpenFOAM/alex-1.6/Breeder_1.6/shipHydrodynamicIG/tutorials/shipFoam/roll$
I am quite an OpenFOAM n00b... all I know is from reading the PDF's from the workshops.... And by what i have seen... I've still got a long way to go... Can this solver be used to simulate aerodynamic rigid motion instability in incompressible flow regime ( a subsonic twin engine aircraft in a single engine failure condition, for instance )...? BRGDS Alex |
|
|
|
|
|
April 12, 2010, 04:51
|
|
#122 |
|
New Member
Julius
Join Date: Mar 2009
Posts: 27
Rep Power: 17 |
Hi Alexvader,
depending on your case there are 2 possibilities: 1. Rename the patch "walls" to "hull" (the one that was requested by the solver for the dynamic mesh, found in dynamicMeshDict) 2. Change the name of the moved patch in the dynamicMeshDict to "walls", then the solver will find it If you already modified all the occurrences of hull in the "0" directory I would prefer the first solution. Good luck, Julius |
|
|
|
|
|
|
April 12, 2010, 09:21
|
|
#123 |
|
Member
Alex
Join Date: Apr 2010
Posts: 32
Rep Power: 16 |
Hi JuliusLein
I will try this... There is something strange though... : checkMesh says Hull_region0 exists in my mesh : Code:
alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ checkMesh
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5.x |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : checkMesh
Date : Apr 12 2010
Time : 14:01:32
Host : iskandhar
PID : 4622
Case : /home/alex/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll
nProcs : 1
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time
Create polyMesh for time = constant
Time = constant
Mesh stats
points: 98496
faces: 282100
internal faces: 269150
cells: 91875
boundary patches: 2
point zones: 0
face zones: 0
cell zones: 1
Number of cells of each type:
hexahedra: 91875
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0
Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
Walls 11725 11796 ok (not multiply connected)
Hull 1225 1296 ok (not multiply connected)
Checking geometry...
Domain bounding box: (-3 -3 -3) (3 3 2)
Boundary openness (7.22458e-16 7.65537e-17 1.10227e-17) OK.
Max cell openness = 1.51788e-16 OK.
Max aspect ratio = 2.57143 OK.
Minumum face area = 0.0114286. Maximum face area = 0.0293878. Face area magnitudes OK.
Min volume = 0.00195918. Max volume = 0.00195918. Total volume = 180. Cell volumes OK.
Mesh non-orthogonality Max: 0 average: 0
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 4.53341e-14 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1
All face flatness OK.
Mesh OK.
--> FOAM Warning :
From function polyMesh::readUpdateState polyMesh::readUpdate()
in file meshes/polyMesh/polyMeshIO.C at line 205
Number of patches has changed. This may have unexpected consequences. Proceed with care.
Time = 0.001
Mesh stats
points: 130500
faces: 357977
internal faces: 338527
cells: 113945
boundary patches: 3
point zones: 0
face zones: 0
cell zones: 1
Number of cells of each type:
hexahedra: 108799
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 5146
Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
Walls 11725 11796 ok (not multiply connected)
Top 1225 1296 ok (not multiply connected)
Hull_region0 6500 6718 ok (closed singly connected surface)
Checking geometry...
Domain bounding box: (-3 -3 -3) (3 3 2)
Boundary openness (-2.24025e-16 -3.70559e-17 2.59882e-17) OK.
Max cell openness = 2.27681e-16 OK.
Max aspect ratio = 2.57175 OK.
Minumum face area = 0.00017856. Maximum face area = 0.0293893. Face area magnitudes OK.
Min volume = 3.82617e-06. Max volume = 0.00195938. Total volume = 179.75. Cell volumes OK.
Mesh non-orthogonality Max: 50.4794 average: 6.82892
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.999953 OK.
All angles in faces OK.
Face flatness (1 = flat, 0 = butterfly) : average = 1 min = 1
All face flatness OK.
Mesh OK.
Time = 0.002
Mesh stats
points: 130500
faces: 357977
internal faces: 338527
cells: 113945
boundary patches: 3
point zones: 0
face zones: 0
cell zones: 1
Number of cells of each type:
hexahedra: 108799
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 5146
Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Topological cell zip-up check OK.
Face vertices OK.
Face-face connectivity OK.
Number of regions: 1 (OK).
Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface
Walls 11725 11796 ok (not multiply connected)
Top 1225 1296 ok (not multiply connected)
Hull_region0 6500 6718 ok (closed singly connected surface)
Checking geometry...
Domain bounding box: (-3 -3 -3) (3 3 2)
Boundary openness (-2.24264e-16 -3.74502e-17 2.63455e-17) OK.
Max cell openness = 2.69775e-16 OK.
Max aspect ratio = 3.53315 OK.
Minumum face area = 5.65956e-05. Maximum face area = 0.0294002. Face area magnitudes OK.
Min volume = 2.68918e-06. Max volume = 0.00196592. Total volume = 179.744. Cell volumes OK.
Mesh non-orthogonality Max: 53.9665 average: 7.07718
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 1.63149 OK.
*There are 4 faces with concave angles between consecutive edges. Max concave angle = 64.1358 degrees.
<<Writing 4 faces with concave angles to set concaveFaces
Face flatness (1 = flat, 0 = butterfly) : average = 0.999942 min = 0.898383
All face flatness OK.
Mesh OK.
End
Yet it cannot be found in any file of directory 0.002, or in any directory of roll testcase whatsoever.... : Code:
alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ grep -nH 'Hull_region0' 0.002/* alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ grep -nH 'Hull' 0.002/* 0.002/cellLevel:113983: Hull 0.002/meshPhi:350294: Hull 0.002/pointLevel:130536: Hull alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ ls 0 0.001 0.002 Hull_motion.txt constant system alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ grep -nH 'Hull' constant/* constant/dynamicMeshDict:35:diffusivity inverseDistance 1(Hull); constant/dynamicMeshDict:37:// diffusivity quadratic inverseDistance 1(Hull); constant/dynamicMeshDict:38:// diffusivity exponential 2000 inverseDistance 1(Hull); constant/shipDict:27: Hull alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ grep -nH 'Hull_region0' constant/* alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ grep -nH 'Hull' system/* system/snappyHexMeshDict:39: Hull system/snappyHexMeshDict:108: Hull system/snappyHexMeshDict:178: Hull alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ grep -nH 'Hull_region0' system/* alex@iskandhar:~/OpenFOAM/alex-1.5.x/openfoam-extend/shipFoamDemos/roll$ Is this sort of Hardcoded in any binary...? I think I will change all my "Hull"s back to "Hull_region0"s in the roll case ... but this is really weird, even using nano and ctrl-W to find this Hull_region0 reports nothing... how come...? BRGDS Alex |
|
|
|
|
|
|
May 4, 2010, 04:28
|
|
#124 |
|
Member
John Wang
Join Date: Mar 2009
Location: Singapore
Posts: 73
Rep Power: 17 |
Hi guys,
I'm wondering if there's a version of the code for shipFoam that would work in OF-1.6? Thanks. John |
|
|
|
|
|
|
November 10, 2010, 06:41
|
|
#125 | |
|
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 15 |
Hello all,
I'm trying to compile shipFoam under OF 1.7.1. I had to change some of the directories but I got an error which I don't understand: EDIT: I got a version of shipFoam from http://openfoam-extend.svn.sourcefor...ydrodynamicIG/ It seems to work under 1.6 (gamma is changed into alpha1) but some work has to be performed to get it right for 1.7.1. After changing p into p_rgh and changing some other files the wmake process goes well. However, when I run shipFoam I got the following error: Quote:
Regards, Ralph Last edited by Ralph M; November 10, 2010 at 11:16. |
||
|
|
|
||
|
November 13, 2010, 07:41
|
|
#126 |
|
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 15 |
I already solved the problem; in some of the files there was still a reference to "p". The solver is now running smoothly, although I have some questions about the shipDict.
1) What is Amax 2) How is the damping defined? is it a number between 0 and 1? 3) What is the "weightFactor" supposed to do? |
|
|
|
|
|
|
November 19, 2010, 09:59
|
|
#127 |
|
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 15 |
Hello all,
1) What is Amax -->maximum acceleration (duh) 2) How is the damping defined? is it a number between 0 and 1? --> this number is multiplied with the specified unit 3) What is the "weightFactor" supposed to do? --> this is used to use the results from earlier calculations. 1st entry is weighting of current force/moment, 2nd of the previous time step and the 3rd entry of two timesteps ago. The updated version of shipFoam can be found in the OF Ship Hydromechanics group http://www.cfd-online.com/Forums/gro...mechanics.html |
|
|
|
|
|
|
March 3, 2011, 12:14
|
|
#128 |
|
New Member
Mark Beal
Join Date: Feb 2011
Posts: 24
Rep Power: 15 |
MarkC (and others?)
I've just read through this thread and have a couple of questions. I'm interested in simulating yacht hulls resistance and lift with 2DoF (pitch and heave). I've only just started but I was using interDyMFoam with dynamicMotionSolverFvMesh and sixDoFRigidBodyMotionConstraint. 1) What does shipFoam do that interDyMFoam does not? 2) Has interDyMFoam been updated since shipFoam was conceived? 3) If I'm trying to simulate a towing tank, what would be the best approach. 4) (Slightly related) Should I be using fixedValue or movingWallVelocity on my hull patch? Kind Regards, Mark |
|
|
|
|
|
|
March 3, 2011, 12:51
|
|
#129 |
|
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 15 |
Hello Mark,
1) The problem with interDyMFoam is that for some cases the pressure explodes. therefore shipFoam is developed; is does the same as interDyMFoam but uses a weighted pressure from (i think) previous time steps. In this way the solver is more stable. ShipFoam was first developed for OF 1.5 (dev?) but I tried to rewrite it for OF 1.7. I think there are still some issues with the pressure but at the moment I don't have time to finish this. 2) I don't know 3) That depends on your problem. If you are looking for steady state trim/sinkage I think it is the easiest way to keep the ship fixed in flow direction and to let the water have a velocity. 4) I would use buoyantPressure for the pressure patch and for the velocity type movingWallVelocity; value uniform (0 0 0); 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 |
|
|
|
|
|
|
March 8, 2011, 06:26
|
problem compiling
|
#130 | |
|
New Member
Francisco Miguel
Join Date: Apr 2010
Posts: 13
Rep Power: 16 |
Hello, I am trying to do some multiphase ship hydrodynamics tests, and would like to try the shipFoam as I think it could save me quite a bit of work. However, I managed to compile the waveBC and the setHydrostatics on 1.6-ext, but I am having problems with the shipFoam solver. I tried to compile the ODE solver from 1.5, as suggested above, but it won't compile, most likely due to the different OF versions. Still, the error seems to be related with the ODE solver. Any clues/suggestions? Below you can find what I am getting:
Quote:
Francisco |
||
|
|
|
||
|
March 19, 2012, 11:14
|
|
#131 |
|
New Member
Jan Löhrmann
Join Date: Sep 2010
Posts: 21
Rep Power: 15 |
Hi,
I hope somebody is still following this thread as I do have exactly the same problem as in the post before by Francisco. Does anybody have an answer on that problem? Here is what I've done: (1) I downloaded all files from this URL: http://openfoam-extend.svn.sourcefor...ydrodynamicIG/ and copied them into the folders in my OpenFoam-1.6-ext directory. (2) Afterwards I went to /applications/solvers/shipFoam1.6.2 and typed wmake, which produced the same error as mentioned before. (3) Then i tried to compile setHydrostatic with wmake which went through. (4) At last I tried to compile waveBC which produced an error as dynamicFvMesh.H couldn't be found. Can I simply copy the dynamicFvMesh from OF-1.5-dev e.g. into my OF-1.6-ext and compile it, so that waveBC has acces to this file? Does anybody know if that also solves the problem with the bodyMotion.o file as mentioned in the previous post when running wmake in the shipFoam-solver? My original intention was to evaluate the shipFoam solver. Does anybody have more experience with this solver? Are there any papers about it (e.g. how does it perform against benchmarks, etc.)? Best Regards Jan |
|
|
|
|
|
|
March 20, 2012, 02:26
|
|
#132 |
|
Senior Member
Ralph Moolenaar
Join Date: Aug 2010
Location: 's-Hertogenbosch, the Netherlands
Posts: 120
Rep Power: 15 |
Dear Jan,
It has been a while since I used shipFoam for the last time. I think that your problems can be solved by changing the wmake files a bit; did you already had a look at the hydro-group (see link below in my signature)? Furthermore I'm quit sure that no benchmarks have been performed for shipFoam. Can you send me a PM what you are planning to do and what your background is? 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 |
|
|
|
|
|
|
October 7, 2012, 22:50
|
|
#133 | |
|
New Member
Bob Zhang
Join Date: Apr 2012
Posts: 1
Rep Power: 0 |
Quote:
i met the same problem, have you resolve it? |
||
|
|
|
||
 |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| hydrostatic pressure | Amanda | FLUENT | 6 | April 20, 2016 11:00 |
| Hydrostatic pressure(rho*g*h) | Pranesh | FloEFD, FloWorks & FloTHERM | 3 | October 17, 2008 06:18 |
| hydrostatic pressure | multiphase | FLUENT | 0 | May 18, 2003 15:16 |
| Hydrostatic pressure in 5.5 | Jens | CFX | 3 | August 21, 2002 11:05 |
| Hydrostatic Pressure | Rhydar | CFX | 3 | March 6, 2002 09:54 |
4Likes
rintStack(Foam::Ostream&) in "/opt/openfoam171/lib/linuxGccDPOpt/libOpenFOAM.so"
perator<< <Foam::error>(Foam::Ostream&, Foam::errorManip<Foam::error>) in "/home/ralphmoolenaar/OpenFOAM/ralphmoolenaar-1.7.1/applications/bin/linuxGccDPOpt/shipFoam"
Linear Mode
