Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Comment: Migration of unmigrated content due to installation of a new plugin
Section
Column
width50%
Panel
titleCommon Case Information:Overview

This circle boundary layer case demonstrates the method of manufactured solutions in which an analytical heat source is derived from an analytical expression of the thermal boundary layer. This derived heat source is applied to a vortex flow around a circular annulus of radius r = [ 1 , 2 ], and the numerical heat flux at the inner boundary is checked with the analytical heat flux to a specified tolerance. All boundary temperatures are set to the analytical thermal boundary layer solution, and the initial condition is set to 1.

Column
width50%
Panel
titleTemperature Distribution:

Image Added

Note that the singularity at (θ=0) is disregarded.

Section
Column
width50%
  • Flow parameters:
    • Vortex Flow
    • Heat source on inner radius boundary
    • Non-dimensional Viscosity: 10^-2^
  • Solver parameters
    Panel
    titleSolver Parameters
    :
    • Equation Set: Scalar 2D
    • Nonlinear Solver: Newton
    • Basis Type: Hierarchical
    • Initial/Max CFL: 1030
    • target Target residual: 10^10-11^11
    • Preconditioning:
      • Type: ILU (Left)
      • Reordering: MDF
      • Outer GMRES iterations (# restarts): 20
      • Inner GMRES iterations (# Krylov vectors): 200

    Tested Outputs

    • L2 L2 error (compared to an exact solution function)
    •  
    • Heat Numerical heat flux on inner boundary
    • Exact heat flux error on inner boundary
      Image Added
    • L Linf of the exact heat flux error
    Column
    width50%
    conditions
    Panel
    titleFlow Conditions:
    • Scalar Convective Diffusion Equation
      Image Added
    • Vortex Flow
      Image Added
      where
      Image Added
    • Boundary
    • Condition:
    • Picture of Source Boundary
    • Solution Order = 2 (non-adaptive)
      • Heat Flux Adjoint Residual (Mesh Adaptive Case)

    adaptation

    • Solution Order =1 
    • Adaptation Iterations = 2
    • Adaptation Method = Fixed Fraction
    • Output Adapted = Heat Flux
    • Anisotropy Method = Hessian

    results or convergence

    • Both cases reach tolerances in 2 iterations
    • Compared to truth files with tolerance of 
      • tol: e-12
      • Adapt tol : e-10

    Case Details

    Scalar2d_CircleBoundaryLayer

    Scalar2d_CircleBoundaryLayer_Adapt

    • Image Added
      where
      Image Added
    • Derived Heat Source:
      Image Added
      where
      Image Added
    • Non-dimensional viscosity:
      Image Added

    Case Specific Details

    Section
    Column
    width50%
    Panel
    titleCase 1 (Non-Adaptive):
    • Solution Order: 2
    • Iterations to solve: 2
    • All outputs tested to a tolerance of 10-12

    Job File:
    Scalar2d_CircleBoundaryLayer

    Column
    width50%
    Panel
    titleCase 1 Mesh

    Image Added

    Section
    Column
    width50%
    Panel
    titleCase 2 (Adaptive):
    • Solution Order = 1
    • Additional tested outputs: Heat flux adjoint residual
    • Adaptation parameters:
      • Adaptation Iterations: 2
      • Adaptation Method: Fixed fraction
      • Output Adapted: Heat flux
      • Anisotropy Method: Hessian
    • Iterations to solve on each mesh: 2
    • All outputs tested to a tolerance of 10-10

    Job File:
    Scalar2d_CircleBoundaryLayer_Adapt

    Column
    width50%
    Panel
    titleCase 2 Adapted Meshes:

    Image Added
    Image Added

    ...