UNIVERSIDAD DE COLIMA
FACULTAD DE CIENCIAS QUÍMICAS
INGENIERO QUÍMICO EN ALIMENTOS
8vo. Semestre. Grupo C. IQA
SIMULACIÓN DE PROCESOS
Practica 2. Centro de computo
ERNESTO MEDINA AHUMADA
Profesor:
DRA. Julia Verónica Hernández Madrigal
1 de Marzo del 2024 . Coquimatlán, Colima
OBJETIVO GENERAL DE LA PRACTICA:
Elaborar la resolución de un problema en Excel mediante tres diferentes métodos,
el método general, el método del trapecio y el Simpson 1/3 mediante la resolución
de dos diferentes problemas.
INTRODUCCION:
Método general:
T −TR t T −TR
dt 10
t
F = 10 z z
t
0
0
Método del trapecio:
Método de Simpson 1/3
MATERIALES:
-computadora
-software Excel
Resolución del primer problema
Calculation of thermal death time F using the
General method
with numerical
integration
Input Calculated
data data
F°
z= 18
F°
T R= 250
Lethal Lethality
Time Temperature Δt rate F
LR (LR
(min) (Fo) LR average average*Δt) (T-250)/18
-
x0 0 140 0.0000 6.11111
-
x1 5 140 5.0 0.0000 0.000 0.000 6.11111
-
x2 10 140 5.0 0.0000 0.000 0.000 6.11111
-
x3 15 140 5.0 0.0000 0.000 0.000 6.11111
-
x4 20 163 5.0 0.0000 0.000 0.000 4.83333
-
x5 25 185 5.0 0.0002 0.000 0.001 3.61111
-
x6 30 201 5.0 0.0019 0.001 0.005 2.72222
-
x7 35 213 5.0 0.0088 0.005 0.027 2.05556
-
x8 40 224 5.0 0.0359 0.022 0.112 1.44444
-
x9 45 229.4 5.0 0.0717 0.054 0.269 1.14444
-
x10 50 234.5 5.0 0.1377 0.105 0.523 0.86111
-
x11 55 238 5.0 0.2154 0.177 0.883 0.66667
x12 60 241 5.0 0.3162 0.266 1.329 -0.5
-
x13 65 243.5 5.0 0.4354 0.376 1.879 0.36111
-
x14 70 245 5.0 0.5275 0.481 2.407 0.27778
-
x15 75 246.3 5.0 0.6229 0.575 2.876 0.20556
x16 80 247.3 5.0 0.7079 0.665 3.327 -0.15
-
x0 85 247 5.0 0.6813 0.695 3.473 0.16667
-
x1 90 245.2 5.0 0.5412 0.611 3.056 0.26667
-
x2 95 223.5 5.0 0.0337 0.287 1.437 1.47222
-
x3 100 175 5.0 0.0001 0.017 0.084 4.16667
-
x4 105 153.000 5.0 0.0000 0.000 0.000 5.38889
Integral
≈ 21.690
Result
F= 21.6899 min
Aquí los datos se calcularon de la siguiente manera
-Tanto el tiempo como la temperatura fueron dados por el problema además del TR
YZ
-Posteriormente se calculo el Δt que fue de 0.5 para cada casilla correspondiente
-el LR fue calculado mediante la formula del método general 10^(( T-TR)/z)
-EL LR average se calculó ejemplo :(F13+F14)/2
-La letalidad F se calculó de esta manera : (LR average*Δt) Como resultado fue 21.6899
minutos
- en la última columna fue (T-TR)/18
DANDO COMO RESULTADO LAS SIGUIENTES GRAFICAS
300
250
Temperature, oF
200
150
100
50
0
0 50 100 150
Time
TEMPERATURA VS TIEMPO
0.80
0.70
0.60
Lethal rate
0.50
0.40
0.30
0.20
0.10
0.00
0 50 100 150
Time
LETALIDAD VS TIEMPO
En la siguiente hoja de calculo se utilizo el método del trapesio y Simpson 1/3
Time Temperature x
(min) (Fo) (T-250)/18 f(X)=10^E
-
x0 0 140 CALENTAMIENTO 6.11111111 7.74264E-07
-
x1 5 140 CALENTAMIENTO 6.11111111 7.74264E-07
-
x2 10 140 CALENTAMIENTO 6.11111111 7.74264E-07
-
x3 15 140 CALENTAMIENTO 6.11111111 7.74264E-07
-
x4 20 163 CALENTAMIENTO 4.83333333 1.4678E-05
-
x5 25 185 CALENTAMIENTO 3.61111111 0.000244844
-
x6 30 201 CALENTAMIENTO 2.72222222 0.001895736
-
x7 35 213 CALENTAMIENTO 2.05555556 0.008799225
-
x8 40 224 CALENTAMIENTO 1.44444444 0.035938137
-
x9 45 229.4 CALENTAMIENTO 1.14444444 0.07170601
-
x10 50 234.5 CALENTAMIENTO 0.86111111 0.137685716
-
x11 55 238 CALENTAMIENTO 0.66666667 0.215443469
x12 60 241 CALENTAMIENTO -0.5 0.316227766
-
x13 65 243.5 CALENTAMIENTO 0.36111111 0.435400465
-
x14 70 245 CALENTAMIENTO 0.27777778 0.527499706
-
x15 75 246.3 CALENTAMIENTO 0.20555556 0.622937455
x16 80 247.3 CALENTAMIENTO -0.15 0.707945784
x0 80 247.3 ENFRIAMIENTO -0.15 0.707945784
-
x1 85 247 ENFRIAMIENTO 0.16666667 0.681292069
-
x2 90 245.2 ENFRIAMIENTO 0.26666667 0.541169527
-
x3 95 223.5 ENFRIAMIENTO 1.47222222 0.033711477
-
x4 100 175 ENFRIAMIENTO 4.16666667 6.81292E-05
-
x5 105 153.000 ENFRIAMIENTO 5.38888889 4.08424E-06
LETALIDAD 250
Se puede apreciar que se dieron 16 puntos de temperatura de calentamiento y los 5
posteriores fueron de enfriamiento como se puede apreciar en la primer grafica el declive
de la curva de calor en base al tiempo.
Los cálculos de las columnas consecuentes fue (T-TR)/18 y f(X)=10^E
Las formulas usadas para resolver son las siguientes
𝑡𝑛−𝑡0
trapecio h= y
𝑛
𝑡𝑛−𝑡0
simpson 1/3 dt= 𝑛
RESULTADOS
CALENTAMIENTO ENFRIAMIENTO
METODO DE TRAPECIO METODO DE TRAPECIO
H 5 H 5
H/2 2.5 H/2 2.5
∑(X1….X15) 2.373796305 ∑(X1….X4) 1.2562412
2*∑(X1….X15) 4.747592609 2*∑(X1….X15) 2.5124824
F0 13.63884792 F0 8.05108068
SIMPSON SIMPSON
H 5 H 5
H/3 1.666666667 H/3 1.66666667
∑(X1….X15)*4 5.418132068 ∑(X1….X3)*4 2.86001418
∑(X2….X14)*2 2.038525027 ∑(X2….X4)*2 1.08247531
F0 13.60767276 F0 7.75073227
TRAPECIO
FTOT 21.6899286
SIMPSON
FTOT 21.358405
CONCLUSION:
En la elaboración de esta practica pudimos comprender e interpretar los tres
métodos utilizados además de la creación de la hoja de calculo sirve como plantilla
para otros ejercicios que deseemos utilizar en base a estos métodos pues ya se
tiene la formula en cada casilla solamente seris cambiar los datos para que nos
arroje los cambios en las casillas consecuentes que ya cuentan con este código.
Bibliografía
Ibartz A. y Barbosa-Cánovas. G. (2005) Operaciones unitarias en ingeniería de
alimentos. Madrid españa: Mundi Prensa.
Singh P. y Heldman D., (2009) Introduction to Food Engineering; Academic Press.
London: 4th edition. Elsevier Publishing
Rao, Sankara (2007). «7.6 Newton-Cotes integration formulae». Numerical Methods
For Scientists And Engineers (en inglés) (3ª edición). New Delhi (India): Prentice-
Hall of India Learning Private. pp. 151-159
Regla de Simpson
https://es.wikipedia.org/wiki/Regla_de_Simpson recopilado el 29 de febrero del
2024
