CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   ReactingFoam & Dynamic Mesh (https://www.cfd-online.com/Forums/openfoam-solving/120119-reactingfoam-dynamic-mesh.html)

yaqb July 1, 2013 05:57

ReactingFoam & Dynamic Mesh
 
Hello All !


I am trying to run reactingFoam with dynamic Mesh.
My basic idea is very simple - I would like to simulate the regression rate of the fuel during the burning process.

I set dynamicMeshDict

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ( "libfvMotionSolvers.so" );

solver velocityLaplacian;

diffusivity uniform;
// ************************************************** *********************** //


,and run it, but all the time get the same error:

--> FOAM FATAL ERROR:
Different dimensions for +=
dimensions : [1 0 -1 0 0 0 0] = [0 3 -1 0 0 0 0]


From function dimensionSet::operator+=(const dimensionSet&) const
in file dimensionSet/dimensionSet.C at line 241.

FOAM aborting



Do you have any idea where can I search for the mistake?
In BC there is no parameter with dimensions [1 0 -1 0 0 0 0] , so I am not sure what does the solver wants from me?

nimasam July 1, 2013 11:28

OpenFOAM analyses for dimension compatibility, so before it solves an equation, it investigates whether parameters have the same dimension of not ;) ,

so check it with out dynamic mesh, maybe your input dimension for p, u r something else is wrong, i think its not related to dynamic mesh

yaqb July 1, 2013 12:06

Thanks for the replay!

In fact if I just run pure reactingFoam case than it works perfect however when I try to use reactingFoamDynamic it crashes.
There is only one additional BC added - pointMotionU - responsible for the mesh motion.

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class pointVectorField;
object pointMotionU;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);



boundaryField
{
fuel
{
type fixedValue;
value uniform (0.1 0 0);
}
air
{
type fixedValue;
value uniform (0 0 0);
}
outlet
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;
}
Up
{
type fixedValue;
value uniform (0 0 0);
}
Bottom
{
type fixedValue;
value uniform (0.1 0 0);
}
}


And here it is m/s so everything seems to be ok. Any idea where I should search for the mistake?

nimasam July 1, 2013 12:17

post your test case here ;)

yaqb July 1, 2013 12:30

1 Attachment(s)
Here it is


All times are GMT -4. The time now is 01:06.