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

errors in curvilinear grid Jacobian matrix

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 11, 2007, 12:08
Default errors in curvilinear grid Jacobian matrix
  #1
DSS
Guest
 
Posts: n/a
Hi,

I am working on higher order methods to solve euler equations. I have recently extended the code to curvilinear coordinates. The code works fine. But, I noticed that, the calculation of jacobian transformation matrix is not very accurate and this leads to spurious oscillations in the solution.

As of now, I dont know how to calculate the Jacobians accurately so that the error is reduced during transformations...

When the grid lines are straight, the error is quite negligible. But, when I have a grid like x=u, y=v+sin(u) (u,v are from 0 to 1), I have some error terms. I cant reduce the error even if I increase the order of accuracy to 8. If anyone is aware of such problems, please guide me.

Regards, DSS
  Reply With Quote

Old   February 11, 2007, 22:46
Default Re: errors in curvilinear grid Jacobian matrix
  #2
VSZ
Guest
 
Posts: n/a
If you don't mind, can I see yours in detail? I'd like to take a look at it and let's see how to help you.
  Reply With Quote

Old   February 11, 2007, 22:59
Default Re: errors in curvilinear grid Jacobian matrix
  #3
DSS
Guest
 
Posts: n/a
I will try to explain what I did to quantify the errors.

I generated the grid points according to the equations, or by solving elliptic equations.

Now, I calculate the Jacobian matrix ( d(xi)/d(x) ) where xi is the computational coordinate and x is the physical coordinate.

With the available transformations, I solved the standard equation

d(Q)/d(t) + d(F)/d(x) = 0

in the transformed sense. I used the conservative form so that the fluxes become

Ft(xi) = 1/J ( F_x d(xi)/d(x) + F_y d(xi)/d(y) + F_z d(xi)/d(z) )

Now, ideally, when I solve this with F=1 for all values, I should get d(Q)/d(t) = 0. This, I am able to get to the value of 1e-13 if the grid lines are linear. However, when the grid is skewed, this error increases to the order of 1e-3.

I calculate the derivatives using a 5 point FD stencil, which transforms to backward and forward differencing near the boundaries. I witnessed the same with 3 point and 9 point stencils.

As mentioned before, I used the grid of the sort

x = u y = v + sin(u)

where (u,v) range from 0 to 1.

Hope this gives enough information.
  Reply With Quote

Old   February 12, 2007, 01:37
Default Re: errors in curvilinear grid Jacobian matrix
  #4
Harish
Guest
 
Posts: n/a
What is the grid spacing you use for the computation.
  Reply With Quote

Old   February 12, 2007, 02:27
Default Re: errors in curvilinear grid Jacobian matrix
  #5
DSS
Guest
 
Posts: n/a
Used 51 points for each direction for the equation

x = u y = v + sin(u)

u,v from 0 to 1.

The error is more as you move along x axis.
  Reply With Quote

Old   February 12, 2007, 02:38
Default Re: errors in curvilinear grid Jacobian matrix
  #6
Harish
Guest
 
Posts: n/a
Did you check the value of jacobian at each point.does any point have jacobians close to 0.Also how does the derivatives (dxi/dx) vary.Do they have some sort of steep rise in some regions.Abrupt variations in metrics can produce parasite waves.
  Reply With Quote

Old   February 12, 2007, 02:47
Default Re: errors in curvilinear grid Jacobian matrix
  #7
DSS
Guest
 
Posts: n/a
actually, the problem is not abrupt variation of the metric values. But, the inability to calculate the derivatives d(x)/d(xi) accurately, leading to error in calculation of d(xi)/d(x) and the Jacobian. This results in loss of conservation in the transformed equations.

So, when you start the iteration, there are some parasitic waves which are visible if you do a low mach no. case.

  Reply With Quote

Old   February 12, 2007, 02:50
Default Re: errors in curvilinear grid Jacobian matrix
  #8
DSS
Guest
 
Posts: n/a
actually, the problem is not abrupt variation of the metric values. But, the inability to calculate the derivatives d(x)/d(xi) accurately, leading to error in calculation of d(xi)/d(x) and the Jacobian. This results in loss of conservation in the transformed equations. So, when you start the iteration, there are some parasitic waves which are visible if you do a low mach no. case.
  Reply With Quote

Old   February 12, 2007, 07:44
Default Re: errors in curvilinear grid Jacobian matrix
  #9
andy
Guest
 
Posts: n/a
Generally you do not want to calculate the metrics accurately but to arrange the differencing to numerically conserve what is important in your simulation. You can test this by using analytical metrics and comparing with your optimum arrangement of differencing terms.

There were quite a few papers on this subject in the 1970s and early 1980s when the use of curvilinear coordinates started to be heavily used for CFD but I am afraid I have no references to hand.

  Reply With Quote

Old   February 12, 2007, 08:04
Default Re: errors in curvilinear grid Jacobian matrix *NM*
  #10
Jose
Guest
 
Posts: n/a
  Reply With Quote

Old   February 12, 2007, 08:10
Default Re: errors in curvilinear grid Jacobian matrix
  #11
Melky
Guest
 
Posts: n/a
I think Andy is correct. I don't have any references that deal specifically with this issue but I know that you can find them in Yusuf Ozyoruks Ph.D. thesis (The Pennsylvania State University). He implemented a higher order method to solve the Euler/NS equations for the purpose of predicting fan noise. You might be able to download it from Lyle Longs Website. (It appears to be down right now though).

  Reply With Quote

Old   February 12, 2007, 09:44
Default Re: errors in curvilinear grid Jacobian matrix
  #12
ag
Guest
 
Posts: n/a
The classic reference is an AIAA paper by Thomas and Lombard on geometric conservation.
  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
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
OpenFOAM version 1.6 details lakeat OpenFOAM Running, Solving & CFD 42 August 26, 2009 21:47
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07
Non-uniform grid calculation Aspens Main CFD Forum 1 February 23, 2000 15:15
Unstructured Grid in Curvilinear Sakti Main CFD Forum 0 November 23, 1999 06:47


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