UNIVERSIDAD NACIONAL DE

TRUJILLO
ESCUELA DE POSGRADO

Maestría en Ingeniería Química

Mención en Ingeniería Ambiental

DOCENTE:
• Dr. Guillermo Evangelista Benites

Estudiante:
• Urquiaga Ríos, Javier Francisco
• Zelada Sangay, Manuel Alfredo.

CURSO:
SIMULACIÓN Y CÁLCULOS EN INGENIERÍA
QUÍMICA

TRUJILLO – PERU
2022

ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR

1 DEFINICIÓN DE ECUACION DIFERENCIAL
Se llama ecuación diferencial aquella ecuación que contiene una variable
dependiente y sus derivadas con respecto a una o más variables
independientes. Muchas de las leyes generales de la naturaleza se expresan en el
lenguaje de las ecuaciones diferenciales; abundan también las aplicaciones en
ingeniería, economía, en las mismas matemáticas y en muchos otros campos de la
ciencia aplicada. (Nieves Hurtado & Domínguez Sánchez, 2014, p. 536).
Por ejemplo, la siguiente ecuación:
dv c
=g− v (1)
dt m
Donde g es la constante gravitacional, m es la masa y c es el coeficiente de arrastre.
En la ecuación, la cantidad que se está derivando, v, se conoce como variable
dependiente. La cantidad con respecto a la cual v se está derivando, t, se conoce
como variable independiente.
- Cuando la función tiene una variable independiente, la ecuación se llama
ecuación diferencial ordinaria (o EDO).
- Esto contrasta con una ecuación diferencial parcial (o EDP) que involucra dos o
más variables independientes.

2 ECUACIONES DIFERENCIALES ORDINARIAS (EDO)

Por otro lado, la ecuación diferencial ordinaria se representa de las siguientes
maneras:
dh
A =G (2)
dt
2
d y
m 2
=ky (3)
dt
2
d y dx
m 2
+ c + kx=0 (4)
dt dt
La variable dependiente en cada una de estas ecuaciones es y, y la variable
independiente es x o t. Las letras k, m y c representan constantes. Una ecuación
diferencial es ordinaria si sólo tiene una variable independiente, por lo que todas
las derivadas que tiene son ordinarias o totales.

3 ECUACIONES DIFERENCIALES ORDINARIAS (EDO) DE ORDEN

SUPERIOR
Las ecuaciones diferenciales se clasifican también en cuanto a su orden. Por
ejemplo, la ecuación (2) se denomina EDO de primer orden, ya que la derivada mayor
es una primera derivada. Una EDO de segundo orden tiene una segunda derivada,
como la mayor. Por ejemplo, la ecuación (4) que describe la posición x de un sistema
masa-resorte con amortiguamiento es la EDO de segundo orden.
Forma general EDO
n
d y ˊ ˊˊ ( n−1)
n
=f ( x , y , y , y , … . , y )
dx
( n−1 )
ˊ ˊ
y ( x 0 )= y 0 , y ( x 0 )= y 0 , … , y ( x 0 )= y (n−1)
0

y ( x f ) =?

2
d y dx
m 2
+ c + kx=0
dt dt
Donde c es un coeficiente de amortiguamiento y k es una constante del resorte. De
manera similar, una ecuación de n-ésimo orden tiene una n-ésima derivada, como la
mayor.
Las ecuaciones de orden superior pueden reducirse a un sistema de ecuaciones
de primer orden. Para la ecuación (4), esto se logra al definir una nueva variable y,
donde:
dx
y= (5)
dt
que al derivar con respecto a t obtiene
2
dx d x
= (6)
dt dt 2
Las ecuaciones (5) y (6) se sustituyen después en la ecuación (4) para llegar a
dy
m +cy + kx=0 (7)
dt
ó
dy −cy +kx
= (8)
dt m

Así, las ecuaciones (5) y (8) son un sistema de ecuaciones diferenciales de primer
orden, equivalentes a la ecuación de segundo orden original. Como otras ecuaciones
diferenciales de n-ésimo orden pueden reducirse en forma similar, esta parte de
nuestro libro se concentra en la solución de ecuaciones diferenciales de primer orden.
Aplicación
2
d x dx dx (0)
2
+ 3 + 2 x=0 ; condiciones iniciales x ( 0 )=5 ; =10
dt dt dt
 2
d x dx
2
=−3 −2 x
dt dt

y = dsolve("D2x = -3*D1x - 2*x","x(0) = 5, Dx(0) = 10", "t")

