[AB08]

Hi people......I need your help ,,,its a bit urgent,,,,,,,

I have a TEqn defined as --

volScalarField mu1=0.62*1.82*pow(a,2)*g*pow(T,2)*pow(0.68-T,-1);
volScalarField mu2=0.41*pow(a,2)*g*T;
volScalarField mu3=0.41*pow(a,2)*T*T;
Info<< "Solving for T " << nl << endl;
solve
(
fvm::div(phi, T)
//- fvm::laplacian(nu, T)
- fvm::laplacian(mu1, T)
- fvm::laplacian(mu2, T)
==
fvc::laplacian(mu3, g)
);

My UEqn is---

tmp<fvVectorMatrix> UEqn
(
fvm::div(phi, U)
- fvm::laplacian(nu*pow((1-T/0.68),(-1.82)), U)
);

UEqn().relax();

fvOptions.constrain(UEqn());

solve(UEqn() == -fvc::grad(p));

So the value of T has to be between 0 and 0.68......How can I do this in OpenFoam..???

[nimasam]

its some how like this

dimensionedScalar MaxT ("MaxT",dimT, 0.68);
dimensionedScalar MinT ("MinT",dimT, 0.0);

T = max (min(T,MaxT), MinT);

[AB08]

Thanx .. but plz tell me where exactly do i put this.......... plz dont mind i m very new to openFoam...