¿Alguna vez has pensado cómo los ingenieros de Medellín predicen la contaminación del aire en la ciudad? ¿O cómo los diseñadores de buses articulados mejoran su aerodinámica para gastar menos combustible? La respuesta está en la Dinámica de Fluidos Computacional (CFD), una herramienta que usa tu mismo computador para simular el comportamiento de fluidos como el aire, el agua o incluso el café que tomas en el bus. En Colombia, donde la minería, el transporte masivo y la gestión del agua son clave, dominar CFD te abre puertas en empresas de ingeniería, ambientalismo y hasta en el diseño de canoas para el Amazonas. Pero ojo: no basta con saber teoría, hay que aplicarla. Por eso hoy resolverás ejercicios basados en problemas reales de Bogotá, Medellín y la Costa Caribe. ¡Vamos a simular!
Fuga de agua en el acueducto de Medellín: ¿Cuánto se pierde?
En el barrio Laureles de Medellín, se detectó una fuga en una tubería principal que abastece a 500 hogares. El diámetro de la tubería es de 30 cm y el agua sale a una velocidad de 5 m/s. Calcula el caudal volumétrico que se está perdiendo. Si el costo del agua en Medellín es de 3 500 COP por metro cúbico, ¿cuánto dinero se pierde por hora si la fuga no se repara?
Datos
| d | diámetro de la tubería | 30 | cm |
| v | velocidad del agua | 5 | m/s |
| costo_agua | costo del agua | 3500 | COP/m³ |
| n_hogares | número de hogares afectados | 500 |
Se busca
- Q — caudal volumétrico perdido (m³/s)
- pérdida_diaria — pérdida económica diaria (COP)
Pistas progresivas
Pista 1
Recuerda que el caudal volumétrico se calcula como , donde es el área transversal de la tubería.
Pista 2
El área de un círculo es . Convierte primero el diámetro a radio en metros.
Pista 3
Para la pérdida económica, multiplica el caudal por el tiempo (1 hora) y luego por el costo por metro cúbico.
Solución completa
- Cálculo del área transversal — Primero calculamos el área de la tubería usando el diámetro dado. Convertimos 30 cm a metros: m. Luego aplicamos la fórmula del área de un círculo.
- Cálculo del caudal volumétrico — Ahora aplicamos la ecuación de continuidad para fluidos incompresibles: el caudal volumétrico es el producto del área por la velocidad.
- Cálculo de la pérdida económica por hora — Multiplicamos el caudal por el tiempo de 1 hora (3600 segundos) para obtener el volumen total perdido en ese tiempo. Luego multiplicamos por el costo del agua.
- Cálculo final de la pérdida económica — Multiplicamos el volumen perdido por el costo por metro cúbico para obtener la pérdida en pesos colombianos.
→ El caudal perdido es 0.177 m³/s y la pérdida económica por hora es 2 233 200 COP.
Smog en Bogotá: ¿Se dispersa o se queda atrapado?
En el centro de Bogotá, durante un día sin viento, se libera un contaminante (CO₂) desde una chimenea de 20 m de altura. La velocidad del viento en la superficie es de 1 m/s y la temperatura ambiente es de 20°C. Usando un modelo simplificado de CFD, determina la concentración del contaminante a 100 m de distancia horizontal desde la chimenea si el coeficiente de difusión turbulenta es m²/s. Considera que la emisión es continua y uniforme.
Datos
| h | altura de la chimenea | 20 | m |
| v | velocidad del viento | 1 | m/s |
| T | temperatura ambiente | 20 | °C |
| D | coeficiente de difusión turbulenta | 0.5 | m²/s |
| x | distancia horizontal | 100 | m |
Se busca
- C — concentración del contaminante a 100 m (kg/m³)
Pistas progresivas
Pista 1
Modela el problema como una pluma de contaminante en un flujo unidimensional. Usa la ecuación de advección-difusión simplificada: donde es la tasa de emisión (supón kg/s para este ejercicio).
Pista 2
Simplifica la expresión exponencial: .
Pista 3
Asegúrate de que las unidades sean consistentes: la concentración resultante estará en kg/m³.
Solución completa
- Ecuación simplificada de la pluma — Para un contaminante emitido de forma continua, la concentración en un punto aguas abajo se puede aproximar con la ecuación de la pluma de Gauss simplificada. Sustituimos los valores conocidos.
- Sustitución de valores — Reemplazamos los valores numéricos en la ecuación. Usamos kg/s como tasa de emisión típica para una chimenea industrial pequeña.
- Cálculo del denominador — Calculamos primero el término dentro de la raíz cuadrada y luego la raíz misma.
- Cálculo del exponente — Evaluamos la parte exponencial de la ecuación.
- Resultado final — Combinamos todos los términos para obtener la concentración final.
→ La concentración del contaminante a 100 m de distancia es aproximadamente kg/m³, lo que indica una dispersión casi total.
Diseño aerodinámico del Metro de Medellín: ¿Cuánto ahorra el nuevo vagón?
El Metro de Medellín está evaluando el diseño de un nuevo vagón articulado para reducir la resistencia al aire. Actualmente, un vagón tiene una resistencia de 5 000 N a 60 km/h. Si se rediseña la forma frontal con un coeficiente de arrastre que se reduce de 0.8 a 0.6, calcula la nueva fuerza de resistencia. ¿Cuánto combustible se ahorra por viaje redondo (ida y vuelta) entre las estaciones Niquía y La Estrella (18 km en total)? Considera que el motor tiene una eficiencia del 30% y el combustible cuesta 10 000 COP por litro con un poder calorífico de 35 MJ/L.
Datos
| F_actual | fuerza de resistencia actual | 5000 | N |
| C_d_actual | coeficiente de arrastre actual | 0.8 | |
| C_d_nuevo | coeficiente de arrastre nuevo | 0.6 | |
| v | velocidad del vagón | 60 | km/h |
| distancia | distancia viaje redondo | 18 | km |
| eficiencia | eficiencia del motor | 30 | % |
| costo_combustible | costo del combustible | 10000 | COP/L |
| poder_calorífico | poder calorífico del combustible | 35 | MJ/L |
Se busca
- F_nuevo — fuerza de resistencia nueva (N)
- ahorro_combustible — ahorro de combustible por viaje (L)
Pistas progresivas
Pista 1
La fuerza de resistencia se calcula como , donde es la densidad del aire (1.225 kg/m³) y es el área frontal (supón m² para este vagón).
Pista 2
Como la velocidad y el área frontal se mantienen constantes, la fuerza es directamente proporcional al coeficiente de arrastre.
Pista 3
El trabajo realizado por la fuerza de resistencia es . La energía necesaria proviene del combustible, considerando la eficiencia del motor.
Solución completa
- Cálculo de la fuerza nueva — La fuerza de resistencia es proporcional al coeficiente de arrastre. Calculamos la nueva fuerza usando la relación entre los coeficientes.
- Cálculo del trabajo ahorrado por viaje — Calculamos el trabajo ahorrado por viaje redondo. Primero convertimos la distancia a metros: 18 km = 18 000 m.
- Cálculo de la energía necesaria del combustible — La energía requerida se obtiene dividiendo el trabajo por la eficiencia del motor (expresada como decimal).
- Cálculo del ahorro de combustible — Dividimos la energía necesaria por el poder calorífico del combustible para obtener los litros ahorrados.
- Cálculo del ahorro económico — Multiplicamos el volumen ahorrado por el costo del combustible para obtener el valor en pesos.
→ La nueva fuerza de resistencia es 3 750 N. El ahorro de combustible por viaje redondo es aproximadamente 0.076 litros, lo que equivale a 760 COP.
Contaminación en el río Magdalena: ¿Hasta dónde llega el mercurio?
Una empresa minera en el corregimiento de La Loma (Cesar) vierte mercurio al río Magdalena a razón de 0.5 kg/h. El río tiene un caudal promedio de 5 000 m³/s y una velocidad de 2 m/s. Usando un modelo de advección-difusión en estado estacionario, determina la concentración de mercurio a 50 km aguas abajo del punto de vertimiento. El coeficiente de dispersión longitudinal es m²/s. Considera que el mercurio no se degrada y que el río tiene sección rectangular de 100 m de ancho y 5 m de profundidad.
Datos
| Q_vertido | caudal de vertido de mercurio | 0.5 | kg/h |
| Q_rio | caudal del río Magdalena | 5000 | m³/s |
| v | velocidad del río | 2 | m/s |
| K | coeficiente de dispersión longitudinal | 100 | m²/s |
| x | distancia aguas abajo | 50 | km |
| ancho | ancho del río | 100 | m |
| profundidad | profundidad del río | 5 | m |
Se busca
- C — concentración de mercurio a 50 km (µg/L)
Pistas progresivas
Pista 1
La concentración en estado estacionario para un vertido continuo se calcula con . Convierte primero el caudal de vertido a kg/s.
Pista 2
La concentración inicial () es . Luego aplica la atenuación exponencial.
Pista 3
1 µg = 10^{-9} kg. Convierte el resultado final a microgramos por litro para compararlo con estándares ambientales.
Solución completa
- Conversión de unidades del vertido — Convertimos el caudal de mercurio de kg/h a kg/s para que sea consistente con las unidades del río.
- Cálculo de la concentración inicial — La concentración justo después del vertido (sin dispersión) es el cociente entre el caudal de vertido y el caudal del río.
- Cálculo del término exponencial — Evaluamos el exponente de la ecuación de dispersión. Primero convertimos 50 km a metros: 50 000 m.
- Cálculo final de la concentración — Multiplicamos la concentración inicial por el término exponencial para obtener la concentración a 50 km.
- Conversión a µg/L — Convertimos la concentración de kg/m³ a µg/L. Recuerda que 1 kg/m³ = 10^9 µg/L.
→ La concentración de mercurio a 50 km aguas abajo es aproximadamente µg/L, un valor extremadamente bajo que indica una gran dilución en el río Magdalena.
Ventilación en la mina de carbón de La Guajira: ¿Sobresale el metano?
En la mina de carbón de El Cerrejón (La Guajira), se libera metano (CH₄) a un ritmo de 10 m³/min a través de una grieta en la pared. La mina tiene un volumen de 50 000 m³ y se ventila con un flujo de aire fresco de 200 m³/min. Calcula la concentración de metano en la mina después de 30 minutos si inicialmente no había metano. ¿Se supera el límite seguro de 1% en volumen?
Datos
| Q_metano | caudal de metano liberado | 10 | m³/min |
| V_mina | volumen de la mina | 50000 | m³ |
| Q_aire | caudal de ventilación | 200 | m³/min |
| t | tiempo de acumulación | 30 | min |
Se busca
- C_metano — concentración de metano después de 30 min (%)
Pistas progresivas
Pista 1
Modela el problema como un balance de masas. El metano se acumula mientras el aire fresco lo diluye.
Pista 2
La concentración final es .
Pista 3
El volumen total de metano es el caudal por el tiempo. El volumen total de aire es el volumen de la mina más el aire ventilado.
Solución completa
- Cálculo del volumen de metano acumulado — Multiplicamos el caudal de metano por el tiempo para obtener el volumen total liberado.
- Cálculo del volumen total de aire en la mina — El volumen total de aire es el volumen de la mina más el volumen de aire ventilado en 30 minutos.
- Cálculo de la concentración en porcentaje — La concentración es el cociente entre el volumen de metano y el volumen total de aire, multiplicado por 100.
→ La concentración de metano después de 30 minutos es 0.57%, un valor seguro por debajo del límite del 1%.
Simulación de tsunamis en el Pacífico colombiano: ¿Qué altura alcanza la ola en Buenaventura?
Un terremoto de magnitud 8.0 ocurre frente a la costa de Tumaco (Nariño). Usando un modelo simplificado de propagación de tsunamis, estima la altura de la ola cuando llega a Buenaventura (150 km de distancia). Considera que la profundidad promedio del océano en esta ruta es de 3 000 m, la velocidad inicial de la ola es donde m/s², y la altura inicial de la ola en el epicentro es de 2 m. Usa la ley de Green para la atenuación: donde es la distancia recorrida.
Datos
| H_0 | altura inicial de la ola | 2 | m |
| x_0 | distancia de referencia | 1 | km |
| x | distancia a Buenaventura | 150 | km |
| h | profundidad promedio del océano | 3000 | m |
| g | aceleración gravitacional | 9.81 | m/s² |
Se busca
- H — altura de la ola en Buenaventura (m)
Pistas progresivas
Pista 1
Primero calcula la velocidad inicial de la ola usando . Luego usa la ley de Green para la atenuación con la distancia.
Pista 2
La ley de Green es una aproximación para tsunamis en aguas profundas: . Aquí km como referencia.
Pista 3
No necesitas calcular el tiempo de propagación, solo la atenuación de la altura.
Solución completa
- Cálculo de la velocidad inicial de la ola — Usamos la fórmula de velocidad de onda en aguas profundas.
- Cálculo del factor de atenuación — Aplicamos la ley de Green para calcular la altura de la ola después de recorrer 150 km.
- Evaluación numérica — Calculamos el valor del exponente y luego la altura final.
→ La altura estimada de la ola en Buenaventura es aproximadamente 0.63 m, lo que indica una atenuación significativa pero aún peligrosa para zonas costeras bajas.
Diseño de canoas en el Amazonas: ¿Qué forma resiste mejor la corriente?
Dos comunidades indígenas en el Amazonas están comparando el diseño de sus canoas tradicionales. La canoa A tiene un área frontal de 0.4 m² y un coeficiente de arrastre , mientras que la canoa B tiene un área frontal de 0.35 m² y . Si la velocidad de la corriente es de 3 m/s y la densidad del agua es 1 000 kg/m³, ¿qué canoa requiere menos fuerza para avanzar contra la corriente? Calcula la diferencia en energía necesaria para recorrer 1 km.
Datos
| A_A | área frontal de la canoa A | 0.4 | m² |
| C_dA | coeficiente de arrastre de la canoa A | 0.9 | |
| A_B | área frontal de la canoa B | 0.35 | m² |
| C_dB | coeficiente de arrastre de la canoa B | 0.7 | |
| v | velocidad de la corriente | 3 | m/s |
| rho | densidad del agua | 1000 | kg/m³ |
| distancia | distancia a recorrer | 1 | km |
Se busca
- F_A — fuerza requerida para la canoa A (N)
- F_B — fuerza requerida para la canoa B (N)
- E_ahorro — energía ahorrada con la canoa B (J)
Pistas progresivas
Pista 1
La fuerza de arrastre se calcula con . Calcula primero la fuerza para cada canoa.
Pista 2
La energía necesaria es el trabajo realizado: . La distancia debe estar en metros.
Pista 3
La diferencia en energía te dirá cuánto más eficiente es una canoa respecto a la otra.
Solución completa
- Cálculo de la fuerza para la canoa A — Aplicamos la fórmula de arrastre con los valores de la canoa A.
- Cálculo de la fuerza para la canoa B — Aplicamos la fórmula de arrastre con los valores de la canoa B.
- Cálculo del trabajo para cada canoa — Convertimos la distancia a metros (1 km = 1 000 m) y calculamos el trabajo para cada canoa.
- Cálculo de la energía ahorrada — Restamos los trabajos para obtener la diferencia en energía necesaria.
→ La canoa A requiere 1 620 N de fuerza y la canoa B requiere 1 102.5 N. La energía ahorrada al usar la canoa B para recorrer 1 km es 517 500 J.
Contaminación por plaguicidas en Boyacá: ¿Cuánto tarda en llegar al embalse de La Esmeralda?
En una finca de papa en Tunja (Boyacá), se aplica un plaguicida que se filtra al suelo con una velocidad de infiltración de 0.5 cm/h. El embalse de La Esmeralda está a 20 km de distancia y el acuífero tiene una pendiente promedio de 0.1%. Usando un modelo simplificado de transporte de contaminantes, estima el tiempo que tarda el plaguicida en llegar al embalse. Considera que el contaminante viaja con la velocidad del agua subterránea, que es proporcional a la pendiente y la permeabilidad (supón una velocidad efectiva de 0.2 m/h).
Datos
| v_infiltración | velocidad de infiltración | 0.5 | cm/h |
| distancia | distancia al embalse | 20 | km |
| pendiente | pendiente del acuífero | 0.1 | % |
| v_efectiva | velocidad efectiva del agua subterránea | 0.2 | m/h |
Se busca
- t — tiempo de llegada al embalse (días)
Pistas progresivas
Pista 1
Convierte primero la distancia a metros: 20 km = 20 000 m.
Pista 2
Usa la fórmula donde es la distancia y es la velocidad efectiva del agua subterránea.
Pista 3
Convierte el tiempo de horas a días dividiendo entre 24.
Solución completa
- Conversión de unidades — Convertimos la distancia de kilómetros a metros para que sea consistente con la velocidad en m/h.
- Cálculo del tiempo en horas — Aplicamos la fórmula del tiempo usando la velocidad efectiva del agua subterránea.
- Conversión a días — Convertimos el tiempo de horas a días para obtener un valor más comprensible.
- Conversión a años — Para mayor claridad, convertimos los días a años.
→ El plaguicida tarda aproximadamente 4 167 días (11.4 años) en llegar al embalse de La Esmeralda desde la finca en Tunja.
Optimización de la forma de un bus en Cali: ¿Reducir el consumo de gasolina?
Un bus de transporte masivo en Cali tiene una resistencia al avance de 8 000 N a 70 km/h. Si se modifica la forma del bus para reducir el coeficiente de arrastre de 0.9 a 0.5, calcula la nueva fuerza de resistencia. ¿Cuánto combustible se ahorra en un día de operación (16 horas) si el bus recorre 200 km diarios? Considera que el motor tiene una eficiencia del 25%, el poder calorífico de la gasolina es 32 MJ/L y su costo es 11 000 COP/L. Usa la ecuación de arrastre donde kg/m³ y m².
Datos
| F_actual | fuerza de resistencia actual | 8000 | N |
| C_d_actual | coeficiente de arrastre actual | 0.9 | |
| C_d_nuevo | coeficiente de arrastre nuevo | 0.5 | |
| v | velocidad del bus | 70 | km/h |
| rho | densidad del aire | 1.225 | kg/m³ |
| A | área frontal del bus | 8 | m² |
| horas_diarias | horas de operación diarias | 16 | h |
| distancia_diaria | distancia diaria recorrida | 200 | km |
| eficiencia | eficiencia del motor | 25 | % |
| poder_calorífico | poder calorífico de la gasolina | 32 | MJ/L |
| costo_combustible | costo de la gasolina | 11000 | COP/L |
Se busca
- F_nuevo — fuerza de resistencia nueva (N)
- ahorro_combustible — ahorro de combustible diario (L)
- ahorro_económico — ahorro económico diario (COP)
Pistas progresivas
Pista 1
Primero calcula la fuerza nueva usando la relación entre coeficientes de arrastre: .
Pista 2
Calcula el trabajo ahorrado por día: , donde es la distancia diaria en metros (200 km = 200 000 m).
Pista 3
La energía necesaria del combustible se obtiene dividiendo el trabajo por la eficiencia del motor. Luego divide por el poder calorífico para obtener litros ahorrados.
Pista 4
Multiplica el volumen ahorrado por el costo del combustible para obtener el ahorro económico.
Solución completa
- Cálculo de la fuerza nueva — Aplicamos la relación de proporcionalidad entre la fuerza y el coeficiente de arrastre.
- Cálculo del trabajo ahorrado por día — Convertimos la distancia a metros y calculamos el trabajo ahorrado.
- Cálculo de la energía necesaria del combustible — Dividimos el trabajo ahorrado por la eficiencia del motor (expresada como decimal).
- Cálculo del volumen de combustible ahorrado — Dividimos la energía por el poder calorífico del combustible para obtener litros.
- Cálculo del ahorro económico diario — Multiplicamos el volumen ahorrado por el costo del combustible.
→ La nueva fuerza de resistencia es 4 444 N. El ahorro de combustible diario es aproximadamente 1.93 litros, lo que equivale a 21 230 COP por día.
Simulación de viento en Cartagena: ¿Qué presión ejerce el viento en el Castillo San Felipe?
Durante el verano, en Cartagena sopla un viento constante de 50 km/h desde el norte. Calcula la presión dinámica que ejerce este viento sobre las paredes del Castillo San Felipe, que tiene una superficie frontal de 200 m². Usa la fórmula de presión dinámica donde kg/m³ es la densidad del aire. ¿Qué fuerza total ejerce el viento sobre el castillo?
Datos
| v | velocidad del viento | 50 | km/h |
| rho | densidad del aire | 1.225 | kg/m³ |
| A | superficie frontal del castillo | 200 | m² |
Se busca
- P — presión dinámica del viento (Pa)
- F — fuerza total ejercida por el viento (N)
Pistas progresivas
Pista 1
Convierte primero la velocidad del viento de km/h a m/s: m/s.
Pista 2
Calcula la presión dinámica usando .
Pista 3
La fuerza total es el producto de la presión por el área frontal: .
Solución completa
- Conversión de velocidad — Convertimos la velocidad del viento a unidades del SI (m/s).
- Cálculo de la presión dinámica — Aplicamos la fórmula de presión dinámica con los valores convertidos.
- Cálculo de la fuerza total — Multiplicamos la presión por el área frontal del castillo para obtener la fuerza total.
→ La presión dinámica del viento es aproximadamente 117 Pa y la fuerza total ejercida sobre el castillo es 23 400 N.