Open Source CAE Toolchain 11 - Salome>Elmer

Salome_Meca를 Linux에 설치하고,

여기서 Mesh를 생성한 후,

그것을 Elmer에서 읽어들여 테스트 해 보았다. 

Salome 자체는 큰 문제없이 실행되는데,

아직 여기 포함된 Code_Aster가

제대로 작동되도록 하는데는 성공하지 못하고 있다.

Solver Input 설정에서 뭔가 잘못되어 있거나

또는 Mesh에서 문제가 있는 것 같은데 

구체적인 원인은 못 찾고 있다.

아무튼 대신 Elmer로 넘겨서 해석해 보기 위해

*.unv 파일로 Export한 후, 이걸 Elmer에서 Open 하였다.

빠르게 읽어들여진다.

Gmesh에서 Mesh를 만들어내는 것도 좋지만,

Boundary Condition을 만들기 위한 Group 작업을 미리 해 두거나

Multi-body를 위한 다중 Mesh를 만들어내기 위해서는

Salome를 사용하는 것이 더 낫다고 생각된다.

또 Salome로는

기존 Netgen 기반의 Tetrahdron 분만 아니라,

i,j,k Hexahedron을 생성해 낼 수도 있다.

물론 Hexa로 만들어낼 때는

실패 확률이 Tetra 보다는 높기 때문에,

형상을 많이 정돈해 주는 것이 좋겠다.

아무튼 간단히 Heat Eqation만 적용해서 Solve 해 본다.

case.sif 파일은 아래와 같다.

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

Simulation
  Max Output Level = 5
  Coordinate System = Cartesian
  Coordinate Mapping(3) = 1 2 3
  Simulation Type = Steady state
  Steady State Max Iterations = 1
  Output Intervals = 1
  Timestepping Method = BDF
  BDF Order = 1
  Solver Input File = case.sif
  Post File = case.ep
End

Constants
  Gravity(4) = 0 -1 0 9.82
  Stefan Boltzmann = 5.67e-08
  Permittivity of Vacuum = 8.8542e-12
  Boltzmann Constant = 1.3807e-23
  Unit Charge = 1.602e-19
End

Body 1
  Target Bodies(1) = 2
  Name = "Body 1"
  Equation = 1
  Material = 1
  Body Force = 1
End

Solver 1
  Equation = Heat Equation
  Procedure = "HeatSolve" "HeatSolver"
  Variable = Temperature
  Exec Solver = Always
  Stabilize = True
  Bubbles = False
  Lumped Mass Matrix = False
  Optimize Bandwidth = True
  Steady State Convergence Tolerance = 1.0e-5
  Nonlinear System Convergence Tolerance = 1.0e-7
  Nonlinear System Max Iterations = 20
  Nonlinear System Newton After Iterations = 3
  Nonlinear System Newton After Tolerance = 1.0e-3
  Nonlinear System Relaxation Factor = 1
  Linear System Solver = Iterative
  Linear System Iterative Method = BiCGStab
  Linear System Max Iterations = 500
  Linear System Convergence Tolerance = 1.0e-10
  BiCGstabl polynomial degree = 2
  Linear System Preconditioning = Diagonal
  Linear System ILUT Tolerance = 1.0e-3
  Linear System Abort Not Converged = False
  Linear System Residual Output = 1
  Linear System Precondition Recompute = 1
End

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

Material 1
  Name = "Aluminium (generic)"
  Heat Conductivity = 237.0
  Youngs modulus = 70.0e9
  Mesh Poisson ratio = 0.35
  Heat Capacity = 897.0
  Density = 2700.0
  Poisson ratio = 0.35
  Sound speed = 5000.0
  Heat expansion Coefficient = 23.1e-6
End

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

Boundary Condition 1
  Target Boundaries(6) = 64 103 143 208 243 278 
  Name = "RoomTemp"
  Temperature = 293
End

나사 구멍으로 냉각수가 지나간다고 가정하고
상온으로 경계조건을 주었다.
열원은 예제에 따라 단순히 BodyForce로 주었다.