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

Improvement in makeDeltaCoeffs

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 15, 2007, 12:43
Default Hi, this is not a bug but a w
  #1
Member
 
Rolando Maier
Join Date: Mar 2009
Posts: 89
Rep Power: 17
rolando is on a distinguished road
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
rolando is offline   Reply With Quote

Old   May 15, 2007, 12:55
Default Your proposal looks fine and I
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
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
henry is offline   Reply With Quote

Old   May 15, 2007, 13:07
Default I could not say that 0.087 is
  #3
Member
 
Rolando Maier
Join Date: Mar 2009
Posts: 89
Rep Power: 17
rolando is on a distinguished road
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
rolando is offline   Reply With Quote

Old   May 15, 2007, 13:49
Default nice picture and explanation.
  #4
Senior Member
 
Billy
Join Date: Mar 2009
Posts: 167
Rep Power: 17
billy is on a distinguished road
nice picture and explanation.
billy is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
SurfaceTransformPoints rotate suggestion for improvement markc OpenFOAM Pre-Processing 1 August 13, 2011 09:47
SurfaceTransformPoints rotate suggestion for improvement markc OpenFOAM Pre-Processing 0 January 9, 2009 03:28
Improvement to the update of T_ et al in hMixtureThermoC and hhuMixtureThermoC richpaj OpenFOAM Bugs 2 October 9, 2008 23:14
Jet Impingement (results improvement) ... Sri FLUENT 6 July 16, 2007 07:33
sponsership-research-innovative improvement-wind amit Main CFD Forum 0 February 9, 2006 08:51


All times are GMT -4. The time now is 09:04.