2Eva_2024PAOI_T3 EDP Parabólica

2da Evaluación 2024-2025 PAO I. 28/Agosto/2024

Tema 3. (30 puntos)

Para la siguiente Ecuación Diferencial Parcial con b = 2, resuelva usando las condiciones mostradas

\frac{\partial ^2 u}{\partial x^2} = b \frac{\partial u}{\partial t}
0 < x < 1  0 < t < 0.5
Condiciones iniciales:  u(x,0)=0
Condiciones de frontera: u(0,t)=1
u(1,t)= 2

Utilice diferencias finitas centradas y hacia adelante para las variables independientes x,t

a. Plantee las ecuaciones para usar un método numérico en un nodo i,j

b. Realice la gráfica de malla,

c. Desarrolle y obtenga el modelo discreto para u(xi,tj)

d. Realice al menos tres iteraciones en el eje tiempo.

e. Estime el error de u(xi,tj) y adjunte los archivos del algoritmo y resultados.

Rúbrica: Aproximación de las derivadas parciales (5 puntos), construcción de la malla (5), desarrollo de iteraciones (15), literal e (5 puntos)

Referencia: EDP Parabólicas. Chapra & Canale. 5ta Ed. Ejercicio 30.15. P.904

2Eva_2024PAOI_T2 Salto Bungee longitud total de cuerda

2da Evaluación 2024-2025 PAO I. 28/Agosto/2024

Tema 2. (40 puntos) Bungee Jumping 02
Para el salto del Bungee del ejercicio del tema anterior se toman lecturas con un sensor de velocidad sujetado a la persona.

2.1 De la tabla de datos obtenida, se observa que los tamaños de paso en tiempo no siempre son equidistantes.
Se requiere encontrar la distancia recorrida en el intervalo de [0,2.55] usando fórmulas de integración compuestas.

ti vi
0 0.0000
0.25 2.4479
0.5 4.8849
0.75 7.3001
1 9.6832
1.375 13.1763
1.75 16.5451
2.125 19.7641
2.4 22.0193
2.55 23.2075

2.2 Usando los datos de la tabla para el intervalo [2.55, 5.175] donde la velocidad de la caída de la persona al primer salto ha llegado a casi cero, o antes del primer rebote, se ha obtenido un polinomio de interpolación:

v = -3.979t2 + 21.557t – 5.3997

Obtenga la distancia recorrida en el segundo intervalo usando el método de Cuadratura de Gauss.

a. Realice el planteamiento de las ecuaciones para cada sección del ejercicio.

b. Describa el criterio usado para determinar el número de tramos usado en cada caso.

c. Desarrolle las expresiones completas del ejercicio para cada sección.

d. Encuentre la distancia total (profundidad) alcanzada por la persona al dar el salto.

Rúbrica: literal a 2.1 (5 puntos), a 2.2 (5 puntos) literal b (5 puntos), literal c 2.1 (10 puntos), c 2.2 (10 puntos), literal d (5 puntos)

Referencia:[1] Chapra. capítulo 28. Ejercicio 28.41 p852.

[2] Extreme Bungy Jumping with Cliff Jump Shenanigans! Play On in New Zealand! 4K! – devinsupertramp. 23 mar 2015.

2Eva_2024PAOI_T1 EDO Salto Bungee tiempo extiende cuerda

2da Evaluación 2024-2025 PAO I. 28/Agosto/2024

Tema 1. (30 puntos) Salto Bungee 01
La ecuación diferencial para la velocidad de alguien que practica el salto del bungee es diferente según si el saltador ha caído una distancia en la que la cuerda está extendida por completo y comienza a estirarse o encogerse.

Si la distancia recorrida es menor que la longitud de la cuerda, el saltador sólo está sujeto a las fuerzas gravitacional y de arrastre de la cuerda.

\frac{d^2y}{dt^2} = g - signo(v) \frac{C_d}{m}\Big( \frac{dy}{dt}\Big)^2 y \leq L

Salto Bungee 02Suponga que las condiciones iniciales son:

y(0) =0

\frac{dy(0)}{dt} = 0

Una vez que la cuerda comienza a estirarse, también deben incluirse las fuerzas del resorte y del amortiguamiento de la cuerda.

