I get the idea but from the numerical standpoint this is an even worse idea than the way we do fluid-fluid interfaces. Additive decomposition of elliptic equations is not effective at all. Are you looking to model contact between multiple solid materials?

Yes, I'm looking to model the transfer of heat between multiple solid zones with different thermal properties. Eventually I'd like to model contact resistance of such interfaces as well, where the user could specify the contact resistance value of each solid-to-solid interface. This is something that me and Praharsh will be working on.

Okay, then allow me to save you guys some time, implement something similar to what we use for the FEA solver.
See the use of element_properties in CFEASolver.cpp.

Okay, then allow me to save you guys some time, implement something similar to what we use for the FEA solver. See the use of element_properties in CFEASolver.cpp.

Thanks for the suggestion! We'll have a look.

Good job @EvertBunschoten. If you want we can schedule a call together and plan the development a bit.

Was able to run a case with fluid and two solids with contact resistance between the solid domains. By setting a non-zero contact resistance between the solid domains there appears a temperature discontinuity between the solids, but by setting the contact resistance value to zero on the fluid-to-solid interfaces, the temperature trends are continuous.
Further testing still required on more complicated cases, but the core functionality seems to work so far.

I get the idea but from the numerical standpoint this is an even worse idea than the way we do fluid-fluid interfaces. Additive decomposition of elliptic equations is not effective at all. Are you looking to model contact between multiple solid materials?

Hi Pedro. What do you mean by additive decomposition? And what do you mean by "effective" in this context? In terms of stability, the model converges rapidly, even when initializing one zone at a temperature far above that of the neighboring zone.

https://en.wikipedia.org/wiki/Additive_Schwarz_method
How large are the meshes in your examples? I expect convergence will slow down as you go to bigger problems.

You had I comment but I don't see it here.
In any case, it's the total time to solution (same residual drop) that I expect will get worse. At constant CFL.

You had I comment but I don't see it here. In any case, it's the total time to solution (same residual drop) that I expect will get worse. At constant CFL.

Sorry for replying so late, I had to put this project on the back burner for a bit. I made a mistake in my previous analysis, which is why I removed the message. I now ran five simulations with only solid to solid CHT heat transfer with contact resistance. With each simulation I quadrupled the number of mesh elements and ran with a constant CFL.

The number of iterations required until convergence seems to scale linearly with the number of elements. However, this was not the case for total simulation time, which scales quadratically.

@EvertBunschoten even if the performance isn't great this is a pretty small change so I think we can merge it.
Do you have time to wrap it up and add a simple regression test?

@EvertBunschoten even if the performance isn't great this is a pretty small change so I think we can merge it. Do you have time to wrap it up and add a simple regression test?

I just added a test case and tutorial under the TestCases and Tutorials repo under the same branch name. I'm also writing a short tutorial on the SU2 website repo.