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

Floating object tutorial case in OpenFOAM 1.7

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 27, 2010, 03:17
Default Floating object tutorial case in OpenFOAM 1.7
  #1
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Hi.

Just installed OpenFOAM 1.7 on my Ubuntu 10.04 LTS.
Works well - easiest installation so far !

I'm just running the new floating object tutorial case.
Can anybody tell me how the case is structured?

It looks really interesting and I would like to know more about it.
Especially about the field pointDisplacement and the boundary condition sixDoFRigidBodyDisplacement.

Happy brainstorming.
sega
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   June 28, 2010, 15:10
Default
  #2
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
Hello,

Sorry, but have you running the tutorial example without any problem (Floating point exception).

I would like to run the example, but I obtain an error.

Best regards
anmartin is offline   Reply With Quote

Old   June 28, 2010, 16:14
Default
  #3
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Quote:
Originally Posted by anmartin View Post
Hello,

Sorry, but have you running the tutorial example without any problem (Floating point exception).

I would like to run the example, but I obtain an error.

Best regards
Everything is fine here.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   June 29, 2010, 12:00
Default
  #4
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
Sorry to trouble you,

I have installed this version of linux,
Distributor ID: Ubuntu
Description: Ubuntu 10.04 LTS
Release: 10.04
Codename: lucid

I have installed OF170 as is described in webpage (
Ubuntu/Debian pack installation), I try to run the floatingObject tutorial example copy it from multiphase to my run directory and I obtain the following error.
Time = 1.2731

Centre of mass: (0.373595449653 0.351882909424 0.482406400186)
Linear velocity: (-0.357961120598 -0.419740647472 -0.322074163028)
Angular velocity: (1.07106504184 -0.903173170713 -0.0110838047329)
GAMG: Solving for cellDisplacementx, Initial residual = 0.0056044966347, Final residual = 5.8395571122e-06, No Iterations 5
GAMG: Solving for cellDisplacementy, Initial residual = 0.00549399203972, Final residual = 5.71456061284e-06, No Iterations 5
GAMG: Solving for cellDisplacementz, Initial residual = 0.00997063228984, Final residual = 3.85256246304e-06, No Iterations 6
Execution time for mesh.update() = 0.68 s
time step continuity errors : sum local = 3.84344042594e-10, global = 2.49126712936e-11, cumulative = -0.000105905828067
GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 2.81917055336e-06, No Iterations 7
time step continuity errors : sum local = 1.08353248665e-15, global = 2.55464254723e-16, cumulative = -0.000105905828067
MULES: Solving for alpha1
Liquid phase volume fraction = 0.533377598497 Min(alpha1) = -3.16898215629e+294 Max(alpha1) = 2.40723995432
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/opt/openfoam170/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 in "/lib/libc.so.6"
#3 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#4 Foam::fv::gaussGrad<double>::grad(Foam::GeometricF ield<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::tmp<Foam::GeometricField<Foam:uterProduct< Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#6 Foam::tmp<Foam::GeometricField<Foam:uterProduct< Foam::Vector<double>, double>::type, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::grad<double>(Foam::GeometricField<doubl e, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#7 Foam::LimitedScheme<double, Foam::vanLeerLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#8 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#9 Foam::surfaceInterpolationScheme<double>::interpol ate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libinterfaceProperties.so"
#10 Foam::fv::gaussConvectionScheme<double>::interpola te(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#11 Foam::fv::gaussConvectionScheme<double>::flux(Foam ::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/opt/openfoam170/lib/linux64GccDPOpt/libfiniteVolume.so"
#12
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#13
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#14
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
#15 __libc_start_main in "/lib/libc.so.6"
#16
in "/opt/openfoam170/applications/bin/linux64GccDPOpt/interDyMFoam"
Floating point exception

If you could help me I will be very grateful

Best regards,


anmartin is offline   Reply With Quote

Old   June 29, 2010, 12:27
Default
  #5
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
And you have a 64bit machine and operating system?
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   June 29, 2010, 14:12
Default
  #6
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
I have the following system and cpu

Ubuntu version 10.04 64 bits
anmartin@anmartin2:~/OpenFOAM/OpenFOAM-1.7.0/bin$ uname -ar
Linux anmartin2 2.6.32-22-generic #36-Ubuntu SMP Thu Jun 3 19:31:57 UTC 2010 x86_64 GNU/Linux

machine 64 bits too
anmartin@anmartin2:~/OpenFOAM/OpenFOAM-1.7.0/bin$ sudo lshw -C cpu
*-cpu
description: CPU
product: Intel(R) Core(TM) i5 CPU 750 @ 2.67GHz
vendor: Intel Corp.
physical id: 4
bus info: cpu@0
version: Intel(R) Core(TM) i5 CPU 750 @ 2.67GHz
serial: To Be Filled By O.E.M.
slot: CPU 1
size: 1200MHz
capacity: 1200MHz
width: 64 bits
clock: 133MHz
capabilities: fpu fpu_exception wp vme de pse tsc msr pae mce cx8 apic sep mtrr pge mca cmov pat pse36 clflush dts acpi mmx fxsr sse sse2 ss ht tm pbe syscall nx rdtscp x86-64 constant_tsc arch_perfmon pebs bts rep_good xtopology nonstop_tsc aperfmperf pni dtes64 monitor ds_cpl vmx smx est tm2 ssse3 cx16 xtpr pdcm sse4_1 sse4_2 popcnt lahf_lm ida tpr_shadow vnmi flexpriority ept vpid cpufreq

