[fvallejog]

Hi OpenFOAMERS

Would you help me with a tips of C++ plus OpenFOAM????

My problem is:

I created a new Solver based on "scalarTransportFoam": The file ".C" is:

while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
solve
(
fvm::ddt(CT)
+ fvm::div(phi, CT)
- fvm::laplacian(DT, CT)

==
fvOptions(CT)
);

The solver compile without problems.

But, my model do not finish yet!

I want to put a new set of equations, like this: The difficult part is that the equations is only true in the region where U=0, it means in the surface region named "fixedWalls".

if (U.component(0))=0;
solve
(
fvm::ddt(CT)
+ fvm::div(phi, CT)
- fvm::laplacian(DT, CT)
- umax*X/Yp*CT/(Ks+CT)
==
fvOptions(CT)
);

Please, I wait any comment.

Fidel