tiempoinicial=0
tiempofinal=10
for t=tiempoinicial:0.1:tiempofinal
ysalida = eval(y);
tsalida = t;
figure (1)
plot(tsalida,ysalida,".k","LineWidth",1,"MarkerSize",16)
xlabel("tsalida","FontSize",10,"Color","b")
ylabel("ysalida","FontSize",10,"color","c")
title("Ejemplo de EDO de orden superior","FontSize",10,"color","b")
hold on;
grid on;

4 SISTEMAS DE ECUACIONES DIFERENCIALES ORDINARIAS

En un sistema de ecuaciones diferenciales ordinarias de cualquier orden, puede ser
reducido a un sistema equivalente de primer orden, si se introducen nuevas variables y
ecuaciones.
Ejemplo:
El siguiente es un modelo dinámico propuesto por Aiken y Lapidus para un reactor de
lecho fluidizado:
' 4
y 1=1.30 ( y 3 − y 1 )+ 1.04∗10 k y 2
' 3
y 2=1.88∗10 ( y 4− y 2 (1+k ))
 '
y 3=1752+266.7 y 1 −269.3 y 3
'
y 4 =0.1+320 y 2−321 y 4
Condiciones iniciales:

y 1 ( 0 )=759.167

y 2 ( 0 )=0

y 3 ( 0 )=600

y 4 ( 0 )=0.1
Donde:

k =0.006 e
( 20.7 −
y )
15000
1

y 1 es la temperatura de la fase sólida.

y 3 es la temperatura de fase líquida.

y 2 , y 4 son las presiones parciales de cada una de las fases.

Determine los perfiles de las temperaturas y de las presiones y 1, y 2 / y 3 , y 4

respectivamente.

Desarrollamos el Matlab

clearall,clf,clc
M=eye(4,4);M(4,4)=0;
options=odeset('Mass',M,'MassSingular','yes');
%tspan=linspace(0,1000,100); tspan=[0 4*logspace(-6,6)];
y0=[759.167 0 600 0.1];
[t,y]=ode15s(@Ej_mio,tspan,y0,options);
for i=1:4
subplot(2,3,i)
plot(t,y(:,i),'-o','lineWidth',2),grid on
title(['y_',int2str(i)]),xlabel('t'),xlim([0 1000])
end
%disp([t y]) para obtener los datos numéricos
% Escriba la función anidada
function rh=Ej_mio(t,y) 7
if y(2)< 0, error('y(2) es negativo en la función DAE')
 end
k=6e-4*exp(20.7-15000/y(1));
rh=[1.3*(y(3)-y(1))+1.04e+4*k*y(2);
1.88e+3*(y(4)-y(2)*(1+k));
1752+266.7*y(1)-269.3*y(3);
0.1+320*y(2)-321*y(4)];
end

5 FORMULACIÓN DEL PROBLEMA DE VALORES EN LA FRONTERA

Un problema de valores en la frontera (PVF), para ecuaciones diferenciales ordinarias,
puede estar dado, por ejemplo, por una EDO de segundo orden dos condiciones de
frontera: CF1 y CF2.
2
d y
=f ( x , y , y )
ˊ
EDO 2
dt
CF 1 y ( x 0 ) = y 0
CF 2 y ( x f )= y f
y ( x )=? para x 0 < x< x f
Obsérvese que ahora se da como información dos puntos distintos por donde pasa la
curva desconocida y, solución de la EDO; es decir, conocemos el valor de y
correspondiente a dos abscisas distintas: xo y xf y queremos conocer el valor de y en el
intervalo (xo y xf). Esto se ilustra gráficamente en la figura. (Chapra & Canale, 2007,
pág. 794)

Desde luego también contamos para encontrar a y con su segunda derivada, esto es

f ( x , y , y ˊ ). Este tipo de problemas surgen, por ejemplo, cuando se resuelven

ecuaciones diferenciales parciales analíticamente.

Por ejemplo:

Se puede utilizar la conservación del calor para desarrollar un balance de calor para
una barra larga y delgada (figura). Si la barra no está aislada en toda su longitud y el
sistema se encuentra en estado estacionario, la ecuación resultante es:
 2
d T ˊ
2
+ h ( T a−T )=0
dx

Donde, h′ es un coeficiente de transferencia de calor (m –2) que parametriza la

