# Bioheat Numberical problem

Biomedical Heat Transfer

Numerical problem

Due 04/07/16

Budget :20 USD only

Obtain a numerical solution by the method of finite volumes to the heating due to electric current circulation in the block shown in figure below. Solve the thermal problem for cartesian coordinates. The governing equation of the electrical part of this problem is, according to the theory developed in class (equation 2.21):

∇ · ∇(σV ) = 0

or, for cartesian coordinates and a constant electrical conductivity σ

The section to be analyzed has dimensions L =0.1 m by L =0.1 m, and has electrical conductivity σ = 0.188 S/m. The coordinate system is fixed in the upper right corner, as shown in figure. The shaded boundaries are electrically and thermally insulated, while the arrowed boundary represents a uniform current being injected into the section while the temperature there is kept at chilled water temperature. The right and bottom boundaries have a voltage equal to zero and temperatures equal to the basal temperature. The thermal properties of the tissue are k = 0.512 W/m·K, ρ = ρa = 1000kg/m3, C = Ca = 3600 J/kg·K and ω = 0.0005 s−1, and ˙qm=700 W/m3.

The boundary conditions for this electrical problem can be stated as:

V (0,y) = 0 for right face

V (x,L) = 0 for bottom boundary

= 0 for upper insulated boundary

For the left hand side boundary insulated sections

= 0 if y ≤ L/3 or y ≥ 2L/3

For left hand side boundary current injecting section

Amperes if L/3 ≤ y ≤ 2L/3

Use the solution for the electric problem obtained by the method of separation of variables to specify the external heat generation term of Pennes equation, and solve numerically the steady state thermal problem using the method of finite volumes. For the thermal problem the governing equation is Pennes bioheat equation for 2-D steady state cartesian coordinate problems

where the term ˙ql can be taken from the solution to HW 5, by taking

use the following boundary conditions for the thermal problem:

T(0,y) = Ta + q˙m/(ρaωCa) for right face

T(x,L) = Ta + q˙m/(ρaωCa) for bottom boundary

= 0 for upper insulated boundary

For the left hand side boundary insulated sections

= 0 if y ≤ L/3 or y ≥ 2L/3

For left hand side boundary current injecting section

T(0,y) = 5◦C if L/3 ≤ y ≤ 2L/3

Change the current density J leaving from the central part of the left hand side until you find a temperature of 100◦C at some point inside the tissue. Plot your results using matlab contour plots.

For that case, determine the percent of region of tissue where the temperature is above 60◦C.

Figure 1: Geometry of problem is in the attached file

also I'll attach related documents

( 12 comentários ) Erie, United States

ID do Projeto: #10127923

## 4 freelancers estão ofertando em média \$38 para esse trabalho

richginkgo

Wonderful Project! We are pro in Matlab... As we are masters of mathematics and control engineering, we will give a successful result in time. Plz keep in touch with me! Thanks..

\$40 USD em 1 dia
(18 Comentários)
3.8
VirtualBrainInc

I have written more than 3000 academic papers in the last 7 years. Please give me this excellent opportunity to deliver a high quality content on time. I value client satisfaction and that is the main reason as to why Mais

\$60 USD em 1 dia
(1 Comentário)
2.4

DEAR HIRING MANAGER I am an industrial engineer. I am very good at mechanical, electrical, control and automation, fluid mechanical,...I am expert at system modeling,data base creating,data quantitative analysis,sta Mais

\$30 USD em 1 dia
(2 Comentários)
0.8
umecka

Hi, although I am an article writer., I do mathematical analysis and give solutions to problems. I assure to deliver on time. Thanks. Umecka

\$25 USD em 1 dia
(0 Comentários)
0.0
midhunMechE

I am currently working in industry in solving PDEs . This is a steady state problem and is easy to solve. We may do a skype chat if you wish. Cheers Midhun

\$25 USD in 40 dias
(0 Comentários)
0.0