0% encontró este documento útil (0 votos)
135 vistas10 páginas

Métodos Numéricos en Problemas de Ingeniería

El primer problema trata sobre el drenaje de un estanque a través de un tubo. Se presenta la ecuación diferencial que describe el cambio en la profundidad con el tiempo y se pide resolverla numéricamente usando Runge-Kutta de cuarto orden para estimar el tiempo de vaciado. El segundo problema modela la deflexión de una placa sometida a fuerzas, requiriendo resolver la ecuación diferencial correspondiente usando Runge-Kutta con el método de disparo.

Cargado por

Xim Berrones
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como XLSX, PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
135 vistas10 páginas

Métodos Numéricos en Problemas de Ingeniería

El primer problema trata sobre el drenaje de un estanque a través de un tubo. Se presenta la ecuación diferencial que describe el cambio en la profundidad con el tiempo y se pide resolverla numéricamente usando Runge-Kutta de cuarto orden para estimar el tiempo de vaciado. El segundo problema modela la deflexión de una placa sometida a fuerzas, requiriendo resolver la ecuación diferencial correspondiente usando Runge-Kutta con el método de disparo.

Cargado por

Xim Berrones
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como XLSX, PDF, TXT o lee en línea desde Scribd

Problema 1 (35%). La potencia de salida de una celda solar varía con el voltaje aplicado.

El voltaje Vmp a la cual lapotencia de salida es máxima está dada por la ecuación:

donde el Voc es el voltaje de circuito abierto, T la temperatura en K, q la carga de un electrón igual a 1.6022E-19
y kB la constante de Boltzmann igual a 1.3806E-23 J/K.
Para un valor Voc = 0.5 V determine el voltaje Vmp para potencia de salida máxima de la celda, calcule el valor pa
a) Realice una gráfica de la función
b) Resuelva el problema utilizando el método de Newton-Raphson con una tolerancia de 0.005%

Problema 3 (30%). La población de bacterias de un cultivo se puede calcular como

𝑃_𝑡=𝑃_0 [1+0.054𝑡/1!+(0.054𝑡)^2/2!+…+(0.054t)^n/n!]

Donde Pt es la población de bacterias a un tiempo t, P0 es la población inicial y t es el tiempo en horas.


Escriba un programa en Scilab que calcule la población de bacterias de 1 a 5 horas (tome en cuenta solo números
Se puede calcular un error aproximado al ir agregando un término nuevo a la serie
Utilice los primeros 11 términos de la serie o el cálculo del término entre corchetes puede terminar si el error ap
Problema 2 (35%). Un recubrimiento en la superficie de un
La temperatura en el recubrimiento es determinada por pro
Si se supone que la radiación es difusa y gris, se obtiene el s

5.67×10^(−8) 𝑇_𝑐^4+17.41
𝐽_𝑐−0.71𝐽_ℎ+7.46𝑇_𝑐=2352.
un electrón igual a 1.6022E-19 C,
5.67×10^(−8) 𝑇_ℎ^4+1.865
𝐽_ℎ−0.71𝐽_𝑐+7.46𝑇_ℎ=1109
a de la celda, calcule el valor para una temperatura de 350 K.

ncia de 0.005%

es el tiempo en horas.
s (tome en cuenta solo números enteros) para P0 = 100. Figura 1: Panel recubierto expuesto a radiación

es puede terminar si el error aproximado es menor que 0.005 donde Jh y Jc son las radiosidades del calentador y de la sup
Resuelva el sistema para calcular Jh, Jc, Th, y Tc utilizando el
Considere una tolerancia de 0.005%
Nota: Se recomienda utilizar los siguientes valores iniciales:
Si lo desea, puede simplificar el sistema a uno de tres variab
brimiento en la superficie de un panel es curada por energía radiante de un calentador.
imiento es determinada por procesos de transferencia de calor de radiación y convección.
n es difusa y gris, se obtiene el siguiente sistema de ecuaciones no lineales:

0^(−8) 𝑇_𝑐^4+17.41𝑇_𝑐−𝐽_𝑐=5188.18
1𝐽_ℎ+7.46𝑇_𝑐=2352.71
0^(−8) 𝑇_ℎ^4+1.865𝑇_ℎ−𝐽_ℎ=2250
1𝐽_𝑐+7.46𝑇_ℎ=11093

recubierto expuesto a radiación

dades del calentador y de la superficie, respectivamente; Th y Tc son sus respectivas temperaturas.


cular Jh, Jc, Th, y Tc utilizando el método de Newton.

los siguientes valores iniciales: T h = Tc = 298 K, Jh = 5000 W/m^2 y Jc = 3000 W/m^2.


r el sistema a uno de tres variables con tres ecuaciones
Problema 1 (35%). Una forma de evaluar la función de confiabilidad R(t) es a través de la distribución Weibull, la
R(t)=e^(-1/α)^β
donde α y β son los parámetros del modelo y t la variable independiente.
a) Realice una transformación lineal de R(t) para ajustar por mínimos cuadrados los siguientes datos tabulados.
b) Calcule los párametros del modelo
c) ¿Es bueno el ajuste del modelo? Explique con base los valores del factor de determinación y el error relativo p

