# Need some details in SimpleFOAM

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 December 16, 2010, 12:21 Need some details in SimpleFOAM #1 Member   Vishal Jambhekar Join Date: Mar 2009 Location: University Stuttgart, Stuttgart Germany Posts: 90 Blog Entries: 1 Rep Power: 9 Hi All, I have set up problem to simulate flow over a city using exp prof at inlet. Unfortunately I dnot know the Char. Length for Epsilon calculation. However I have calculated it using values of using ratio of (nut/nu). but in tranport properties i see following : nu0, nuInf m and n . I dont understand which value means what.....??? And as I run the simulation I cant even even notice the solver solving for nuTIlda. can someone tell me does the solver solves for nuTilda in SimpleFoam... __________________ Cheers, Vishal Jambhekar... "Simulate the way ahead......!!!"

 December 19, 2010, 09:39 #2 Senior Member   David Boger Join Date: Mar 2009 Location: Penn State Applied Research Laboratory Posts: 146 Rep Power: 9 It sounds to me like the coefficients you are referring to belong to one of the non-Newtonian viscosity models. They should be listed within a dictionary with a name such as CrossPowerLaw, in which case the corresponding source code shows how the parameters are used: Code: ```Foam::viscosityModels::CrossPowerLaw::calcNu() const { return (nu0_ - nuInf_)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_; }``` But whether you are even using these or not depends on which transport model you have selected in transportProperties; e.g. if you are using Code: `transportModel Newtonian;` then the coefficients you are asking about are irrelevant. As far as the solution of nuTilda, the only models I know of that solve for nuTilda are the RAS and LES versions of the Spalart-Allmaras turbulence model, so the code will only solve for that variable if you have selected Spalart-Allmaras as the turbulence model. For example, in your RASProperties file, you might specify Code: `RASModel SpalartAllmaras;` __________________ David A. Boger

 December 23, 2010, 05:48 #3 Member   Vishal Jambhekar Join Date: Mar 2009 Location: University Stuttgart, Stuttgart Germany Posts: 90 Blog Entries: 1 Rep Power: 9 Thanks Bogar, But I am using "Newtonian" Fluid and "Standard k-epsilon model". but In that case i dont know how to calculate epsilon as i dont know the char. length for the domain. I was trying to calculate it using formula specified at following link and give it as inlet condition.(List of values). There they use "mut" for calculation of "epsilon". [IMG]file:///tmp/moz-screenshot-2.png[/IMG]http://www.cfd-online.com/Wiki/Standard_k-epsilon_model However I am trying to validate results for Fluent against OF, where we have following condotions for fluent. Terbulent Intensity : 3 Standard k-epsilon model. Note : I am giving list of values for U, K and Epsilon at inlet I guess terbulent intensity here belongs to the ration of mut/mu. which should be same for incompressible fluid. How can I represtne the same inlet conditions in OF...??? following is the error i am getting : [1] #0 Foam::error:rintStack(Foam::Ostream&) in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #1 Foam::sigFpe::sigFpeHandler(int) in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #2 __restore_rt at sigaction.c:0 [1] #3 Foam::PBiCG::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #4 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field&, Foam::Field const&) const in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #5 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList >&, Foam::PtrList >&, unsigned char) const in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #6 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libOpenFOAM.so" [1] #7 Foam::fvMatrix::solve(Foam::dictionary const&) in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libfiniteVolume.so" [1] #8 Foam::lduMatrix::solverPerformance Foam::solve(Foam::tmp > const&) in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libincompressibleRASModels.so" [1] #9 Foam::incompressible::RASModels::kEpsilon::correct () in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/lib/linux64GccDPOpt/libincompressibleRASModels.so" [1] #10 main in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/simpleFoam" [1] #11 __libc_start_main in "/lib64/libc.so.6" [1] #12 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/vj/OpenFOAM/OpenFOAM-1.7.1/applications/bin/linux64GccDPOpt/simpleFoam" __________________ Cheers, Vishal Jambhekar... "Simulate the way ahead......!!!"

 December 23, 2010, 11:11 #4 Senior Member   David Boger Join Date: Mar 2009 Location: Penn State Applied Research Laboratory Posts: 146 Rep Power: 9 There are other threads on the forum that already discuss aspects of your question. See this one for example, which also discusses a way in OpenFOAM to specify the turbulence intensity as a boundary condition. I don't think there's any way to get around the fact that you'll need to "know" or make some assumptions about the turbulence length scale in your domain. Presumably, Fluent is making some assumption about it if you're not specifying it some way; you'll have to check the Fluent documentation to find out. As far as the definition of turbulence intensity, it is not mut/mu but rather can be related to something proportional to the square root of k. You should be able to easily find the definition on-line. __________________ David A. Boger

 March 10, 2011, 05:15 #5 Member   Vishal Jambhekar Join Date: Mar 2009 Location: University Stuttgart, Stuttgart Germany Posts: 90 Blog Entries: 1 Rep Power: 9 Hi Boger, Thanks a lot, can you please look into this. I have complete details for my case here Help with k epsilon values of turbulence __________________ Cheers, Vishal Jambhekar... "Simulate the way ahead......!!!"

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post brooksmoses OpenFOAM Running, Solving & CFD 9 March 16, 2014 16:19 herenger OpenFOAM Running, Solving & CFD 7 July 11, 2013 06:27 CFDUSER CFX 2 December 9, 2006 07:31 cfduser CFX 0 April 29, 2006 10:58 confused FLUENT 2 July 29, 2001 21:58

All times are GMT -4. The time now is 18:13.

 Contact Us - CFD Online - Top