NACA Regression Tests

Let's test some of the ProjectX solvers for the case of a NACA0012 (? check ?) airfoil (2D) in compressible, inciscid flow:

  • 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

Tests on linear solvers:

  • 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
  • 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
  • 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
  • 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
  • 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
  • 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
  • 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
  • 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
  • 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
  • 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