Skip to content
Este artículo tiene fines educativos. Te animamos a verificar con fuentes oficiales.

¿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?

facileapplication

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

ddiámetro de la tubería30cm
vvelocidad del agua5m/s
costo_aguacosto del agua3500COP/m³
n_hogaresnúmero de hogares afectados500

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 Q=A×v, donde A es el área transversal de la tubería.

Pista 2

El área de un círculo es A=πr2. 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
  1. Cálculo del área transversal — Primero calculamos el área de la tubería usando el diámetro dado. Convertimos 30 cm a metros: r=0.15 m. Luego aplicamos la fórmula del área de un círculo.
    A=πr2=π(0.15)2
  2. Cálculo del caudal volumétrico — Ahora aplicamos la ecuación de continuidad para fluidos incompresibles: el caudal volumétrico Q es el producto del área por la velocidad.
    Q=A×v
  3. 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.
    Vpe´rdido=Q×3600 s
  4. 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.
    Pérdida=Vpe´rdido×3500 COP/m3

Q=0.177 m3/sy2 233 200 COP/hora

→ 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?

moyenmodeling

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 D=0.5 m²/s. Considera que la emisión es continua y uniforme.

Datos

haltura de la chimenea20m
vvelocidad del viento1m/s
Ttemperatura ambiente20°C
Dcoeficiente de difusión turbulenta0.5m²/s
xdistancia horizontal100m

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: C(x)=Qv4πDxexp(vx24Dx) donde Q es la tasa de emisión (supón Q=1 kg/s para este ejercicio).

Pista 2

Simplifica la expresión exponencial: exp(vx4D).

Pista 3

Asegúrate de que las unidades sean consistentes: la concentración resultante estará en kg/m³.

Solución completa
  1. 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.
    C(x)=Qv4πDxexp(vx4D)
  2. Sustitución de valores — Reemplazamos los valores numéricos en la ecuación. Usamos Q=1 kg/s como tasa de emisión típica para una chimenea industrial pequeña.
    C(100)=11×4π×0.5×100exp(1×1004×0.5)
  3. Cálculo del denominador — Calculamos primero el término dentro de la raíz cuadrada y luego la raíz misma.
    4π×0.5×100=200π
  4. Cálculo del exponente — Evaluamos la parte exponencial de la ecuación.
    exp(50)1.93×1022
  5. Resultado final — Combinamos todos los términos para obtener la concentración final.
    C(100)=1200π×1.93×1022

C4.5×1024 kg/m3

→ La concentración del contaminante a 100 m de distancia es aproximadamente 4.5×1024 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?

moyenoptimization

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 Cd 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_actualfuerza de resistencia actual5000N
C_d_actualcoeficiente de arrastre actual0.8
C_d_nuevocoeficiente de arrastre nuevo0.6
vvelocidad del vagón60km/h
distanciadistancia viaje redondo18km
eficienciaeficiencia del motor30%
costo_combustiblecosto del combustible10000COP/L
poder_caloríficopoder calorífico del combustible35MJ/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 F=12ρv2CdA, donde ρ es la densidad del aire (1.225 kg/m³) y A es el área frontal (supón A=10 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 W=F×d. La energía necesaria proviene del combustible, considerando la eficiencia del motor.

Solución completa
  1. 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.
    Fnuevo=Factual×Cd_nuevoCd_actual
  2. 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.
    Wahorro=(FactualFnuevo)×18 000 m
  3. 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).
    Ecombustible=Wahorro0.3
  4. Cálculo del ahorro de combustible — Dividimos la energía necesaria por el poder calorífico del combustible para obtener los litros ahorrados.
    Vahorro=Ecombustible35 MJ/L
  5. Cálculo del ahorro económico — Multiplicamos el volumen ahorrado por el costo del combustible para obtener el valor en pesos.
    Ahorro=Vahorro×10 000 COP/L

Fnuevo=3 750 Ny760 COP/ viaje

→ 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?

