CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Programming & Development (https://www.cfd-online.com/Forums/openfoam-programming-development/)
-   -   adding temperature to simpleFoam (https://www.cfd-online.com/Forums/openfoam-programming-development/84480-adding-temperature-simplefoam.html)

Goutam March 1, 2012 05:32

my_SimpleFoam solver problem
 
I am using OpenFOAM-2.1.0. I have tried to use my_simpleFoam. I did wclean and wmake and got 1 error.
readSIMPLEControls.H: no such file or directory.

I checked in src/finitevloume/lnInclude/...

There is no such file.

Could you please tell me how to solve this problem ...

Thanks

akidess March 1, 2012 05:53

Goutam, I think your question has already been answered. See post #22 and onward in this thread.

Goutam March 1, 2012 06:05

Quote:

Originally Posted by akidess (Post 347051)
Goutam, I think your question has already been answered. See post #22 and onward in this thread.


Thanks Akidess...

I have seen this after I reply to the post. Sorry.

Thanks again.

libia87 June 19, 2012 04:13

Hi!
I'm trying to add concentration following the tutorial "how to add temperature to icoFoam", but I'm using buoyantboussinesqPimpleFoam and it doesn't work, anyone knows why??
when I do the WMAKE appears that error:

Making dependency list for source file my_buoyantBoussinesqPimpleFoam.C
could not open file readTransportProperties.H for source file my_buoyantBoussinesqPimpleFoam.C
SOURCE=my_buoyantBoussinesqPimpleFoam.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I../buoyantBoussinesqSimpleFoam -I/opt/openfoam210/src/finiteVolume/lnInclude -I/opt/openfoam210/src/turbulenceModels -I/opt/openfoam210/src/turbulenceModels/incompressible/RAS/lnInclude -I/opt/openfoam210/src/transportModels -I/opt/openfoam210/src/transportModels/incompressible/singlePhaseTransportModel -IlnInclude -I. -I/opt/openfoam210/src/OpenFOAM/lnInclude -I/opt/openfoam210/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/my_buoyantBoussinesqPimpleFoam.o
In file included from my_buoyantBoussinesqPimpleFoam.C:61:
createFields.H:47:41: error: readTransportProperties.H: No existe el fichero o el directorio
In file included from my_buoyantBoussinesqPimpleFoam.C:61:
createFields.H: In function ‘int main(int, char**)’:
createFields.H:52: error: ‘laminarTransport’ was not declared in this scope
createFields.H:64: error: ‘beta’ was not declared in this scope
createFields.H:64: error: ‘TRef’ was not declared in this scope
In file included from my_buoyantBoussinesqPimpleFoam.C:86:
TEqn.H:2: error: ‘Prt’ was not declared in this scope
TEqn.H:5: error: ‘Pr’ was not declared in this scope
/opt/openfoam210/src/finiteVolume/lnInclude/readTimeControls.H:38: warning: unused variable ‘maxDeltaT’
make: *** [Make/linux64GccDPOpt/my_buoyantBoussinesqPimpleFoam.o] Error 1

msrinath80 August 7, 2012 14:27

Hello Martin,

Thanks very much for the clear explanation and the modified sources. Your mySimpleFoam solver is indeed very useful. I wonder however, if you can point me in the right direction for my task (described below):

Basically your solver is a steady-state conjugate heat transfer solver, in that it solved for the flow field first, then advects temperature using that flow field. The energy equation has a diffusion term and a source term (viscous dissipation).

Is it possible to modify the solver to do the following instead. I want to make the viscosity of the flowing material a function of the scalar T. Basically, if the scalar T is within a certain range, I want the material to have viscosity X, and if it is within another range, I want it to have viscosity Y. I of course assume that X will be roughly the same order as Y (i.e. no sharp discontinuities).

I look forward to your response.

Thanks in advance for your help!

Best Regards,
Srinath

Quote:

Originally Posted by MartinB (Post 329943)
Hi Nickolas,

you can try the attached solver (for OpenFOAM 2.0.x) and case.

Comments are added to the code and the case.

Increase of temperature is calculated via controlDict (4.585 K) and nearly matches to the appropriate formula for high viscous newtonian fluids (4.73 K).

Ah, and you can access the viscosity for your formula with "laminarTransport.nu()".

Martin


MartinB August 7, 2012 14:41

Hi Srinath,

you should modify the viscosity model to be dependent on your scalar T. You can find the viscosity models here:
OpenFOAM-2.1.x/src/transportModels/incompressible/viscosityModels/
If your scalar T is temperature, then there is Arrhenius shift or WLF shift to describe temperature dependent viscosities and which can be easily implemented in OpenFOAM.
Can you post your specific function for nu(T) or nu(T, shear rate)?

Martin

msrinath80 August 7, 2012 14:47

Hello again Martin,

Thanks very much for the really prompt reply. Here is my exact situation. I need the viscosity in the flow problem to change depending on my value of my scalar T. Think of 'T' as being some kind of species concentration as opposed to temperature.

If a<=T<=b, then nu = X

else if

b<=T<=c then nu = Y

Is this feasible?

Thanks once again for your help!

Best Regards,
Srinath

Quote:

Originally Posted by MartinB (Post 375879)
Hi Srinath,

you should modify the viscosity model to be dependent on your scalar T. You can find the viscosity models here:
OpenFOAM-2.1.x/src/transportModels/incompressible/viscosityModels/
If your scalar T is temperature, then there is Arrhenius shift or WLF shift to describe temperature dependent viscosities and which can be easily implemented in OpenFOAM.
Can you post your specific function for nu(T) or nu(T, shear rate)?

Martin


MartinB August 7, 2012 16:46

2 Attachment(s)
Hi Srinath,

here is a template for fine tuning.

In srinathSimpleFoam.tar.gz you find a solver based on simpleFoam with scalar transport of T. Furthermore there is a viscosity model with name "scalarDependentViscosity" which must be compiled with "wmake libso" (or use the Allwmake script). The results are a user library for your new viscosity model and the solver itself. The srinathSimpleFoam solver already links against the viscosity model. To use it in other solvers or utilities you can include it in the controlDict with
Code:

libs ("libuserscalarDependentViscosity.so" "libOpenFOAM.so");
There is a test case attached, too. Run it with the standard sequence
Code:

blockMesh
srinathSimpleFoam

Right now only nu1 and nu2, as well as b in constant/transportProperties are used for the viscosity model. You should customize the formula in "srinathSimpleFoam/scalarDependentViscosity/scalarDependentViscosity.C" line 57ff.

To keep the solvers convergence stable you might want to change the viscosity law a bit, so that the transition from nu1 to nu2 is a bit smoother, and not "binary" as it is now.

Have fun

Martin

msrinath80 August 7, 2012 17:40

Wow Martin! That's awesome. Let me look into the code and understand it and modify it (ultimately I want this to work with BirdCarreau). I will follow up with you on this thread with my developments. Thanks once again for the really quick response!

And yes, I will use some kind of a smoothed heaviside transition between the values :-)

Thanks again for your help!

Best Regards,
Srinath

Quote:

Originally Posted by MartinB (Post 375899)
Hi Srinath,

here is a template for fine tuning.

In srinathSimpleFoam.tar.gz you find a solver based on simpleFoam with scalar transport of T. Furthermore there is a viscosity model with name "scalarDependentViscosity" which must be compiled with "wmake libso" (or use the Allwmake script). The results are a user library for your new viscosity model and the solver itself. The srinathSimpleFoam solver already links against the viscosity model. To use it in other solvers or utilities you can include it in the controlDict with
Code:

libs ("libuserscalarDependentViscosity.so" "libOpenFOAM.so");
There is a test case attached, too. Run it with the standard sequence
Code:

blockMesh
srinathSimpleFoam

Right now only nu1 and nu2, as well as b in constant/transportProperties are used for the viscosity model. You should customize the formula in "srinathSimpleFoam/scalarDependentViscosity/scalarDependentViscosity.C" line 57ff.

To keep the solvers convergence stable you might want to change the viscosity law a bit, so that the transition from nu1 to nu2 is a bit smoother, and not "binary" as it is now.

Have fun

Martin


sdharmar October 21, 2012 04:50

Hi Martin,
Thank you for sharing your mySimpleFom solver and the case. I tried to run in one of my cases but I keep on getting error massage

keyword SIMPLE is undefined in dictionary "/home/naren/OpenFOAM/naren-2.1.1/run/filmCoolingCourse1/system/fvSolution"

Any suggestion to get rid of this.

Best regards,
Suranga.

MartinB October 21, 2012 04:58

Hi Suranga,

in your fvSolution the dictionary entry for "SIMPLE" is missing. I think you run your case with a PIMPLE or PISO algorithm otherwise, so to use the mySimpleFoam solver you must add parameters for the SIMPLE algorithm.

You can have a look at the attached fvSolution file at the case attached to the mySimpleFoam solvers.

You might need to edit your fvSchemes, too: you must switch the ddtSchemes to steadyState, and probably you might want to change controlDict to use delta 1 as an iteration counter.

Martin

sdharmar October 21, 2012 09:47

Hi Martin,

Thank you for the reply. I think I have SIMPLE written in my fvSolution dictionary. Please be kind enough to go through the listing.

solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
}
T
{
solver BICCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0.01

U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;

residualControl
{
p 1e-5;
U 1e-5;
T 1e-5;
"(k|epsilon|omega)" 1e-3;
}
}