\frac{d^2y}{dt^2} = g - signo(v) \frac{C_d}{m}\Big( \frac{dy}{dt}\Big)^2 -\frac{k}{m}(y-L) - \frac{\gamma}{m}( \frac{dy}{dt}\Big) y \gt L
dy/dt m/s velocidad (v)
t s tiempo
g 9.8 m/s2 gravedad
cd 0.25 kg/m coeficiente de arrastre
m 68.1 Kg masa
L 30 m Longitud de la cuerda
k 40 N/m constante de resorte de la cuerda
γ 8 N s/m coeficiente de amortiguamiento de la cuerda
signo(v) función que devuelve –1, 0 y 1, para v negativa, cero y positiva, respectivamente

Encuentre el tiempo tc y la velocidad de la persona cuando se alcanza la longitud de la cuerda extendida y sin estirar (30 m), es decir y<L, aún se entra cayendo signo(v)=1. (solo primera ecuación)

a. Realice el planteamiento del ejercicio usando Runge-Kutta de 2do Orden

b. Desarrolle tres iteraciones para y(t) con tamaño de paso h=0.5

c. Usando el algoritmo, aproxime la solución para y en el intervalo entre [0,tc], adjunte sus resultados.txt

d. Indique el valor de tc, muestre cómo mejorar la precisión y realice sus observaciones sobre los resultados.

Observación: dy/dt = v, función signo(v) en Numpy:np.sign(v)

Rúbrica: literal a (5 puntos), literal b (15 puntos), literal c resultados.txt y grafica.png (5 puntos), literal d (5 puntos),

Referencia: [1] Chapra. capítulo 28. Ejercicio 28.41 p852.

[2] Extreme Bungy Jumping with Cliff Jump Shenanigans! Play On in New Zealand! 4K! – devinsupertramp. 23 mar 2015.

 

2Eva_2023PAOII_T3 EDP desarrolle expresión

2ra Evaluación 2023-2024 PAO II. 30/Enero/2024

Tema 3 (30 puntos) Para la siguiente Ecuación Diferencial Parcial con b = 2, resuelva usando las condiciones mostradas

\frac{\partial ^2 u}{\partial x^2} + b\frac{\partial u}{\partial x} = \frac{\partial u}{\partial dt}
0 < x < 1

0 < t < 0.5

Condiciones de frontera:
u(0,t)=0, u(1,t)= 1, 0≤t≤0.5
Condiciones iniciales:
u(x,0)=0, 0≤x≤1

Utilice diferencias finitas centradas y hacia adelante para las variables independientes x,t

a. Plantee las ecuaciones para usar un método numérico en un nodo i,j

b. Realice la gráfica de malla,

c. Desarrolle y obtenga el modelo discreto para u(xi,tj)

d. Realice al menos tres iteraciones en el eje tiempo.

e. Estime el error de u(xi,tj) y adjunte los archivos del algoritmo y resultados.

f. Con el algoritmo, estime la solución para b = 0 y b=-4. Realice las observaciones de resultados para cada caso.

Rúbrica: Aproximación de las derivadas parciales (5 puntos), construcción de la malla (5), desarrollo de iteraciones (10), literal e (10 puntos), literal f (5 puntos extra)

Referencia: EDP Parabólicas. Chapra & Canale. 5ta Ed. Ejercicio 30.15. P.904

2Eva_2023PAOII_T2 Cable cuelga entre apoyos A y B

2ra Evaluación 2023-2024 PAO II. 30/Enero/2024

Tema 2 (40 puntos) Un cable cuelga de dos apoyos en A y B. cable colgante entre apoyos

El cable sostiene una carga distribuida cuya magnitud varía con x según la ecuación

w = w_0 \Big[ 1+ \sin \Big(\frac{\pi x}{2l_B} \Big) \Big]

donde w0 = 1 000 lbs/ft y T0. = 0.588345×106.
La pendiente del cable (dy/dx) = 0 en x = 0, que es el punto más bajo del cable.

También es el punto donde la tensión del cable alcanza un mínimo de T0. La ecuación diferencial que gobierna el cable es

\frac{d^2y}{dx^2} = \frac{w_0}{T_0} \Big[ 1+ \sin \Big(\frac{\pi x}{2l_B} \Big) \Big]

