[SOLVED] Diffusion for tetrahedral mesh

[jjcasmar]

Hello,

is there any component to solve the heat diffusion equation over a tetrahedral mesh? If not, is there anyway to add the temperature as a value of the mechanical state and solve it as other state values (position, velocities) or should I define a container for the temperature and solve it using my own solver?

I develop a diffusion solver for tetrahedral meshes using FEM discretization, but I have lost the code and I want to check first if SOFA can do this.

Thanks!

[Hugo]

Hi Juan Jo,

That’s a good point: today the heat equation isn’t available in SOFA. However, I worked during 3 years (PhD) around this topic. I intend to release a plugin for heat diffusion in tetrahedra. How urgent is that for you (this would need few days maybe some weeks max to officially release it with Inria)? I also implemented different solvers dedicated to heat transfer, would this be of interest for you too?

Best regards,

Hugo

[jjcasmar]

Hi Hugo,

sounds really interesting. Please, let me know when it is released!

[Hugo]

Dear Juan Jo,

The diffusion forcefield has been added in SofaSimpleFEM module in the master branch. You can find an example for this in: examples/Components/forcefield/TetrahedronDiffusionFEMForceField.scn. Enjoy!

Cheers,

Hugo

[Author]

Posts