halbux / sparselizard

C++ FEM library | user-friendly | multi-physics | hp-adaptive | HPC
http://www.sparselizard.org
Other
334 stars 62 forks source link

Thermal radiation #14

Closed hbadi closed 4 years ago

hbadi commented 4 years ago

Hello Alex,

I want to try Sparcelizard for simulating coupled conduction-radiation 2D problems. Does Sparselizard handle nonlinear Stefan-Boltzmann thermal radiation model (gray model body) ?

Best regards.

halbux commented 4 years ago

You mean the heat radiation proportional to epsilonsigmaT^4? It definitely does.

Simply write your thermal conduction equation (see e.g. "CONDUCTOR HEATING DUE TO DC CURRENT" online example) and add your radiation term so something like this (check the signs etc.):

field T("h1); ... formulation heatequation; heatequation += integral(volume, grad(tf(T))kgrad(dof(T))); // Add the radiation heat going out: heatequation += integral(volumeskin, -epsilonsigmapow(T,4) );

And solve with a nonlinear loop as done in the conductor heating example.

Alex