Wall function help!

 Register Blogs Members List Search Today's Posts Mark Forums Read

February 21, 2011, 15:52
Wall function help!
#1
Member

Ammar Tareen
Join Date: Jan 2011
Location: Boston University
Posts: 61
Rep Power: 8
Dear Foamers,

Basic setup: I am passing fluid through a laval nozzle and I expect the velocity to exceeded Mach 1, so I decided to use the K-epsilon model for my simulation; my solver is sonicFoam. The thermophysical properties for my fluid are as follows:

Quote:
 thermoType ePsiThermo>>>>; mixture KrHe // 20%helium-80%krypton mixture. 1 //nMoles (??) 69.04 //molecular weight [kg/kmol] 1041.8 // heat capacity [K/kgK] 0.0 //standard formation enthalpy [J/kg] muPolynomial (1.27906e-005 4.61187e-008) //dynamic viscosity [kg/ms] kappaPolynomial (0.0406885 0.000146747); //heat conductivity [W/mK]
(here the viscosity and heat conductivity are linear temperature-dependant polynomials)
Confusion: As you can see from the mupolynomial, I want the magnitude of the viscosity to be somewhere around the 10e-5 range. However, because I am using the k-e model there is a mut file in the 0 directory andI am afraid that this mut file is causing the viscosity to be higher than I actually want in my simulation. The value of mut is some around 10e-2. The contents of 0/mut are as follows:

Quote:
 FoamFile { version 2.0; format ascii; class volScalarField; location "0"; object mut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [1 -1 -1 0 0 0 0]; internalField uniform 0; boundaryField { inle2 { type calculated; value uniform 0; } outl3 { type calculated; value uniform 0; } wall1 { type mutWallFunction; Cmu 0.09; kappa 0.41; E 9.8; value uniform 0; } frontAndBack { type empty; } }

I don't quite understand what the mut file is for? Shouldn't muPolynomial in my thermophysical properties be sufficient for specifying the viscosity?

Ideally I want the viscosity to be specified by my thermoPhysical properties dictionary only and I don't want a mut file, in fact I don't even understand what it is for or what a wall function means. Can some body please please please help me or give me some advice?

Best regards,
Ammar.

 February 22, 2011, 13:46 #2 Senior Member     Marco A. Turcios Join Date: Mar 2009 Location: Vancouver, BC, Canada Posts: 735 Rep Power: 21 Read up on law of the wall and wall models/functions on the CFD Wiki: http://www.cfd-online.com/Wiki/Law_of_the_wall http://www.cfd-online.com/Wiki/Wall_functions If you have a wall in your domain and you are running turbulent simulation, you need to specify wall functions or else your simulation isn't going to make sense. Your simulation should be using whatever is in your thermoPhysical properties file. Calculated boundary conditions mean nothing; they simply show what the calculated value for a variable is for the particular timestep. Only boundary conditions like fixedValue, zeroGradient (basically anthing in finiteVolume/fields/fvPatchFields and similar) will affect your simulation. Hope this helps.

 February 22, 2011, 13:48 #3 Member   Ammar Tareen Join Date: Jan 2011 Location: Boston University Posts: 61 Rep Power: 8 Thanks soo much!

 Tags k-e model, viscosity, wallfunction

 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 AlmostSurelyRob OpenFOAM 3 June 24, 2011 13:06 Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42 yka8150 Main CFD Forum 0 September 21, 2009 23:08 Christian Main CFD Forum 2 February 27, 2007 07:27 tod Phoenics 1 May 19, 2003 05:05

All times are GMT -4. The time now is 06:40.