Discontinuous Galerkin for convective heat transfer problems

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

Discontinuous Galerkin for convective heat transfer problems

Post by vencels »

Hi,

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

Result (temperature dstribution):
Image

Expected result:
Image

Any ideas why this could happen?
Attachments
hexInflation2D.zip
Mesh
(21.25 KiB) Downloaded 227 times
heatSolver-Elmer.sif
Case file
(3.04 KiB) Downloaded 253 times
kevinarden
Posts: 2237
Joined: 25 Jan 2019, 01:28
Antispam: Yes

Re: Discontinuous Galerkin for convective heat transfer problems

Post by kevinarden »

temperature.png
temperature.png (84.98 KiB) Viewed 3172 times
case.sif
(3.36 KiB) Downloaded 286 times
vencels
Posts: 63
Joined: 20 Sep 2016, 17:05
Antispam: Yes
Location: Latvia
Contact:

Re: Discontinuous Galerkin for convective heat transfer problems

Post by vencels »

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: 1168
Joined: 27 Aug 2013, 13:51
Antispam: Yes

Re: Discontinuous Galerkin for convective heat transfer problems

Post by annier »

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
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice
Post Reply