# Bugs in LES filter codes

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

 February 17, 2008, 04:12 Bug in LES Filter Codes: Desc #1 New Member   YekaniYasari Join Date: Mar 2009 Location: Göteborg Posts: 6 Rep Power: 8 Bug in LES Filter Codes: Description: 1-The problem is about the laplace filter and anisotropic filter: with running a case in LES solvers (oodles, channelOodles …) and using "laplace" filter or "anisotropic" filter with "dynamic or mix models" an error about laplace or anisotropic filtercoeffs occurred .for solving this problem we should add the following line in "turbulenceproperties" file in "constant" directory of relevant case. For example: mixedSmagorinskyCoeffs { ce ce [0 0 0 0 0 0 0] 1.05; ck ck [0 0 0 0 0 0 0] 0.07; filter anisotropic; anisotropicCoeffs { widthCoeff 24; or 6 (for box or Gaussian filter) } } 2.The second problem was seen about definition of "widthcoeff" vector in "anisotropicfilter.C" file that was placed in"/src/lesmodel/lesfilter" . In order to modify this problem, instead of the statement below: " for (direction d=0; d tmpFilteredField ( new volTensorField ( IOobject ( "anisotropicFilteredTensorField", mesh().time().timeName(), mesh(), IOobject::MUST_READ, IOobject::AUTO_WRITE ), mesh(), dimLength ) ); With this form : tmp tmpFilteredField ( new volTensorField ( IOobject ( "anisotropicFilteredTensorField", mesh().time().timeName(), mesh(), IOobject::NO_READ, IOobject::AUTO_WRITE ), mesh(), dimLength*dimLength/(dimTime*dimTime) ) ); the problem is solved. You can download the correct file from here. anisotropicFilter.C Testcase: All case for example: channel395 Version: OpenFOAM version 1.3

 February 17, 2008, 04:41 in number 2 change " "widthcoe #2 New Member   YekaniYasari Join Date: Mar 2009 Location: Göteborg Posts: 6 Rep Power: 8 in number 2 change " "widthcoeff" vector " to " "coeff_" vector "

 February 17, 2008, 07:10 I think it would be better to #3 Senior Member   Join Date: Mar 2009 Posts: 854 Rep Power: 13 I think it would be better to use unFilteredField().dimensions() rather than hard-coding dimLength*dimLength/(dimTime*dimTime) as the dimensions of the tensor field resulting from the filtering of given field. Henry

 February 18, 2008, 02:59 Dr Henry Weller, Thanks for #4 New Member   YekaniYasari Join Date: Mar 2009 Location: Göteborg Posts: 6 Rep Power: 8 Dr Henry Weller, Thanks for your attention. Could you kindly comment about a problem in : http://www.cfd-online.com/OpenFOAM_D...ages/1/40.html By ehsan yasari on Tuesday, January 01, 2008 - 12:59 am Best Regards

 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 su_junwei OpenFOAM Bugs 4 July 7, 2013 10:29 luca OpenFOAM Bugs 2 January 27, 2009 16:40 Andy R Main CFD Forum 0 July 25, 2008 13:15 whitemelon CD-adapco 0 July 11, 2008 06:26 whitemelon CD-adapco 6 May 16, 2008 10:44

All times are GMT -4. The time now is 05:34.