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

floating object

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   June 15, 2022, 03:50
Default floating object
  #1
New Member
 
leon
Join Date: Feb 2019
Posts: 6
Rep Power: 7
nang19960916 is on a distinguished road
I am doing a simulation related to a thin floating object. I modified the floating object tutorial from OpenFoam. I changed the size of the object to 2.6m*0.6m*0.1m (Lx*Ly*Lz) and the simulation diverged.
I realized that when the shape of the object becomes thinner, the simulation will diverge easily.

blockMesh:
PHP Code:
scale   1;

vertices
(
    (
0 0 0)
    (
5 0 0)
    (
5 3 0)
    (
0 3 0)
    (
0 0 1)
    (
5 0 1)
    (
5 3 1)
    (
0 3 1)
);

blocks
(
    
hex (0 1 2 3 4 5 6 7) (100 60 20simpleGrading (1 1 1)
); 

DynamicMesh:
PHP Code:
dynamicFvMesh       dynamicMotionSolverFvMesh;

motionSolverLibs    (sixDoFRigidBodyMotion);

motionSolver        sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{
    
patches         (floatingObject);
    
innerDistance   0.05;
    
outerDistance   1.0;

    
centreOfMass    (2.5 1.5 0.55);

    
// Cuboid dimensions
    
Lx              2.6;
    
Ly              0.6;
    
Lz              0.1;

    
// Density of the solid
    
rhoSolid        160;

    
// Cuboid mass
    
mass            #eval{ $rhoSolid*$Lx*$Ly*$Lz };

    // Cuboid moment of inertia about the centre of mass
    
momentOfInertia #codeStream
    
{
        
codeInclude
        
#{
            #include "diagTensor.H"
        #};

        
code
        
#{
            
scalar sqrLx sqr($Lx);
            
scalar sqrLy sqr($Ly);
            
scalar sqrLz sqr($Lz);
            
os  <<
                
$mass
               
*diagTensor(sqrLy sqrLzsqrLx sqrLzsqrLx sqrLy)/12.0;
        
#};
    
};

    
report          on;
    
accelerationRelaxation 0.7;
    
//accelerationDamping 0;

    
solver
    
{
        
type Newmark;
    }

    
constraints
    
{

    }



interFoam log:
PHP Code:
--------------------------------------------------------------------------


Applying solid body motion to entire mesh
Selecting sixDoFSolver Newmark
Translational constraint tensor 
(1 0 0 0 1 0 0 0 1)
Rotational constraint tensor (1 0 0 0 1 0 0 0 1)

PIMPLEno residual control data foundCalculations will employ 3 corrector loops

Reading field p_rgh

Reading field U

Reading
/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    
RASModel        kEpsilon;
    
turbulence      on;
    
printCoeffs     on;
    
Cmu             0.09;
    
C1              1.44;
    
C2              1.92;
    
C3              0;
    
sigmak          1;
    
sigmaEps        1.3;
}


Reading g

Reading hRef
Calculating field g
.h

No MRF models present

No finite volume options present
GAMGPCG
:  Solving for pcorrInitial residual 0, Final residual 0No Iterations 0
time step continuity errors 
sum local 0, global = 0cumulative 0
Constructing face velocity Uf

Courant Number mean
0 max0

Starting time loop

Courant Number mean
0 max0
Interface Courant Number mean0 max0
deltaT 
0.0125
Time 
0.0125

PIMPLE
iteration 1
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.5 0.541657)
    
Centre of mass: (2.5 1.5 0.541657)
    
Orientation: (1 4.03281e-19 -1.90322e-18 -4.52272e-19 0.999667 -0.0258116 1.89218e-18 0.0258116 0.999667)
    
Linear velocity: (-1.06752e-17 -2.13504e-17 -1.33482)
    
Angular velocity: (4.13031 -3.02663e-16 -7.19431e-17)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 3.45605e-06No Iterations 6
time step continuity errors 
sum local 5.34885e-09, global = 5.79173e-10cumulative 5.79173e-10
smoothSolver
:  Solving for alpha.waterInitial residual 0.00443175, Final residual 1.96762e-09No Iterations 3
Phase
-1 volume fraction 0.550752  Min(alpha.water) = 0  Max(alpha.water) = 1
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -2.14784e-08  Max(alpha.water) = 1
GAMGPCG
:  Solving for p_rghInitial residual 1, Final residual 7.21983e-09No Iterations 7
time step continuity errors 
sum local 4.44523e-09, global = 4.61743e-10cumulative 1.04092e-09
PIMPLE
iteration 2
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.49977 0.561033)
    
