|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Rolando Maier
Join Date: Mar 2009
Posts: 89
Rep Power: 17 ![]() |
Hi,
this is not a bug but a way to overcome some trouble resulting from surfaceInterpolation::deltaCoeffs(). The problem arises, if there are some concave cells in the mesh. (I know, thatīs not fine, but they are outcome of a mesh moving process and do not affect the quality of my calculation.) ![]() As there is: DeltaCoeffs = 1.0/(mag(unitArea & delta) + VSMALL); deltaCoeffs tends to infinity in the above case. This causes the problems in my calculation. To overcome this problem I use a limiter: DeltaCoeffs = 1.0/mag(max(unitArea & delta, 0.087 * mag(delta))); The parameter 0.087 could be set smaller (~0.2 should also do fine) to approximate the original implementation. The parameter 0.087 is the output of some test calculations. Rolando |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 22 ![]() |
Your proposal looks fine and I will consider it further. How did you find that 0.087 is the optimal coefficient and how general do you think it is?
Henry |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Rolando Maier
Join Date: Mar 2009
Posts: 89
Rep Power: 17 ![]() |
I could not say that 0.087 is the optimal value, I think that a somewhat smaller value would be more general. With 0.03 or 0.04 you should be save.
The value 0.087 corresponds to a maximum angle of 85° between delta and unitArea. By choosing this value there was no more influence of the concave cells in my calculation. Rolando |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Billy
Join Date: Mar 2009
Posts: 167
Rep Power: 17 ![]() |
nice picture and explanation.
|
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
SurfaceTransformPoints rotate suggestion for improvement | markc | OpenFOAM Pre-Processing | 1 | August 13, 2011 08:47 |
SurfaceTransformPoints rotate suggestion for improvement | markc | OpenFOAM Pre-Processing | 0 | January 9, 2009 02:28 |
Improvement to the update of T_ et al in hMixtureThermoC and hhuMixtureThermoC | richpaj | OpenFOAM Bugs | 2 | October 9, 2008 22:14 |
Jet Impingement (results improvement) ... | Sri | FLUENT | 6 | July 16, 2007 06:33 |
sponsership-research-innovative improvement-wind | amit | Main CFD Forum | 0 | February 9, 2006 07:51 |