摘要
An implicit non-linear lower-upper symmetric Gauss-Seidel(LU-SGS)solution algorithm has been developed for a high-order spectral difference Navier-Stokes solver on unstructured hexahedral grids.The non-linear LU-SGS solver is preconditioned by a block element matrix,and the system of equations is then solved with the LU decomposition approach.The large sparse Jacobian matrix is computed numerically,resulting in extremely simple operations for arbitrarily complex residual operators.Several inviscid and viscous test cases were performed to evaluate the performance.The implicit solver has shown speedup of 1 to 2 orders of magnitude over the multi-stage Runge-Kutta time integration scheme.
基金
The study was partially funded by Rockwell Scientific/DARPA contract W911NF-04-C-0102,AFOSR grant FA9550-06-1-0146,and DOE grant DE-FG02-05ER25677
The views and conclusions contained herein are those of the authors and should not be interpreted as necessarily representing the official policies or endorsements,either expressed or im-plied,of DARPA,AFOSR,DOE,or the U.S.Government.
