Code: Select all

```
Simulation
Max Output Level = 3
Coordinate System = "Cartesian 2D"
Coordinate Mapping(3) = 1 2 3
Simulation Type = "transient"
Timestep Sizes(1) = 1e-3
Timestep Intervals(1) = 60
Timestepping Method = BDF
BDF Order = 1
Steady State Max Iterations = 10
Solver Input File = TBHSModel1_melting.sif
End
Solver 1
Exec Solver = "Always"
Equation = "Heat Equation"
Procedure = "HeatSolve" "HeatSolver"
Variable = "Temperature"
Variable Dofs = 1
Optimize Bandwidth = True
Linear System Solver = "iterative"
Linear System Iterative Method = "BiCGStabl"
Linear System Max Iterations = 1000
Linear System Convergence Tolerance = 1.0e-08
Linear System Abort Not Converged = True
Linear System Preconditioning = "ILU1"
Linear System Residual Output = 10
Steady State Convergence Tolerance = 1.0e-6
Stabilize = False
Bubbles = False
Lumped Mass Matrix = False
Nonlinear System Convergence Tolerance = 1.0e-5
Nonlinear System Max Iterations = 5
Nonlinear System Newton After Iterations = 3
Nonlinear System Newton After Tolerance = 1.0e-3
Nonlinear System Relaxation Factor = 1.0
End
```

Code: Select all

```
HeatSolve:
HeatSolve: Starting Assembly...
HeatSolve: Assembly done
10 0.4729E-07 0.4729E-07
16 0.9798E-08 0.9798E-08
ComputeChange: NS (ITER=5) (NRM,RELC): ( 59.312385 0.22720930E-04 ) :: heat equation
```