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

attachmentPt for floating body at each time step for overInterDyMFoam solver

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 24, 2022, 11:13
Question 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
mahsankhan is on a distinguished road
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
mahsankhan is offline   Reply With Quote

Reply

Tags
attachmentpt, floating body, linear spring, openfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
bash script for pseudo-parallel usage of reconstructPar kwardle OpenFOAM Post-Processing 42 May 7, 2024 23:17
[solidMechanics] Support thread for "Solid Mechanics Solvers added to OpenFOAM Extend" bigphil OpenFOAM CC Toolkits for Fluid-Structure Interaction 686 December 22, 2022 09:10
fluent divergence for no reason sufjanst FLUENT 2 March 23, 2016 16:08
AMI interDyMFoam for mixer nu problem danny123 OpenFOAM Programming & Development 8 September 6, 2013 02:34
dynamic Mesh is faster than MRF???? sharonyue OpenFOAM Running, Solving & CFD 14 August 26, 2013 07:47


All times are GMT -4. The time now is 21:35.