CFD Online Logo CFD Online URL
Home > Forums > General Forums > Main CFD Forum


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

LinkBack Thread Tools Display Modes
Old   October 1, 2007, 07:22
Default r-refinement
Posts: n/a
Dear Friends,

I am currently looking at r-refinement, wherein the nodes are moved according to the solution strengths. Currently, I am looking at a modified laplacian strategy employed by Marchant and Weatherill. In essence they make use of the following expression.

r_o^(n+1) = r_o^(n) + delta

delta = sum(w_i*(r_i - r_o))/sum(w_i)

where w_i is the weights computed as w_i = k1 + k2 * dabs((phi_i - phi_o)/(phi_i + phi_o)).

Here, phi stands for the solution variable, r the position vector and the summation is over all neighbours i of the point o. k1 and k2 are constants for noise damping and gradient amplification. Also, a slight modification of the aforementioned strategy is employed for boundary nodes. I have two questions.

1. What are the commonly used values of k1 and k2, and what is the basic idea behind selecting those ?

2. On the boundaries, how is the procedure modified, so as to also take care of the fact that it always lies on the boundary itself.

I will be grateful if anyone with experience on these issues is able to help me.

Thanks and Regards,


  Reply With Quote

Old   October 1, 2007, 10:12
Default Re: r-refinement
Praveen. C
Posts: n/a
In constant or nearly-constant regions, you would expect the grid to be isotropic. The node movement should reduce to an averaging formula, all weights should be equal. But there will be usually noise in the solution, may be due to limiters. Adding a k1 term which is larger than the magnitude of the second term due to noise would give nearly constant weights.

k2 seems to be just a scaling term for the gradients. I guess one needs to do some trial and error to determine good values for these.
  Reply With Quote


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
Number of cells in mesh don't match with size of cellLevel colinB OpenFOAM Native Meshers: snappyHexMesh and Others 13 July 4, 2017 16:27
Boundary layer in a pipe Clementhuon OpenFOAM Native Meshers: snappyHexMesh and Others 6 March 12, 2012 13:41
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
snappyHexMesh aborting Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 0 November 10, 2010 04:23
external flow with snappyHexMesh chelvistero OpenFOAM 11 January 15, 2010 20:43

All times are GMT -4. The time now is 02:27.