attachmentPt for floating body at each time step for overInterDyMFoam solver

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

 January 24, 2022, 11:13 attachmentPt for floating body at each time step for overInterDyMFoam solver #1 New Member   Ahsan Join Date: Nov 2019 Location: Bologna, Italy Posts: 27 Rep Power: 6 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

 Tags attachmentpt, floating body, linear spring, openfoam