A longitudinaly cooled helix with HDG method

In this test case we will simulate a longitudinaly cooled helix with the HDG method. We consider only the thermoelectric behaviour of the helix.

1. Running the case

mpirun -np 32 feelpp_toolbox_hdg_thermoelectric --config-file  HL-31_H1.cfg

If you use the IBC feature to set the total current:

mpirun -np 32 feelpp_toolbox_hdg_thermoelectric --config-file  HL-31_H1.cfg

The thermolelectric HDG method is detailed in this section.

If you want to change the number of procs used, you need to create a proper partition mesh. See the next sections for a detailed guide.

1.1. Convert the mesh

The mesh may be converted into an another format using gmsh compiled with proper med support

gmsh -0 -bin HL-31_H1.med -o HL-31_H1.msh

The resulting mesh is in mm. For convergence sake, HDG method requires the mesh to be left in mm. As a consequence, remember that the physical properties of materials shall be scaled to "mm".

1.2. Partitionning the mesh:

feelpp_mesh_partitioner --ifile HL-31_H1.msh [--ofile HL-31_H1_p32] --part 32

For calculation with IBC, it is mandatory to use the following command to properly partion the mesh:

feelpp_mesh_partitioner --ifile HL-31_H1[_shapes].med --ofile HL-31_H1_p  --part 32 --sc.ibc_partitioning=true --sc.ibc_partitioning.marker=<marker>

where <marker> is the marker of the Surface marker where we apply the IBC.

2. Data files

The data files may be retreived from this repository. The mesh files are stored in collection hifimanget/cases/v0.108/HL-31_H1-HDG collection on Unistra Girder

2.1. Mesh file

  • HL-31_H1-woglue.msh

  • HL-31_H1-woglue.med

  • HL-31_H1-woglue_p.json, HL-31_H1-woglue_p.j.h5

2.2. Simulation cfg file

  • HL-31_H1.cfg,

  • HL-31_H1_ibc.cfg with IBC.

2.3. Model files

  • HL-31_H1_nonlinear.json

The data are given to be consistant with the mesh in mm. If you ever want to use a mesh in m:

  • Change the physical values,

  • Change the Cooling Boundary Conditions values,

3. Input parameters

Name Description Value Unit

\(V_D\)

electrical potential

9

\(V\)

4. Outputs

Name Description Value Unit

\(min(T)\)

min of Temperature Field

293.32

\(K\)

\(max(T)\)

max of Temperature Field

400.781

\(K\)

\(I\)

current

28.771 (-28.768)

\(kA\)