Por qué las unidades forman parte del modelo

Con este capítulo comienza la segunda parte del tomo, dedicada a las herramientas matemáticas que usaremos a lo largo del libro. La primera parece elemental y consiste en distinguir variables, unidades y escalas. Muchos errores de modelamiento aparecen, sin embargo, antes de cualquier cálculo difícil. Una ecuación puede manipularse algebraicamente y carecer de sentido porque suma magnitudes incompatibles, atribuye unidades incorrectas a un parámetro o compara sistemas mediante números que dependen de una elección arbitraria de escala.

Consideremos una sala que recibe pacientes y los da de alta a lo largo del día. Sea N(t)N(t) el número de personas presentes en el instante tt. Si las admisiones ocurren a una tasa aa y las altas a una tasa proporcional al número de pacientes, una ecuación mínima es

dNdt=a−k N.\frac{dN}{dt} = a - k\,N.

Antes de resolverla ya puede extraerse información. La variable NN es un número de pacientes. Si tt se mide en días, dN/dtdN/dt tiene unidades de pacientes por día, de modo que aa también debe medirse en pacientes por día. Para que kNkN tenga esas mismas unidades, kk debe tener unidades de 1/día. Si alguien expresara kk en pacientes por día, la ecuación dejaría de ser dimensionalmente consistente, porque los términos que se restan representarían magnitudes de naturaleza diferente.

La consistencia dimensional no demuestra que el modelo sea verdadero, porque muchas ecuaciones falsas tienen unidades perfectamente coherentes. Su función es más modesta y muy útil. Descarta representaciones que no pueden ser correctas tal como están escritas y hace explícita la interpretación física, biológica u operacional de los parámetros. A lo largo del libro aplicaremos este control antes de interpretar cualquier resultado numérico.

Magnitud, dimensión y unidad

Conviene separar tres niveles. Una magnitud es una propiedad susceptible de comparación cuantitativa dentro de un esquema de medición. La dimensión expresa el tipo de magnitud en términos de dimensiones básicas, mientras que la unidad fija una convención concreta para asignarle números. La longitud tiene dimensión L y puede expresarse en metros, centímetros o kilómetros. El valor numérico cambia con la unidad y la dimensión se mantiene.

El Sistema Internacional de Unidades, que mantiene la Oficina Internacional de Pesas y Medidas (BIPM, 2019), define siete magnitudes básicas, que son la longitud, la masa, el tiempo, la corriente eléctrica, la temperatura termodinámica, la cantidad de sustancia y la intensidad luminosa. Muchas variables se construyen a partir de ellas. Una velocidad tiene dimensión L T−1\mathrm{L\,T^{-1}}, una aceleración L T−2\mathrm{L\,T^{-2}} y una constante de decaimiento exponencial T−1\mathrm{T^{-1}}. En otros campos las unidades son operacionales y no se reducen de manera útil a dimensiones físicas. Un conteo de pacientes, una probabilidad o una fracción de ocupación requieren igualmente una definición explícita, aunque sean adimensionales o se traten como conteos.

Que una cantidad sea adimensional no significa que carezca de interpretación. Una proporción, una probabilidad y el cociente entre dos longitudes son adimensionales y representan objetos conceptualmente distintos. El número 0,5 puede designar una probabilidad, una prevalencia relativa, una concentración normalizada o la mitad de una longitud de referencia, de modo que el valor numérico no determina qué representa la variable.

También hay que distinguir variables de parámetros. En el modelo dN/dt=a−kNdN/dt = a - kN, NN cambia durante la evolución y es una variable de estado. Los parámetros aa y kk se consideran constantes mientras se analiza una ejecución concreta. Nada impide construir después un modelo en que aa cambie con la hora o kk dependa de la dotación disponible. En ese caso dejarían de desempeñar el mismo papel matemático, porque la clasificación depende de la representación elegida.

Consistencia dimensional

Una ecuación físicamente interpretable debe respetar la homogeneidad dimensional. Los términos que se suman o restan deben compartir dimensiones y ambos lados de una igualdad deben representar magnitudes compatibles. La regla permite detectar errores sin conocer la solución.

Supongamos un modelo para una cantidad x(t)x(t) con dimensión de masa M,

dxdt=r x−b x2.\frac{dx}{dt} = r\,x - b\,x^2 .

El lado izquierdo tiene dimensión M T−1\mathrm{M\,T^{-1}}. El término rxrx debe tener la misma dimensión, de modo que rr tiene dimensión T−1\mathrm{T^{-1}}. Para que bx2bx^2 tenga dimensión M T−1\mathrm{M\,T^{-1}}, bb debe tener dimensión M−1 T−1\mathrm{M^{-1}\,T^{-1}}. Asignar a rr y bb las mismas unidades porque ambos aparecen como coeficientes sería un error. La dimensión de un parámetro depende de la magnitud a la que multiplica y de la exigencia de que todos los términos sean compatibles, con independencia de su posición algebraica.

