|
[Sponsors] |
Force Oscillations in sloshing simulation with dynamicMeshDict (6DoF) |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
February 9, 2024, 08:47 |
Force Oscillations in sloshing simulation with dynamicMeshDict (6DoF)
|
#1 |
New Member
thr
Join Date: Dec 2023
Posts: 2
Rep Power: 0 |
Hello,
im try to simulate sloshing in a closed tank. I use OF11 and based on the tutorial sloshingTank6DoF i built up my own simulation. In the end i want the forces acting on the tank due to the liquid sloshing. When i read the data from the postProcessing folder i get really confusing Oscillations. I have no clue why this is happening. I have tried several different options in fvSchemes and fvSolution. controlDict: application foamRun; solver incompressibleVoF; startFrom latestTime; startTime 0; stopAt endTime; endTime 1; deltaT 0.0001; writeControl adjustableRunTime; writeInterval 0.1; purgeWrite 0; writeFormat binary; writePrecision 9; writeCompression off; timeFormat general; timePrecision 9; runTimeModifiable yes; adjustTimeStep yes; maxCo 0.4; maxAlphaCo 0.4; maxDeltaT 0.01; functions { forces { type forces; libs ("libforces.so"); patches ("sloshingCylinder"); log on; writeControl adjustableRunTime; writeInterval 0.001; rhoInf 1000; rho rhoInf; CofR (0.4 -0.4 0.1); // Point for calculation of moments } } fvSolution solvers { alpha.water { nAlphaCorr 1; nAlphaSubCycles 3; } "pcorr.*" { solver PCG; preconditioner { preconditioner GAMG; tolerance 1e-05; relTol 0; smoother DICGaussSeidel; cacheAgglomeration no; } tolerance 1e-05; relTol 0; maxIter 100; } p_rgh { solver GAMG; tolerance 1e-08; relTol 0.01; smoother DIC; } p_rghFinal { solver PCG; preconditioner { preconditioner GAMG; tolerance 2e-09; relTol 0; nVcycles 2; smoother DICGaussSeidel; nPreSweeps 2; } tolerance 2e-09; relTol 0; maxIter 20; } "(U|k|epsilon|omega).*" { solver smoothSolver; smoother GaussSeidel; tolerance 1e-06; relTol 0; nSweeps 1; } } PIMPLE { momentumPredictor no; nCorrectors 5; nNonOrthogonalCorrectors 1; correctPhi no; pRefPoint (0.4 -0.4 0.1); pRefValue 1e5; } relaxationFactors { equations { ".*" 1; } } dynamicMeshDict mover { type motionSolver; libs ("libfvMeshMovers.so" "libfvMotionSolvers.so"); motionSolver solidBody; select all; solidBodyMotionFunction sixDoFMotion; CofG (0.4 -0.4 0.1); translationRotation table; //interpolationScheme spline; file "$FOAM_CASE/constant/6DoF.dat"; } 6DoF.dat is a table with a frequency of 800hz. simple acceleration in x direction with 9.81m/s^2. I have attached a plot of the force. Its the same for Fx,Fy,Fz. Every 10nth data point drops. The sloshing in Paraview look good. I think this could be a problem with the dynamictMeshDict, the forcesFunction or the adjustableTimeStep. Best, T |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Drag Force Ratio for Flat Plate | Rob Wilk | Main CFD Forum | 40 | May 10, 2020 04:47 |
Z-dimension for 2D 6DoF simulation | Fieschi | OpenFOAM Pre-Processing | 1 | March 29, 2020 23:01 |
Incorrect lift force in simpleFoam airfoil simulation | olov | OpenFOAM Running, Solving & CFD | 2 | May 18, 2019 06:51 |
Massless 6DOF Simulation | Mojtaba.a | OpenFOAM Running, Solving & CFD | 1 | February 14, 2016 07:54 |
VOF Sloshing - beating pattern in force output? | Brian | FLUENT | 0 | May 18, 2006 10:59 |