my_SimpleFoam solver problem
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 
Goutam, I think your question has already been answered. See post #22 and onward in this thread.

Quote:
Thanks Akidess... I have seen this after I reply to the post. Sorry. Thanks again. 
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 
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:

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 
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:

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"); Code:
blockMesh 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 
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:

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. 
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 
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. 
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 
2 Attachment(s)
Thanks Martin,
Herewith I have attached 0 and system directories. Thanks. 
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; } // ************************************************** *********************** // 
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 
Hello,
may a } be missing?? Quote:

Hi Maddalena,
you are right, the bracket is missing ;) oh, and a semicolon directly before the missing bracket... Thanks Martin 
Hi Maddalena,
Yes it worked. Thank you very much for your concern. Best regards, Suranga. 
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 18:16. 