El principio restringe también los argumentos de ciertas funciones. Expresiones como eze^{z}, sin⁡z\sin z o log⁡z\log z requieren un argumento adimensional. En e−kte^{-kt}, el producto ktkt debe ser adimensional, lo que exige que kk tenga dimensión T−1\mathrm{T^{-1}} cuando tt tiene dimensión T. El razonamiento reaparecerá al estudiar el decaimiento, la relajación, las oscilaciones y las distribuciones de probabilidad.

La consistencia dimensional sirve también para auditar el código. Si un programa recibe tiempos en horas y usa sin convertirla una tasa estimada por día, la sintaxis es válida y el computador devuelve números. Con una tasa de 0,2 por día interpretada como 0,2 por hora, el decaimiento simulado es 24 veces más rápido que el previsto. El error reside en el modelo implementado. Conviene acompañar cada parámetro con su unidad y mantener explícitas las conversiones, porque una simulación numéricamente estable puede ser científicamente incorrecta si mezcla escalas incompatibles.

Ecuaciones de balance

Muchos modelos comienzan con una idea de balance. Para una cantidad acumulada XX, el cambio temporal puede escribirse como

dXdt=entradas−salidas+produccioˊn−consumo.\frac{dX}{dt} = \text{entradas} - \text{salidas} + \text{producción} - \text{consumo}.

La expresión no especifica todavía ningún mecanismo. Es una contabilidad dinámica. La física la usa para la masa, la energía o la cantidad de movimiento, la biología para poblaciones y sustancias, la farmacocinética para cantidades de fármaco y los sistemas sanitarios para personas en espera o camas ocupadas (Sterman, 2000). La semejanza formal no implica que los mecanismos sean equivalentes, porque cada término debe derivarse o justificarse según el sistema real.

Consideremos un depósito que recibe material a una tasa constante qq y lo pierde a una tasa proporcional a la cantidad presente,

dXdt=q−k X.\frac{dX}{dt} = q - k\,X.

Si XX se mide en miligramos y tt en horas, qq tiene unidades de mg/h y kk de h−1\mathrm{h^{-1}}, de modo que la salida kXkX vuelve a expresarse en mg/h. El modelo supone mezcla homogénea y una pérdida proporcional al contenido. Son hipótesis mecanicistas que la ecuación de balance no impone.

El equilibrio X∗X^* se obtiene con dX/dt=0dX/dt = 0, es decir, q−kX∗=0q - kX^* = 0 y

X∗=qk.X^* = \frac{q}{k}.

Las unidades confirman el resultado, porque (mg/h)/(h−1)=mg(\mathrm{mg/h})/(\mathrm{h^{-1}}) = \mathrm{mg}. El cociente q/kq/k define una escala natural de cantidad para este modelo. Si qq se duplica con kk fijo, el equilibrio se duplica. Si kk se duplica con qq fijo, el equilibrio se reduce a la mitad. Sin resolver todavía la trayectoria, el análisis de unidades y de balances ya identificó relaciones importantes.

Escalas características

Una escala característica es una magnitud de referencia que representa un tamaño, un tiempo, una velocidad u otra cantidad pertinente para el comportamiento del modelo. No siempre hay una única elección y escogerla forma parte del análisis. Una buena escala permite comparar términos, reconocer qué parámetros dominan y expresar las variables en forma adimensional.

En dX/dt=q−kXdX/dt = q - kX aparece de inmediato una escala temporal. Como kk tiene dimensión T−1\mathrm{T^{-1}}, su inverso tiene dimensión de tiempo, T=1/kT = 1/k. Esta duración describe el ritmo con que el sistema responde. Si qq se interrumpe, la ecuación queda dX/dt=−kXdX/dt = -kX y la solución decae como e−kte^{-kt}. Después de un tiempo T=1/kT = 1/k la cantidad se ha multiplicado por e−1≈0,368e^{-1} \approx 0{,}368. Después de varios múltiplos de TT, el efecto de la condición inicial se vuelve pequeño.

La escala de estado puede tomarse como el equilibrio X∗=q/kX^* = q/k, una referencia natural para medir cuánto vale XX respecto de su nivel estacionario. La variable normalizada vale uno cuando X=X∗X = X^* y dos cuando X=2X∗X = 2X^*, lo que facilita comparar sistemas con valores absolutos muy diferentes.

Una escala característica no coincide necesariamente con un promedio observado ni con un parámetro medido directamente. Puede surgir de una combinación de parámetros, como q/kq/k. Su utilidad se juzga por la estructura del problema y por la pregunta que se intenta responder.

Adimensionalización paso a paso

