[T.D.]

Hi Foamers,

I have this complicated boundary condition:

In my solver i calculates a flux vector filed J=function(c)
where c is a scalarField to be later solved in the same solver.

lets say this looks my solver:
.
.
.
and here i calculate the volVectorField J
J=something complicated as function of c;here J=div(pow(1-c/0.5)*Q);
where Q is a defined 2nd order tensor of constants.
.
.
after that in same solver i solve c equation
.
solve
(
fvm::ddt(c) + fvm::div(phi,c) == fvc::div(J)
);
.
.
etc


since i solve the "c" scalar equation, so of course i must define the BC for it

and here comes the problem since:
when i go to define the Patches BC for "c", i need in my case to set one Patch BC where the value of Flux J=0, so BC here for "c" must be an expression "Jy=0" (value of the component of vector J in direction of y is zero, that is normal to the surface xz).

where J is a vectorField calculated as function of c previously in the solver
as i mentioned.

How to do that?
Any Ideas?

Help please?