# convergence analysis for viscous problems

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

 September 22, 2009, 04:19 convergence analysis for viscous problems #1 New Member   Shyam Join Date: Apr 2009 Posts: 29 Rep Power: 10 Hi, I am trying to benchmark my higher order DG code for viscous problems. I devised a new method of representing the viscous flux at cell interface and am required to establish the order of accuracy of the proposed scheme. I took the standard test case of isentropic vortex evolution problem as the initial setup and am trying to estimate the error as the difference between the solutions (at t=5) of the current grid and the refined grid. For example, if I have grids G1, G2, G3, the error for G1 would be |Q1-Q2| where Q is the solution on the grid. The error for G2 would be |Q2-Q3| and so on. Is this the right way to estimate the grid convergence if we do not know the analytical solution? When I do the above test, I get acceptable grid convergence results. I am just wondering if it is the right way to do... Shyam

September 25, 2009, 13:54
#2
New Member

Josh
Join Date: Sep 2009
Posts: 14
Rep Power: 9
Quote:
 Originally Posted by shyamdsundar For example, if I have grids G1, G2, G3, the error for G1 would be |Q1-Q2| where Q is the solution on the grid. The error for G2 would be |Q2-Q3| and so on. Is this the right way to estimate the grid convergence if we do not know the analytical solution? When I do the above test, I get acceptable grid convergence results. I am just wondering if it is the right way to do...
That approach should give you the right order of convergence estimate (slope of log(err) v log(dx)), but won't give you an estimate of the constant multiplier out front (intercept of log(err) v log(dx)). To get that you would want to do something like: |Q1-Q4|, |Q2-Q4|, |Q3-Q4|, where the grid for Q4 is much finer than 1,2 and 3, but 1,2 and 3 are still in the asymptotic range. So you're basically treating Q4 as your analytic solution.

Not many people care about estimating the constant on the truncation error though.

 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 Vidya FLUENT 24 January 1, 2016 14:01 HSeldon FLUENT 1 March 15, 2007 02:04 Carlos FLUENT 0 March 12, 2007 03:44 Simone Siemens 5 June 29, 2005 10:48 Olivier Macchion CFX 4 January 22, 2003 02:57

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