CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   How to bound fields... (http://www.cfd-online.com/Forums/openfoam-solving/120106-how-bound-fields.html)

AB08 July 1, 2013 02:29

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

nimasam July 1, 2013 11:04

its some how like this ;)
Code:

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

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


AB08 July 1, 2013 13:01

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


All times are GMT -4. The time now is 07:32.