I am writing a code in C++ to find out the temperature profile in a 2d conduction problem.

suppose if i want to vary thermal conductivity which is the function of Coordinates, like k = x + y.

So where this function wll be added in the discritizes equation.

if k varies with temperature then i thnk ,it should as ap = an+aw+ae+as+(k);

ap = an+aw+ae+as+(aT+b).

is this right. please let me know

i am solving this problem using Finite volume method, Guass Seidal method.

In your PDE, the second order term should be d/dx(K*dT/dx)

ya i know that but if my k varies with temperature then can u tll me the algorithm to write the code.

take care of the flux continuity at the interface (harmonic mean). see patankar for details - very well described.

ok thanks