velocidad con que se disipa el calor en el medio ambiente, y T a es la temperatura del
medio ambiente (°C). Las temperaturas en los extremos de la barra son: T (0) = T 1 y
T(L) = T2. La barra es de 10 metros con Ta = 20, T1 = 40, T2 = 200 y h′ = 0.01.

Determina la ecuación para encontrar el valor de T :

Respuesta: T =73.4523 e 0.1 x −53.4523 e−0.1 x + 20

6 ECUACIÓN DIFERENCIAL ORDINARIA RÍGIDA

La rigidez es un problema especial que puede surgir en la solución de ecuaciones

diferenciales ordinarias. Un sistema rígido es uno que tiene componentes que cambian
rápidamente, junto con componentes de cambio lento. En muchos casos, los
componentes de variación rápida son transitorios, desaparecen rápidamente, después
de lo cual los componentes de variación lenta son los que dominan la solución. Pero,
aunque el fenómeno transitorio existe en sólo un período pequeño, puede determinar
el paso del tiempo en la aplicación de un método numérico. (Conceptos básicos. (s/f).
Edu.ar. Recuperado el 2 de julio de 2022, s.f.)

Ejemplo:
Considere las siguientes reacciones irreversibles simples en serie:
A k1 B k2 C
→ →

El conjunto de Ecuaciones diferenciales para un reactor de mezcla completa es:

dCA
=−r 1=−k 1 C A
dt
 d CB
=r 1−r 2=k 1 C A−k 2 C B
dt
d CC
=r 2=−k 2 C B
dt

C A (0)=1 C B (0)=0 C C (0)=0

Las soluciones analíticas son:

C A (t ) −k B t
=e
C A0
C B(t ) kB
= [ e−k t −e−k
B C
t
]
CA 0 k C −k B

CC (t ) C () C
=1− A t − B(t )
CA 0 C A 0 CA 0

Donde: k 1=2 k 2=1

Hallar las concentraciones de C A, C By C c para un tiempo de 2 segundos.

Desarrollamos el Matlab

Clear all, clf, clc

Eq=inline('[-2*x(1);2*x(1)-x(2);x(2)]','t','x');

[t,x]=ode15s(Eq,[0:0.01:2],[1 0 0]);

disp([t,x])

plot(t,x(:,1),'-r',t,x(:,2),'-b',t,x(:,3),'-g','LineWidth',2),grid on

title('Solución de un sistema de ecuaciones reactivas','Fontsize',15);

text(0.2,0.8,'CA');text(0.8,0.1,'CB');text(0.3,0.3,'CC');

xlabel('Tiempo');ylabel('Concentración');

A la Salida del programa

T CA CB Cc

0 1.0000 0 0

0.0100 0.9802 0.0197 0.0001

0.0200 0.9608 0.0388 0.0004

0.0300 0.9418 0.0574 0.0009

0.0400 0.9231 0.0754 0.0015

0.0500 0.9048 0.0928 0.0024

0.0600 0.8869 0.1097 0.0034

0.0700 0.8694 0.1261 0.0046

0.0800 0.8521 0.1420 0.0059

0.0900 0.8353 0.1573 0.0074

0.1000 0.8187 0.1722 0.0091

0.1100 0.8025 0.1866 0.0108

0.1200 0.7866 0.2006 0.0128

0.1300 0.7711 0.2141 0.0149

0.1400 0.7558 0.2271 0.0171

0.1500 0.7408 0.2398 0.0194

0.1600 0.7261 0.2520 0.0219

0.1700 0.7118 0.2638 0.0244

0.1800 0.6977 0.2752 0.0271

0.1900 0.6839 0.2862 0.0299

0.2000 0.6703 0.2968 0.0329

0.2100 0.6570 0.3071 0.0359

0.2200 0.6440 0.3170 0.0390

0.2300 0.6313 0.3265 0.0422

0.2400 0.6188 0.3357 0.0455

0.2500 0.6065 0.3445 0.0489

0.2600 0.5945 0.3531 0.0524

0.2700 0.5827 0.3613 0.0560

0.2800 0.5712 0.3692 0.0596

0.2900 0.5599 0.3767 0.0634

0.3000 0.5488 0.3840 0.0672

0.3100 0.5379 0.3910 0.0710

0.3200 0.5273 0.3977 0.0750

0.3300 0.5168 0.4042 0.0790

0.3400 0.5066 0.4103 0.0831

0.3500 0.4966 0.4162 0.0872

0.3600 0.4867 0.4219 0.0914

0.3700 0.4771 0.4273 0.0956

