Ecuación de calor transitorio en un tubo (coordenadas polares) con temperatura inicial uniforme y temperatura límite uniforme

Traté de resolver este problema:

tengo un tubo de radio interior y exterior R 1 y R 2 (como tubo de cobre para frigoríficos). La temperatura en el exterior del tubo ( ρ = R 2 ) se mantiene a T mi . denotamos T 1 la temperatura dentro del tubo y T 2 la temperatura entre R 1 y R 2 . Los coeficientes de difusividad del material son k 1 (como el fluido dentro del tubo) y k 2 (coeficiente de cobre).Tubo de cobre, fluido en el interior, temperatura mantenida a <span class= T mi fuera del tubo" />

Entonces las ecuaciones de calor son:

[ k i 1 ρ ρ ( ρ ρ ) t ] T i = 0

entonces las condiciones de contorno son:

T 1 ( R 1 , t ) = T 2 ( R 1 , t )
yo 1 ρ T 1 ( ρ = R 1 , t ) = yo 2 ρ T 2 ( ρ = R 1 , t )
T 2 ( R 2 , t ) = T mi

con yo i la conductividad térmica, y la condición inicial es:

T 1 ( ρ , 0 ) = T 2 ( ρ , 0 ) = T 0

Por separación de variables encontramos soluciones generales:

T 1 = A + norte = 1 a norte mi λ norte t j 0 ( λ norte / k 1 ρ )
T 2 = B + norte = 1 b norte mi m norte t j 0 ( m norte / k 2 ρ ) + C norte mi v norte t Y 0 ( v norte / k 2 ρ )

como t + se alcanza el estado estacionario por lo que lógicamente tenemos: A = B = T mi pero tengo dificultades para encontrar los valores propios λ norte , m norte y v norte y los coeficientes de la serie. ¿Alguien sabe la solución de este problema?

Lo que yo sé

Para un caso más simple, un disco de radio R , con temperatura inicial dentro del disco T 0 y temperatura límite T ( R , t ) = T mi encontramos la solución:

T ( ρ , t ) = T mi + norte = 1 + a norte j 0 ( λ norte / k 1 ρ ) mi λ norte t

con el λ norte resolver (vinculado a los ceros de j 0 función de Bessel):

j 0 ( λ norte / k 1 R ) = 0

y coeficientes a norte son (usando producto escalar ponderado):

a norte = ( T 0 T mi ) 0 R 1 ρ j 0 ( λ norte / k 1 ρ ) d ρ 0 R 1 ρ j 0 2 ( λ norte / k 1 ρ ) d ρ

Tengo una idea, usando el 3 BC y asumiendo v norte = λ norte = m norte para deshacerme de los términos de tiempo exponencial, encontré una ecuación determinante en λ norte que da el λ norte (que implica, por supuesto, cálculos numéricos). Para tener la temperatura adentro, T 1 (entonces los coeficientes a norte ), puedo integrar en el radio R norte , que corresponden a j 0 ( R norte ( λ norte / k 1 ) ) = 0 y me deshago de la suma infinita debido a la ortogonalidad. ¿Qué piensan ustedes, amigos?

Respuestas (1)

Finalmente lo encontré usando este papel (complicado):

http://verl.npre.illinois.edu/Documentos/J-08-01.pdf

Como dije en el comentario, necesitamos aplicar BC y encontramos una matriz igual a cero con coeficientes, lo que da una ecuación determinante cero, capaz de encontrar el λ norte valores propios.