Modelos Continuos · Tema 12

Leyes naturales

La física como fundamento del modelo: balances que escriben la ecuación por nosotros.

01 · Punto de partida

Cuando la regla ya existe

La taza no «elige» cómo enfriarse: el flujo de calor va del cuerpo caliente al ambiente según la diferencia de temperaturas. El tanque no negocia: lo que entra menos lo que sale se acumula. Estas leyes naturales —conservación de masa y energía, Newton, Fourier— escriben f antes de que abramos el editor.

Modelar desde leyes es partir con ventaja: la estructura es conocida y solo quedan pocos parámetros por calibrar.

  • ¿Qué se conserva: masa, energía, población?
  • ¿Qué flujo cruza la frontera y según qué diferencia?
  • ¿Qué parámetro resume el material o el ambiente (k)?
  • ¿Dónde termina el sistema si nadie lo toca?

02 · Definición

¿Qué es un modelo basado en leyes?

Tres leyes cubren casi todo el curso:

◈

Balance de masa

Lo que entra menos lo que sale se acumula: dN/dt = q_in − q_out. Ejemplo: tanques, SIR, poblaciones.

〜

Flujo por diferencia

El flujo es k por la diferencia: q = k·(T − T_amb). Ejemplo: enfriamiento, calor, difusión.

●

Newton (movimiento)

Fuerza = masa × aceleración: mx'' = −kx − cv + F. Ejemplo: resorte, caída con aire.

La ley fija la forma (lineal en la diferencia); los datos solo ajustan k y las condiciones. Por eso estos modelos extrapolan mejor que los puramente ajustados.

Cada ley trae su EDO bajo el brazo.
LeyEDOEquilibrio
Balance de tanquedN/dt = q − k·NN* = q/k.
Enfriamiento de NewtondT/dt = −k·(T − T_amb)T* = T_amb.
Resorte + frenomx'' = −kx − cvReposo en x = 0.
DecaimientodM/dt = −λ·M0 (se agota).

03 · Anatomía

Balance, flujo y parámetro k

Todo modelo de leyes se arma igual:

Frontera y X
Qué se acumula dentro: litros, joules, individuos. Sin frontera no hay balance.
Flujos con signo
Entradas positivas, salidas negativas, cada una con sus unidades por tiempo.
k del material
Conductividad, desagüe, transferencia: resume geometría y material en un número.
Ambiente impuesto
T_amb, q_in, F: el entorno que el sistema no decide pero que lo arrastra.

Desde leyes

Estructura conocida, pocos parámetros

Enfriamiento: forma −k·(T − T_amb) garantizada; solo k se calibra. Predice fuera del rango medido.

Desde datos (Tema 13)

Forma propuesta, varios parámetros

Polinomio o exponencial ajustada: flexible dentro de los datos, frágil fuera de ellos.

04 · Leyes en modelos clásicos

La misma física, cinco disfraces

Cada modelo del curso es un balance con su diferencia:

Ley + frontera→Flujo k·Δ→EDO
  1. 1
    Tanque.

    Balance de masa con salida k·N: dN/dt = q − k·N.

  2. 2
    Taza.

    Energía con flujo k·(T − T_amb): dT/dt = −k·(T − T_amb).

  3. 3
    Regulador.

    Dos cuerpos que intercambian k·(T₁ − T₂): convergen al medio.

  4. 4
    Caída con aire.

    Newton con freno −c·v: velocidad límite mg/c.

  5. 5
    Mezcla salina.

    Balance de sal con caudales: concentración que converge.

05 · Representación en Python

Del balance a f(estado, t)

En Python, la ley es una función de flujos con nombre: cada término conserva sus unidades y su signo. El balance los suma.

Python en tu navegador. Estos bloques se ejecutan localmente con Pyodide, sin instalar nada. Cambiá t_amb y observá que el equilibrio se muda con él.

def flujo_calor(t, t_amb, k):
    return k * (t - t_amb)


def f_enfriamiento(t, t_amb=20, k=0.1):
    return -flujo_calor(t, t_amb, k)


for t in [80, 50, 30]:
    print(t, "→", round(f_enfriamiento(t), 2))

−6, −3, −1 °C/min: el flujo muere con la diferencia. La ley está en el código.

Solución con la ley adentro

La exacta hereda la forma exponencial hacia el ambiente:

