QuarterTurn

In this example, we will estimate the rise in temperature due to Joules losses in a stranded conductor. An electrical potential V0 is applied to the entry/exit of the conductor which is also cooled by a force flow.
The geometry of the conductor is choosen as to have an analytical expression for the temperature.

1. Running the case

The command line to run this case in 2D is

mpirun -np 4 feelpp_toolbox_thermoelectric --case "github:{path:toolboxes/thermoelectric/cases/cvg}" --case.config-file 2d.cfg
mpirun
Case option
--case "github:{path:toolboxes/thermoelectric/cases/cvg}"
Case config file option
--case.config-file 2d.cfg

The command line to run this case in 3D is

mpirun -np 4 feelpp_toolbox_thermoelectric --case "github:{path:toolboxes/thermoelectric/cases/cvg}" --case.config-file 3d.cfg
mpirun
Case option
--case "github:{path:toolboxes/thermoelectric/cases/cvg}"
Case config file option
--case.config-file 3d.cfg

2. Geometry

The conductor consists in a rectangular cross section torus which is somehow "cut" to allow for applying electrical potential. The conductor is cooled with a force flow along its cylindrical faces.+ In 2D, the geometry is as follow: geometry
In 3D, this is the same geometry, but extruded along the z axis.

3. Input parameters

Name Description Value Unit

ri

internal radius

1

m

re

external radius

2

m

δ

angle

π/2

rad

VD

electrical potential

9

V

hi

internal transfer coefficient

60e3

Wm2K1

Twi

internal water temperature

303

K

he

external transfer coefficient

58e3

Wm2K1

Twe

external water temperature

293

K

3.1. Model & Toolbox

  • This problem is fully described by a Thermo-Electric model, namely a poisson equation for the electrical potential V and a standard heat equation for the temperature field T with Joules losses as a source term.

  • toolbox: thermoelectric

3.2. Materials

Name Description Marker Value Unit

σ

electric conductivity

omega

4.8e7

S.m1

3.3. Boundary conditions

The boundary conditions for the electrical probleme are introduced as simple Dirichlet boundary conditions for the electric potential on the entry/exit of the conductor. For the remaining faces, as no current is flowing througth these faces, we add Homogeneous Neumann conditions.

Marker Type Value

V0

Dirichlet

0

V1

Dirichlet

VD

Rint, Rext, top*, bottom*

Neumann

0

As for the heat equation, the forced water cooling is modeled by robin boundary condition with Tw the temperature of the coolant and h an heat exchange coefficient.

Marker Type Value

Rint

Robin

hi(TTwi)

Rext

Robin

he(TTwe)

V0, V1, top*, bottom*

Neumann

0

*: only in 3D

4. Outputs

The main fields of concern are the electric potential V, the temperature T and the current density j or the electric field E. // presented in the following figure.

5. Verification Benchmark

The analytical solutions are given by:

V=VDδθ=VDδatan2(y,x)E=(VDδyx2+y2,VDδxx2+y2)T=Alog(r)2+Blog(r)+C=Alog(x2+y2)2+Blog(x2+y2)+CA=σ2k(VDδ)2B=BeBiDC=CeCiDBe=2Tweδ2hehikreri+V2Dhehireriσlog(re)2+V2Dhikriσlog(r2e)Bi=2Twiδ2hehikreri+V2Dhehireriσlog(ri)2V2Dhekreσlog(r2i)Ce=(herelog(re)+k)(2Twiδ2hikri+V2Dhiriσlog(ri)2V2Dkσlog(r2i))Ci=(hirilog(ri)k)(2Tweδ2hekre+V2Dhereσlog(re)2+V2Dkσlog(r2e))D=2δ2k(hehirerilog(re)hehirerilog(ri)+hekre+hikri)

We will check if the approximations converge at the appropriate rate:

  • k+1 for the L2 norm for V and T

  • k for the H1 norm for V and T

  • k for the L2 norm for E and j

  • k-1 for the H1 norm for E and j

Table 1. Temperature 2D and 3D

temperature 2D

potential 2D

Table 2. Electric potential 2D and 3D

potential 2D

potential 3D

Table 3. Electric field 2D and 3D

electric field 2D

electric field 3D