# MUSCL scheme implementation

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

 June 26, 2010, 10:28 MUSCL scheme implementation #1 New Member   J Kyle Miller Join Date: Jun 2010 Posts: 2 Rep Power: 0 This is a very basic question, regarding MUSCL schemes as described in http://en.wikipedia.org/wiki/MUSCL_scheme 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 -ANPDEstudent AshwaniAssam likes this.

 June 26, 2010, 23:57 #2 Super Moderator     Praveen. C Join Date: Mar 2009 Location: Bangalore Posts: 256 Blog Entries: 6 Rep Power: 10 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.

 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 Aurélien FLUENT 2 January 23, 2013 08:42 Romuald Skoda Main CFD Forum 10 August 3, 2010 12:23 Chandra Main CFD Forum 6 February 14, 2007 12:21 Markus Main CFD Forum 4 January 6, 2007 02:53 B. Nilsson Main CFD Forum 1 October 30, 1998 03:11

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