import math


def temperatura(t0, t_amb, k, t):
    return t_amb + (t0 - t_amb) * math.exp(-k * t)


print(round(temperatura(80, 20, 0.1, 10), 1))

06 · Exploración

Laboratorio: la taza que obedece a Newton

Enfriamiento dT/dt = −k·(T − T_amb) desde T₀ = 80 °C. La punteada es T_amb: el atractor impuesto por el ambiente. Movete k y T_amb y observá que la forma es siempre la misma exponencial, solo cambian destino y velocidad.

EXPERIMENTO 12

Convergencia al ambiente

dT/dt = −k · (T − T_amb) · T(0) = 80

Los resultados numéricos aparecen debajo.
T final (t=60)20,15
Flujo final−0,01
Vida media ln2/k6,93
LecturaLlegó al ambiente

Con k 0,1 y ambiente 20, en t = 60 queda 20,15.

Todas las curvas son la misma exponencial trasladada: la ley fija la forma, el ambiente fija el destino, k fija la velocidad.

Preguntas para explorar

  1. Con T_amb = 0 y T(0) = 80, ¿a dónde converge? ¿Cambia la forma?
  2. Duplicá k. ¿La vida media se duplica o se parte a la mitad?
  3. Si T(0) = 10 con ambiente 20, ¿la curva sube o baja? ¿La ley sigue valiendo?
Ver respuestas sugeridas
  1. A 0: la forma es idéntica, solo baja el destino. La ley no distingue calentar de enfriar.
  2. Se parte a la mitad (ln2/k): más k, convergencia más rápida.
  3. Sube hacia 20: el flujo cambia de signo con la diferencia. La misma f calienta y enfría.

07 · Comprensión

Confusiones frecuentes

«k es universal»

k resume material, geometría y ventilación de tu taza. Otra taza, otro k: se calibra, no se recita.

«El ambiente es condición inicial»

No: T_amb es entrada impuesta que dura todo el experimento. La condición inicial es T(0).

«Si ajusta, es ley»

Ajustar no es derivar: un polinomio puede ajustar y no conservar nada. La ley trae balance y unidades; el ajuste solo trae números.

«Más física siempre gana»

Solo en su rango: fuera de él, una forma empírica bien calibrada predice mejor que una ley forzada.

08 · Práctica guiada

Ejercicios con Python

Ejercicio 1: flujo con signo

Con k = 0,1 y T_amb = 20, calculá el flujo k·(T − T_amb) y dT/dt en T = 80 y T = 10.

def flujo(t, t_amb=20, k=0.1):
    return k * (t - t_amb)


for t in [80, 10]:
    print(t, "flujo:", round(flujo(t), 2), "dT/dt:", round(-flujo(t), 2))
Ver solución

En 80: flujo +6, dT/dt −6. En 10: flujo −1, dT/dt +1: el signo se da vuelta con la diferencia.

Ejercicio 2: vida media

¿En cuánto se reduce a la mitad la diferencia T − T_amb con k = 0,1? ¿Y con k = 0,05?

import math


def vida_media(k):
    return round(math.log(2) / k, 2)


print(vida_media(0.1), vida_media(0.05))
Ver solución

6,93 y 13,86 min: T₀ y T_amb no aparecen. La velocidad la manda k.

Ejercicio 3: balance completo

Agregá una fuente (microondas de 2 °C/min) al enfriamiento y hallá el nuevo equilibrio. Probalo con k = 0,1 y T_amb = 20.

Ver una posible respuesta
def equilibrio_fuente(potencia, k=0.1, t_amb=20):
    return t_amb + potencia / k


print(equilibrio_fuente(2))

Da 40 °C: la fuente traslada el equilibrio en potencia/k. Balance con dos flujos, un solo cero.

09 · Síntesis

Ideas para recordar

  • Ley = balance + flujo por diferencia: la forma de f viene de la física.
  • k resume material y geometría; el ambiente impone el destino (T_amb, q).
  • La misma exponencial sirve para enfriar y calentar: el signo lo pone la diferencia.
  • En Python: función por flujo + función por balance; exacta T_amb + (T₀ − T_amb)e^(−kt).
  • Las leyes extrapolan mejor, pero solo en su rango: fuera de él, otra f.

En el próximo tema construiremos modelos al revés: desde las observaciones hacia la ecuación.