Centre of mass: (2.5 1.49977 0.561033)
    
Orientation: (1 1.20291e-09 -3.3882e-08 1.58164e-10 0.999194 0.0401426 3.3903e-08 -0.0401426 0.999194)
    
Linear velocity: (-3.9926e-08 -0.0372947 1.76526)
    
Angular velocity: (-6.42454 -5.42066e-06 1.42289e-08)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 3.07093e-06No Iterations 6
time step continuity errors 
sum local 1.26647e-08, global = -6.36036e-11cumulative 9.77313e-10
smoothSolver
:  Solving for alpha.waterInitial residual 0.011091, Final residual 4.94102e-09No Iterations 5
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -5.67741e-32  Max(alpha.water) = 1
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -0.000973394  Max(alpha.water) = 1
GAMGPCG
:  Solving for p_rghInitial residual 0.0295466, Final residual 8.63573e-09No Iterations 5
time step continuity errors 
sum local 3.08501e-07, global = 7.08797e-08cumulative 7.1857e-08
PIMPLE
iteration 3
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.49875 0.523842)
    
Centre of mass: (2.5 1.49875 0.523842)
    
Orientation: (-1.45676e-09 -1.1496e-08 3.54242e-10 0.995446 -0.0953277 1.15825e-08 0.0953277 0.995446)
    
Linear velocity: (1.31346e-07 -0.20049 -4.18524)
    
Angular velocity: (15.2756 -1.84582e-06 6.55606e-08)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 9.35698e-06No Iterations 5
time step continuity errors 
sum local 9.44162e-08, global = 4.15212e-09cumulative 7.60091e-08
smoothSolver
:  Solving for alpha.waterInitial residual 0.0220165, Final residual 1.04293e-09No Iterations 6
Phase
-1 volume fraction 0.550752  Min(alpha.water) = 0  Max(alpha.water) = 1
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -0.000953931  Max(alpha.water) = 1.00341
GAMGPCG
:  Solving for p_rghInitial residual 0.0176407, Final residual 1.08902e-09No Iterations 6
time step continuity errors 
sum local 1.7405e-08, global = -1.80868e-10cumulative 7.58282e-08
smoothSolver
:  Solving for epsilonInitial residual 0.124924, Final residual 4.03731e-07No Iterations 12
smoothSolver
:  Solving for kInitial residual 1, Final residual 5.21168e-07No Iterations 14
ExecutionTime 
1.13 s  ClockTime 3 s

Courant Number mean
0.413601 max3.64401
Interface Courant Number mean0.0758947 max2.528
deltaT 
0.00273438
Time 
0.0152344

PIMPLE
iteration 1
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.49789 0.51534)
    
Centre of mass: (2.5 1.49789 0.51534)
    
Orientation: (-3.07488e-09 -2.06507e-08 4.62783e-10 0.992117 -0.125315 2.08733e-08 0.125315 0.992117)
    
Linear velocity: (9.47395e-08 -0.42316 -2.03328)
    
Angular velocity: (6.79387 -4.94599e-06 1.93534e-07)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 1.13549e-06No Iterations 6
time step continuity errors 
sum local 1.30284e-09, global = 6.11995e-11cumulative 7.58894e-08
smoothSolver
:  Solving for alpha.waterInitial residual 0.00648737, Final residual 9.76138e-09No Iterations 4
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -0.000841302  Max(alpha.water) = 1.00555
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -0.000841302  Max(alpha.water) = 1.47813
GAMGPCG
:  Solving for p_rghInitial residual 0.00238556, Final residual 0.000380392No Iterations 20
time step continuity errors 
sum local 0.00100119, global = -0.000423867cumulative = -0.000423791
PIMPLE
iteration 2
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.49717 0.526951)
    
Centre of mass: (2.5 1.49717 0.526951)
    
Orientation: (1 3.63539e-06 -2.53848e-06 -3.82961e-06 0.996694 -0.0812428 2.23474e-06 0.0812428 0.996694)
    
Linear velocity: (-0.000568451 -0.951659 6.45878)
    
