# Discussion Forum

## Boundary conditions in heat transfer

RSS feed   |   Turn on email notifications   |   3 Replies   Last post: January 4, 2014 2:55pm UTC

Yael Ruschkewitz

December 21, 2009 7:36pm UTC

Boundary conditions in heat transfer

Hi,
I want to model a muscle tissue containing a capillary blood vessel. Now I’m using a simple model for start… My model geometry is in the attached file.
As I want to model heat convection between the tissue and blood I want to use on interior boundaries between the blood vessel the “heat flux discontinuity” boundary condition.

My main question is regarding Tinf: according to comsol help files: “Tinf is the external bulk temperature” and this is from the section of interior boundaries (heat flux discontinuity- interior boundaries only)… I don’t understand which subdomain is considered as having external bulk temperature? The blood vessel or the tissue surrounding it? This is a critical question as the temperature of the blood and the tissue are not equal in my model… I’m modeling the tissue as having 37C and the blood as 35C… Does Tinf belong to the subdomain of the higher temperature? in another case of the model I’m representing both blood and tissue as having 37C, but when the tissue has a higher heat production rate from metabolism…
Also, in my model I am representing heat production from metabolism using the subdomain settings of the bioheat modulo, so I thought of setting q=0… I’m not sure whether I best set q=0 or maybe I should set it’s value as the same value I set the metabolic heat of blood? What do you think?

Ivar Kjelberg

December 27, 2009 10:35am UTC in response to Yael Ruschkewitz

Re: Boundary conditions in heat transfer

Hi Yeal

You might have found out by now, I hope so, but in anycase, in the boundary condition Heat flux you have 2 variables Q and h, you might either choose one or the other, or both. Now to understand how COMSOL isusing these and how they are related (Tinf is only linked to "h") take a look at the equations:

Once you have selected 2D "htgh" BC's "Heat transfer" and the boundary of interest go to :
Physics - Equation Systems - Boundary Settings

and take alook at i.e. 4th row (for "T" as the 4th dependent variable u,v,p,T,J), you should see something like:

-h_htgh*T+q0_htgh+h_htgh*Tinf_htgh

that once rearranged this means:

h_htgh*(Tinf_htgh-T) + q0_htgh

Hence the Q variable (or q0_htgh internally) enters directly thermal power onto the boundary, while the "h" is proportional to the temperature difference (in K) of (Tinf-T) and T is the temperature you are solving for, so Tinf yould be your constant blood temperature, then h should be scaled to take into account the liquid solid exchange, but perhaps next step to take into account the blood flow, or total volume passing through, as the blood temperature might not manage to take/give all the heat, unaffected, if you have a sever feber, but I would propose not to use non-linear terms to start with.

you should have a look, it's excellent for your type of activities (I just got it a week ago)

Good luck
Ivar

Abhishek Kumar Singh

September 12, 2013 3:48pm UTC in response to Ivar Kjelberg

Re: Boundary conditions in heat transfer

In the Joule heating module (2 D), I want to know what are the boundary conditions required for the temperature rise in nanowire, when 1 mA current is passed through the wire having cross-sectional area150 x 150 nm^2 and 5 micron length. I have tried all the possibilities but still I am not able to get the solution. please help in this regard.

Waleed Salman

January 4, 2014 2:55pm UTC in response to Ivar Kjelberg

Re: Boundary conditions in heat transfer

Hi all
How to find the boundary equation setting in the new versions 4.4
I knew that in v3.5 it was physics>equation system>boundary setting

Rules and guidelines