# Elliptic Grid Generation

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

 May 22, 2009, 16:59 Elliptic Grid Generation #1 New Member   Join Date: Apr 2009 Posts: 13 Rep Power: 8 Hi; I need your help to solve my problem. My problem is grid generation.I have a geometry.İt isn`t rectangular. I attach image. Problem.JPG İn the geometry, There is a radius. I have to use elliptical grid generation method to solve problem.Firstly, I convert physical domain to computational domain. How can ı make this? I wait your help......

May 23, 2009, 11:12
#2
Senior Member

John Chawner
Join Date: Mar 2009
Location: Fort Worth, Texas, USA
Posts: 221
Rep Power: 9
Hello student_58.

It's not clear exactly what kind of help you're seeking (but not knowing the question never stopped me from offering an answer).

Since you're interested in elliptic grid generation, I'll assume your grid is structured quads (also known as a mapped mesh). So the first issue is one of topology; how will you map the region in the image to a rectangle.

I propose that the easiest way to do this is to locate the 4 corners of the computational domain as follows.
1. upper left corner
2. upper right corner
3. lower right corner
4. top of the quarter circle on the left
See the topology in the attached picture.

Assuming you can create the geometry of the bounding edges, the next step is to distribute grid points along the edges. In the attached picture the points are simply equally spaced. Of course, you have to use the same number of points on the top and bottom edges and the left and right edges.

Before applying the elliptic PDE technique you have to first initialize the points on the interior of the domain. While I suppose you could just initialize them all to (0,0) a better approach would be to apply a transfinite interpolation (TFI) method. TFI methods are well documented in most every CFD text that discusses grid generation.

With the interior mesh points initialized, you then iteratively apply the elliptic PDE method to move the points to improve smoothness, clustering, orthogonality, or any other geometric property of the grid. I assume that your professor has described how the elliptic PDE method works so I won't go into that here.

Hope this helps.
Attached Images
 student_58.png (16.3 KB, 65 views)
__________________
John Chawner / jrc@pointwise.com / www.pointwise.com
Blog: http://blog.pointwise.com/

 May 23, 2009, 12:36 #3 New Member   Join Date: Apr 2009 Posts: 13 Rep Power: 8 Thanks for your interest,Jchawner. I read your answer and ı see your grid.But it isn`t right answer.Because ı want to convert phsyical domain to computational domain.For this, ı use numerical methods and elliptic grid generator method.You probably make this grid at a packet program(gambit,gridgen or otherwise).But ı have to make a program for this, using MATLAB. My computational domain has to be rectangular.I attach domain image which ı want to make.Problem.JPG

 May 23, 2009, 13:14 #4 Senior Member   John Chawner Join Date: Mar 2009 Location: Fort Worth, Texas, USA Posts: 221 Rep Power: 9 I guess that I really don't understand your question. I don't know anything about MATLAB so I can't help you with that. You write that your computational domain "has to be rectangular". The computational domain for the grid in the picture I attached is rectangular. Generating a computational domain is trivial. You could create the computational domain in Microsoft Word by simply typing in the coordinates. What exactly is the problem? __________________ John Chawner / jrc@pointwise.com / www.pointwise.com Blog: http://blog.pointwise.com/ on Twitter: @jchawner

 May 24, 2009, 03:17 #5 Senior Member   N/A Join Date: Mar 2009 Posts: 188 Rep Power: 8 You are confusing two things. The way you do it is generate the grid using elliptic grid generator and then solve the Navier-Stokes equations in curvilinear coordinates. So even though you have a curvy grid, your solution is performed in a rectangular coordinate. Check the grid generation book of Thompson available online to understand how it should be done. When you post on a forum, ensure that you had already searched using google and do not just post a question and expect someone else to solve it fuly for you.

 May 24, 2009, 06:53 #6 New Member   Join Date: Apr 2009 Posts: 13 Rep Power: 8 Sure that, ı search a lot for this,harishg. My question is about radius area.My physical cordinate are x,y.Comptutational cordinate ksi,eta.Radius area prevent me to solve the problem. x=x(ksi,eta) y=y(ksi,eta) I haven`t any transforming equation, like ksi=x or eta=x/y^2. How do ı convert x,y cordinate to ksi and eta? Your grid cell are rectangular.But ı want to move my geometry to a rectangular domain.Radius area will be a rectangular.Then ı will solve the problem with mapped mesh.

 May 24, 2009, 12:25 #7 Senior Member   N/A Join Date: Mar 2009 Posts: 188 Rep Power: 8 You generate an elliptic grid using some elliptic grid generation method. You can assume deta, dpsi = 1 respectively. With the values of x and y known, you calculate dphi/dx = - J dx/dsi ... , you have four such formulae. Check some grid generation or CFD book. Then say your equation is df/dt + df/dx =0 You transform it to df/dt + df/dphi*dphi/dx + df/dsi*dsi/dx =0 , and since u know the grid changes, you calculate df/dphi in the rectangular grid and proceed with your computation. Take a look at any CFD book and you would have the whole procedure explained clearly.

 May 24, 2009, 14:57 #8 Senior Member   John Chawner Join Date: Mar 2009 Location: Fort Worth, Texas, USA Posts: 221 Rep Power: 9 Hello Again student_58. Perhaps the best help I can offer is a link to Thompson's classic book on grid generration. http://www.hpc.msstate.edu/publicati...book/index.php __________________ John Chawner / jrc@pointwise.com / www.pointwise.com Blog: http://blog.pointwise.com/ on Twitter: @jchawner

 June 9, 2010, 02:31 #9 New Member   Uygur Join Date: Jun 2010 Posts: 1 Rep Power: 0 hi student_58 my problem is same, http://www2.itu.edu.tr/~guneshasa/cfd/homework/hw5.pdf if you solve your problem could you help me

 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 vasu Main CFD Forum 6 March 24, 2011 07:22 prapanj Main CFD Forum 1 November 19, 2007 17:53 Art Stretton Phoenics 5 April 2, 2002 05:59 R.V.Ramkumar Main CFD Forum 2 July 5, 1999 11:06 A.Rajani Kumar Main CFD Forum 1 March 15, 1999 13:16

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