relaxationFactors
{
fields
{
p 0.3;
T 0.7;
}
equations
{
U 0.7;
k 0.7;
epsilon 0.7;
R 0.7;
nuTilda 0.7;
}
}


// ************************************************** *********************** //

Still trying to figure out what's wrong with this. Thanks for your time in advance.

MartinB October 21, 2012 09:56

I can't see any obvious mistake... can you upload the case?
Or at least: 0/*, system/*, constant/transportProperties, altogether in tar.gz archive...

Martin

sdharmar October 21, 2012 10:35

2 Attachment(s)
Thanks Martin,

Herewith I have attached 0 and system directories.

Thanks.

sdharmar October 21, 2012 10:41

Sorry I couldn't attach the transportProperties to the previous reply. Here it is.


/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

transportModel Newtonian;

nu nu [ 0 2 -1 0 0 0 0 ] 1e-05;
DT DT [ 0 2 -1 0 0 0 0 ] 2e-05;

CrossPowerLawCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
m m [ 0 0 1 0 0 0 0 ] 1;
n n [ 0 0 0 0 0 0 0 ] 1;
}

BirdCarreauCoeffs
{
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
k k [ 0 0 1 0 0 0 0 ] 0;
n n [ 0 0 0 0 0 0 0 ] 1;
}


// ************************************************** *********************** //

MartinB October 21, 2012 12:47

I can't find any problem in your files... they seem to be fine for me.

You can send me your case, if you like, so I can try to run it here... I'll send you my eMail adress via boardmail...

Martin

maddalena October 23, 2012 02:48

Hello,
may a } be missing??

Quote:

Originally Posted by sdharmar (Post 387758)
solvers
{
p
{
solver PCG;
preconditioner DIC;
tolerance 1e-06;
relTol 0.1;
}
T
{
solver BICCG;
preconditioner DILU;
tolerance 1e-07;
relTol 0.01
}
U
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

k
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

epsilon
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

R
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}

nuTilda
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.01;
}
}

SIMPLE
{
nNonOrthogonalCorrectors 0;

residualControl
{
p 1e-5;
U 1e-5;
T 1e-5;
"(k|epsilon|omega)" 1e-3;
}
}

relaxationFactors
{
fields
{
p 0.3;
T 0.7;
}
equations
{
U 0.7;
k 0.7;
epsilon 0.7;
R 0.7;
nuTilda 0.7;
}
}


MartinB October 23, 2012 03:53

Hi Maddalena,

you are right, the bracket is missing ;-)
oh, and a semicolon directly before the missing bracket...

Thanks

Martin

sdharmar October 24, 2012 23:17

Hi Maddalena,

Yes it worked. Thank you very much for your concern.

Best regards,
Suranga.

wlamont November 7, 2012 15:40

Hey Martin B,

I want to implement a simpler version of your srinathSimpleFoam example. I want nu to be a function of T and p (T for me is temperature) and I want to implement something like 0.0000171*pow(T/273,0.7)/(p/(8314/28.96*T). Can you modify your code to do this? I tried modifying your code but I get errors. I think it is my lack of familiarity.

Thanks,

Warren


All times are GMT -4. The time now is 01:14.