|
[Sponsors] |
![]() |
![]() |
#1 |
New Member
|
Hi there,
I am using TEACH code for my project. To simulate a conducting wall inside a close cavity I faced a strange problem. What I did: 1- I wrote the code for natural convection in a cavity 2- Defined a vertical wall in the cavity and assigned zero velocity for those cells 3- For energy equation I changed diffusion coefficient from (1/(Ra.Pr)^0.5) to (Kr/(Ra.Pr)^0.5) which Kr is conductivity coefficient of wall over fluid. (dimensionless form of governing equation are shown in the picture.) I expected code to solve fluid part and when reaches to solid section conductivity changes and only conduction should be in the wall. However, interestingly the coefficient I add does not change thermal conductivity of wall and results show that conductivitis of wall and fluid are same. Any body knows how to fix it? Patankar explains it in his book (page 149) but I do not understand it. Any help appreciated. ![]() |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
I do not understand your question, you solve for the fluid and the solid zone two different temperature equations, each one with different diffusion coefficient (if Kr is not 1), isn't that? Now what you are checking in this two zones? could you plot the fields?
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
New Member
|
I am studying partitioned cavity as showed in the following figure:
![]() As you can see in the above figure, although I changed thermal conductivity of the partition in the middle of the cavity, it shows same thermal conductivity as fluid. If thermal conductivity was different, we should see breaking point for isotherms at the interface of partition-fluid. Now I am not should what I did wrong. Should I change more or just changing diffusion coefficient and making velocities zero in the partition is enough? Since usually most of people are familiar with TEACH code, I was hopping people help me with the items I should change to be able to have functioning conducting partition in the cavity. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
At the wall you prescribe the continuity of the heat flux, that is q_wall = K grad T, right? Now, what about the value of Kr? if it is of O(1) the difference is disregardable.
|
|
![]() |
![]() |
![]() |
![]() |
#5 |
New Member
|
Oh, thanks god you are here...
We should apply the continuity of flux at the interface but Patankar says instead we can increase viscosity of solid section to make velocity zero. For the energy equation, he says we can set solid property in the energy equation and solve whole domain without worrying about continuity. Now I did that but I am not getting big difference for partition section as I assumed Kr=1000 !! Should I also apply that constant heat flux for interface ? or other thing ? How ? ![]() |
|
![]() |
![]() |
![]() |
![]() |
#6 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
Quote:
yes, you must apply the continuity of the flux at the wall and, simultaneously, setting zero velocity for the solid zone |
||
![]() |
![]() |
![]() |
![]() |
#7 |
New Member
|
||
![]() |
![]() |
![]() |
![]() |
#8 |
New Member
|
I did as follow but my code stopped working and got error:
Tw=T(w-1) - (Kf/Kw)(dxf/dxw)(T(w-2)-T(w-1)) Should it be something like above or I am totally in wrong direction ? |
|
![]() |
![]() |
![]() |
![]() |
#9 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
what kind of error? just after one time step?
|
|
![]() |
![]() |
![]() |
![]() |
#10 |
New Member
|
It stopped working. Even not one step. I am guessing, I do not know correct way of making sure of flux continuity.
Could you please give me some details about how I can write between partition and fluid cells ? Much appreciated |
|
![]() |
![]() |
![]() |
![]() |
#11 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
q=-k dTdn, you discretize the RHS with a backwar/forward second order derivative and consider that q = h* Delta T|0
|
|
![]() |
![]() |
![]() |
![]() |
#12 | |
New Member
|
Quote:
(k dTdn )fluid = (k dTdn )solid so I have two cells from fluid side for left equation and two cells from solid side for right equation. The question is how can I write it for code to be able to solve it ? ![]() ![]() |
||
![]() |
![]() |
![]() |
![]() |
#13 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
Quote:
no, the flux in the solid is q = h *(T-Tr), not the Fourier flux |
||
![]() |
![]() |
![]() |
![]() |
#14 | |
New Member
|
Quote:
and I have checked all the papers in this regard and they all say this : ![]() as a boundary condition for conducting partition. Should I write this continuity flux for both of partition ? |
||
![]() |
![]() |
![]() |
![]() |
#15 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
the value of h depends on the material of the solid zone
|
|
![]() |
![]() |
![]() |
![]() |
#16 |
New Member
|
Do you mean something like the picture ?
![]() What is Tr in your note ? (Ts - ? ) Should I do this calculation for both sides ? For the left side I am going to calculate Ts (as in picture) what about right side ? T(E+1) ? |
|
![]() |
![]() |
![]() |
![]() |
#17 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,830
Rep Power: 73 ![]() ![]() ![]() |
yes, correct, the expression of q is well known, see for example here
https://en.wikipedia.org/wiki/Convective_heat_transfer you have the difference between fluid and wall temperature |
|
![]() |
![]() |
![]() |
Tags |
coding, conjugate heat transfer, partition |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Problem with conjugate heat transfer in Fluent | JianT88 | FLUENT | 0 | March 16, 2015 08:17 |
problem in conjugate heat transfer | reza67 | FLUENT | 2 | December 2, 2014 03:16 |
Transient conjugate heat transfer problem | troyker | FLUENT | 0 | June 21, 2013 03:37 |
Error finding variable "THERMX" | sunilpatil | CFX | 8 | April 26, 2013 07:00 |
Conjugate heat transfer problem with porous media | piko | Siemens | 1 | April 17, 2009 15:41 |