La adimensionalización expresa las variables como cocientes respecto de escalas de referencia (Strogatz, 2015). No elimina la física ni la interpretación del modelo. Reorganiza la ecuación para separar lo que depende de unidades particulares de lo que determina su forma relativa.

Partimos de dX/dt=q−kXdX/dt = q - kX con q>0q > 0 y k>0k > 0, de modo que Xc=q/kX_c = q/k y Tc=1/kT_c = 1/k son escalas positivas. Si q=0q = 0, la escala q/kq/k se anula y hay que elegir otra escala de cantidad, por ejemplo una condición inicial positiva. Si k=0k = 0, desaparece la escala de relajación y se requiere otra referencia temporal. Definimos las variables adimensionales

y=XXc,τ=tTc=k t.y = \frac{X}{X_c}, \qquad \tau = \frac{t}{T_c} = k\,t .

Como X=Xc yX = X_c\,y y dτ/dt=kd\tau/dt = k, la regla de la cadena da

dXdt=Xc dydτ dτdt=Xc k dydτ.\frac{dX}{dt} = X_c\,\frac{dy}{d\tau}\,\frac{d\tau}{dt} = X_c\,k\,\frac{dy}{d\tau}.

Al sustituir en la ecuación original se obtiene Xck dy/dτ=q−kXc yX_c k\, dy/d\tau = q - k X_c\, y. Como Xck=qX_c k = q, basta dividir por qq para llegar a

dydτ=1−y.\frac{dy}{d\tau} = 1 - y .

Los parámetros qq y kk desaparecieron de la ecuación adimensional, aunque siguen importando en el sistema real, porque determinan cómo convertir yy y τ\tau de vuelta a unidades físicas. Lo que se mostró es que, dentro de este modelo, todos los sistemas con q>0q > 0 y k>0k > 0 comparten la misma evolución cuando el estado y el tiempo se expresan en sus propias escalas características. La solución adimensional es

y(τ)=1+[y(0)−1] e−τ.y(\tau) = 1 + \big[y(0) - 1\big]\,e^{-\tau}.

Un depósito con k=1 h−1k = 1\ \mathrm{h^{-1}} tiene Tc=1T_c = 1 h y otro con k=0,1 h−1k = 0{,}1\ \mathrm{h^{-1}} tiene Tc=10T_c = 10 h. Ambos alcanzan la misma fracción de aproximación al equilibrio después del mismo número de tiempos característicos, una comparación más informativa que decir que tardaron mucho o poco sin especificar la escala.

Dos sistemas distintos con una misma dinámica normalizada

El sistema A recibe q=20q = 20 mg/h y tiene k=0,5 h−1k = 0{,}5\ \mathrm{h^{-1}}, de modo que su equilibrio es X∗=40X^* = 40 mg y su escala temporal, 2 h. El sistema B recibe q=300q = 300 mg/día y tiene k=0,1 dıˊa−1k = 0{,}1\ \text{día}^{-1}, con equilibrio X∗=3000X^* = 3000 mg y escala temporal de 10 días. Los valores absolutos difieren en órdenes de magnitud y las unidades de tiempo no son las mismas.

Si ambos parten vacíos, la solución es X(t)=X∗(1−e−kt)X(t) = X^*\big(1 - e^{-kt}\big). Después de una escala temporal, 2 horas en A y 10 días en B, se tiene kt=1kt = 1 en los dos y X/X∗=1−e−1≈0,632X/X^* = 1 - e^{-1} \approx 0{,}632. Después de tres escalas temporales, kt=3kt = 3 y la fracción es 1−e−3≈0,9501 - e^{-3} \approx 0{,}950, lo que A logra en 6 horas y B en 30 días. La diferencia absoluta es enorme, mientras la dinámica relativa es idéntica, como muestra la figura 11.1.

Dos sistemas de acumulación con equilibrios y escalas temporales muy distintos.
Figura 11.1. Dos sistemas de acumulación con equilibrios y escalas temporales muy distintos. Los paneles a y b muestran sus trayectorias en unidades físicas, en horas y en días. En el panel c, expresadas como fracción del equilibrio y en múltiplos de la escala temporal, las dos trayectorias coinciden. Ampliar figura ↗

Supongamos que alguien compara solo los tiempos absolutos y concluye que B es menos resiliente porque tarda más en acercarse a su equilibrio. La afirmación exigiría justificar que el tiempo absoluto es el criterio pertinente y que ambos sistemas deben compararse en la misma unidad operacional. En el modelo normalizado ambos tienen la misma forma de relajación. La elección de la escala puede modificar, así, la interpretación de una comparación.

El resultado depende de la estructura lineal propuesta. Si la salida se saturara, dependiera de un umbral o cambiara con el tiempo, la adimensionalización podría dejar uno o varios parámetros adimensionales en la ecuación. Esos parámetros suelen ser especialmente informativos, porque expresan relaciones entre procesos que ningún cambio de unidades elimina.