a. Realice el planteamiento del ejercicio usando Runge-Kutta de 2do Orden

b. Desarrolle tres iteraciones para y(x) con tamaño de paso h=0.5

c. Usando el algoritmo, aproxime la solución para x en el intervalo entre [0,200], adjunte sus resultados.txt en la evaluación.

d. Realice sus observaciones sobre los resultados obtenidos sobre la altura y(200) alcanzada en el extremo derecho del cable y lo indicado en la gráfica del enunciado.

Rúbrica: literal a (5 puntos), literal b (15 puntos), literal c resultados.txt (10 puntos), grafica.png (5 puntos), literal d (5 puntos),

Referencia: Cable entre dos apoyos con carga distribuida. Chapra & Canale. 5ta Ed. Ejercicio 28.21. P849

2Eva_2023PAOII_T1 Volumen por solido de revolución

2da Evaluación 2023-2024 PAO II. 30/Enero/2024

Tema 1 (30 puntos) Los sólidos de revolución se generan al girar una región plana alrededor de un eje. solido de revolucion 1

V = \int_{a}^{b} \pi (f(x))^2 dx

El volumen generado al girar la región de la función f(x) en el intervalo [a,b], se puede calcular como el volumen del disco de radio f(x) y anchura dx.

f(x) = \sqrt{\sin (x/2)} g(x) = e^{x/3} - 1

Calcule el volumen de revolución generado por la región sombreada de la gráfica que ese encuentra entre: f(x) y g(x).
Las funciones se usan en el intervalo [0.1 , 1.8]:

Realice el planteamiento de las ecuaciones para el ejercicio, considerando que

a. Para el integral con f(x), use formulas de Simpson con al menos 3 tramos, mientras que

b. Para el integral con g(x) use Cuadratura de Gauss de dos puntos con al menos 2 tramos.

c. Desarrolle las expresiones completas del ejercicio para cada función.

d. Indique el resultado obtenido para el área requerida y la cota de error.

e. Determine el volumen de revolución generado por la región sombreada

Rúbrica: literal a (5 puntos), literal b (5 puntos), literal c (10 puntos), literal d (5 puntos), literal e (5 puntos)

Referencia:  [1] Volúmenes de sólidos de revolución. Moisés Villena Muñoz.Capítulo 4 p78. https://www.dspace.espol.edu.ec/bitstream/123456789/4800/4/7417.pdf
[2] 8.2.2 Gráficas en 3D en Python, sólidos de revolución. http://blog.espol.edu.ec/ccpg1001/graficas-en-3d-en-python-sistema-de-ecuaciones-y-planos/
[3] Volumes: Washer Method Animation 2. Stacey Roshan. 24 Abril 2016.

2Eva_2023PAOI_T3 EDP elíptica, placa rectangular con frontera variable

2da Evaluación 2023-2024 PAO I. 29/Agosto/2023

Tema 3 (35 puntos) Aproxime la solución de la Ecuación Diferencial Parcial

\frac{\partial ^2 u}{\partial x^2} + \frac{\partial ^2 u}{\partial y^2} = \Big( x^2 + y^2 \Big) e^{xy} 0 \lt x \lt 1 0 \lt y \lt 0.5

Con las condiciones de frontera:

u(0,y)=1, u(1,y)= y, 0≤y≤0.5
u(x,0)=1, u(x,0.5)=x/2, 0≤x≤1

Aproxime la solución con tamaños de paso Δx = 0.25, Δy = 0.25
Utilice diferencias finitas centradas para las variables independientes x,y

a. Plantee las ecuaciones para usar un método numérico en un nodo i,j

b. Realice la gráfica de malla,

c. desarrolle y obtenga el modelo discreto para u(xi,tj)

d. Realice al menos tres iteraciones en el eje tiempo.

e. Estime el error de u(xi,tj) y adjunte los archivos del algoritmo y resultados.

Rúbrica: Aproximación de las derivadas parciales (5 puntos), construcción de la malla (10), construcción del sistema lineal (15), resolución del sistema (5 puntos).

Referencia: 2Eva_IT2012_T3 EDP elíptica, placa rectangular

2Eva_2023PAOI_T2 Péndulo vertical amortiguado

2da Evaluación 2023-2024 PAO I. 29/Agosto/2023

