$ prefix = "ISOsert" $ visc = 1.0e-10 Header CHECK KEYWORDS Warn Mesh DB "." "Mesh_1" Include Path "" Results Directory "" End Simulation Max Output Level = 7 Coordinate System = Cartesian Coordinate Mapping(3) = 1 2 3 Coordinate Scaling = 0.001 Simulation Type = Steady State Steady State Max Iterations = 5 Output Intervals = 0 ! Post File = $prefix$.ep Mesh Levels = 2 Partition Mesh = Logical True Partitioning Method = String "Zoltan" End Constants Gravity(4) = 0 -1 0 9.82 Stefan Boltzmann = 5.67e-08 Boltzmann Constant = 1.3807e-23 End Body 1 Target Bodies(1) = 1 Name = "WallExtLayer" Equation = 1 Material = 1 End Body 2 Target Bodies(1) = 2 Name = "Insulation" Equation = 1 Material = 2 End Body 3 Target Bodies(1) = 3 Name = "HorizontalConcrete" Equation = 1 Material = 3 End Body 4 Target Bodies(1) = 4 Name = "AlphaWallIntLayer" Equation = 1 Material = 4 End Body 5 Target Bodies(1) = 5 Name = "BetaWallIntLayer" Equation = 1 Material = 4 End Body 6 Target Bodies(1) = 6 Name = "FloorFinishing" Equation = 1 Material = 5 End ! Heat equation for the iterative solution of natural convection !---------------------------------------------------------------- Solver 1 Equation = "Heat Equation" Procedure = "HeatSolve" "HeatSolver" Variable = -dofs 1 Temperature Exec Solver = Always Stabilize = True Bubbles = False Lumped Mass Matrix = False Optimize Bandwidth = True Steady State Convergence Tolerance = 1.0e-8 Nonlinear System Convergence Tolerance = 1.0e-8 Nonlinear System Max Iterations = 20 Nonlinear System Newton After Iterations = 3 Nonlinear System Newton After Tolerance = 1.0e-5 Nonlinear System Relaxation Factor = 0.2 Linear System Solver = Iterative Linear System Iterative Method = BiCGStabL BiCgStabl Polynomial Degree = 4 Linear System Max Iterations = 500 Linear System Convergence Tolerance = 1.0e-8 Linear System Preconditioning = ILU1 Linear System ILUT Tolerance = 1.0e-4 Linear System Abort Not Converged = False Linear System Residual Output = 10 Linear System Precondition Recompute = 1 Calculate Loads = Logical True End ! This solver calculates the flux vector in the whole domain ! using the Galerkin method from equation q=-\kappa \grad T !--------------------------------------------------------------- Solver 2 Exec Solver = after all Equation = "flux compute" Procedure = "FluxSolver" "FluxSolver" Calculate Flux = Logical True Flux Variable = String Temperature Flux Coefficient = String "Heat Conductivity" Linear System Solver = "Iterative" Linear System Iterative Method = "cg" Linear System Preconditioning = ILU0 Linear System Residual Output = 10 Linear System Max Iterations = Integer 500 Linear System Convergence Tolerance = 1.0e-10 Linear System Abort Not Converged = False End ! This solver saves the results in other formats ! Perhaps most importantly, the vtu format for Paraview !------------------------------------------------------------ Solver 3 Exec Solver = after all Equation = "result output" Procedure = "ResultOutputSolve" "ResultOutputSolver" Output File Name = $prefix$ ! Enforce Parallel Format = Logical True No Fileindex = Logical True Binary Output = Logical True Single Precision = Logical True Vtu Format = Logical True ! ElmerPost Format = Logical True ! GiD Format = Logical True ! Gmsh Format = Logical True Scalar Field 1 = String Temperature ! Scalar Field 2 = String Pressure ! Vector Field 1 = String Velocity End ! This solver computes the boundary flux and saves them ! to external files, to be read by Matlab, for example. !----------------------------------------------------- !Solver 4 ! Exec Solver =after all ! Equation = "save line" ! Procedure = "SaveData" "SaveLine" ! Filename = $prefix$_line.dat ! Polyline Coordinates(2,3) = 0.0 0.0 0.0 1.0 1.0 1.0 ! Save Flux = Logical True ! Flux Variable = String Temperature ! Flux Coefficient = String "Heat Conductivity" !End ! This solver computes the total flux in two different ways ! 1) using integration points of the boundary ! 2) summing up the nodal heat loads from the residual, r=Ax-b ! Of these the latter is more accurate. !------------------------------------------------------------- Solver 4 Exec Solver = after all Equation = "save scalars" Procedure = "SaveData" "SaveScalars" Filename = $prefix$_tot.dat Operator 1 = "diffusive flux" Variable 1 = "Temperature" Coefficient 1 = "Heat Conductivity" Operator 2 = "boundary sum" Variable 2 = "Temperature Loads" End Equation 1 Name = "WallHeatTransferCalculation" Active Solvers(2) = 1 2 End Equation 1 Name = "Equation 1" ! Convection = None ! Active Solvers(7) = 1 2 3 4 5 6 7 Active Solvers(1) = 1 End Material 1 Name = "WallExtLayer" Heat Conductivity = 1 Heat Capacity = 1000 Density = 2000 End Material 2 Name = "Insulation" Heat Conductivity = 0.04 Heat Capacity = 1000 Density = 200 End Material 3 Name = "HorizontalConcrete" Heat Conductivity = 2.5 Heat Capacity = 600 Density = 5000 End Material 4 Name = "WallIntLayer" Heat Conductivity = 0.7 Heat Capacity = 800 Density = 1700 End Material 5 Name = "FloorFinishing" Heat Conductivity = 1 Heat Capacity = 800 Density = 1000 End Boundary Condition 1 Target Boundaries(3) = 39 40 41 Name = "INdorUpperFloorBeta" Heat Flux BC = True Save Line = Logical True Save Scalars = Logical True Heat Transfer Coefficient = 5 External Temperature = 288.15 End Boundary Condition 2 Target Boundaries(3) = 33 34 35 Name = "INdorLowerFloorAlpha" Heat Flux BC = True Save Line = Logical True Save Scalars = Logical True Heat Transfer Coefficient = 5 External Temperature = 293.15 End Boundary Condition 3 Target Boundaries(6) = 1 2 11 12 13 14 Name = "OutdoorGamma" Heat Flux BC = True Save Line = Logical True Save Scalars = Logical True Heat Transfer Coefficient = 20 External Temperature = 273.15 End