Changes between Version 28 and Version 29 of u/ehansen/buildcode

Timestamp:

11/14/11 13:30:20 (13 years ago)

Author:

ehansen

Comment:

—

u/ehansen/buildcode

@@ -123 +123 @@
 
 [[BR]]
+== The MUSCL-Hancock Scheme ==
+
+[[BR]]
 == Program Outline ==
-These are the basic steps in my program SEEQUOD (Solver for Euler EQUations in One Dimension).  SEEQUOD currently uses either an exact Riemann solver, an HLLC solver, or a Roe solver.  SEEQUOD currently uses a first order Godunov upwind scheme, but will later be updated to use higher order update schemes.  These steps should be general enough to apply to many codes that are used to solve the Euler Equations.
+These are the basic steps in my program SEEQUOD (Solver for Euler EQUations in One Dimension).  SEEQUOD currently uses either an exact Riemann solver, an HLLC solver, or a Roe solver.  SEEQUOD can be easily modified to use either the aforementioned first order Godunov scheme or the second order MUSCL-Hancock scheme.  These steps should be general enough to apply to many codes that are used to solve the Euler Equations.
 
 * Read all input data necessary for problem
@@ -157 +160 @@
 Test 1 was also run with the Roe solver in two different ways.  First as a "barebones" Roe solver, and second with an entropy fix.  This entropy fix is required for transonic rarefactions, which is the case for the left rarefaction in Test 1. You can see that without this entropy fix, the gap inside that left rarefaction is very large and completely inaccurate.
 
+Test 1 was also run with the MUSCL-Hancock scheme with the HLLC solver.  A comparison with the first order plots shows that the MUSCL-Hancock method does indeed give more accurate results. 
+
 For all the plots, the numerical results (points) are alongside the exact solution (line).  
 
@@ -190 +195 @@
 || [[Image(roe_nofix_test1_density.png, width=500)]] || [[Image(roe_withfix_test1_density.png, width=500)]] ||
 [[CollapsibleEnd]]
+[[CollapsibleStart(Test 1 (Godunov vs. MUSCL-Hancock)]]
+||= Godunov =||= MUSCL-Hancock =||
+||= [[Image(hllc_test1_density.png, width = 500)]] =||= [[Image(muscl_test1_density.png, width = 500)]] =||
+||= [[Image(hllc_test1_velocity.png, width = 500)]] =||= [[Image(muscl_test1_velocity.png, width = 500)]] =||
+||= [[Image(hllc_test1_pressure.png, width = 500)]] =||= [[Image(muscl_test1_pressure.png, width = 500)]] =||
+||= [[Image(hllc_test1_intenergy.png, width = 500)]] =||= [[Image(muscl_test1_intenergy.png, width = 500)]] =||
+[[CollapsibleEnd]]