0.3800 0.4677 0.4324 0.0999

0.3900 0.4584 0.4373 0.1043

0.4000 0.4493 0.4420 0.1087

0.4100 0.4404 0.4465 0.1131

0.4200 0.4317 0.4507 0.1176

0.4300 0.4232 0.4547 0.1221

0.4400 0.4148 0.4585 0.1267

0.4500 0.4066 0.4621 0.1313

0.4600 0.3985 0.4655 0.1359

0.4700 0.3906 0.4688 0.1406

0.4800 0.3829 0.4718 0.1453

0.4900 0.3753 0.4747 0.1500

0.5000 0.3679 0.4773 0.1548

0.5100 0.3606 0.4798 0.1596

0.5200 0.3534 0.4822 0.1644

0.5300 0.3464 0.4843 0.1692

0.5400 0.3396 0.4863 0.1741

0.5500 0.3329 0.4882 0.1790

0.5600 0.3263 0.4899 0.1839

0.5700 0.3198 0.4914 0.1888

0.5800 0.3135 0.4928 0.1937

0.5900 0.3073 0.4941 0.1986

0.6000 0.3012 0.4953 0.2036

0.6100 0.2952 0.4963 0.2085

0.6200 0.2894 0.4971 0.2135

0.6300 0.2836 0.4979 0.2185

0.6400 0.2780 0.4985 0.2234

0.6500 0.2725 0.4991 0.2284

0.6600 0.2671 0.4995 0.2334

0.6700 0.2618 0.4998 0.2384

0.6800 0.2566 0.5000 0.2434

0.6900 0.2515 0.5000 0.2484

0.7000 0.2466 0.5000 0.2534

0.7100 0.2417 0.4999 0.2584

0.7200 0.2369 0.4997 0.2634

0.7300 0.2322 0.4994 0.2684

0.7400 0.2276 0.4990 0.2734

0.7500 0.2231 0.4986 0.2784

0.7600 0.2187 0.4980 0.2833

0.7700 0.2143 0.4974 0.2883

0.7800 0.2101 0.4966 0.2933

0.7900 0.2059 0.4958 0.2982

0.8000 0.2018 0.4950 0.3032

0.8100 0.1978 0.4940 0.3081

0.8200 0.1939 0.4930 0.3131

0.8300 0.1901 0.4920 0.3180

0.8400 0.1863 0.4908 0.3229

0.8500 0.1826 0.4896 0.3278

0.8600 0.1790 0.4884 0.3327

0.8700 0.1754 0.4871 0.3375

0.8800 0.1719 0.4857 0.3424

0.8900 0.1685 0.4842 0.3472

0.9000 0.1652 0.4828 0.3521

0.9100 0.1619 0.4812 0.3569

0.9200 0.1587 0.4796 0.3617

0.9300 0.1556 0.4780 0.3665

0.9400 0.1525 0.4763 0.3712

0.9500 0.1494 0.4746 0.3760

0.9600 0.1465 0.4728 0.3807

0.9700 0.1436 0.4710 0.3854

0.9800 0.1407 0.4692 0.3901

0.9900 0.1379 0.4673 0.3948

1.0000 0.1352 0.4654 0.3995

1.0100 0.1325 0.4634 0.4041

1.0200 0.1299 0.4614 0.4087

1.0300 0.1273 0.4594 0.4133

1.0400 0.1248 0.4573 0.4179

1.0500 0.1223 0.4553 0.4224

1.0600 0.1199 0.4531 0.4270

1.0700 0.1175 0.4510 0.4315

1.0800 0.1152 0.4488 0.4360

1.0900 0.1129 0.4466 0.4405

1.1000 0.1107 0.4444 0.4449

1.1100 0.1085 0.4422 0.4494

1.1200 0.1063 0.4399 0.4538

1.1300 0.1042 0.4376 0.4582

1.1400 0.1021 0.4354 0.4625

1.1500 0.1001 0.4330 0.4669

1.1600 0.0981 0.4307 0.4712

1.1700 0.0962 0.4284 0.4755

1.1800 0.0943 0.4260 0.4797

1.1900 0.0924 0.4236 0.4840

1.2000 0.0906 0.4212 0.4882

1.2100 0.0888 0.4188 0.4924

1.2200 0.0870 0.4164 0.4966

1.2300 0.0853 0.4140 0.5007

1.2400 0.0836 0.4115 0.5049

1.2500 0.0819 0.4091 0.5090