Parámetros adimensionales y competencia entre procesos

En el ejemplo anterior, la elección de Xc=q/kX_c = q/k y Tc=1/kT_c = 1/k eliminó todos los parámetros, porque el modelo tiene una sola escala temporal independiente y una escala de estado construida con qq y kk. En sistemas más ricos, la adimensionalización deja parámetros sin dimensión que condensan la competencia entre procesos.

Supongamos que el sistema de relajación recibe una entrada periódica,

dXdt=q0[1+Acos⁡(ωt)]−k X,\frac{dX}{dt} = q_0\big[1 + A\cos(\omega t)\big] - k\,X,

donde AA es una amplitud relativa y ω\omega tiene dimensión T−1\mathrm{T^{-1}}. Para interpretar la entrada como un aporte no negativo, supondremos q0>0q_0 > 0 y ∣A∣≤1|A| \le 1. Con Tc=1/kT_c = 1/k y Xc=q0/kX_c = q_0/k, el argumento ωt\omega t se convierte en (ω/k) τ(\omega/k)\,\tau y la ecuación queda

dydτ=1+Acos⁡(Ωτ)−y,Ω=ωk.\frac{dy}{d\tau} = 1 + A\cos(\Omega\tau) - y, \qquad \Omega = \frac{\omega}{k}.

Quedan la amplitud relativa AA y el cociente Ω\Omega, además de la condición inicial normalizada. El cociente Ω\Omega compara la frecuencia del forzamiento con la tasa interna de relajación. Como ambas tienen dimensión de tiempo inversa, no cambia al pasar de horas a días.

Pasado el transitorio, la respuesta oscila alrededor de 1 con amplitud A/1+Ω2A/\sqrt{1 + \Omega^2} y un retraso de fase arctan⁡Ω\arctan\Omega, el mismo resultado del ejercicio 9.4. Si Ω\Omega es pequeño, el forzamiento varía lentamente respecto del ajuste interno y el estado lo sigue casi por completo. Si es grande, el sistema no alcanza a responder y filtra las oscilaciones. Con Ω=10\Omega = 10, solo llega al estado un 10 % de la amplitud del forzamiento. La figura 11.2 muestra ambos regímenes.

La misma lógica aparece en muchos campos. Un parámetro adimensional puede comparar transporte con reacción, inercia con viscosidad, producción con eliminación o una escala de aprendizaje con otra de cambio ambiental. Los nombres varían según la disciplina, aunque la idea es la misma. Cuando dos procesos tienen dimensiones compatibles, su razón indica cuál domina o si ambos compiten en una región intermedia.

Decir que un proceso es rápido carece de contenido mientras no se indique respecto de qué otro proceso o de qué escala. Una respuesta de una hora es lenta frente a una señal que cambia cada minuto y rápida frente a otra que cambia cada semana. Los cocientes adimensionales convierten esa comparación en parte explícita del modelo.

El principio puede usarse en aplicaciones clínicas u organizacionales aunque no todos los procesos admitan una medición precisa. Supongamos que una variable clínica responde a una intervención con una constante temporal de ocho semanas y que la intervención alterna cada semana entre dos intensidades, con un período de dos semanas. Entonces ω=2π/2=π\omega = 2\pi/2 = \pi por semana, k=1/8k = 1/8 por semana y Ω≈25\Omega \approx 25. Según el modelo, a la variable llega menos del 5 % de la amplitud de esa alternancia. Cuantificar la separación exige datos y una definición del estado, aunque la intuición sobre las escalas puede orientar la pregunta.

Relajación con una entrada periódica en variables adimensionales.
Figura 11.2. Relajación con una entrada periódica en variables adimensionales. En el panel a, la respuesta estacionaria con amplitud relativa A = 0,5 sigue casi por completo un forzamiento lento (Ω = 0,1) y filtra uno rápido (Ω = 10). En el panel b, la fracción de la amplitud que llega al estado, 1/√(1 + Ω²), disminuye con Ω. El cuadrado marca el ejemplo clínico del texto, con Ω ≈ 25. Ampliar figura ↗

El teorema π de Buckingham y los grupos adimensionales

El teorema π de Buckingham (1914) sistematiza la búsqueda de combinaciones adimensionales. Se escribe una matriz con los exponentes dimensionales de las nn cantidades consideradas. Si esa matriz tiene rango rr, existen n−rn - r grupos adimensionales independientes. El número rr es el rango de la matriz, que no coincide necesariamente con el número de dimensiones fundamentales disponibles. Una relación dimensionalmente homogénea puede reformularse como una relación entre esos grupos (Bridgman, 1922; Barenblatt, 2003). El teorema reduce el número de argumentos sin determinar la función que los relaciona. Las cantidades que ya son adimensionales, como una amplitud angular, cuentan también entre los argumentos del problema.

