# 2nd order TVD scheme

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

 March 13, 2013, 09:03 2nd order TVD scheme #1 New Member   Join Date: Mar 2013 Posts: 14 Rep Power: 6 Hi, I am learning the second order TVD scheme which is listed in the attachment, and I want to verify the order of the scheme by the linear advection equation u_t+u_x=0, periodic BC, initial condition u0=sin(pi*x), x in [0, 2]. But the order I got is about 1.3, and the wave shape becomes flat at valley and peak. Why is the order less than 2 ? My code is in the attachment too. Thank you very much for your kind reply.TVD scheme.png order.doc solution.pdf

 March 13, 2013, 17:01 #2 Senior Member   cfdnewbie Join Date: Mar 2010 Posts: 557 Rep Power: 13 looks like your TVD limiter works as expected and cuts off your peaks when a grid point coincides with the maxima/minina. Your initial projection looks ok, you should approach second order when refining the mesh. Order of convergence is a concept of h to zero, from your plot it looks like h is 2/6.... Wenjun likes this.

 March 13, 2013, 21:45 #3 New Member   Join Date: Mar 2013 Posts: 14 Rep Power: 6 Thank you for your consideration. But even if I choose the step h=2/6,2/12..., the order is almost the same as before. I don't know why. I have already checked my code for several times.

 March 13, 2013, 21:55 #4 New Member   Join Date: Mar 2013 Posts: 14 Rep Power: 6 I use the L1 error for order calculation, it confirms the second order accuracy. However, L_\infty does not work. Normally the two different measure both can confirm the order. Do you have some idea about this problem?

 March 14, 2013, 03:06 #5 Senior Member   cfdnewbie Join Date: Mar 2010 Posts: 557 Rep Power: 13 The L infty will always be dominated by the limiter at the maxima, so since you are doing a first order approximation at the maxima (see your plot, the maxima are approximated by constants, not linears), the L infty will not give you second order. The L1 and L2 will converge to second order, since the influence of the "bad region" will decrease with increasing h!

 March 14, 2013, 03:24 #6 New Member   kian Join Date: Feb 2013 Posts: 14 Rep Power: 6 Hi, I've got a fortran code for solving flow over an airfoil with the method of first order van leer flux vector splitting . I have to change it to second order . Do you know how I can do that ? can you help me with it ?

 March 14, 2013, 03:25 #7 New Member   Join Date: Mar 2013 Posts: 14 Rep Power: 6 Yes, I understand now. I tried L1, it works. Thank you very much. It is my pleasure to meet you in the forum.

 March 14, 2013, 03:51 #8 New Member   Join Date: Mar 2013 Posts: 14 Rep Power: 6 Sorry, I only code by Matlab. I am a fresh guy in CFD and my problem is rather simple.

 March 14, 2013, 13:31 #9 New Member   kian Join Date: Feb 2013 Posts: 14 Rep Power: 6 Thanks anyway

 Tags tvd; second order; shape

 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 nottern Main CFD Forum 1 May 11, 2012 00:30 quarkz Main CFD Forum 30 December 26, 2011 08:12 makaveli_lcf OpenFOAM Running, Solving & CFD 0 November 12, 2009 08:15 Pavel O FLUENT 0 September 28, 2009 09:15 Valle OpenFOAM 0 August 26, 2009 08:12

All times are GMT -4. The time now is 06:21.

 Contact Us - CFD Online - Privacy Statement - Top