case.sif

!----------------------------------------------------------------
! Test case for automated computation of boundary loads for each 
! boundary. This is the most efficient way to compute fluxes when
! it is possible. There are limitations related to cases where
! nodes are shared between different boundaries. In such a case the 
! procedure cannot separate between the fluxes on that node from 
! different boundaries.
! 
! Similar procedure would be available through "Calculate Loads" and
! summing of nodal loads through "boundary sum" operator of SaveScalars.
! This one sweep operation is somewhat easier to use. 
!
! P.R. / 26.11.2014
!-----------------------------------------------------------------

Header
  CHECK KEYWORDS Warn
  Mesh DB "." "angle"
  Include Path ""
  Results Directory ""
End

Simulation
  Max Output Level = 20

  Coordinate System = "Cartesian"
  Coordinate Mapping(3) = 1 2 3

  Simulation Type = "Steady State"
  Steady State Max Iterations = 1
  Output Intervals = 1

!  Post File = "case.vtu"
End

Constants
  Gravity(4) = 0 -1 0 9.82
  Stefan Boltzmann = 5.67e-08
End

Body 1
  Name = "Body"
  Body Force = 1
  Equation = 1
  Material = 1
End

Equation 1
  Name = "Equations"
  Active Solvers(2) = 1 2
End

Solver 1
  Exec Solver = "Always"
  Equation = "Heat Equation"
  Variable = "Temperature"
  Variable Dofs = 1
  Linear System Solver = "Iterative"
  Linear System Iterative Method = "BiCGStab"
  Linear System Max Iterations = 350
  Linear System Convergence Tolerance = 1.0e-08
  Linear System Abort Not Converged = True
  Linear System Preconditioning = "ILU0"
  Linear System Residual Output = 1
  Steady State Convergence Tolerance = 1.0e-05
  Stabilize = True
  Nonlinear System Convergence Tolerance = 1.0e-05
  Nonlinear System Max Iterations = 1
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-02
  Nonlinear System Relaxation Factor = 1.0

  Calculate Boundary Fluxes = Logical True
End


Solver 2
  Equation = SaveScalars
  Procedure = "SaveData" "SaveScalars"

! This will make the heat flux through boundary 2 the value 
! to monitor in the consistency test. SaveScalars is only needed
! to implement the test, not to compute the flux. 
  Show Norm Index = 2
!  Filename = f.dat
End

Material 1
  Name = "Material"
  Density = 1
  Heat Conductivity = 1.0
End

Body Force 1
  Name = "BodyForce"
  Heat Source = 0
End

Boundary Condition 1
  Name = "Int"
  Target Boundaries(1) = 1 
  Temperature = 0.0
End

Boundary Condition 2
  Name = "Ext"
  Target Boundaries(1) = 2 
  Temperature = 1.0
End



Solver 2 :: Reference Norm = Real 2.5655331
RUN
!End Of File