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

MUSCL scheme implementation

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

Like Tree2Likes
  • 1 Post By ANPDEstudent
  • 1 Post By praveen

LinkBack Thread Tools Display Modes
Old   June 26, 2010, 10:28
Default MUSCL scheme implementation
New Member
J Kyle Miller
Join Date: Jun 2010
Posts: 2
Rep Power: 0
ANPDEstudent is on a distinguished road
This is a very basic question, regarding MUSCL schemes as described in

When solving a conservation law u_t + d/dx[f(u)] = 0 using a MUSCL scheme,
we arrive at the semi-discrete formula d/dt[u_i] + (1/deltax)*[F(ustar_i+half) - F(ustar_i-half)] = 0.

ustar_i+/- is calculated through the use of a flux limiter phi(r), where r_i = (u_i - u_i-1) / (u_i+1 - u_i). Since r is evaluated at grid point i-1, i, and i+1, this results in the use of five spatial grid points (i-2 through i+2).

My question is the following. Given that this scheme is five points wide, what is typically done for the grid points just inside the left and right boundary? That is, if the grid goes from i=1 to N, then i=2 and i=N-1 require information from points outside the grid. How are these points treated? How is treatment different for Dirichlet vs. Neumann conditions?

Thank you very much for all your help
AshwaniAssam likes this.
ANPDEstudent is offline   Reply With Quote

Old   June 26, 2010, 23:57
Super Moderator
praveen's Avatar
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 277
Blog Entries: 6
Rep Power: 11
praveen is on a distinguished road
You can drop the accuracy on the boundary, use a lower order scheme. Usually, dropping the accuracy by one for boundary points does not affect the overall accuracy. But this is a difficult thing to prove in practice, maybe you can do numerical studies.

You can introduce ghost cells and find the solution in them using boundary conditions.
SangeethCFD likes this.
praveen is offline   Reply With Quote


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
Implementation of QUICK scheme Romuald Skoda Main CFD Forum 11 November 6, 2017 22:20
Trouble when using MUSCL scheme on ONERA M6 wing Aurélien FLUENT 2 January 23, 2013 08:42
extrapolation in MUSCL scheme Chandra Main CFD Forum 6 February 14, 2007 12:21
FEM Implementation of pressure-correction scheme Markus Main CFD Forum 4 January 6, 2007 02:53
MUSCL scheme B. Nilsson Main CFD Forum 1 October 30, 1998 03:11

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