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?
Discontinuous Galerkin for convective heat transfer problems
Discontinuous Galerkin for convective heat transfer problems
- Attachments
-
- hexInflation2D.zip
- Mesh
- (21.25 KiB) Downloaded 230 times
-
- heatSolver-Elmer.sif
- Case file
- (3.04 KiB) Downloaded 257 times
-
- Posts: 2312
- Joined: 25 Jan 2019, 01:28
- Antispam: Yes
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!
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
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
Faculty of Mechanical Engineering, Silesian University of Technology, Gliwice