# r-refinement

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 October 1, 2007, 07:22 r-refinement #1 ganesh Guest   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, Ganesh

 October 1, 2007, 10:12 Re: r-refinement #2 Praveen. C Guest   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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post colinB OpenFOAM Native Meshers: snappyHexMesh and Others 13 July 4, 2017 16:27 Clementhuon OpenFOAM Native Meshers: snappyHexMesh and Others 6 March 12, 2012 13:41 sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11 Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 0 November 10, 2010 04:23 chelvistero OpenFOAM 11 January 15, 2010 20:43

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

 Contact Us - CFD Online - Privacy Statement - Top