new_at_this 
April 6, 2012 11:09 
Implement constant heat flux boundary condition
I have written a matlab code that will solve for the developing temperature field given a velocity field. At the moment I have already implemented the constant wall temperature boundary condition and the results look correct.
Now I am trying to implement the constant heat flux boundary condition. I have the equation below
where T_1 is the first interior point and T_w is the wall temperature.
I convert to non dimensional temp
and calculate the appropriate so that it satisfies the flux and implement it as a dirichlet condition.
When I calculate Nu, defined as
it doesn't ever reach steady state. When it is supposed to reach steady state, the wall temperature increases more slowly than the bulk temperature and I get a decreasing Nu. Anyone seen this type of behavior before or know whats going on? Am I implementing the constant flux boundary condition properly?
