Changes between Version 3 and Version 4 of u/erica/ApproximateRS
- Timestamp:
- 05/13/13 13:35:49 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/ApproximateRS
v3 v4 3 3 The Godunov method relies on the solution of pressure and velocity inside of the star region set up within local Riemann problems, which before now, has been done using an iterative approach inside of the Exact Riemann Solver. The iterative approach is computationally costly, and so if we could devise schemes which approximate pstar and ustar, we can make a more efficient time marching procedure. 4 4 5 = Linearization of the Euler equations =5 = Linearization of the Euler equations / Primitive Variable Riemann Solver = 6 6 7 7 There are 2 methods which make use of approximating the Euler equations themselves in linearized form, to get simple algebraic expressions for pstar and ustar. The first makes a substitution in the coefficient matrix of the Euler equations, such that all fluid variables are replaced with constants, making the system of equations a constant coefficient system, which can be solved exactly algebraically. A obvious choice for a constant, may be the mean of the left and right state's variable to be substituted, (i.e. rho = mean(rho_L, rho_R)). The 2nd method for this type of approximation is matching states using the characteristic equations, to solve for the star regions variables (p, u, and rho). After the star region is specified, the variables are fed into the sampling routine of the Exact Riemann Solver, slightly modified as it no longer has to solve for rho in the star region, and the Godunov method continues as before. 8 8 waves generated 9 9 These methods for approximation are the cheapest computationally, but also the least robust. You would not want to use such a scheme over the entire computational domain for most practical problems. 10 10 11 11 = Non-linear Wave Assumption Approximations = 12 12 13 If one assumes the generation of either 2 shocks or 2 rarefaction from the Riemann problem, then there are 2 possibilities for approximating pstar, from which follow ustar and rho. 14 13 15 == Two-Rarefaction == 16 17 This method is computationally costly (in terms of cost it goes - TRRS > TSRS > PVRS, where TRRS = two-rarefaction Riemann solver, TSRS = two-shock Riemann solver, PVRS = primitive variable Riemann solver) 14 18 15 19 == Two-Shock ==