Heat conduction in rectangular domain with periodic heat transfer coefficient

1 reply · 2K views · last reply
bancux
Messages
11
Reaction score
0
Dear all, I have difficulty in solving this problem (see the figure in attached thumbnails)

I have a rectangular shape with length/height of L and the thickness/width of [tex]\delta[/tex]
Within the rectangular area, a heat conduction occurred. I would like to determine the temperature profile within the rectangular.

Based on the right hand side figure, I take a small cell with size of dx times dy. Applying the energy balance within the cell.

I have attached all the equation that I can think of, and also the boundary conditions. I have difficulty in solving the differential equation.

Thanks in advance.

note:

[tex]\lambda[/tex] is thermal conductivity
[tex]\alpha[/tex] is heat transfer coefficient, it is suppose to be dependent on x-location and time. I am going to simulate the heat transfer coefficient as a periodic function.
 
Attachments
  • Left: a rectangular domain with x-axis pointing down and y-axis pointing right, width labelled delta, height labelled L, with a heat transfer coefficient alpha(x,t) on the right boundary. Right: a control volume box with heat flows q''x in from top, q''y in from left, q''x+dx and q''y+dy out, and internal generation phi'''.
    schema1.jpg
    12.5 KB · Views: 400
  • Handwritten-style derivation of 2D heat conduction energy balance: Fourier flux equations in x and y, second-derivative fluxes at x+dx and y+dy, and the final balance lambda d2T/dx2 dy + lambda d2T/dy2 dx + phi dx dy = rho cp dT/dt dx dy.
    equation1.jpg
    22.1 KB · Views: 404
  • Boundary condition equations: x = 0 and x = L give -λ ∂T/∂x = 0; y = 0 gives -λ ∂T/∂y = 0; y = δ gives -λ ∂T/∂y = α(x, t)(T_w − T_fl).
    BC.jpg
    14.7 KB · Views: 456