Angular velocity: (-25.6181 -0.00156541 -0.00283657)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 2.98597e-06No Iterations 16
time step continuity errors 
sum local 2.79038e-08, global = 2.03909e-10cumulative = -0.000423791
smoothSolver
:  Solving for alpha.waterInitial residual 0.00945137, Final residual 8.82505e-09No Iterations 29
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -6.15717e-05  Max(alpha.water) = 1.00584
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550752  Min(alpha.water) = -5.69898e-05  Max(alpha.water) = 1.00655
GAMGPCG
:  Solving for p_rghInitial residual 0.0024599, Final residual 8.10289e-09No Iterations 5
time step continuity errors 
sum local 4.83906e-08, global = 1.5697e-09cumulative = -0.000423789
PIMPLE
iteration 3
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.50046 0.470671)
    
Centre of mass: (2.5 1.50046 0.470671)
    
Orientation: (-1.18391e-05 -6.7268e-05 -7.96725e-06 0.957921 -0.287033 6.78357e-05 0.287033 0.957921)
    
Linear velocity: (-0.000952581 1.45159 -34.7061)
    
Angular velocity: (127.833 -0.0487111 -0.00622385)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 5.09782e-06No Iterations 7
time step continuity errors 
sum local 8.83202e-08, global = -1.73598e-08cumulative = -0.000423806
smoothSolver
:  Solving for alpha.waterInitial residual 0.0346697, Final residual 4.07643e-09No Iterations 12
Phase
-1 volume fraction 0.550738  Min(alpha.water) = -0.000587794  Max(alpha.water) = 1.00445
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.00251015  Max(alpha.water) = 1.07466
[CFI-KO8UC1K6BP6:091288 more processes have sent help message help-btl-vader.txt cma-permission-denied
[CFI-KO8UC1K6BP6:09128Set MCA parameter "orte_base_help_aggregate" to 0 to see all help error messages
GAMGPCG
:  Solving for p_rghInitial residual 0.0336004, Final residual 0.0336005No Iterations 20
time step continuity errors 
sum local = -0.318064, global = -7.71198e-05cumulative = -0.000500926
smoothSolver
:  Solving for epsilonInitial residual 0.988594, Final residual 3.03644e-07No Iterations 11
bounding epsilon
min: -0.028384 max19526.6 average6.29484
smoothSolver
:  Solving for kInitial residual 0.98574, Final residual 4.49156e-07No Iterations 11
bounding k
min: -0.0654643 max179.829 average0.176057
ExecutionTime 
2.55 s  ClockTime 5 s

Courant Number mean
1.35532 max731.368
Interface Courant Number mean0.608158 max731.368
deltaT 
2.99092e-06
Time 
0.0152374

PIMPLE
iteration 1
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.50046 0.470567)
    
Centre of mass: (2.5 1.50046 0.470567)
    
Orientation: (-1.18632e-05 -6.74131e-05 -8.01186e-06 0.957811 -0.2874 6.79785e-05 0.2874 0.957811)
    
Linear velocity: (-0.000952433 1.46358 -34.783)
    
Angular velocity: (128.122 -0.0488329 -0.00624071)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 0.999994No Iterations 100
time step continuity errors 
sum local = -0.000342755, global = -8.57687e-08cumulative = -0.000501012
smoothSolver
:  Solving for alpha.waterInitial residual 0.000482363, Final residual 2.91048e-12No Iterations 3
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -173.6  Max(alpha.water) = 179.819
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -88.7178  Max(alpha.water) = 97.2301
GAMGPCG
:  Solving for p_rghInitial residual 0.934422, Final residual 0.934452No Iterations 20
time step continuity errors 
sum local = -6.07187e-05, global = -8.58974e-08cumulative = -0.000501098
PIMPLE
iteration 2
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.50046 0.470567)
    
Centre of mass: (2.5 1.50046 0.470567)
    
Orientation: (-1.18632e-05 -6.74131e-05 -8.01186e-06 0.957811 -0.2874 6.79785e-05 0.2874 0.957811)
    
Linear velocity: (-0.000951759 1.46607 -34.7867)
    
Angular velocity: (128.136 -0.0488398 -0.00624143)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 0.999995No Iterations 100
time step continuity errors 
sum local = -6.0717e-05, global = -8.59274e-08cumulative = -0.000501184
smoothSolver
:  Solving for alpha.waterInitial residual 0.00172738, Final residual 1.65336e-12No Iterations 3
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.607241  Max(alpha.water) = 2.53212
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.00189122  Max(alpha.water) = 1.40816
GAMGPCG
:  Solving for p_rghInitial residual 0.928015, Final residual 0.929433No Iterations 20
time step continuity errors 
sum local = -5.14401e-05, global = -8.58966e-08cumulative = -0.00050127
PIMPLE
iteration 3
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.50046 0.470567)
    
