CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

2D linear advection equation

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

LinkBack Thread Tools Display Modes
Old   February 11, 2011, 11:00
Default 2D linear advection equation
New Member
Miguel Caro
Join Date: Apr 2010
Posts: 26
Rep Power: 9
mcaro is on a distinguished road
Well, trying to solve a 2D linear advection equation
u_t + au_x + bu_y = 0;
u_0(x,y,0) = sin( 2pi* x ) sin( 2pi y), (x,y) 0,1) x (0x1) , periodic boundary conditions

with exact solutions u(x,y,t) = sin (2pi (x-t) ) sin (2pi (y-t) )

i implemented this discretization :

u_i,j^{n+1} = u_i,j^n - dt/dx(Fi+1/2 - Fi-1/2) - dt/dy(Gi+1/2 - Gi-1/2);

CFL = max(dt/dx;dt/dy);

But the solution did not agree with the exact solution at a time t=2.0

Does anybody has a suggestion of what kind of algorithm to implement?

I would like to imlement a new-high resoluton scheme and i woud like to evaluate the convergence order. For this i would lik e to implement a technique including a flux limiter.

I realyy appreciate any help

mcaro is offline   Reply With Quote

Old   February 11, 2011, 12:31
Senior Member
Join Date: Nov 2009
Posts: 411
Rep Power: 12
DoHander is on a distinguished road
Try to increase the order of your temporal discretization by using a Runge-Kutta method (order 4 should do). This will allow you to use a reasonable time step and to obtain a more precise solution. Or, you may consider using an implicit temporal discretization.

DoHander is offline   Reply With Quote

Old   February 11, 2011, 13:30
New Member
Miguel Caro
Join Date: Apr 2010
Posts: 26
Rep Power: 9
mcaro is on a distinguished road
Thanks Do, well i implemented yet a third order temporal RK, but with a more refined grid the numerical solution has an acceptable approximation to the reference solution. But I think i have to implement a fourth order temporal RK (usually it is called Strong Stability Prerserving Runge-Kutta of fourth order, for example. Gottlieb, Shu and others) to obtain better approximation.

mcaro is offline   Reply With Quote


2d linear advection

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
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
Singular Linear System - Pressure Poisson Equation Allan Main CFD Forum 11 March 17, 2011 10:39
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56
solution diverges when linear upwind interpolation scheme is used subash OpenFOAM 0 May 29, 2010 01:23
non linear continuity equation S.Venkat Main CFD Forum 0 May 30, 2000 07:43

All times are GMT -4. The time now is 17:03.