
[Sponsors] 
March 1, 2012, 06:32 
my_SimpleFoam solver problem

#41 
Senior Member
Goutam Saha
Join Date: Dec 2011
Location: UK
Posts: 131
Rep Power: 6 
I am using OpenFOAM2.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 

March 1, 2012, 06:53 

#42 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,130
Rep Power: 20 
Goutam, I think your question has already been answered. See post #22 and onward in this thread.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. *Join the OpenFOAM stackexchange Q&A site: http://area51.stackexchange.com/prop...oHPxcPqde7HtA2 

June 19, 2012, 04:13 

#44 
Member
anonymous
Join Date: Mar 2012
Posts: 45
Rep Power: 6 
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 Wnounusedparameter Woldstylecast Wnonvirtualdtor O3 DNoRepository ftemplatedepth100 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 

August 7, 2012, 14:27 

#45  
Senior Member
Srinath Madhavan (a.k.a pUl)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13 
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 steadystate 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:


August 7, 2012, 14:41 

#46 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 12 
Hi Srinath,
you should modify the viscosity model to be dependent on your scalar T. You can find the viscosity models here: OpenFOAM2.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 

August 7, 2012, 14:47 

#47  
Senior Member
Srinath Madhavan (a.k.a pUl)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13 
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:


August 7, 2012, 16:46 

#48 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 12 
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"); Code:
blockMesh srinathSimpleFoam 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 

August 7, 2012, 17:40 

#49  
Senior Member
Srinath Madhavan (a.k.a pUl)
Join Date: Mar 2009
Location: Edmonton, AB, Canada
Posts: 703
Rep Power: 13 
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:


October 21, 2012, 04:50 

#50 
Member
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 7 
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/naren2.1.1/run/filmCoolingCourse1/system/fvSolution" Any suggestion to get rid of this. Best regards, Suranga. 

October 21, 2012, 04:58 

#51 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 12 
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 

October 21, 2012, 09:47 

#52 
Member
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 7 
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 1e06; relTol 0.1; } T { solver BICCG; preconditioner DILU; tolerance 1e07; relTol 0.01 U { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.01; } k { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.01; } epsilon { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.01; } R { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.01; } nuTilda { solver PBiCG; preconditioner DILU; tolerance 1e05; relTol 0.01; } } SIMPLE { nNonOrthogonalCorrectors 0; residualControl { p 1e5; U 1e5; T 1e5; "(kepsilonomega)" 1e3; } } 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. 

October 21, 2012, 09:56 

#53 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 12 
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 

October 21, 2012, 10:35 

#54 
Member
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 7 
Thanks Martin,
Herewith I have attached 0 and system directories. Thanks. 

October 21, 2012, 10:41 

#55 
Member
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 7 
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 ] 1e05; DT DT [ 0 2 1 0 0 0 0 ] 2e05; CrossPowerLawCoeffs { nu0 nu0 [ 0 2 1 0 0 0 0 ] 1e06; nuInf nuInf [ 0 2 1 0 0 0 0 ] 1e06; 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 ] 1e06; nuInf nuInf [ 0 2 1 0 0 0 0 ] 1e06; k k [ 0 0 1 0 0 0 0 ] 0; n n [ 0 0 0 0 0 0 0 ] 1; } // ************************************************** *********************** // 

October 21, 2012, 12:47 

#56 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 12 
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 

October 23, 2012, 02:48 

#57  
Senior Member
maddalena
Join Date: Mar 2009
Posts: 436
Rep Power: 14 
Hello,
may a } be missing?? Quote:


October 23, 2012, 03:53 

#58 
Senior Member
Martin
Join Date: Oct 2009
Location: Aachen, Germany
Posts: 252
Rep Power: 12 
Hi Maddalena,
you are right, the bracket is missing ;) oh, and a semicolon directly before the missing bracket... Thanks Martin 

October 24, 2012, 23:17 

#59 
Member
Suranga Dharmarathne
Join Date: Jan 2011
Location: TX, USA
Posts: 39
Rep Power: 7 
Hi Maddalena,
Yes it worked. Thank you very much for your concern. Best regards, Suranga. 

November 7, 2012, 16:40 

#60 
New Member
Warren Lamont
Join Date: Nov 2012
Posts: 5
Rep Power: 5 
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 

Tags 
simplefoam, temperature 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Calculation of the Governing Equations  Mihail  CFX  7  September 7, 2014 06:27 
Adding a new temperature dependent viscositymodel?  dgadensg  OpenFOAM Programming & Development  10  May 22, 2010 05:47 
Adding temperature equation in settlingFoam  sachin  OpenFOAM Running, Solving & CFD  2  March 31, 2010 03:21 
Adding temperature field to InterFoam  yapalparvi  OpenFOAM Running, Solving & CFD  8  October 14, 2009 20:18 
Adding coriolis forces in simplefoam  Xabi  OpenFOAM Running, Solving & CFD  1  April 24, 2009 04:43 