Centre of mass: (2.5 1.50046 0.470567)
    
Orientation: (-1.18632e-05 -6.74131e-05 -8.01186e-06 0.957811 -0.2874 6.79785e-05 0.2874 0.957811)
    
Linear velocity: (-0.000951557 1.46682 -34.7878)
    
Angular velocity: (128.141 -0.0488414 -0.00624151)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 1.00149No Iterations 100
time step continuity errors 
sum local = -5.14565e-05, global = -8.59266e-08cumulative = -0.000501356
smoothSolver
:  Solving for alpha.waterInitial residual 0.000137159, Final residual 2.75897e-11No Iterations 3
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.00193996  Max(alpha.water) = 1.91297
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.00189376  Max(alpha.water) = 1.67964
GAMGPCG
:  Solving for p_rghInitial residual 0.923641, Final residual 0.925003No Iterations 20
time step continuity errors 
sum local = -0.000142742, global = -8.58965e-08cumulative = -0.000501442
smoothSolver
:  Solving for epsilonInitial residual 0.845553, Final residual 7.83306e-08No Iterations 4
bounding epsilon
min3.32231e-22 max3.37541e+06 average146.892
smoothSolver
:  Solving for kInitial residual 0.855809, Final residual 9.43268e-07No Iterations 26
ExecutionTime 
6.39 s  ClockTime 9 s

Courant Number mean
0.0019287 max8.88122
Interface Courant Number mean0.00077699 max8.88122
deltaT 
2.69415e-07
Time 
0.0152376

PIMPLE
iteration 1
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.50046 0.470557)
    
Centre of mass: (2.5 1.50046 0.470557)
    
Orientation: (-1.18653e-05 -6.74262e-05 -8.01589e-06 0.957801 -0.287433 6.79914e-05 0.287433 0.957801)
    
Linear velocity: (-0.000951175 1.46925 -34.7967)
    
Angular velocity: (128.175 -0.0488532 -0.00624253)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 1.00431No Iterations 100
time step continuity errors 
sum local = -1.31756e-05, global = -7.73992e-09cumulative = -0.000501449
smoothSolver
:  Solving for alpha.waterInitial residual 3.85187e-05, Final residual 1.7671e-11No Iterations 2
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.300998  Max(alpha.water) = 2.63512
Applying the previous iteration compression flux
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
MULES
Correcting alpha.water
Phase
-1 volume fraction 0.550745  Min(alpha.water) = -0.299621  Max(alpha.water) = 2.5644
GAMGPCG
:  Solving for p_rghInitial residual 0.990967, Final residual 23.9154No Iterations 20
time step continuity errors 
sum local = -0.000825312, global = -2.40862e-08cumulative = -0.000501473
PIMPLE
iteration 2
forces forces
:
    
rhorho
    Not including porosity effects

6
-DoF rigid body motion
    Centre of rotation
: (2.5 1.50046 0.470569)
    
Centre of mass: (2.5 1.50046 0.470569)
    
Orientation: (-1.18945e-05 -6.8939e-05 -8.41955e-06 0.957814 -0.287389 6.94491e-05 0.287389 0.957814)
    
Linear velocity: (-0.182999 -20.8448 54.1237)
    
Angular velocity: (-213.045 -10.8675 -3.02622)
GAMGPCG:  Solving for pcorrInitial residual 1, Final residual 29.8837No Iterations 100
time step continuity errors 
sum local = -0.0104674, global = -5.69035e-08cumulative = -0.00050153
-------------------------------------------------------
Primary job  terminated normallybut 1 process returned
a non
-zero exit code.. Per user-directionthe job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero statusthus causing
the job to be terminated
The first process to do so was:

  
Process name: [[58911,1],8]
  Exit 
code:    145
-------------------------------------------------------------------------- 
nang19960916 is offline   Reply With Quote

Reply


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
engineFoam with layers - pressure problems when adding layers mturcios777 OpenFOAM Running, Solving & CFD 23 January 4, 2023 21:56
floating object shaiashe OpenFOAM Running, Solving & CFD 0 May 18, 2022 02:12
[OpenFOAM] Floating Object Tutorial. No fluid shown when I click on "Play" Ivy Zhou ParaView 6 October 9, 2013 03:02
Floating object with 6DoF in InterDyMFoam wavytracy OpenFOAM Programming & Development 2 November 5, 2010 10:28
Floating object in a fluid benjaminm CFX 7 November 14, 2009 05:46


All times are GMT -4. The time now is 17:59.