NACA Regression Tests
Let's test some of the ProjectX solvers for the case of a NACA0012 (? check ?) airfoil (2D) in compressible, inciscid flow:
Inputs:
- Common flow parameters:
- Freestream Mach number: M ∞ = 0.5
- Angle of attack: α = 0 °
- Total pressure: Pt =
- Total temperature: Tt =
- Static temperature:
- CFL: 1030
- Common solver parameters:
- Preconditioning:
- Outer GMRES iterations (no. restarts): 20
- Inner GMRES iterations (no. Krylov vectors): 200
- Solution order: 2
- PSequencing: False
- Preconditioning:
Tests on linear solvers:
Newton-Jacobi (Left):
- Solvers:
- Nonlinear: Newton
- Preconditioner: Block Jacobi (left)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .53 counts
- Lift coefficient: CL = .96 counts
Newton-Jacobi (Right):
- Solvers:
- Nonlinear: Newton
- Preconditioner: Block Jacobi (right)
- Results:
- Residual norm converges to 4.35·10^-13^ after 5 iterations
- Drag coefficient: CD = .53 counts
- Lift coefficient: CL = .96 counts
Newton-ILU
- Solvers:
- Nonlinear: Newton
- Preconditioner: Incomplete LU (0) (left)
- Results:
- Residual norm converges to 4.357·10^-13^ after 5 iterations
- Drag coefficient: CD = .5318 counts
- Lift coefficient: CL = .9638 counts
Newton-ILU (Parallel)
- Solvers:
- Nonlinear: Newton
- Preconditioner: Jacobi (right)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .5318 counts
- Lift coefficient: CL = .9638 counts
Newton-ILUT
- Solvers:
- Nonlinear: Newton
- Preconditioner: Threshold-based Incomplete LU (ILUT) (right)
- Reordering type: Block Minimum Discarded Fill (MDF)
- Fill level: 2
- Drop threshold: 10^-3^
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .53 counts
- Lift coefficient: CL = .96 counts
Newton-UMFPACK
- Solvers:
- Nonlinear: Newton
- Preconditioner: UMFPACK (right)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .5318 counts
- Lift coefficient: CL = .9638 counts
Newton-ILU-ASM-MG (Parallel)
- Solvers:
- Nonlinear: Newton
- Preconditioner: Jacobi (left)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .53 counts
- Lift coefficient: CL = .96 counts
Newton-ILU-ASM (Parallel)
- Solvers:
- Nonlinear: Newton
- Preconditioner: Jacobi (right)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .5318 counts
- Lift coefficient: CL = .9638 counts
Lean Jacobi
- Solvers:
- Nonlinear: Newton
- Preconditioner: Jacobi (left)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .53 counts
- Lift coefficient: CL = .96 counts
Newton-Osher
- Solvers:
- Nonlinear: Newton
- Preconditioner: Jacobi (right)
- Results:
- Residual norm converges to 4.36·10^-13^ after 5 iterations
- Drag coefficient: CD = .5318 counts
- Lift coefficient: CL = .9638 counts