Problemas de Ecuaciones No Lineales
Problemas de Ecuaciones No Lineales
Problemas de Ecuaciones No Lineales
. (1)
Donde:
P
T
R
V
a, b
Calcular el volumen molar V a 80C para una presin P=10 atm, para el gas CO 2, cuyas
constantes a=3,599 y b=0,04267.
Solucin del problema
La ecuacin (1) de Van der Walls para un gas real al multiplicarla por V 2 y reordenando sus
trminos se puede escribir como:
3 2 2 + = 0
. (2)
3 ( + ) 2 + = 0
. (3)
De esta manera la ecuacin (3) es una ecuacin no lineal de la forma f(V)=0, la cul puede ser
resuelta usando los mtodos de punto fijo, newton raphson, secante, posicin falsa, biseccin,
etc.
En primera instancia usaremos el mtodo de newton raphson donde F representa la funcin y
D su derivada:
= 3 ( + ) 2 +
. (4)
= 2 2 2( + ) +
. (5)
. (6)
De esta forma establecemos una formulacin en hoja de clculo Excel para resolver el proceso
iterativo hasta lograr el valor V deseado a un nivel de error preestablecido.
La siguiente hoja se muestra haciendo referencia al contenido de las celdas, es decir se pueden
observar las formulas ingresadas para lograr luego de 5 iteraciones el valor requerido.
De esta forma podemos concluir que el valor del volumen molar de CO 2, V para una presin de
10 atm y 80 C de temperatura es de 2,81475 L/gmol. Cabe destacar que solo cambiando los
contenidos de las celdas B2 a B6, podemos hacer clculos adicionales variando presin,
temperatura o tipo de gas.
Ahora plantearemos una solucin diferente, establecemos un cdigo computacional en
MATLAB, para la solucin del algoritmo de newton raphson.
clc, clear;
%Datos de entrada
P=10;
T=353.2;
R=0.08205;
a=3.599;
b=0.04267;
E=0.0001; %Proceso de clculo e iniciacin
V=R*T/P;
d=0.01;
while abs(d)>E
F=P*V^3-(P*b+R*T)*V^2+a*V-a*b;
D=2*P*V^2-2*(P*b+R*T)*V+a;
d=F/D;
V1=V-d;
V=V1;
end
Vol=V
Este programa codificado en MATLAB, arroja como resultado: Vol = 2.8147e+000
ESCUELA DE INGENIERIA GEOLOGICA-GEOTECNIA, UNJBG TACNA
TELEFONOS: 952298638, 241595, #701907 E-MAIL: [email protected]
5 + 22 3 4
120
. (1)
Usar el mtodo de Newton Raphson para determinar el punto de mxima deflexin (es decir el
valor de x donde dy/dx=0). Usar L=450 cm, E=50 000 kN/cm2, I=30 000 cm4, wo=1,75 kN/cm.
Solucin del problema
En primer lugar, para determinar la mxima deflexin debemos plantear que dy/dx=0, esto
garantizara que obtengamos un mximo. Por tanto, la funcin a tratar ser:
=
=
5 + 22 3 4
120
5 4 + 62 2 4
120
. (2)
. (3)
=
5 4 + 62 2 4
120
20 3 + 122
120
. (4)
. (5)
. (6)
De esta forma establecemos una formulacin en hoja de clculo Excel para resolver el proceso
iterativo hasta lograr el valor x deseado a un nivel de error preestablecido.
La siguiente hoja se muestra haciendo referencia al contenido de las celdas, es decir se pueden
observar las formulas ingresadas para lograr luego de 5 iteraciones el valor requerido.
De esta forma podemos concluir que el valor de LA mxima deflexin en la viga ser de 0,1141
cm hacia abajo. Ahora plantearemos una solucin diferente, establecemos un cdigo
computacional en MATLAB, para la solucin del algoritmo de newton raphson.
ESCUELA DE INGENIERIA GEOLOGICA-GEOTECNIA, UNJBG TACNA
TELEFONOS: 952298638, 241595, #701907 E-MAIL: [email protected]
1 2/3 1/2
. (1)
Donde:
Q
n
A
R
S
: Caudal en m3/s.
: Coeficiente de rugosidad de manning adimensional.
: rea Hidrulica en m.
: Radio hidrulico.
: Pendiente longitudinal de fondo del canal.
Usar el mtodo de Newton Raphson para determinar el tirante de agua para una seccin
trapezoidal. Si se sabe que el caudal es Q=2,3 m3/s, ancho de solera b=1,5 m., talud Z=1,5,
rugosidad n=0,014.
Solucin del problema
En primer lugar, debemos establecer la forma de la funcin. Por tanto, la funcin a tratar ser:
= 2/3
1/2
= 1/2
. (2)
3
. (3)
Como R=A/P
=
1/2
2
. (4)
4
5 4
3
. (5)
Donde:
T=b+2*Z*h
L=(1+Z2)1/2
A=b*h+Z*h2
P=b+2*h*L
El algoritmo de newton raphson puede ser escrito de la siguiente manera:
+1 = ( )/( )
. (6)
De esta forma establecemos una formulacin en hoja de clculo Excel para resolver el proceso
iterativo hasta lograr el valor h deseado a un nivel de error preestablecido.
La siguiente hoja se muestra haciendo referencia al contenido de las celdas, es decir se pueden
observar las formulas ingresadas para lograr luego de 5 iteraciones el valor requerido.
De esta forma podemos concluir que el valor del tirante de agua ser de 0,8274 m. Ahora
plantearemos una solucin diferente, establecemos un cdigo computacional en MATLAB,
para la solucin del algoritmo de newton raphson.
clc, clear; %Datos de entrada
b=1.5;
Q=2.3;
n=0.014;
z=1.5;
So=0.0005; %Proceso de clculo e iniciacin
y=b/2;
precis=0.00001;
d=0.01;
while abs(d)>precis
A=(b*y+z*y^2);
P=(b+2*y*sqrt(1+z^2));
R=A/P;
dA=b+2*z*y;
dP=2*sqrt(1+z^2);
f=A*R^(2/3)-Q*n/sqrt(So);
df=(2/3)*A*R^(-1/3)*((P*dA-A*dP)/P^2)+R^(2/3)*dA;
d=f/df;
y1=y-d;
y=y1;
end
h=y
ESCUELA DE INGENIERIA GEOLOGICA-GEOTECNIA, UNJBG TACNA
TELEFONOS: 952298638, 241595, #701907 E-MAIL: [email protected]
. (1)
... (2)
108
0,00012378 1 = 0
0,00012378
. (1)
Producto
Alimentacin
Reactor de flujo tipo tapn
Reciclante
Suponga que se est procesando una sustancia qumica A, para generar un producto B. Para el
caso que A forma a B para una reaccin autocataltica (es decir que uno de los productos acta
como catalizador o estimulante en la reaccin), es posible demostrar que una razn optima de
la recirculacin debe satisfacer:
1 + (1 )
+1
=
(1 )
1 + (1 )
. (2)
Donde XAf, es la fraccin del reactante A que se convierte en el producto B. La razn ptima de
recirculacin corresponde a un reactor de tamao mnimo necesario para alcanzar el nivel
deseado de conversin.
Utilice al menos dos mtodos numricos (biseccin y newton raphson) para determinar la
razn de recirculacin necesaria, de manera que se minimice el tamao del reactor para una
conversin fraccional XAf igual a 0,9. La solucin debe incluir:
a) Formulacin en hoja de clculo para los dos mtodos.
b) Cdigo MATLAB para cada caso, e incluir una solucin grafica.
. (1)
. (1)
=
1+
. (2)