Thank you for your help

Best regards
anmartin is offline   Reply With Quote

Old   July 13, 2010, 11:19
Default
  #7
New Member
 
yannH
Join Date: Feb 2010
Posts: 26
Rep Power: 7
yannH is on a distinguished road
hi Anmartin,

Have you succeed with the floating point error? I have the same trouble but for me it works two times before this problem..

best regards,

yann
yannH is offline   Reply With Quote

Old   July 14, 2010, 04:05
Default
  #8
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
Hi Yann,

Yes, Henry Weller solved my problem, (it is an issue with the initialization (or lack of it) of the pressure (p_rgh) which imparts an unphysical impulse on the object which then moves further than the mesh-motion can cope with).
You only need to pull the latest OpenFOAM-1.7.x recompile it and try it out, it works for me.
OpenFoam.1.7.x floatingObject tutorial case & MULES:alpha1 greater than 1
best regards,
anmartin is offline   Reply With Quote

Old   January 17, 2011, 08:09
Default 2D case of floatingObject
  #9
Member
 
Stefan
Join Date: Jan 2010
Location: Kiel, Germany
Posts: 72
Rep Power: 7
SD@TUB is on a distinguished road
Hello FOAMers,

Before posting a new thread regarding the floatingObject tutorial i am continuing within this discussion.
I want to simualte free floating geometries in 2D and i am insecure how to proceed with setting appropriate boundary conditions especially for the pointDisplacement dictionary due to the dimensions expected from this input file. The values reffering to 3D, e.g. mass, moment of inertia, etc.
Is it possible to simulate 2D cases with the mesh motion solver provided in version 1.7.x, e.g. 2D wave tank (structures floating in long crested waves)?
Maybe someone could give me advice to perform such simulations!

Thanks in advance!
/Stefan
SD@TUB is offline   Reply With Quote

Old   November 7, 2011, 11:02
Default 2D floating object
  #10
New Member
 
Rose Lambert
Join Date: Jun 2011
Posts: 6
Rep Power: 6
Rose2011 is on a distinguished road
Hi Stefan (or anyone else),

Did you manage to simulate a 2D case of a floating object? I would like to do the same but am also having some difficulty in determining the correct parameters for the pointDisplacement file.

Thanks,
Rose
Rose2011 is offline   Reply With Quote

Old   November 10, 2011, 09:26
Default 2D floating object
  #11
New Member
 
Rose Lambert
Join Date: Jun 2011
Posts: 6
Rep Power: 6
Rose2011 is on a distinguished road
No one?


My problem is this: I have a floating box that should respond freely to waves generated at the inlet. I need three degrees of freedom, up/down, left/right and tilt side to side. My case is 2D in the xz plane. To me the obvious choice of the available constraints would be fixedPlane with a normal of (0 1 0). However this case only runs for 1 or 2 time steps and crashes with the following error:
--> FOAM FATAL ERROR:

Maximum number of sixDoFRigidBodyMotion constraint iterations (500) exceeded.

From function Foam::sixDoFRigidBodyMotion::applyConstraints(scal ar deltaT)
in file pointPatchFields/derived/sixDoFRigidBodyMotion/sixDoFRigidBodyMotion.C at line 134.

FOAM exiting

The simulation does run for a few seconds using fixedPlane with normal (0 0 1) but does not allow the box to rise or sink. The simulation crashes when my box twists the mesh.

Am I using the wrong combination of constraints? I am hoping that someone who has modelled a 2D floating object with 3DoF can give me some tips.

Thanks,
Rose
Rose2011 is offline   Reply With Quote

Old   February 9, 2012, 05:29
Default fixedAxis?
  #12
CKH
New Member
 
Kie Hian
Join Date: Aug 2011
Location: Singapore
Posts: 19
Rep Power: 6
CKH is on a distinguished road
Have you tried fixedAxis?

For example, rolling (y-axis rotation) only:

Code:
constraints
        {
            maxIterations   500;

            fixedAxis1
            {
                sixDoFRigidBodyMotionConstraint fixedAxis;
                tolerance       1e-04;
                relaxationFactor 0.7;
                fixedAxisCoeffs
                {
                    axis            ( 0 1 0 );
                }
            }
        }
CKH 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
Modified OpenFOAM Forum Structure and New Mailing-List pete Site News & Announcements 0 June 29, 2009 05:56
Request for pipe bifurcation fluent tutorial case Giorgos Momferratos FLUENT 0 October 20, 2008 15:10
CheckMesh error using a tutorial from OpenFOAM 114 with openFOAM 13 martapajon OpenFOAM Native Meshers: blockMesh 7 January 21, 2008 13:52
case and dat file of FLUENT tutorial Sajad Ranjbaran FLUENT 0 November 8, 2006 07:55
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24


All times are GMT -4. The time now is 00:48.