CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Bugs

Bugs in LES filter codes

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 17, 2008, 03:12
Default Bug in LES Filter Codes: Desc
  #1
New Member
 
YekaniYasari
Join Date: Mar 2009
Location: Göteborg
Posts: 6
Rep Power: 7
yekaniyasari is on a distinguished road
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<vector::nComponents; d++)
{
coeff_.internalField().replace
(
d,
(2.0/widthCoeff_)*mesh.V()
/fvc::surfaceSum(mag(mesh.Sf().component(d)))().int ernalField()
);
}
"
We should use the following code:
"
for (direction d=0; d<vector::nComponents; d++)
{
coeff_.internalField().replace
(
d,
(4.0/widthCoeff_)*(mesh.V()/fvc::surfaceSum(mag(mesh.Sf().component(d)))().int ernalField())*(mesh.V()/fvc::surfaceSum(mag(mesh.Sf().component(d)))().int ernalField())
);
}
"
3. The third one , was about dimension problem : in "anisotropicfilter.C" file , in definition of filtration for tensor field , with substituting this part of code :
tmp<voltensorfield> tmpFilteredField
(
new volTensorField
(
IOobject
(
"anisotropicFilteredTensorField",
mesh().time().timeName(),
mesh(),
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh(),
dimLength
)
);

With this form :
tmp<voltensorfield> 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
yekaniyasari is offline   Reply With Quote

Old   February 17, 2008, 03:41
Default in number 2 change " "widthcoe
  #2
New Member
 
YekaniYasari
Join Date: Mar 2009
Location: Göteborg
Posts: 6
Rep Power: 7
yekaniyasari is on a distinguished road
in number 2 change " "widthcoeff" vector " to " "coeff_" vector "
yekaniyasari is offline   Reply With Quote

Old   February 17, 2008, 06:10
Default I think it would be better to
  #3
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 12
henry is on a distinguished road
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
henry is offline   Reply With Quote

Old   February 18, 2008, 01:59
Default Dr Henry Weller, Thanks for
  #4
New Member
 
YekaniYasari
Join Date: Mar 2009
Location: Göteborg
Posts: 6
Rep Power: 7
yekaniyasari is on a distinguished road
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
yekaniyasari is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Bugs in fvcsurfaceIntegrate su_junwei OpenFOAM Bugs 4 July 7, 2013 10:29
Bugs in FFT luca OpenFOAM Bugs 2 January 27, 2009 15:40
Wiki - Bugs Andy R Main CFD Forum 0 July 25, 2008 13:15
bugs in starcd 4.06 whitemelon CD-adapco 0 July 11, 2008 06:26
starcd 3.24 has bugs? whitemelon CD-adapco 6 May 16, 2008 10:44


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