difficilemodeling

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 K=100 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_vertidocaudal de vertido de mercurio0.5kg/h
Q_riocaudal del río Magdalena5000m³/s
vvelocidad del río2m/s
Kcoeficiente de dispersión longitudinal100m²/s
xdistancia aguas abajo50km
anchoancho del río100m
profundidadprofundidad del río5m

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 C(x)=QvertidoQrioexp(vx2K). Convierte primero el caudal de vertido a kg/s.

Pista 2

La concentración inicial (C0) es Qvertido/Qrio. 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
  1. 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.
    Qvertido=0.5 kg/h=0.53600 kg/s
  2. 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.
    C0=QvertidoQrio
  3. 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.
    exp(vx2K)=exp(2×50 0002×100)=exp(500)
  4. 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.
    C(50 000)=C0×exp(500)
  5. Conversión a µg/L — Convertimos la concentración de kg/m³ a µg/L. Recuerda que 1 kg/m³ = 10^9 µg/L.
    C(μg/L)=C(50 000)×109

C2.78×1012 μg/L

→ La concentración de mercurio a 50 km aguas abajo es aproximadamente 2.78×1012 µ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?

facileapplication

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_metanocaudal de metano liberado10m³/min
V_minavolumen de la mina50000
Q_airecaudal de ventilación200m³/min
ttiempo de acumulación30min

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 C=Volumen total de metano acumuladoVolumen total de aire en la mina×100.

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
  1. Cálculo del volumen de metano acumulado — Multiplicamos el caudal de metano por el tiempo para obtener el volumen total liberado.
    Vmetano=Qmetano×t
  2. 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.
    Vtotal=Vmina+Qaire×t
  3. 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.
    C=VmetanoVtotal×100

C=0.57 %

→ 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?

difficilemodeling

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 v=gh donde g=9.81 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: H=H0(x0x)1/4 donde x es la distancia recorrida.

Datos

H_0altura inicial de la ola2m
x_0distancia de referencia1km
xdistancia a Buenaventura150km
hprofundidad promedio del océano3000m
gaceleración gravitacional9.81m/s²

Se busca

  • H — altura de la ola en Buenaventura (m)

Pistas progresivas

Pista 1

Primero calcula la velocidad inicial de la ola usando v=gh. 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: H=H0(x0x)1/4. Aquí x0=1 km como referencia.

Pista 3

No necesitas calcular el tiempo de propagación, solo la atenuación de la altura.

Solución completa
  1. Cálculo de la velocidad inicial de la ola — Usamos la fórmula de velocidad de onda en aguas profundas.
    v=gh=9.81×3000
  2. 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.
    H=H0(1150)1/4
  3. Evaluación numérica — Calculamos el valor del exponente y luego la altura final.
    H=2×(0.00667)0.25

H0.63 m

→ 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?

moyenmodeling

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 CdA=0.9, mientras que la canoa B tiene un área frontal de 0.35 m² y CdB=0.7. 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 A0.4
C_dAcoeficiente de arrastre de la canoa A0.9
A_Bárea frontal de la canoa B0.35
C_dBcoeficiente de arrastre de la canoa B0.7
vvelocidad de la corriente3m/s
rhodensidad del agua1000kg/m³
distanciadistancia a recorrer1km

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 F=12ρv2CdA. Calcula primero la fuerza para cada canoa.

Pista 2

La energía necesaria es el trabajo realizado: W=F×d. 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
  1. Cálculo de la fuerza para la canoa A — Aplicamos la fórmula de arrastre con los valores de la canoa A.
    FA=12×1000×32×0.9×0.4
  2. Cálculo de la fuerza para la canoa B — Aplicamos la fórmula de arrastre con los valores de la canoa B.
    FB=12×1000×32×0.7×0.35
  3. Cálculo del trabajo para cada canoa — Convertimos la distancia a metros (1 km = 1 000 m) y calculamos el trabajo para cada canoa.
    WA=FA×1000yWB=FB×1000
  4. Cálculo de la energía ahorrada — Restamos los trabajos para obtener la diferencia en energía necesaria.
    Eahorro=WAWB

FA=1 620 N, FB=1 102.5 NyEahorro=517 500 J

→ 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?