Un ejemplo clásico es el período PP de un péndulo simple en el régimen de oscilaciones pequeñas. Supongamos que depende de la longitud LL, de la aceleración gravitatoria gg y, quizá, de la masa mm. Las dimensiones son [P]=T[P] = \mathrm{T}, [L]=L[L] = \mathrm{L}, [g]=L T−2[g] = \mathrm{L\,T^{-2}} y [m]=M[m] = \mathrm{M}. Una combinación adimensional de la forma P La gb mcP\,L^{a}\,g^{b}\,m^{c} exige c=0c = 0 para eliminar la masa, a+b=0a + b = 0 para la longitud y 1−2b=01 - 2b = 0 para el tiempo. Se obtiene b=1/2b = 1/2, a=−1/2a = -1/2 y

Π=Pg/L.\Pi = P\sqrt{g/L}.

El análisis dimensional indica que Π\Pi debe ser una constante si esas son las únicas variables relevantes en el régimen considerado. No da el valor 2π2\pi, que requiere resolver la dinámica. Tampoco demuestra que la masa nunca importe en un péndulo real, porque en otras condiciones pueden intervenir la fricción, la amplitud, la geometría del cuerpo o la resistencia del aire.

Si no se fija el régimen de pequeña amplitud, hay que incluir el ángulo inicial θ0\theta_0 como argumento adimensional. La conclusión pasa a ser Pg/L=F(θ0)P\sqrt{g/L} = F(\theta_0), una función de la amplitud. El límite F(θ0)→2πF(\theta_0) \to 2\pi cuando θ0→0\theta_0 \to 0 proviene de la dinámica del péndulo ideal. Esa misma dinámica muestra que con una amplitud de 30° el período es cerca de un 1,7 % mayor que el de las oscilaciones pequeñas y con una de 90°, cerca de un 18 % mayor. La selección previa de variables es una hipótesis física indispensable para aplicar el teorema.

La utilidad epistemológica del procedimiento es precisa. Permite descartar dependencias dimensionalmente imposibles, identificar combinaciones de variables y diseñar comparaciones entre sistemas. No sustituye la teoría del mecanismo ni los datos empíricos, porque una ecuación dimensionalmente correcta puede ser físicamente equivocada.

Normalización, estandarización y cambio de unidades

En investigación se usan varios procedimientos que cambian los números de una variable y que no deben confundirse. Convertir 120 minutos en 2 horas es un cambio de unidad. Dividir XX por una escala XcX_c para obtener y=X/Xcy = X/X_c es una normalización adimensional. Restar la media y dividir por la desviación estándar es una estandarización estadística. Las tres operaciones pueden producir números más convenientes y responden a preguntas diferentes.

La normalización física expresa una variable respecto de una escala con significado en el modelo. Si una concentración se divide por la de equilibrio, el valor uno tiene una interpretación dinámica. En una puntuación estandarizada, en cambio, el cero representa la media del conjunto de referencia y el uno, una desviación estándar por encima. Cambiar la población de referencia cambia esos valores aunque el individuo no haya cambiado. Un puntaje de 30 en una escala de síntomas equivale a z=2z = 2 en una muestra comunitaria con media 20 y desviación estándar 5. En una muestra clínica con media 32 y desviación estándar 8 equivale a z=−0,25z = -0{,}25.

La distinción será importante con datos clínicos. Convertir síntomas o señales en puntuaciones zz puede servir para comparar magnitudes estadísticas, aunque no crea una variable de estado adimensional con significado dinámico. Para afirmarlo hay que especificar qué escala caracteriza el proceso y por qué.

También hay que conservar las transformaciones usadas, porque una curva en variables normalizadas puede ocultar diferencias absolutas relevantes. Dos servicios pueden mostrar la misma ocupación relativa del 90 %, con nueve camas ocupadas de diez en uno y novecientas de mil en el otro. Ante cinco pacientes adicionales, el primero se desborda y el segundo los absorbe. Si interesa la vulnerabilidad ante una demanda adicional fija, la escala absoluta vuelve a importar. Una representación útil para una pregunta puede ser insuficiente para otra.

Escalas en sistemas biológicos, clínicos y organizacionales

Los sistemas que estudiaremos contienen varias escalas simultáneas. En neurociencia coexisten milisegundos de actividad neuronal, segundos de integración perceptiva y ritmos de horas. En farmacología, la absorción, la distribución, la unión a receptores y la adaptación fisiológica pueden tener tiempos distintos. En un sistema asistencial, una decisión se toma en minutos mientras una modificación de la dotación requiere meses. Que estos procesos compartan un sistema no autoriza a representarlos con una única escala temporal.

