Versions Compared

Old Version 12

changes.mady.by.user Philip Claude Caplan

Saved on

compared with

New Version 13

changes.mady.by.user Philip Claude Caplan

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

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

Section
Column
width50%
Panel
titleInputs:
  • 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
Column
width50%
Wiki Markup

{center}
{panel:title=Boundary conditions:}
!nacabcs.png|border=1,scale=.6,width=388,height=292!
{panel}
{center}

Tests on linear solvers:

Section
Column
width50%
Panel
titleNewton-Jacobi (Left):
  • Nonlinear solver:
    • Newton
  • Preconditioner:
    • Block Jacobi (left)
  • Serial results:
    • Residual norm converges to 4.3571·103571·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
Column
width50%
Panel
titleNewton-Jacobi (Right):
  • Nonlinear solver:
    • Newton
  • Preconditioner:
    • Block Jacobi (right)
  • Results:
    • Residual norm converges to 4.3532·103532·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
Section
Column
width50%
Panel
titleNewton-ILU
  • Nonlinear solver:
    • Newton
  • Preconditioner:
    • Incomplete LU (0) (left)
    • Reordering type: Minimum discarded fill
    • Side: Left
  • Serial results:
    • Residual norm converges to 4.3566·103566·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
  • Parallel results:
    • Residual norm converges to 4.3501·103501·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
Column
width50%
Panel
titleNewton-ILUK
  • Nonlinear solver:
    • Newton
  • Preconditioner:
    ** ?? Incomplete LU (ILUK)
    • Reordering type: Minimum discarded fill
    • Side: Left
  • Serial results:
    • Residual norm converges to 4.??·10·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .53 counts
    • Lift coefficient: CL = .96 counts
  • Parallel results:
    • Residual norm converges to 4.3545·103545·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
Section
Column
width50%
Panel
titleNewton-UMFPACK
  • Nonlinear solver:
    • Newton
  • Preconditioner:
    • Threshold-based Incomplete LU
    • Side: Right
    • Reordering type: Minimum discarded fill
  • Serial results:
    • Residual norm converges to 4.3559·103559·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
  • Parallel results:
    • Residual norm converges to 4.3506·103506·10^-1313^ after 5 iterations
    • Drag coefficient: CD = .5318 counts
    • Lift coefficient: CL = .9638 counts
Column
width50%
Panel
titleLean Jacobi
  • Nonlinear solver:
    • Lean Jacobi
  • Preconditioner:
    • None
  • Serial results:
    • Residual norm converges to 4.7345·107345·10^-44^ after 200 iterations
    • Drag coefficient: CD = .5098 counts
    • Lift coefficient: CL = .9478 counts
  • Parallel results:
    • Residual norm converges to 4.7345·107345·10^-44^ after 200 iterations
    • Drag coefficient: CD = .5098 counts
    • Lift coefficient: CL = .9478 counts