Scribd
Upload
4 views

CLASE 18

Uploaded by

Cruz Paniagua Rene Angel
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
4 views

CLASE 18

Uploaded by

Cruz Paniagua Rene Angel
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 15

CLASE 17: RESOLUCIÓN DE EDO PVF METODO

DEL DISPARO NO LINEAL


AUXILIAR: DIEGO FLORES URIBE
Consider a tubular reactor of length L = 0.1 (m) and radius R = 0.005 (m). The
volumetric heat source generated due to the exothermic reaction is Q = 1.4108
W/m3. The thermal conductivity of fuel is kf = 1.9W/(m°C), and the outer surface
temperature of the reactor is Ts = 650 °C.
Account the heat transfer only in the radial direction of the tubular reactor and
steady state.
𝑇 𝑟 = 𝑅 = 𝑇𝑠 𝑑𝑇
𝑞𝑟 = −𝑘
𝑑𝑟
𝑟 𝑑𝑇
Condiciones 𝑇 𝑟 = 0 =? ; 𝑞𝑟 = 0 0 = −𝑘
𝑧 𝑑𝑟
Frontera.
𝑑𝑇
𝑟=0 =0
𝑑𝑟

1. No conozco la 2. La
𝑇 = (𝑇(𝑟) , 𝑇(𝜃) , 𝑇(𝑧) ) velocidad del transferencia de
𝑣Ԧ = (𝑣𝑟 , 𝑣𝜃 , 𝑣𝑧 ) fluido por lo que calor solo es en
despreciara su sentido radial.
efecto.
𝑇 ≠ 𝑓(𝜃, 𝑧)
𝐸. 𝐸.

+ 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛

𝑣𝑟 = 0, 𝑣𝜃 = 0, 𝑣𝑧 = 0
+ 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛
1𝜕 𝜕𝑇
0=𝑘 𝑟 + 𝑔𝑒𝑛𝑒𝑟𝑎𝑐𝑖ó𝑛
𝑟 𝜕𝑟 𝜕𝑟
1𝑑 𝑑𝑇
0=𝑘 𝑟 +𝑄
𝑟 𝑑𝑟 𝑑𝑟
1 𝑑 𝑑𝑇 𝑑 𝑑𝑇
0=𝑘 𝑟 +𝑟 +𝑄
𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟

1 𝑑 𝑑𝑇 𝑑 𝑑𝑇 𝑄
0= 𝑟 +𝑟 +
𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟 𝑑𝑟 𝑘
1 𝑑𝑇 𝑑2 𝑇 𝑄
0= 1 +𝑟 2 +
𝑟 𝑑𝑟 𝑑𝑟 𝑘
𝑑2 𝑇 1 𝑑𝑇 𝑄
0= 2+ +
𝑑𝑟 𝑟 𝑑𝑟 𝑘
𝑑2 𝑇 1 𝑑𝑇 𝑄
0= 2+ +
𝑑𝑟 𝑟 𝑑𝑟 𝑘
𝑑𝑇
𝑟=0 =0
𝐵. 𝐶. ቐ 𝑑𝑟
𝑇 𝑟 = 𝑅 = 𝑇𝑠

𝑑𝑇
= 𝑧 … (1)
𝑑𝑟
𝑑𝑧 1 𝑄
=− 𝑧− … (2) MATLAB
𝑑𝑟 𝑟 𝑘
𝑧 𝑟=0 =0
𝐵. 𝐶. ቊ
𝑇 𝑟 = 𝑅 = 𝑇𝑠
Solve the following boundary value problem using the shooting method, in this
problem, we will look at the convergence of the solution to a boundary value
problem as a function of the grid spacing.
Where the first order spatially dependent reaction term is
𝒌(𝒙) = 𝟏𝟎[𝟏 + 𝐬𝐢𝐧(𝝅𝒙)]
The boundary conditions for the problem are:
𝑪(−𝟏) = 1 𝒚 𝑪(𝟏) = 0. 5
Mass diffusivity: D = 1.
𝑥𝑐 𝑥 = 1 = 0.5

Condiciones
Frontera.

𝑐 𝑥 = −1 = 1

1. No existe 2. La
𝑐 = (𝑐(𝑥) , 𝑐(𝑦) , 𝑐(𝑧) ) movimiento. concentración
𝑣Ԧ = (𝑣𝑥 , 𝑣𝑦 , 𝑣𝑧 ) es función de x.
𝐸. 𝐸. 𝑐 ≠ 𝑓(𝑦, 𝑧)
𝑣𝑟 = 0, 𝑣𝜃 = 0, 𝑣𝑧 = 0
𝜕2𝐶
0=𝐷 2 +𝑅
𝜕𝑥
𝜕2𝐶
0 = 𝐷 2 − 𝑘𝐶
𝜕𝑥
𝑑2 𝐶 10 1 + sin 𝜋𝑥
2
= 𝐶
𝑑𝑥 1
𝑑2 𝐶 𝐶 𝑥 = −1 = 1
2
= 10 1 + sin 𝜋𝑥 𝐶 𝐵. 𝐶. ቊ
𝑑𝑥 𝐶 𝑥 = 1 = 0.5
𝑑2 𝐶
2
= 10 1 + sin 𝜋𝑥 𝐶
𝑑𝑥
𝐶 𝑥 = −1 = 1
𝐵. 𝐶. ቊ
𝐶 𝑥 = 1 = 0.5

𝑑𝐶
= 𝑧 … (1)
𝑑𝑥
𝑑𝑧 MATLAB
= 10 1 + sin 𝜋𝑥 𝐶 … (2)
𝑑𝑥
𝐶 𝑥 = −1 = 1
𝐵. 𝐶. ቊ
𝑇 𝑥 = 1 = 0.5
Gracias
…

You might also like