Cuando las escalas están muy separadas aparecen aproximaciones útiles. Un proceso rápido puede alcanzar casi su régimen estacionario mientras otro cambia lentamente. Más adelante estudiaremos los sistemas lento-rápidos y la eliminación de variables. Por ahora basta reconocer que la razón entre escalas puede informar más que cada valor por separado. Si Traˊpido/TlentoT_{\text{rápido}}/T_{\text{lento}} es muy pequeño, como el cociente de unos 0,002 entre un proceso de diez minutos y otro de tres días, esa relación puede justificar ciertas simplificaciones. Si es cercano a uno, los procesos pueden requerir una representación acoplada.

En psiquiatría conviene ser especialmente cuidadoso. Decir que el cerebro opera en una escala y la experiencia en otra no establece cómo se conectan ambas descripciones. Tampoco basta observar patrones de duraciones diferentes para inferir una jerarquía causal. Las escalas deben definirse a partir de variables y mecanismos concretos. Un episodio clínico, una fluctuación afectiva horaria y una trayectoria de años pueden requerir representaciones distintas, cuya relación debe investigarse sin darla por supuesta a partir del vocabulario multiescala.

La misma cautela vale para los indicadores organizacionales. Una tasa mensual puede ocultar variaciones diarias que producen saturaciones locales. Un registro minuto a minuto puede contener, a la inversa, más detalle del necesario para una pregunta estratégica. La resolución adecuada depende de lo que se intenta explicar, predecir o decidir.

Consolidación

Ante una ecuación, el procedimiento de este capítulo empieza por identificar qué representa cada variable y en qué unidades se expresa. Sigue con la comprobación de que los términos sumados o igualados sean dimensionalmente compatibles y con la búsqueda de combinaciones de parámetros que definan escalas naturales de estado o de tiempo. Termina expresando las variables respecto de esas escalas y observando qué parámetros permanecen en la ecuación adimensional.

En el modelo dX/dt=q−kXdX/dt = q - kX, la consistencia dimensional exigió que qq tuviera dimensión de cantidad por tiempo y kk de tiempo inverso. De esas dimensiones surgieron las escalas Xc=q/kX_c = q/k y Tc=1/kT_c = 1/k. La transformación y=X/Xcy = X/X_c, τ=t/Tc\tau = t/T_c redujo la dinámica a dy/dτ=1−ydy/d\tau = 1 - y, de modo que una familia completa de sistemas comparte la misma evolución relativa aunque sus cantidades y tiempos absolutos difieran. Con una entrada periódica quedó un parámetro, Ω\Omega, que decide si el estado sigue o filtra el forzamiento.

El teorema π de Buckingham generalizó la idea, porque las combinaciones adimensionales reducen un problema y revelan qué razones entre magnitudes controlan su comportamiento. El análisis dimensional no selecciona mecanismos ni garantiza que las variables elegidas sean suficientes. Su fuerza reside en imponer estructura antes del cálculo y en facilitar comparaciones que no dependan de una elección arbitraria de unidades. El cuadro 11.1 reúne los conceptos del capítulo.

Cuadro 11.1. Conceptos de unidades y escalas.

ConceptoDefinición operativaEjemploError frecuente
UnidadConvención para expresar numéricamente una magnitudhora, mg, pacientes por díaConfundir un cambio de unidad con un cambio del fenómeno
DimensiónTipo de magnitud expresado mediante dimensiones básicasT⁻¹ para una tasa de decaimientoAsignar unidades según la posición algebraica del parámetro
Escala característicaMagnitud de referencia natural o elegida para comparar estados o tiemposT_c = 1/kTomarla automáticamente como un promedio observado
Variable adimensionalCociente u otra combinación sin dimensión físicay = X/X_cSuponer que adimensional significa sin interpretación
Parámetro adimensional residualRazón entre procesos que ningún cambio de unidades eliminaΩ = ω/kIgnorarlo al comparar sistemas
Grupo πCombinación adimensional de variables dimensionalesP√(g/L)Creer que determina la ley dinámica completa
EstandarizaciónTransformación respecto de la media y la desviación de una poblaciónPuntuación zTratarla como una variable de estado con significado dinámico

En los próximos capítulos representaremos estados mediante vectores, estudiaremos transformaciones lineales y retomaremos valores propios, derivadas y matrices jacobianas. Las unidades seguirán presentes aunque dejen de ocupar el centro de la exposición. Ante cada ecuación nueva, una de las primeras preguntas seguirá siendo qué significa cada término y qué dimensiones debe tener para que la expresión sea coherente.

Ejercicios