1.2600 0.0803 0.4066 0.5130

1.2700 0.0787 0.4042 0.5171

1.2800 0.0772 0.4017 0.5211

1.2900 0.0756 0.3992 0.5251

1.3000 0.0741 0.3968 0.5291

1.3100 0.0727 0.3943 0.5331

1.3200 0.0712 0.3918 0.5370

1.3300 0.0698 0.3893 0.5409

1.3400 0.0684 0.3868 0.5448

1.3500 0.0671 0.3843 0.5486

1.3600 0.0657 0.3818 0.5525

1.3700 0.0644 0.3793 0.5563

1.3800 0.0632 0.3768 0.5600

1.3900 0.0619 0.3743 0.5638

1.4000 0.0607 0.3718 0.5675

1.4100 0.0595 0.3693 0.5712

1.4200 0.0583 0.3668 0.5749

1.4300 0.0572 0.3643 0.5786

1.4400 0.0560 0.3618 0.5822

1.4500 0.0549 0.3593 0.5858

1.4600 0.0538 0.3568 0.5894

1.4700 0.0528 0.3543 0.5929

1.4800 0.0517 0.3518 0.5965

1.4900 0.0507 0.3494 0.6000

1.5000 0.0497 0.3469 0.6034

1.5100 0.0487 0.3444 0.6069

1.5200 0.0477 0.3419 0.6103

1.5300 0.0468 0.3395 0.6137

1.5400 0.0459 0.3370 0.6171

1.5500 0.0449 0.3346 0.6205

1.5600 0.0441 0.3321 0.6238

1.5700 0.0432 0.3297 0.6271

1.5800 0.0423 0.3273 0.6304

1.5900 0.0415 0.3249 0.6337

1.6000 0.0407 0.3224 0.6369

1.6100 0.0399 0.3200 0.6401

1.6200 0.0391 0.3176 0.6433

1.6300 0.0383 0.3153 0.6465

1.6400 0.0375 0.3129 0.6496

1.6500 0.0368 0.3105 0.6527

1.6600 0.0361 0.3081 0.6558

1.6700 0.0353 0.3058 0.6589

1.6800 0.0346 0.3034 0.6619

1.6900 0.0340 0.3011 0.6649

1.7000 0.0333 0.2988 0.6679

1.7100 0.0326 0.2965 0.6709

1.7200 0.0320 0.2942 0.6739

1.7300 0.0313 0.2919 0.6768

1.7400 0.0307 0.2896 0.6797

1.7500 0.0301 0.2873 0.6826

1.7600 0.0295 0.2850 0.6854

1.7700 0.0289 0.2828 0.6883

1.7800 0.0284 0.2805 0.6911

1.7900 0.0278 0.2783 0.6939

1.8000 0.0272 0.2761 0.6967

1.8100 0.0267 0.2739 0.6994

1.8200 0.0262 0.2717 0.7021

1.8300 0.0257 0.2695 0.7049

1.8400 0.0251 0.2673 0.7075

1.8500 0.0246 0.2652 0.7102

1.8600 0.0242 0.2630 0.7128

1.8700 0.0237 0.2609 0.7155

1.8800 0.0232 0.2587 0.7181

1.8900 0.0227 0.2566 0.7206

1.9000 0.0223 0.2545 0.7232

1.9100 0.0219 0.2524 0.7257

 1.9200 0.0214 0.2504 0.7282

1.9300 0.0210 0.2483 0.7307

1.9400 0.0206 0.2462 0.7332

1.9500 0.0202 0.2442 0.7356

1.9600 0.0198 0.2422 0.7381

1.9700 0.0194 0.2401 0.7405

1.9800 0.0190 0.2381 0.7429

1.9900 0.0186 0.2361 0.7453

2.0000 0.0182 0.2341 0.7476

7 REFERENCIAS
Chapra, S., & Canale, R. (2007). Métodos numéricos para ingenieros. México:
McGRAW-HILL/INTERAMERICANA EDITORES, S.A. DE C.V.
Conceptos básicos. (s/f). Edu.ar. Recuperado el 2 de julio de 2022, d. (s.f.). Obtenido
de http://www.frsn.utn.edu.ar/gie/an/mnedo/36_rigidas.html
Nieves Hurtado, A., & Domínguez Sánchez, F. (2014). Métodos Numéricos: Aplicados
a la Ingeniería. México, D.F.: GRUPO PATRIA CULTURAL, S.A. DE C.V.