|June 26, 2010, 10:28||
MUSCL scheme implementation
J Kyle Miller
Join Date: Jun 2010
Posts: 2Rep Power: 0
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
|June 26, 2010, 23:57||
Join Date: Mar 2009
Blog Entries: 6Rep Power: 9
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.
|Thread||Thread Starter||Forum||Replies||Last Post|
|Trouble when using MUSCL scheme on ONERA M6 wing||Aurélien||FLUENT||2||January 23, 2013 08:42|
|Implementation of QUICK scheme||Romuald Skoda||Main CFD Forum||10||August 3, 2010 12:23|
|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|