Changes between Version 4 and Version 5 of u/erica/ApproximateRS
 Timestamp:
 05/13/13 13:55:54 (11 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

u/erica/ApproximateRS
v4 v5 1 1 = Approximating pstar and ustar = 2 2 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.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 (ERS). 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 5 = Linearization of the Euler equations / Primitive Variable Riemann Solver = 6 6 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.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 (ERS), 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. … … 15 15 == TwoRarefaction == 16 16 17 This method is computationally costly (in terms of cost it goes  TRRS > TSRS > PVRS, where TRRS = tworarefaction Riemann solver, TSRS = twoshock Riemann solver, PVRS = primitive variable Riemann solver) 17 This method is computationally costly (in terms of cost it goes  TRRS > TSRS > PVRS, where TRRS = tworarefaction Riemann solver, TSRS = twoshock Riemann solver, PVRS = primitive variable Riemann solver), although extremely robust. If you assume 2 rarefactions, then you can derive closedform expressions for pstar and ustar from the pressure function formalism for rarefactions, given in chapter 4 of Toro. From there, you would feed these values of pstar and ustar into the sampling routine of the ERS, and continue with Godunov method as before. Assuming there are actually 2 rarefaction waves, this method would produce the exact solution. Even if there are not 2 rarefactions, this method produces sufficient results. 18 18 19 19 == TwoShock == 20 20 21 Unlike the TRRS, there does not exist an exact closed form expression for pstar and ustar, derivable from the pressure functions for the case of a twoshock scenario. Instead, one could approximate one of the unknowns in the equations equal to a constant, and then the system would be solvable. This leads to robust results that are not too computationally costly. 22 21 23 = Hybrid Schemes = 22 24 23 == Iterative == 24 25 == NonIterative == 25 One could make use of various approximation methods in different regions of the flow as needed. In regions of smooth flow, without sharp gradients, one could use the computationally cheap (and also least robust) PVRS methods. In regions with sharp gradients (near reflective boundaries, or contacts, shocks, and shears), one would opt for the more robust methods of TRRS or TSRS. Such adaptive, or hybrid schemes, can either be iterative or not, meaning near sharp gradients, pstar can either be solved for using the iteration method of the ERS, or estimated using one of the TSRS or TRRS, whereas both iterative and noniterative schemes would use the same simple PVRS in smoother flow regions. It is interesting to note that if one were to choose an iterative hybrid scheme, that most times only 1 iteration step produces sufficient results, provided there is a good enough initial estimate for pEst. 26 26 27 27 = Results and Discussion = 28  Same as God 29  Most problems use a simple over 90% of grid, and a more sophisticated over the disc 30  Can further improve the results using approximation methods to the God method (later chapters) 28 Numerical results are presented at the end of chapter 9. These results show that the 2 methods tested, 1) TSRS everywhere, and 2) a noniterative hybrid scheme using PVRS in smooth flow, and TSRS in sharp gradient flow, both produce results that exactly match up to the Godunov method that uses the ERS. This shows that 1) you can get the same order of accuracy of the Godunov method, using computationally cheaper methods, and 2) using a noniterative TSRS method can further reduce the computational cost without reducing accuracy of the solution. It is interesting to note that according to Toro, most problems use are able to use a simple PVRS over 90% of grid, and a more sophisticated scheme over the remaining domain. 29 30 = Further directions = 31 These schemes can be further improved by adding approximation methods to the Godunov method, to be discussed in later chapters such as 15.