Tema 2 (35 puntos) Una mejor aproximación a un péndulo oscilante con un ángulo θ más amplio y con un coeficiente de amortiguamiento μ se expresa con una ecuación diferencial ordinaria de segundo orden.

\frac{d^2 \theta}{dt^2} = -\mu \frac{d\theta}{ dt}-\frac{g}{L}\sin (\theta)

g = 9.81 m/s2
L = 2 m
θ(0) = π/4 rad
θ’ (0) = 0 rad/s

El péndulo se suelta desde el reposo, desde un ángulo de π/4 respecto al eje vertical. El coeficiente de amortiguamiento μ=0.5 es proporcional a la velocidad angular.

a. Realice el planteamiento del ejercicio usando Runge-Kutta de 2do Orden

b. Desarrolle tres iteraciones para θ(t) con tamaño de paso h=0.2

c. Usando el algoritmo, aproxime la solución entre t=0 a t=10 s, adjunte sus resultados en la evaluación.

d. Realice una observación sobre el movimiento estimado del péndulo a lo largo del tiempo.

Rúbrica: literal a (5 puntos), literal b (15 puntos), literal c (10 puntos), literal d (5 puntos)

Referencia: 2Eva_IT2019_T2 Péndulo vertical

Vista general de ecuaciones diferenciales I Capítulo 1, 6min 54s. 3Blue1Brown 31-Marzo-2023.

2Eva_2023PAOI_T1 Material para medalla de academia

2da Evaluación 2023-2024 PAO I. 29/Agosto/2023

Tema 1 (30 puntos) medalla area con integral numerico
Una academia encarga a un joyero un modelo de medalla cuyo costo unitario se determina por el área descrita entre las funciones f(x) y g(x) presentadas.

Se considera que el grosor de la medalla es único e independiente de la forma de la medalla.

f(x) = 2-8\Big( \frac{1}{2} - x \Big)^2 0 \le x \lt 1 g(x) = -\Big( 1-x\Big)\ln \Big( 1- x \Big)

Para el desarrollo numérico, use diferentes métodos de Simpson para cada función.

a. Realice el planteamiento de las ecuaciones para el ejercicio.

b. Describa el criterio usado para determinar el número de tramos usado en cada caso.

c. Desarrolle las expresiones completas del ejercicio para cada función.

d. Indique el resultado obtenido para el área requerida y la cota de error.

e. Encuentre el valor del tamaño de paso si se requiere una cota de error de 0.00032

Nota: en Python ln(x) se escribe np.log(x).

Rúbrica: literal a (5 puntos), literal b (5 puntos), literal c (10 puntos), literal d (5 puntos), literal e (5 puntos)

Referencia: Star Trek https://intl.startrek.com/
¿A quien se le ocurrió crear la moneda? | Discovery en Español Youtube.com 8 nov 2016.

2Eva_2022PAOII_T3 EDP Parabólica con coseno 3/4π

2da Evaluación 2022-2023 PAO II. 24/Enero/2023

Tema 3. (35 puntos) Aproxime la solución a la siguiente ecuación diferencial parcial parabólica

\frac{\partial^2 u}{\partial x^2} = b \frac{\partial u}{\partial t}

2Eva2022PAOII_T3 EDP ParabolicaCon las siguientes condiciones de frontera:
u(0,t)=1
u(1,t)=0

Y las condiciones iniciales
u(x,0) = \cos \Big( \frac{3π}{2}x\Big)

Utilice diferencias finitas centradas para x, para t hacia adelante.

a. Plantee las ecuaciones para usar un método numérico en un nodo i,j
b. Realice la gráfica de malla,
c. desarrolle y obtenga el modelo discreto para u(xi,tj)

Suponga que b = 2, Aproxime la solución con Δx = 0.2, Δt = Δx/100.

d. Realice al menos tres iteraciones en el eje tiempo.
e. Estime el error de u(xi,tj), y presente observaciones sobre la convergencia del método.

Rúbrica: literal a (5 puntos), literal b (5 puntos), literal c (5 puntos), literal d (15 puntos), literal e (5 puntos).

Referencia: Chapra & R. Canale (2010). Métodos Numéricos para Ingenieros. Ejercicio 30.15 p904,
Solving the heat equation | DE3. 3Blue1Brown 16 Junio 2019.