## Discontinuous Galerkin for convective heat transfer problems

Numerical methods and mathematical models of Elmer
vencels
Posts: 26
Joined: 20 Sep 2016, 17:05
Antispam: Yes
Location: Latvia
Contact:

### Discontinuous Galerkin for convective heat transfer problems

Hi,

I am trying DG for simple convective heat transfer problem - curved pipe, 2D

Result (temperature dstribution):

Expected result:

Any ideas why this could happen?
Attachments
hexInflation2D.zip
Mesh
heatSolver-Elmer.sif
Case file

kevinarden
Posts: 319
Joined: 25 Jan 2019, 01:28
Antispam: Yes

### Re: Discontinuous Galerkin for convective heat transfer problems

temperature.png (84.98 KiB) Viewed 452 times
case.sif

vencels
Posts: 26
Joined: 20 Sep 2016, 17:05
Antispam: Yes
Location: Latvia
Contact:

### Re: Discontinuous Galerkin for convective heat transfer problems

Although it is not Discontinuous Galerkin solution, your solution helped me to find a bug in my viscosity value. I was using kinematic instead of dynamic. Thanks!

annier
Posts: 1123
Joined: 27 Aug 2013, 13:51
Antispam: Yes

### Re: Discontinuous Galerkin for convective heat transfer problems

Juris,
I would say it is difficult to get convergent solutions for such a low value of viscosity for transient simulations (asssuming the original kinematic viscosity magnitude in the order of 10^{-7} supplied into the N-S equation is interpreted as dynamic viscosity by ElmerSolver ). In your case you have performed steady state simulations, and so the solutions were obtained without difficulty.
Yes, the viscosity in Elmer is always supplied to solver input file as dynamic viscosity (Pa s) , and not as kinematic viscosity (m^{2}/s).
The unit of mu (Div(velocity)) must return the unit of Div(stress), and this is true for the unit of viscosity (mu) in [Pa S] or [ Ns/m^{2}].

Yours Sincerely,
Anil Kunwar
Anil Kunwar
Department of Materials Engineering, KU Leuven, Belgium