[STutexas]

Hi

I am using OpenFoam to solve the magnetization decay problem and the boundary condition that I have looks like this:

D.n.delta(M)+b.M = 0

where,
D is scalar constant,
n is the surface normal to the boundary face
delta is the gradient operator
b is a constant
M is a vector field

I tried using:

boundaryField
{
fixedWalls
{
type groovyBC;
refValue uniform (0 0 1);
refGradient uniform (0 0 0);
valueFraction uniform 1;
value uniform (0 0 1);
valueExpression "Tinf";
gradientExpression "0";
fractionExpression "f";
evaluateDuringConstruction 0;
variables
4
(
"D=3e-09;"
"rho=10e-06;"
"Tinf=D/rho*snGrad(M);"
"f=I;"
)
;
timelines (
);
lookuptables (
);
}
}
but it gives me an error. Pleas let me know how I can code my boundary condition.

Thank you.