[spenceme]

I am trying to simulate a porous region that is rotating at the same omega as my rotating region using a rotating motion and AMI. The simulation runs fine, but does not give the behavior though the porous region as I would expect/like. I want the porous region to 'hold' the fluid it contains. I.e. with a very low porosity it would act like a pseudo wall and approach being non-existent with a very high porosity.

I have a similar steady state case running using MRF and the branch feature-porous-mrf. MRF model presents the anticipated velocities.

Two pictures attached: 1) the 'incorrect' behavior from pimplefoam currently, 2) a more correct expected velocity field coming from the 'feature-porous-mrf' branch running porousSimpleFoam.

This can be thought of as a simple annulus of a porous region that is rotated in a fluid. The porosity should pull the fluid along and eventually cause a flow similar to an impeller due to centrifugal forces.

I am currently using pimpleFoam with fvOptions to establish the porous region. Any suggestions if openfoam (v2112 currently) supports my needs, or how to change the setup/solver to reflect the porous region rotating?

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

porosity
{
	type            explicitPorositySource;

	explicitPorositySourceCoeffs
	{
		selectionMode   cellZone;
		cellZone        porusBlockage;
		type            DarcyForchheimer;

	//	D	100;
	//	D	200;
		D	5e3;
	//	D	1000;
	
		d	($D $D $D);
//		f	(0 0 0);
		f	($D $D $D);

		coordinateSystem
		{
			type	cartesian;
			origin  (0 0 0);
			coordinateRotation
			{
				type axesRotation; //local cartesian coordinates
				e1      (1 0 0);
				e2      (0 1 0);
			}
		}
	}
}



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

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

dynamicFvMesh   dynamicMotionSolverFvMesh;

motionSolverLibs (fvMotionSolvers);

motionSolver    solidBody;

cellZone        innerCylinderSmall;

solidBodyMotionFunction  rotatingMotion;

origin      (0 0 0);
axis        (0 1 0);
omega       52; // rad/s (496.56 rpm)


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

[Tempest]

Did you solve this issue?

[spenceme]

Originally Posted by Tempest

Did you solve this issue?

No, i have not. It appears the porous region is assumed static in position. i.e. as if the 'brush' was stationary in my pictures.