New Member
leon
Join Date: Feb 2019
Posts: 6
Rep Power: 7
|
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 20) simpleGrading (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 + sqrLz, sqrLx + sqrLz, sqrLx + 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)
PIMPLE: no residual control data found. Calculations 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 pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Constructing face velocity Uf
Courant Number mean: 0 max: 0
Starting time loop
Courant Number mean: 0 max: 0
Interface Courant Number mean: 0 max: 0
deltaT = 0.0125
Time = 0.0125
PIMPLE: iteration 1
forces forces:
rho: rho
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 pcorr, Initial residual = 1, Final residual = 3.45605e-06, No Iterations 6
time step continuity errors : sum local = 5.34885e-09, global = 5.79173e-10, cumulative = 5.79173e-10
smoothSolver: Solving for alpha.water, Initial residual = 0.00443175, Final residual = 1.96762e-09, No 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_rgh, Initial residual = 1, Final residual = 7.21983e-09, No Iterations 7
time step continuity errors : sum local = 4.44523e-09, global = 4.61743e-10, cumulative = 1.04092e-09
PIMPLE: iteration 2
forces forces:
rho: rho
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 pcorr, Initial residual = 1, Final residual = 3.07093e-06, No Iterations 6
time step continuity errors : sum local = 1.26647e-08, global = -6.36036e-11, cumulative = 9.77313e-10
smoothSolver: Solving for alpha.water, Initial residual = 0.011091, Final residual = 4.94102e-09, No 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_rgh, Initial residual = 0.0295466, Final residual = 8.63573e-09, No Iterations 5
time step continuity errors : sum local = 3.08501e-07, global = 7.08797e-08, cumulative = 7.1857e-08
PIMPLE: iteration 3
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 9.35698e-06, No Iterations 5
time step continuity errors : sum local = 9.44162e-08, global = 4.15212e-09, cumulative = 7.60091e-08
smoothSolver: Solving for alpha.water, Initial residual = 0.0220165, Final residual = 1.04293e-09, No 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_rgh, Initial residual = 0.0176407, Final residual = 1.08902e-09, No Iterations 6
time step continuity errors : sum local = 1.7405e-08, global = -1.80868e-10, cumulative = 7.58282e-08
smoothSolver: Solving for epsilon, Initial residual = 0.124924, Final residual = 4.03731e-07, No Iterations 12
smoothSolver: Solving for k, Initial residual = 1, Final residual = 5.21168e-07, No Iterations 14
ExecutionTime = 1.13 s ClockTime = 3 s
Courant Number mean: 0.413601 max: 3.64401
Interface Courant Number mean: 0.0758947 max: 2.528
deltaT = 0.00273438
Time = 0.0152344
PIMPLE: iteration 1
forces forces:
rho: rho
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: (1 -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 pcorr, Initial residual = 1, Final residual = 1.13549e-06, No Iterations 6
time step continuity errors : sum local = 1.30284e-09, global = 6.11995e-11, cumulative = 7.58894e-08
smoothSolver: Solving for alpha.water, Initial residual = 0.00648737, Final residual = 9.76138e-09, No 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_rgh, Initial residual = 0.00238556, Final residual = 0.000380392, No Iterations 20
time step continuity errors : sum local = 0.00100119, global = -0.000423867, cumulative = -0.000423791
PIMPLE: iteration 2
forces forces:
rho: rho
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 pcorr, Initial residual = 1, Final residual = 2.98597e-06, No Iterations 16
time step continuity errors : sum local = 2.79038e-08, global = 2.03909e-10, cumulative = -0.000423791
smoothSolver: Solving for alpha.water, Initial residual = 0.00945137, Final residual = 8.82505e-09, No 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_rgh, Initial residual = 0.0024599, Final residual = 8.10289e-09, No Iterations 5
time step continuity errors : sum local = 4.83906e-08, global = 1.5697e-09, cumulative = -0.000423789
PIMPLE: iteration 3
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 5.09782e-06, No Iterations 7
time step continuity errors : sum local = 8.83202e-08, global = -1.73598e-08, cumulative = -0.000423806
smoothSolver: Solving for alpha.water, Initial residual = 0.0346697, Final residual = 4.07643e-09, No 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:09128] 8 more processes have sent help message help-btl-vader.txt / cma-permission-denied
[CFI-KO8UC1K6BP6:09128] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages
GAMGPCG: Solving for p_rgh, Initial residual = 0.0336004, Final residual = 0.0336005, No Iterations 20
time step continuity errors : sum local = -0.318064, global = -7.71198e-05, cumulative = -0.000500926
smoothSolver: Solving for epsilon, Initial residual = 0.988594, Final residual = 3.03644e-07, No Iterations 11
bounding epsilon, min: -0.028384 max: 19526.6 average: 6.29484
smoothSolver: Solving for k, Initial residual = 0.98574, Final residual = 4.49156e-07, No Iterations 11
bounding k, min: -0.0654643 max: 179.829 average: 0.176057
ExecutionTime = 2.55 s ClockTime = 5 s
Courant Number mean: 1.35532 max: 731.368
Interface Courant Number mean: 0.608158 max: 731.368
deltaT = 2.99092e-06
Time = 0.0152374
PIMPLE: iteration 1
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 0.999994, No Iterations 100
time step continuity errors : sum local = -0.000342755, global = -8.57687e-08, cumulative = -0.000501012
smoothSolver: Solving for alpha.water, Initial residual = 0.000482363, Final residual = 2.91048e-12, No 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_rgh, Initial residual = 0.934422, Final residual = 0.934452, No Iterations 20
time step continuity errors : sum local = -6.07187e-05, global = -8.58974e-08, cumulative = -0.000501098
PIMPLE: iteration 2
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 0.999995, No Iterations 100
time step continuity errors : sum local = -6.0717e-05, global = -8.59274e-08, cumulative = -0.000501184
smoothSolver: Solving for alpha.water, Initial residual = 0.00172738, Final residual = 1.65336e-12, No 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_rgh, Initial residual = 0.928015, Final residual = 0.929433, No Iterations 20
time step continuity errors : sum local = -5.14401e-05, global = -8.58966e-08, cumulative = -0.00050127
PIMPLE: iteration 3
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 1.00149, No Iterations 100
time step continuity errors : sum local = -5.14565e-05, global = -8.59266e-08, cumulative = -0.000501356
smoothSolver: Solving for alpha.water, Initial residual = 0.000137159, Final residual = 2.75897e-11, No 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_rgh, Initial residual = 0.923641, Final residual = 0.925003, No Iterations 20
time step continuity errors : sum local = -0.000142742, global = -8.58965e-08, cumulative = -0.000501442
smoothSolver: Solving for epsilon, Initial residual = 0.845553, Final residual = 7.83306e-08, No Iterations 4
bounding epsilon, min: 3.32231e-22 max: 3.37541e+06 average: 146.892
smoothSolver: Solving for k, Initial residual = 0.855809, Final residual = 9.43268e-07, No Iterations 26
ExecutionTime = 6.39 s ClockTime = 9 s
Courant Number mean: 0.0019287 max: 8.88122
Interface Courant Number mean: 0.00077699 max: 8.88122
deltaT = 2.69415e-07
Time = 0.0152376
PIMPLE: iteration 1
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 1.00431, No Iterations 100
time step continuity errors : sum local = -1.31756e-05, global = -7.73992e-09, cumulative = -0.000501449
smoothSolver: Solving for alpha.water, Initial residual = 3.85187e-05, Final residual = 1.7671e-11, No 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_rgh, Initial residual = 0.990967, Final residual = 23.9154, No Iterations 20
time step continuity errors : sum local = -0.000825312, global = -2.40862e-08, cumulative = -0.000501473
PIMPLE: iteration 2
forces forces:
rho: rho
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 -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 pcorr, Initial residual = 1, Final residual = 29.8837, No Iterations 100
time step continuity errors : sum local = -0.0104674, global = -5.69035e-08, cumulative = -0.00050153
-------------------------------------------------------
Primary job terminated normally, but 1 process returned
a non-zero exit code.. Per user-direction, the job has been aborted.
-------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:
Process name: [[58911,1],8]
Exit code: 145
--------------------------------------------------------------------------
|