[KW95]

Hi, I have a case where I have two rotating bodies in two different dynamic meshes, so for each one I use the AMI method with the dynamicMotionSolver.

According to this site: https://openfoamwiki.net/index.php/P...onSolverFvMesh
This solver has the ability to report for each time step in the uniform file inside each timestep file.

However when I run it with two bodies this no longer works. Anyone here know why this happens & how I can fix it?

Sorry I'm new to OpenFOAM.


Here you have the DynamicMeshDictFile I'm using:


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2306 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicMultiMotionSolverFvMesh;
motionSolverLibs ("libfvMotionSolvers.so");

reportToFile on;

dynamicMultiMotionSolverFvMeshCoeffs
{
generator1
{

dynamicFvMesh dynamicMotionSolverFvMesh;
cellZone innerCylinder1;
motionSolverLibs ("libsixDoFRigidBodyMotion.so");

solver sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{
patches (generator1);
innerDistance 51;
outerDistance 51.5;
rho rhoInf;
rhoInf 1.293;
mass 994523.68572;
momentOfInertia (72280656.68796 124473124.00092 72280656.68568);
centreOfMass (0 -1.395 0);
g ( 0 0 -9.8065 );
// orientation
// (
// 1 0 0
// 0 1 0
// 0 0 1
// );
velocity (0 0 0);
acceleration (0 0 0);
angularMomentum (0 262638291.6419412 0);
torque (0 0 0);
report on;
reportToFile on;
accelerationRelaxation 1.0; //1.0;
accelerationDamping 0.3; //1.0;

solver
{
type Newmark;
}

constraints
{
zAxis
{
sixDoFRigidBodyMotionConstraint axis;
axis (0 1 0);
}
fixedPt
{
sixDoFRigidBodyMotionConstraint point;
point (0 0.1 0);
}
}
restraints
{
sixDoFRigidBodyMotionRestraint sphericalAngularDamper;
coeff 0.02;
}
}
}
generator2
{

dynamicFvMesh dynamicMotionSolverFvMesh;
cellZone innerCylinder2;
motionSolverLibs ("libsixDoFRigidBodyMotion.so");

solver sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{
patches (generator2);
innerDistance 51;
outerDistance 51.5;
rho rhoInf;
rhoInf 1.293;
mass 994523.68572;
momentOfInertia (72280656.68796 124473124.00092 72280656.68568);
centreOfMass (0 -701.395 0);
g ( 0 0 -9.8065 );
// orientation
// (
// 1 0 0
// 0 1 0
// 0 0 1
// );
velocity (0 0 0);
acceleration (0 0 0);
angularMomentum (0 262638291.6419412 0);
torque (0 0 0);
report on;
reportToFile on;
accelerationRelaxation 1.0; //1.0;
accelerationDamping 0.3; //1.0;

solver
{
type Newmark;
}

constraints
{
zAxis
{
sixDoFRigidBodyMotionConstraint axis;
axis (0 1 0);
}
fixedPt
{
sixDoFRigidBodyMotionConstraint point;
point (0 -700.1 0);
}
}
restraints
{
sixDoFRigidBodyMotionRestraint sphericalAngularDamper;
coeff 0.02;
}
}
}
}

// ************************************************** *********************** //