mahsankhan |
January 24, 2022 11:13 |
attachmentPt for floating body at each time step for overInterDyMFoam solver
Hello everyone,
Can someone please help me with how can I get the value of "attachmentPt" at each time step?
I used the linearSpring restraint for my floating body to attach it to the ground.
The "anchor" points are the coordinates of the 4 anchor points for 4 spings on the ground.
In dynamicMeshDict I also provided the "attachmentPt" for each spring on my moving body. I would like to know the position of those points for each time step as well. But as the OpenFOAM solves the case for "overInterDyMFoam", it is not giving me the values for the "attachmentPt".
Here is a part of the simulation output:
Code:
Courant Number mean: 0.000367988206492 max: 1.30042055746
Interface Courant Number mean: 1.41151225072e-05 max: 0.276089110203
deltaT = 0.00137661866959
Time = 0.264208
PIMPLE: iteration 1
forces forces:
rho: rho
Not including porosity effects
Restraint verticalSpring1: attachmentPt - anchor (0.675369668994 0.346996919813 0.527406581676) spring length 0.924493782863 force (-0.0173035203843 -0.00889034339402 -0.0135125857079)
Restraint verticalSpring2: attachmentPt - anchor (0.675368834426 -0.347003080204 0.527408065947) spring length 0.924496332178 force (-0.0173053136288 0.00889143357967 -0.013514040812)
Restraint verticalSpring3: attachmentPt - anchor (-0.67535116558 -0.347002245643 0.527408524159) spring length 0.924483372856 force (0.017295636459 0.00888667258898 -0.0135068488279)
Restraint verticalSpring4: attachmentPt - anchor (-0.675350331011 0.346997754374 0.527407039888) spring length 0.924480230649 force (0.01729337843 -0.00888540836539 -0.0135050641254)
6-DoF rigid body motion
Centre of rotation: (9.31246335781e-06 -2.70002806243e-06 0.11241267961)
Centre of mass: (9.31246335781e-06 -2.70002806243e-06 0.11241267961)
Orientation: (0.999999999937 -9.78173701946e-06 -5.48358517723e-06 9.78164007906e-06 0.999999999796 -1.76780352421e-05 5.483758098e-06 1.76779816025e-05 0.999999999829)
Linear velocity: (4.4234784925e-05 -2.68512259873e-05 0.00375371956529)
Angular velocity: (0.000302560862431 -0.000112926348239 5.63382227438e-05)
cellVolumeWeight : detected 2 mesh regions
zone:0 nCells:117810
zone:1 nCells:93548
Marking patch-cells on zone 0
Marking patch-cells on zone 1
After patch analysis : nCells : 211358
other : 156798
patch : 50666
overset: 3894
Creating mesh-to-mesh addressing for region0 and region0 regions using cellVolumeWeight
meshToMesh: Using AABBTree method
Overlap volume: 0.114317031147
cellVolumeWeight : detected 3 mesh regions after overset
Front : fraction:1 size:5229
Overset analysis : nCells : 211358
calculated : 205772
interpolated : 4578
hole : 1008
Updating StokesV wave model for patch inlet
Updating shallowWaterAbsorption wave model for patch outlet
Execution time for mesh.update() = 4.52 s
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 3.98197792029e-06, Final residual = 4.07322251154e-09, No Iterations 1
Phase-1 volume fraction = 0.670666212252 Min(alpha.water) = -1.8853175244e-09 Max(alpha.water) = 1.00000000011
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.670666211374 Min(alpha.water) = -1.86467766824e-09 Max(alpha.water) = 1.00000000011
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 3.98020495633e-06, Final residual = 4.21676115869e-09, No Iterations 1
Phase-1 volume fraction = 0.670663134233 Min(alpha.water) = -1.76164326384e-09 Max(alpha.water) = 1.00000000011
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.670663133389 Min(alpha.water) = -1.74120946968e-09 Max(alpha.water) = 1.00000000011
Updating StokesV wave model for patch inlet
DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.00244957787231, Final residual = 2.3253173653e-06, No Iterations 1
time step continuity errors : sum local = 2.69188594725e-06, global = 1.1047555752e-07, cumulative = 0.000263033779715
DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.70864207765e-06, Final residual = 9.40763534329e-10, No Iterations 27
time step continuity errors : sum local = 1.16722161447e-09, global = 3.42850852485e-10, cumulative = 0.000263034122566
PIMPLE: iteration 2
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 4.98681607805e-06, Final residual = 4.59646402518e-09, No Iterations 1
Phase-1 volume fraction = 0.670666197477 Min(alpha.water) = -1.88443190357e-09 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.670666198094 Min(alpha.water) = -1.86365959624e-09 Max(alpha.water) = 1.00000000002
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 4.00087672356e-06, Final residual = 3.98868344977e-09, No Iterations 1
Phase-1 volume fraction = 0.670663107707 Min(alpha.water) = -1.76076241918e-09 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.67066310705 Min(alpha.water) = -1.7400820768e-09 Max(alpha.water) = 1.00000000003
Updating StokesV wave model for patch inlet
DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.2371744886e-05, Final residual = 2.07562969402e-07, No Iterations 1
time step continuity errors : sum local = 1.92418107577e-07, global = -2.84068300766e-08, cumulative = 0.000263005715736
DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.89256586301e-07, Final residual = 9.55818249316e-10, No Iterations 20
time step continuity errors : sum local = 1.220193276e-09, global = 9.36248905173e-10, cumulative = 0.000263006651985
ExecutionTime = 1023.67 s ClockTime = 1031 s
Courant Number mean: 0.000410680261873 max: 1.28811016855
Interface Courant Number mean: 1.6589285851e-05 max: 0.286103027656
deltaT = 0.00155617762649
Time = 0.265764
PIMPLE: iteration 1
forces forces:
rho: rho
Not including porosity effects
Restraint verticalSpring1: attachmentPt - anchor (0.675369733081 0.34699687937 0.527411683655) spring length 0.924496725094 force (-0.0173056163364 -0.00889141838932 -0.0135143519195)
Restraint verticalSpring2: attachmentPt - anchor (0.675368891851 -0.347003120647 0.527413203961) spring length 0.924499320461 force (-0.0173074421758 0.00889252750295 -0.0135158335547)
Restraint verticalSpring3: attachmentPt - anchor (-0.675351108154 -0.347002279426 0.527413675564) spring length 0.924486282421 force (0.017297706036 0.00888773757957 -0.0135085981338)
Restraint verticalSpring4: attachmentPt - anchor (-0.675350266925 0.346997720591 0.527412155258) spring length 0.924483089431 force (0.0172954117007 -0.00888645304631 -0.0135067842687)
6-DoF rigid body motion
Centre of rotation: (9.38149026062e-06 -2.74160685292e-06 0.112418573328)
Centre of mass: (9.38149026062e-06 -2.74160685292e-06 0.112418573328)
Orientation: (0.999999999935 -9.86954951546e-06 -5.65881379727e-06 9.86944683233e-06 0.999999999787 -1.81454384919e-05 5.65899288337e-06 1.81453826413e-05 0.999999999819)
Linear velocity: (4.44786133845e-05 -2.65858492943e-05 0.00382088873648)
Angular velocity: (0.000298145326304 -0.000112279537068 5.65143921208e-05)
cellVolumeWeight : detected 2 mesh regions
zone:0 nCells:117810
zone:1 nCells:93548
Marking patch-cells on zone 0
Marking patch-cells on zone 1
After patch analysis : nCells : 211358
other : 156798
patch : 50666
overset: 3894
Creating mesh-to-mesh addressing for region0 and region0 regions using cellVolumeWeight
meshToMesh: Using AABBTree method
Overlap volume: 0.114317063947
cellVolumeWeight : detected 3 mesh regions after overset
Front : fraction:1 size:5229
Overset analysis : nCells : 211358
calculated : 205772
interpolated : 4578
hole : 1008
Updating StokesV wave model for patch inlet
Updating shallowWaterAbsorption wave model for patch outlet
Execution time for mesh.update() = 4.53 s
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 4.52165734097e-06, Final residual = 5.21939646806e-09, No Iterations 1
Phase-1 volume fraction = 0.670659612339 Min(alpha.water) = -1.62853632761e-09 Max(alpha.water) = 1.00000000015
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.670659611448 Min(alpha.water) = -1.60526824823e-09 Max(alpha.water) = 1.00000000015
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 4.51905319739e-06, Final residual = 5.45819623977e-09, No Iterations 1
Phase-1 volume fraction = 0.670656115556 Min(alpha.water) = -1.50043628112e-09 Max(alpha.water) = 1.00000000014
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.67065611471 Min(alpha.water) = -1.4775418421e-09 Max(alpha.water) = 1.00000000014
Updating StokesV wave model for patch inlet
DILUPBiCGStab: Solving for p_rgh, Initial residual = 0.0026856675826, Final residual = 2.62237078736e-06, No Iterations 1
time step continuity errors : sum local = 3.82284177183e-06, global = 1.39543331081e-07, cumulative = 0.000263146195316
DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.04279733626e-06, Final residual = 8.58826776316e-10, No Iterations 30
time step continuity errors : sum local = 1.38199049382e-09, global = 4.41989198394e-10, cumulative = 0.000263146637305
PIMPLE: iteration 2
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 5.68823736792e-06, Final residual = 5.85074040543e-09, No Iterations 1
Phase-1 volume fraction = 0.670659593399 Min(alpha.water) = -1.62732225817e-09 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.670659594239 Min(alpha.water) = -1.6038779044e-09 Max(alpha.water) = 1.00000000003
Updating StokesV wave model for patch inlet
smoothSolver: Solving for alpha.water, Initial residual = 4.54395149033e-06, Final residual = 5.18582133497e-09, No Iterations 1
Phase-1 volume fraction = 0.670656081221 Min(alpha.water) = -1.49911033101e-09 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.670656080587 Min(alpha.water) = -1.47589074005e-09 Max(alpha.water) = 1.00000000004
Updating StokesV wave model for patch inlet
DILUPBiCGStab: Solving for p_rgh, Initial residual = 2.73365601231e-05, Final residual = 2.38666793307e-07, No Iterations 1
time step continuity errors : sum local = 2.8324455281e-07, global = -4.15303920812e-08, cumulative = 0.000263105106913
DILUPBiCGStab: Solving for p_rgh, Initial residual = 3.27184089818e-07, Final residual = 9.95079080361e-10, No Iterations 24
time step continuity errors : sum local = 1.57914158151e-09, global = 1.24724361871e-09, cumulative = 0.000263106354157
ExecutionTime = 1028.81 s ClockTime = 1036 s
Courant Number mean: 0.00046656111542 max: 1.33424010549
Interface Courant Number mean: 1.87452560599e-05 max: 0.389957117035
deltaT = 0.00171179538914
Time = 0.267476
PIMPLE: iteration 1
forces forces:
rho: rho
Not including porosity effects
Restraint verticalSpring1: attachmentPt - anchor (0.675369805884 0.346996834016 0.527417549739) spring length 0.924500107782 force (-0.0173080260177 -0.00889265433381 -0.0135163825708)
Restraint verticalSpring2: attachmentPt - anchor (0.675368957102 -0.347003166002 0.527419110242) spring length 0.92450275461 force (-0.0173098882774 0.00889378460799 -0.0135178938529)
Restraint verticalSpring3: attachmentPt - anchor (-0.675351042903 -0.34700231723 0.527419596916) spring length 0.92448962705 force (0.017300085078 0.00888896178281 -0.0135106090297)
Restraint verticalSpring4: attachmentPt - anchor (-0.675350194122 0.346997682789 0.527418036413) spring length 0.924486377236 force (0.0172977501146 -0.00888765452276 -0.0135087625342)
Thanks for any help.
Regards,
Ahsan
|