𝑅(𝑡)=𝑒^((−1/
t,h R(t)
37.5 0.90
46.0 0.74
48.0
51.5
0.58
0.42 𝛼)^𝛽 )
53.0 0.26
54.5 0.10

Problema 2 (35%). La siguiente Tabla muestra la concentración de salida de un reactor en distintos momentos du

El caudal de salida en m^3 s^-1 se puede calcular con la siguiente ecuación:


Q(t)=20+10sen[2π/24(t-10)]
Utilice la integración numérica para calcular la concentración promedio durante el periodo de 16 horas con la sig
ċ=(∫{0,t} Q(t)c(t)dt)/∫{0,t} Q(t)dt
Nota: Utilice la mayor cantidad de daos de la tabla y la (s) regla (s) de integración adecuadas para la estimación y

Nota del pizzaron: Utilizar t en horas (como está en la tabla) el resultado está en m^3/s

t, h c, mg/L
0 1
1 1.5
5.5 2.3
10 2.1
12 4
14 5
16 5.5
és de la distribución Weibull, la cual se expresa como sigue: Problema 3 (30%): Si se conoce la

os siguientes datos tabulados.

erminación y el error relativo promedio


Escribe un programa en Scilab que c
Se recomienda el us del comando "fo
La función a integrar se puede defini
o se pueden definir los valores de la

actor en distintos momentos durante un periodo de 16 horas.

el periodo de 16 horas con la siguiente ecuación:

adecuadas para la estimación y lograr la mejor precisión. Verifique la consistencia de unidades.


roblema 3 (30%): Si se conoce la distribución de la velocidad de un fluido a través de un tubo circular, el flujo volumétrico (cm^3 s^-1) se calcula como:

donde r(cm) es la distancia radial, R es el radio total (3cm en este c

scribe un programa en Scilab que calcule Q utilizando la regla trapezoidal de segmentos múltiples con 10 segmentos.
e recomienda el us del comando "for" para evualuar la sumatoria de la ecuación.
a función a integrar se puede definir con el comando "function" y utilizar un ciclo "for" para pedirle a Scilab qu calcule la función en cada punto necesario
se pueden definir los valores de la función simplemente con un vector.
étrico (cm^3 s^-1) se calcula como:

al, R es el radio total (3cm en este caso).

e la función en cada punto necesario,


Problema 1 (50%) Un estanque se drena a través de un tubo como se observa en la siguiente figura:
Con algunas suposiciones, la ecuación que describe el cambio de la profundidad on el tiempo es
dh/dt=-πd^2/4A raiz(2g(h+e))
donde h es la profundidad (m), g es la constante gravitacional (9.81 m/s^2), e es la profundidad de l
A es el área de la superficie del estanque (m^2), l cual depende de h, como una aproximación se pu
A=0.114h^2+1.285h
Tomando en cuenta los siguientes valores: h(0)=6m, d=0.25m, e=1m
a) Resuelva la ecuacion diferencial para dar un estimado del tiempo en que se tarda en vaciar el est
Para la solución utilice el método de Runge-Kutta de cuarto orden con un tamaño de paso de 2 s
b) Muestre sus resultados en forma de tabla y gráfica

𝐴=0.114ℎ^

Problema 2 (50%). La siguiente ecuación diferencial permite modelar la deflexión de una placa rectangular de ac
d^2y/dx^2 -S/D*y = -qL/2D*x + q/2D*x^2
Sujeta a las condiciones de frontera: y(0) = 0; y(L) = 0
a) Utilice el método Runge-Kutta con el método del disparo para calcular la deflexión de la placa con
Considere un valor de h igual a 10 cm.
b) Construya una tabla y una gráfica con los resultados obtenidos.
la siguiente figura:
profundidad on el tiempo es

m/s^2), e es la profundidad de la salida del tubo por debajo del fondo del estanque (m).
, como una aproximación se puede establecer la función de A respecto a h como:

en que se tarda en vaciar el estanque.


on un tamaño de paso de 2 s

𝐴=0.114ℎ^2+1.285ℎ

de una placa rectangular de acero larga, con carga uniformemente distribuida sujeta a una fuerza de tension axial.

cular la deflexión de la placa con los siguiente valores de los parámetros: q=15kgcm^-2, S=18kg/cm, D=1.02E8kg/cm y L=130cm
2E8kg/cm y L=130cm

También podría gustarte