[emaportr]

Hi,


I am trying to get a paper or tutorial to simulate a cylinder with rotation and translation inside a pipe using openFoam.


Tobias Holzmann made this and publish this video.

https://www.youtube.com/watch?v=XC-_F42CZjg


Anyone has an idea on the strategy to make this happened.


there is a forum thread on this matter without any response on this matter but has already many years that was posted. I hope there is new info that I could use to model my case,


thanks

[JBeilke]

ACME is not really working as Tobias wrote in the description of the video. Maybe it can be done with overset.

[emaportr]

http://ibb.co/sJ25bpC


I was able to do rotation and translation using oversetMesh. Works pretty well when the direction of rotating and translating are both counterclockwise or clockwise. The issue is that I need to translate counterclockwise and rotate clockwise, and seems that openfoam does not understand that with the dynamicMeshDict I wrote, only rotate or translate (not both).


Any help would be very appreciated



this is the dynamicMeshDict:


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

dynamicFvMesh dynamicOversetFvMesh;

dynamicOversetFvMeshCoeffs
{
}

motionSolverLibs ("libsixDoFRigidBodyMotion.so");


solver multiSolidBodyMotionSolver;

multiSolidBodyMotionSolverCoeffs
{
rotor
{

cellZone rotor;
CofG (0.0045005 0 0.011);
mass 0.1; // Mass of the object in kg
velocity (0 0 0);
g (0 0 -9.8065)

// Translational motion (clockwise) (6DoF table)
solidBodyMotionFunction tabulated6DoFMotion;
timeDataFileName "../pEcc/stator/constant/6DoFCylin.dat";

// Rotational motion (counterclockwise with dynamic origin)
solidBodyMotionFunction rotatingMotion;
rotatingMotionCoeffs
{
origin (0.0045005 0 0.011); // Fixed origin of rotation
axis (0 0 1); // Axis of rotation (z-axis)
omega 69.8133; // Angular velocity in rad/s (counterclockwise)
}

solver
{
type CrankNicolson;
}

restraints
{
// Define any restraints here if necessary
}

// Function to update the center of rotation dynamically
originFunction // centreOfRotationFunction changes while translating
{
type codeStream;
code
#{
// Access to time and angular velocity
const Time& time = mesh().time();
scalar t = time.value();
scalar omega = 209.44; // Angular velocity in rad/s

// Calculate new center of rotation based on the given formula
scalar x = cos(omega * t) * 0.0045005;
scalar y = sin(omega * t) * 0.0045005;
scalar z = 0.0;

point newOrigin = point(x, y, z);

// Output the new center of rotation
os << newOrigin;
#};
}
}
}