CFD Online Discussion Forums - How to bound fields...

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)

- OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)

- - How to bound fields... (https://www.cfd-online.com/Forums/openfoam-solving/120106-how-bound-fields.html)

How to bound fields...

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..???

its some how like this ;)

Code:

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

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



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

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