Changes between Version 80 and Version 81 of u/erica/RoeSolver

Timestamp:

05/20/13 17:26:29 (12 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/RoeSolver

@@ -212 +212 @@
 
 
-= The code = 
-
+= The code - Non-Entropy Fix= 
+To get a 1D version of the method Toro was describing, without an 'entropy fix', I had the intuition to do the following- keep the flux equation the same, but ignore certain vectors and eigenvalues associated with y and z directions. This got rid of alpha's 3 and 4, K's 3 and 4, and lambda's 3 and 4. I changed the remaining  eigenvectors to have 3 components each, again choosing to get rid of the y and z elements. This is the resulting code that I added to the Godunov solver. Unfortunately it did not work, and produced wonky results, see below. 
+ 
 [[Image(CallRoe.png, 50%)]]
 [[Image(rhoSolver.png, 100%)]]
 
 = Results =
+The rho result should look like the result presented above in the shortcomings section. The reason it doesn't may be due to modifying the eigen-structure of the equations, while keeping the form of the numerical flux the same... Another potential problem that I can think of other than typos, which seem to be absent, would be that the numerical flux isn't in terms of conserved quantities. But again, this seems to not be right given the theory that developed the equation, drew heavily upon conservative methods and integration.. I would be curious to add in the full expressions to this solver and see if it fixes things. However, I am planning on putting this in astrobear, so will instead focus on just writing an entropy correcting Roe Solver, which uses a different flux function given in Toro. 
 
 [[Image(RoeRho.png, 35%)]]