|
[Sponsors] |
June 1, 2010, 05:01 |
Anyone using IDDES
|
#1 |
Member
Markus Weinmann
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77
Rep Power: 17 |
Hi all,
I was wondering whether anyone is using IDDES using the IDDESdelta in LESProperties. I have applied IDDES to two different cases and the results were not exactly what I would expect them to be (compared to the case in the original paper). So I started digging into the code in order to check the implementation. In line 77 in src/turbulenceModels/.../SpalartAllmarasIDDES/IDDESdelta/IDDESDelta.C, I found something which does not make sense in my opinion. The variable deltamax= 2*deltamaxTemp is currently defined as a scalar which is obtained as the maximum cell dimension that appears in the whole domain. To my understanding deltamax should be a FIELD variable (deltamax=hmax = max(dx,dy,dz)) and should be put inside of the loop over all the cells in the domain. What do you guys think? Markus |
|
June 1, 2010, 13:16 |
... yes, we are ... and massivley so ... and ... I am shocked...
|
#2 |
Member
Ulf Bunge
Join Date: Mar 2009
Location: Wolfsburg, Germany
Posts: 34
Rep Power: 17 |
Hi Markus,
thanks for the hint!!! deltamax is the maximum grid cell edge length in the orignal formulation for delta for IDDES and it is used in Code:
{ delta_.internalField() = deltaCoeff_ *min ( max(max(cw_*wallDist(mesh()).y(), cw_*deltamax), hwn), deltamax ); } What do others think!? However, do not forget, that the hybrid numerics which is an important ingredient for a DES is not implemented. This might also explain that you do not achieve the results you want to achieve!? We will present this additional implementation with some results in Gothenburg, however, the results have to be checked carefully, now, and a possible influence of this mistake will have to be evaluated. Best regards, Ulf. Last edited by djbungee; June 1, 2010 at 15:13. Reason: forgot remark about hybrid scheme, corrected a typo |
|
June 1, 2010, 13:38 |
Delta in IDDES
|
#3 |
New Member
Charles Mockett
Join Date: Dec 2009
Location: Berlin, Germany
Posts: 9
Rep Power: 16 |
Dear Markus and Ulf,
I can confirm that the h_max in the IDDES Delta formula should be a local measure of the maximum cell length, not the global maximum. I also agree that this won't be a problem for decaying, isotropic turbulence on a uniform mesh. It also won't be a problem when running a channel flow on the usual kind of structured grid. For anything real however, it will usually be a big problem. Thanks very much for bringing this to our attention! It would be very good to plot the value of Delta in the field, then you can be sure what is happening and also sure that your remedy is effective. Out of interest, what cases have you been running to test the IDDES in OF? I'd be very interested in seeing some results if you wish to share them. Best regards, Charlie. |
|
June 1, 2010, 15:49 |
code adjustement available
|
#4 |
Member
Ulf Bunge
Join Date: Mar 2009
Location: Wolfsburg, Germany
Posts: 34
Rep Power: 17 |
Hi Markus,
do you already have an adjustment of the code for IDDESDelta that you can provide to this forum? My last post of code adjustment for solving unsteady flows with arbitrarily large time steps is some weeks ago - so it might be my turn again . However, at the moment other topics are pressing much more, so it might take some days. Nonetheless, I would like to update our running computations as soon as possible.... Best regards, Ulf. |
|
June 2, 2010, 05:12 |
|
#5 |
Member
Markus Weinmann
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77
Rep Power: 17 |
Hi
I am glad to hear that at least some people agree with me. I have submitted a bug report about this issue a long time ago but no one was interrested in fixing that bug at that time. Ulf: You have mentioned that you have also observed some unexpected behviour of IDDES. I would be very interested to hear more details about that. Charles: I was using IDDES for a channel flow and a tandem cylinder configuration. Best regards, Markus |
|
June 9, 2010, 04:25 |
Our observations
|
#6 |
Member
Ulf Bunge
Join Date: Mar 2009
Location: Wolfsburg, Germany
Posts: 34
Rep Power: 17 |
Hi Markus,
what we computed so far is, decay of isotropic turbulence, inviscid isotropic turbulence, and a NACA-profile wake flow for high angle of attack. Based on several comparisons that Charlie has made he found the behavior of our IDDES too diffusive for the wake flow. This, however, was only one of the observations, so we have to perform more computations, especially for a very long time to have statistically significant data. The dissipative behavior can be explained by the wrong delta, but this has to be evaluated in detail. Best regards, Ulf. |
|
June 30, 2010, 03:35 |
|
#7 |
Member
Fabian Peng Karrholm
Join Date: Mar 2009
Posts: 61
Rep Power: 17 |
Here is a small modification in IDDESDelta.C that changes deltamax into a field instead of a scalar. I would be happy to hear from anyone who tries it.
|
|
September 20, 2010, 08:08 |
anything else wrong in IDDES implementation?! hmax in SpalartAllmarasIDDES.C
|
#8 |
Member
Ulf Bunge
Join Date: Mar 2009
Location: Wolfsburg, Germany
Posts: 34
Rep Power: 17 |
Dear Fabian, dear all,
there are some more peculiarities in the IDDES implementation regarding the value hmax in the turbulence model itself that we are now checking in detail. In case anyone has already checked or confirmed that and found some/none/more/less mistakes, please let us know. E.g., in SpalartAllmarasIDDES.C we find for alpha: Code:
tmp<volScalarField> SpalartAllmarasIDDES::alpha() const { return max ( 0.25 - y_/dimensionedScalar("hMax", dimLength, max(cmptMax(delta()))), scalar(-5) ); } Best regards, Ulf. |
|
September 22, 2010, 06:16 |
|
#9 |
Member
Markus Weinmann
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77
Rep Power: 17 |
Hi Ulf,
I have to agree with you - I don't think the present alpha() is correct. In fact, I have checked the IDDES implemantation a while ago and I have changed alpha() such that it uses hmax (deltamax) as computed in IDDESDelta.C instead of dimLength, max(cmptMax(delta()))) Markus |
|
September 22, 2010, 08:38 |
|
#10 |
New Member
Christian Graurock
Join Date: Jul 2010
Location: Berlin, Germany
Posts: 6
Rep Power: 16 |
Dear Markus,
I'm also interested in this topic. As mentioned before the h_max used in the alpha-calculation should be a local value as well? Thanks, Christian. |
|
September 22, 2010, 08:50 |
|
#11 |
Member
Markus Weinmann
Join Date: Mar 2009
Location: Stuttgart, Germany
Posts: 77
Rep Power: 17 |
yes that is correct - hmax should be a local value and not a single value for the whole domain. I am not sure whether that as already been corrected in the file IDDESDelata.C for the latest 1.7.x release.
Markus |
|
November 30, 2010, 17:24 |
|
#12 |
New Member
James
Join Date: Jan 2010
Posts: 6
Rep Power: 16 |
Dear All,
I'm very interested to know how you have got on with IDDES in OpenFOAM since your last posts? I'm currently looking at using it for incompressible flow over the Ahmed body, and am wondering if you have come across any other problems with the current implementation? EDIT: I've had a look at alpha() and it appears it hasn't been changed yet in the latest 1.7.x. Markus please could you tell me what changes you made to use deltamax from IDDESDelta? Regards James Last edited by tookey_1989; December 2, 2010 at 14:03. |
|
January 7, 2012, 17:00 |
IDDES bugs in Openfoam 201
|
#13 |
New Member
Join Date: Nov 2011
Posts: 13
Rep Power: 15 |
Dear all,
I'm using IDDES in OpenFOAM 201, but I'm not sure whether the alpha and delta bugs have been removed - I don't have enough experience with C to go by my opinion alone, so I'd be grateful for your advice. Delta: see attached code from OF 201 Alpha: Lines 45-52 of SpalartAllmarasIDDES.C: Code:
tmp<volScalarField> SpalartAllmarasIDDES::alpha() const { return max ( 0.25 - y_/static_cast<const volScalarField&>(hmax_()), scalar(-5) ); } |
|
June 29, 2012, 11:12 |
|
#14 |
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 16 |
Hey,
what about the bugs, are they fixed in OpenFOAM-2.1.1? Thanks Anne |
|
July 3, 2012, 12:15 |
|
#15 | ||
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 16 |
By the way, the IDDESDelta is missing in the OpenFOAM user guide.
Here what is written inside user guide Quote:
Quote:
I set them with the help of literature to { cw 0.15; deltaCoeff 1; } inside LESProperties. Cheers Anne |
|||
December 14, 2015, 14:27 |
RANS/LES Region.
|
#16 |
Senior Member
Hasan K.J.
Join Date: Dec 2011
Location: Bristol, United Kingdom
Posts: 200
Rep Power: 15 |
Dear All,
As it looks like there are lot of DES users here I have decided to ask my questions here, I have good experience with RANS simulations and reasonable experience with LES simulations, it is due to some time constraints have started using DES recently and I have started with IDDES. - I am using DES for flow around an airfoil ! - Can some one tell me what criteria is that the IDDES looks at before switching to LES and RANS, for my case i feel the LES regions are the wrong regions. - If it is mesh density !! what ratio of the density is it ? and how do I calculate and define it !!. does the LES region have to be highly dense or the RANS region has to be highly Dense ? - because as of now the LES appears to be on the coarse regions of the mesh only !! that doesnt make sense LES requires DENSE mesh right ? - I have been playing with this for more than a month but still unsuccessfull results in getting LES and RANS in the regions I particularly want. eg. Image attached, the LES is in the wake of the airfoil, I want LES everywhere except at the airfoil close to the wall regions. Can some one please help me, Thanks for your time and Effort. Edit: If the post here hasnt given enough information, I have also posted a similar question here http://www.cfd-online.com/Forums/ope...tml#post576089 Regards, Hasan K.J
__________________
"Real knowledge is to know the extent of one's ignorance." - Confucius Last edited by Alhasan; December 15, 2015 at 14:40. |
|
December 16, 2015, 10:37 |
|
#17 |
Senior Member
Anne Gerdes
Join Date: Aug 2010
Location: Hamburg
Posts: 168
Rep Power: 16 |
Hey Hasan,
the IDDES model is quite complex and it has been a while ago since I dealt with it. In general, RANS should be applied next to the wall where the mesh density is high. LES should be applied in regions far away from the wall where the mesh density is lower. The IDDES model does not only take into account geometrical properties but also flow properties (turbulent eddy viscosity) to compute the switch between RANS and LES. For a more detailed explanation, have a look at the Paper of Shur, Spalart, Strelets and Travin in the Journal of Heat and Fluid Flow from 2008. As far as I know, the IDDES model implemented in OpenFOAM is the same as stated in this paper. Hope this helps. Kind Regards Anne |
|
July 6, 2018, 01:30 |
Modification of hmax
|
#18 | |
Member
王莹
Join Date: May 2017
Posts: 51
Rep Power: 9 |
Hello,
I agree with your findings. I am trying to modify the defination of hmax in OF 5.0 and I have found the related code, which is "maxdeltaxyz.C". However, I really cannot find the location of the defination. In my mind, hmax=max(deltax, deltay, deltaz), so that there should be a max function. can anyone figure out what I have missed? Thanks a lot! Alisa. Quote:
|
||
June 19, 2021, 05:06 |
IDDES tutorials
|
#19 |
Member
Join Date: Apr 2021
Posts: 41
Rep Power: 5 |
Hello,
I can see there is a Spalart-Allmaras IDDES tutorial here : http://="https://www.openfoam.com/do...ted-cube.html" However I cannot see any tutorial regarding the k-omega-SST IDDES. Here is my question: is there anything to add to the following lines in its formulation in the turbulenceProperties file ? Code:
simulationType LES; LES { turbulence on; LESModel kOmegaSSTIDDES; ... ? } |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem running IDDES with pimpleFoam | charlie | OpenFOAM Running, Solving & CFD | 7 | July 11, 2013 12:45 |
IDDES delta | cfdmarkus | OpenFOAM Bugs | 1 | July 17, 2012 05:28 |