Ejercicio 11.1. (a) En dN/dt=a−kNdN/dt = a - kN, con NN en pacientes y tt en días, indique las unidades de aa y kk. (b) En dx/dt=rx−bx2dx/dt = rx - bx^2, con xx en gramos y tt en horas, indique las unidades de rr y bb. (c) Un programa usa tt en horas y una tasa de eliminación k=0,2k = 0{,}2 estimada por día. ¿Qué tasa diaria simula en realidad y qué consecuencia tiene? (d) ¿Es dimensionalmente válida la expresión dX/dt=q−etdX/dt = q - e^{t}? ¿Cómo podría corregirse?

Ejercicio 11.2. Un depósito vacío sigue dX/dt=q−kXdX/dt = q - kX con q=12q = 12 mg/h y k=0,3 h−1k = 0{,}3\ \mathrm{h^{-1}}. (a) Calcule el equilibrio y la escala temporal. (b) ¿Cuánto tarda en alcanzar el 90 % del equilibrio? (c) ¿Cómo cambian el equilibrio y la escala temporal si qq se duplica y kk se reduce a la mitad?

Ejercicio 11.3. (a) Adimensionalice el modelo logístico dx/dt=r x (1−x/K)dx/dt = r\,x\,(1 - x/K) con las escalas xc=Kx_c = K y tc=1/rt_c = 1/r. ¿Quedan parámetros? (b) Agregue una extracción constante hh, de modo que dx/dt=r x (1−x/K)−hdx/dt = r\,x\,(1 - x/K) - h. Muestre que la ecuación adimensional contiene un único parámetro, H=h/(rK)H = h/(rK). (c) Muestre que hay equilibrios positivos solo si H≤1/4H \le 1/4 y calcúlelos para H=0,2H = 0{,}2. ¿Qué ocurre con H=0,3H = 0{,}3?

Ejercicio 11.4. Considere dy/dτ=1+Acos⁡(Ωτ)−ydy/d\tau = 1 + A\cos(\Omega\tau) - y con A=0,5A = 0{,}5. (a) Calcule la amplitud estacionaria de la oscilación con Ω\Omega igual a 0,1; 1 y 10. (b) En el ejemplo clínico del texto, compruebe que Ω≈25\Omega \approx 25 y calcule qué fracción de la amplitud llega al estado. (c) ¿Con qué período debería alternar la intervención para que llegara la mitad de la amplitud?

Ejercicio 11.5. (a) Obtenga el grupo adimensional del péndulo simple con el método del texto. (b) Con g=9,81 m/s2g = 9{,}81\ \mathrm{m/s^2}, calcule el período de un péndulo de 1 m con oscilaciones pequeñas y con una amplitud de 90°, usando en este caso el factor 1,18. (c) La distancia dd recorrida en caída libre desde el reposo depende del tiempo tt, de gg y quizá de la masa mm. Obtenga el grupo adimensional y diga qué parte del resultado no puede obtenerse por análisis dimensional. (d) Un problema tiene cinco cantidades cuya matriz dimensional tiene rango 3. ¿Cuántos grupos adimensionales independientes hay?

Ejercicio 11.6. (a) Dos servicios tienen una ocupación del 90 %, uno con 10 camas y otro con 1000. ¿Cuál se desborda si llegan cinco pacientes adicionales? (b) Un puntaje de 30 se compara con una muestra comunitaria de media 20 y desviación estándar 5 y con una clínica de media 32 y desviación estándar 8. Calcule ambas puntuaciones zz. (c) ¿Por qué ninguna de ellas es una variable de estado adimensional en el sentido de este capítulo?

Ejercicio 11.7. Un proceso rápido tiene una constante temporal de 10 minutos y uno lento, de 3 días. (a) Calcule la razón entre ambas escalas. (b) ¿Qué simplificación sugiere? (c) ¿Qué habría que verificar antes de aplicarla?

Consultar las soluciones de los ejercicios

Soluciones de los ejercicios

Solución 11.1. (a) La tasa aa se mide en pacientes por día y kk en 1/día. (b) La tasa rr se mide en 1/h y bb en g−1 h−1\mathrm{g^{-1}\,h^{-1}}. (c) El programa aplica 0,2 por hora, es decir, 4,8 por día, una tasa 24 veces mayor que la estimada. La semivida simulada es de unas 3,5 horas en lugar de 3,5 días. (d) No es válida, porque el exponente debe ser adimensional y ete^{t} tampoco tiene las unidades de dX/dtdX/dt. Una corrección posible es dX/dt=q−c et/TdX/dt = q - c\,e^{t/T}, con TT una escala temporal y cc una constante con unidades de cantidad por tiempo.

Solución 11.2. (a) El equilibrio es X∗=12/0,3=40X^* = 12/0{,}3 = 40 mg y la escala temporal, 1/0,3≈3,331/0{,}3 \approx 3{,}33 h. (b) Se requiere 1−e−kt=0,91 - e^{-kt} = 0{,}9, es decir, t=ln⁡10/k≈7,68t = \ln 10/k \approx 7{,}68 h. (c) El equilibrio pasa a 2q/(k/2)=4 X∗=1602q/(k/2) = 4\,X^* = 160 mg y la escala temporal se duplica, hasta unas 6,67 h.

