CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Elliptic Grid Generation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 22, 2009, 16:59
Default Elliptic Grid Generation
  #1
New Member
 
Join Date: Apr 2009
Posts: 13
Rep Power: 8
student_58 is on a distinguished road
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......
student_58 is offline   Reply With Quote

Old   May 23, 2009, 11:12
Default
  #2
Senior Member
 
John Chawner
Join Date: Mar 2009
Location: Fort Worth, Texas, USA
Posts: 221
Rep Power: 9
jchawner is on a distinguished road
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
File Type: png student_58.png (16.3 KB, 65 views)
__________________
John Chawner / jrc@pointwise.com / www.pointwise.com
Blog: http://blog.pointwise.com/
on Twitter: @jchawner
jchawner is offline   Reply With Quote

Old   May 23, 2009, 12:36
Default
  #3
New Member
 
Join Date: Apr 2009
Posts: 13
Rep Power: 8
student_58 is on a distinguished road
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
student_58 is offline   Reply With Quote

Old   May 23, 2009, 13:14
Default
  #4
Senior Member
 
John Chawner
Join Date: Mar 2009
Location: Fort Worth, Texas, USA
Posts: 221
Rep Power: 9
jchawner is on a distinguished road
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
jchawner is offline   Reply With Quote

Old   May 24, 2009, 03:17
Default
  #5
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 188
Rep Power: 8
harishg is on a distinguished road
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.
harishg is offline   Reply With Quote

Old   May 24, 2009, 06:53
Default
  #6
New Member
 
Join Date: Apr 2009
Posts: 13
Rep Power: 8
student_58 is on a distinguished road
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.
student_58 is offline   Reply With Quote

Old   May 24, 2009, 12:25
Default
  #7
Senior Member
 
N/A
Join Date: Mar 2009
Posts: 188
Rep Power: 8
harishg is on a distinguished road
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.
harishg is offline   Reply With Quote

Old   May 24, 2009, 14:57
Default
  #8
Senior Member
 
John Chawner
Join Date: Mar 2009
Location: Fort Worth, Texas, USA
Posts: 221
Rep Power: 9
jchawner is on a distinguished road
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
jchawner is offline   Reply With Quote

Old   June 9, 2010, 02:31
Default
  #9
New Member
 
Uygur
Join Date: Jun 2010
Posts: 1
Rep Power: 0
makinay is on a distinguished road
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
makinay is offline   Reply With Quote

Reply

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
elliptic grid generation (orthogonal) vasu Main CFD Forum 6 March 24, 2011 07:22
elliptic grid generation prapanj Main CFD Forum 1 November 19, 2007 17:53
Combustion Convergence problems Art Stretton Phoenics 5 April 2, 2002 05:59
Elliptic Grid Generation R.V.Ramkumar Main CFD Forum 2 July 5, 1999 11:06
A quetion about elliptic grid generation A.Rajani Kumar Main CFD Forum 1 March 15, 1999 13:16


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