moyenanalysis

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ónvelocidad de infiltración0.5cm/h
distanciadistancia al embalse20km
pendientependiente del acuífero0.1%
v_efectivavelocidad efectiva del agua subterránea0.2m/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 t=dv donde d es la distancia y v es la velocidad efectiva del agua subterránea.

Pista 3

Convierte el tiempo de horas a días dividiendo entre 24.

Solución completa
  1. Conversión de unidades — Convertimos la distancia de kilómetros a metros para que sea consistente con la velocidad en m/h.
    d=20 000 m
  2. Cálculo del tiempo en horas — Aplicamos la fórmula del tiempo usando la velocidad efectiva del agua subterránea.
    t=20 0000.2=100 000 horas
  3. Conversión a días — Convertimos el tiempo de horas a días para obtener un valor más comprensible.
    t=100 000244 167 días
  4. Conversión a años — Para mayor claridad, convertimos los días a años.
    t=4 16736511.4 años

t4 167 días(11.4 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?

difficileoptimization

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 F=12ρv2CdA donde ρ=1.225 kg/m³ y A=8 m².

Datos

F_actualfuerza de resistencia actual8000N
C_d_actualcoeficiente de arrastre actual0.9
C_d_nuevocoeficiente de arrastre nuevo0.5
vvelocidad del bus70km/h
rhodensidad del aire1.225kg/m³
Aárea frontal del bus8
horas_diariashoras de operación diarias16h
distancia_diariadistancia diaria recorrida200km
eficienciaeficiencia del motor25%
poder_caloríficopoder calorífico de la gasolina32MJ/L
costo_combustiblecosto de la gasolina11000COP/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: Fnuevo=Factual×Cd_nuevoCd_actual.

Pista 2

Calcula el trabajo ahorrado por día: Wahorro=(FactualFnuevo)×d, donde d 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
  1. Cálculo de la fuerza nueva — Aplicamos la relación de proporcionalidad entre la fuerza y el coeficiente de arrastre.
    Fnuevo=8 000×0.50.9
  2. Cálculo del trabajo ahorrado por día — Convertimos la distancia a metros y calculamos el trabajo ahorrado.
    Wahorro=(8 000Fnuevo)×200 000
  3. Cálculo de la energía necesaria del combustible — Dividimos el trabajo ahorrado por la eficiencia del motor (expresada como decimal).
    Ecombustible=Wahorro0.25
  4. Cálculo del volumen de combustible ahorrado — Dividimos la energía por el poder calorífico del combustible para obtener litros.
    Vahorro=Ecombustible32 MJ/L
  5. Cálculo del ahorro económico diario — Multiplicamos el volumen ahorrado por el costo del combustible.
    Ahorro=Vahorro×11 000

Fnuevo=4 444 N, Vahorro=1.93 L/díay21 230 COP/día

→ 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?

facileapplication

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 P=12ρv2 donde ρ=1.225 kg/m³ es la densidad del aire. ¿Qué fuerza total ejerce el viento sobre el castillo?

Datos

vvelocidad del viento50km/h
rhodensidad del aire1.225kg/m³
Asuperficie frontal del castillo200

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: 50 km/h=50×10003600 m/s.

Pista 2

Calcula la presión dinámica usando P=12ρv2.

Pista 3

La fuerza total es el producto de la presión por el área frontal: F=P×A.

Solución completa
  1. Conversión de velocidad — Convertimos la velocidad del viento a unidades del SI (m/s).
    v=50×10003600=13.89 m/s
  2. Cálculo de la presión dinámica — Aplicamos la fórmula de presión dinámica con los valores convertidos.
    P=12×1.225×(13.89)2
  3. Cálculo de la fuerza total — Multiplicamos la presión por el área frontal del castillo para obtener la fuerza total.
    F=P×200

P=117 PayF=23 400 N

→ La presión dinámica del viento es aproximadamente 117 Pa y la fuerza total ejercida sobre el castillo es 23 400 N.

Fuentes

  1. en.wikipedia.org
  2. ui.adsabs.harvard.edu
  3. doi.org
  4. www.lanl.gov
  5. zenodo.org
  6. search.worldcat.org
  7. www.mw.tum.de
  8. hdl.handle.net