# Stable implicit method to solve convection-heat diffusion in 3D

The last couple of hours I have been looking for an unconditionally stable method to solve the convection-diffusion equation within a 3D inhomogeneous material.

Here's the well known diffusion-convection equation: $$\rho(\vec{r})C_p(\vec{r})\frac{\partial T(\vec{r},t)}{\partial t}=\nabla [\kappa(\vec{r})\nabla T(\vec{r},t)]-\rho(\vec{r})C_p(\vec{r})\nabla[\vec{u}T(\vec{r},t)]+g(\vec{r},t)$$ where $\rho$ is the material density, $C_p$ is the heat capacity, $\kappa$ is the thermal conductivity, $\vec{u}$ is the velocity vector and $g$ represents the heat source.

So far I haven't found any method which can be used for that. The articles I came across all made severe simplifications like $\kappa$=const., ignored the convection term (2nd term on the R.H.S.) or reduced the problem to 1 or 2 dimensions. So I'm wondering now if this can be solved at all.

• You would get more answers if you stated the equations you want to solve as a formula. I think it would also useful if you stated what you are looking for. Why is the method described there not sufficient, for example? – Wolfgang Bangerth Sep 26 '16 at 19:10
• Thanks for your comment. I guess this is just not the right place for my question. Is there a way to delete it? – OD IUM Sep 26 '16 at 19:13
• @ODIUM, I disagree and this is likely the proper place for your question. However, it's rather specialized in its current form and would gather more attention if you present it more explicitly (i.e. relevant questions, boundary conditions, other numerical schemes, etc.) – Spencer Bryngelson Sep 27 '16 at 5:29
• Centered differences in space and backward Euler in time should work just fine. – David Ketcheson Sep 27 '16 at 15:51

The solving procedure should go with the deferred correction method, which allows us to correct solution step by step inside an iteration. Upwind scheme blending with the central difference can be applied to the convection term; it also has second order accuracy. The diffusion term with a non-constant $\kappa (\vec{r})$ can be handled by Hrvoje Jasak PhD Thesis, p.83. For a generic FVM, we can find more useful things within this thesis.