Solución 11.3. (a) Con y=x/Ky = x/K y τ=rt\tau = rt se obtiene dy/dτ=y(1−y)dy/d\tau = y(1 - y), sin parámetros. (b) La extracción se divide por rKrK al adimensionalizar, de modo que dy/dτ=y(1−y)−Hdy/d\tau = y(1 - y) - H con H=h/(rK)H = h/(rK). (c) Los equilibrios cumplen y2−y+H=0y^2 - y + H = 0, es decir, y=(1±1−4H)/2y = \big(1 \pm \sqrt{1 - 4H}\big)/2, que son reales solo si H≤1/4H \le 1/4. Con H=0,2H = 0{,}2 se obtienen 0,276 y 0,724. Con H=0,3H = 0{,}3 no hay equilibrios y la población se agota, porque la extracción supera la tasa máxima de crecimiento, rK/4rK/4. El parámetro residual HH decide así un cambio cualitativo que estudiaremos más adelante como bifurcación.

Solución 11.4. (a) La amplitud es 0,5/1+Ω20{,}5/\sqrt{1 + \Omega^2}, es decir, 0,497; 0,354 y 0,050. (b) Con ω=π\omega = \pi por semana y k=1/8k = 1/8 por semana, Ω=8π≈25,1\Omega = 8\pi \approx 25{,}1 y la fracción es 1/1+25,12≈0,041/\sqrt{1 + 25{,}1^2} \approx 0{,}04. (c) Hace falta 1/1+Ω2=0,51/\sqrt{1 + \Omega^2} = 0{,}5, es decir, Ω=3≈1,73\Omega = \sqrt{3} \approx 1{,}73. Entonces ω=1,73/8≈0,217\omega = 1{,}73/8 \approx 0{,}217 por semana y el período es 2π/ω≈292\pi/\omega \approx 29 semanas.

Solución 11.5. (a) Se obtiene Π=Pg/L\Pi = P\sqrt{g/L}, como en el texto. (b) El período de las oscilaciones pequeñas es 2π1/9,81≈2,012\pi\sqrt{1/9{,}81} \approx 2{,}01 s. Con 90° de amplitud es cercano a 1,18⋅2,01≈2,371{,}18 \cdot 2{,}01 \approx 2{,}37 s. (c) En d ga tb mcd\,g^{a}\,t^{b}\,m^{c}, la masa exige c=0c = 0, la longitud 1+a=01 + a = 0 y el tiempo −2a+b=0-2a + b = 0. Se obtiene Π=d/(g t2)\Pi = d/(g\,t^2), de modo que d=C g t2d = C\,g\,t^2. El valor C=1/2C = 1/2 solo se obtiene resolviendo la dinámica. (d) Hay 5−3=25 - 3 = 2 grupos independientes.

Solución 11.6. (a) El servicio pequeño se desborda, porque necesitaría 14 camas y tiene 10. El grande queda con 905 de 1000. (b) Las puntuaciones son z=(30−20)/5=2z = (30 - 20)/5 = 2 y z=(30−32)/8=−0,25z = (30 - 32)/8 = -0{,}25. (c) Ambas dependen de la población de referencia, de modo que cambian aunque el individuo no cambie. No expresan el estado respecto de una escala característica de su propia dinámica, como un nivel basal o un equilibrio.

Solución 11.7. (a) La razón es 10/(3⋅24⋅60)≈0,002310/(3 \cdot 24 \cdot 60) \approx 0{,}0023. (b) Sugiere tratar la variable rápida como si alcanzara al instante su equilibrio para cada valor de la variable lenta. (c) Habría que verificar que el subsistema rápido tiene un equilibrio estable para cada valor de la variable lenta, que ninguna entrada externa relevante actúa en la escala rápida, que la pregunta concierne a la escala lenta y que las estimaciones de ambas escalas son confiables.

Referencias

Barenblatt, G. I. (2003). Scaling. Cambridge University Press.

Bridgman, P. W. (1922). Dimensional analysis. Yale University Press.

Buckingham, E. (1914). On physically similar systems; illustrations of the use of dimensional equations. Physical Review, 4(4), 345–376. https://doi.org/10.1103/PhysRev.4.345

Bureau International des Poids et Mesures. (2019). The International System of Units (SI) (9.ª ed.). https://www.bipm.org/en/publications/si-brochure

Sterman, J. D. (2000). Business dynamics: Systems thinking and modeling for a complex world. Irwin/McGraw-Hill.

Strogatz, S. H. (2015). Nonlinear dynamics and chaos: With applications to physics, biology, chemistry, and engineering (2.ª ed.). Westview Press.