# Sistema de masa variable

En Mecánica, un sistema de masa variable es un conjunto de materia cuya masa varía con el tiempo. La Segunda Ley de Newton no se puede aplicar directamente dado que sólo es válida para sistemas de masa constante. Por lo tanto, la dependencia de la masa m respecto del tiempo se analizará tratando al sistema como si de uno cerrado se tratase; esto significa considerar un sistema que contenga tanto a la masa principal como a la masa variante. Al añadir un término que considera el momento de la masa que se adhiere o desprende, la ecuación general de movimiento de una masa variable (que expresa el cambio del momento respecto al tiempo) puede escribirse como:
    $\mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m{\mathrm {d} \mathbf {v} \over \mathrm {d} t}$
Donde F$_{ext}$ es la fuerza neta externa ejercida en el sistema cerrado, v$_{rel}$ es la velocidad relativa de la masa que está escapando o ingresando con respecto al centro de masa del cuerpo, y v es la velocidad del cuerpo. En la ingeniería aeroespacial, la cual estudia la mecánica de los cohetes, el término v$_{rel}$ se lo llama la velocidad efectiva de escape y se denomina como v$_e$.
Un error de concepto frecuente es que un sistema de masa variable puede describirse como la derivada respecto del tiempo del producto de la masa con la velocidad, ya que la fuerza neta externa es el cambio en el momento respecto del tiempo. Si bien este hecho no cambia, en un sistema abierto la variación del momento ya no puede describirse como el cambio en el producto de la masa con la velocidad, pues al hacerlo se estaría olvidando al impulso asociado a la variación de masa del sistema. Dado que estos sistemas, como el del cohete que pierde combustible y eyecta gases, no son sistemas cerrados y dado que la masa no se debe tratar como una variable en función del tiempo, la siguiente fórmula no es correcta:
    $\mathbf {F} _{\mathrm {net} }={\frac {\mathrm {d} }{\mathrm {d} t}}{\bigm(t)\mathbf {v} (t){\big ]}=m(t){\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}+\mathbf {v} (t){\frac {\mathrm {d} m}{\mathrm {d} t}}\qquad \mathrm {(Incorrecto)}$
La falacia de esta fórmula puede verse en que no respeta la Invariancia galileana la cual sostiene que un objeto de masa variable con F = 0 en un marco de referencia, tendrá F ≠ 0 en otro.

## Deducción de la fórmula

Existen distintas derivaciones para la ecuación de movimiento del sistema de masa variable, que depende de si la masa está escapando o ingresando al cuerpo (en otras palabras, si la masa del cuerpo está aumentando o disminuyendo). Para simplificar los cálculos, tratamos a los cuerpos como puntuales. También se asume que la masa no es capaz de aplicarle una fuerza externa al cuerpo fuera del evento de acreción/ablación.

### Caso acreción de masa

La siguiente derivación es para un cuerpo que está ganando masa. Un cuerpo con masa variable en el tiempo m se mueve con velocidad v a un tiempo inicial t. En el mismo instante, una partícula de masa dm se mueve con una velocidad u. El momento inicial se puede escribir como:
    $\mathbf {p} _{\mathrm {1} }=m\mathbf {v} +\mathbf {u} \mathrm {d} m$
Ahora, a un tiempo t + dt, el cuerpo principal y la masa puntual se acrecientan en un solo cuerpo de velocidad v + dv. Por lo tanto, el nuevo momento del sistema quedaría:
    $\mathbf {p} _{\mathrm {2} }=(m+\mathrm {d} m)(\mathbf {v} +\mathrm {d} \mathbf {v} )=m\mathbf {v} +m\mathrm {d} \mathbf {v} +\mathbf {v} \mathrm {d} m+\mathrm {d} m\mathrm {d} \mathbf {v}$
Como dmdv es el producto de dos valores muy pequeños, este término se puede despreciar. Es decir que para un cierto dt el momento del sistema sería:
    $\mathrm {d} \mathbf {p} =\mathbf {p} _{\mathrm {2} }-\mathbf {p} _{\mathrm {1} }=(m\mathbf {v} +m\mathrm {d} \mathbf {v} +\mathbf {v} \mathrm {d} m)-(m\mathbf {v} +\mathbf {u} \mathrm {d} m)=m\mathrm {d} \mathbf {v} -(\mathbf {u} -\mathbf {v} )\mathrm {d} m$
Entonces, por la segunda ley de Newton:
    $\mathbf {F} _{\mathrm {net} }={\frac {\mathrm {d} \mathbf {p} }{\mathrm {d} t}}={\frac {m\mathrm {d} \mathbf {v} -(\mathbf {u} -\mathbf {v} )\mathrm {d} m}{\mathrm {d} t}}=m{\frac {\mathrm {d} \mathbf {v} }{\mathrm {d} t}}-(\mathbf {u} -\mathbf {v} ){\frac {\mathrm {d} m}{\mathrm {d} t}}$
Fijarse en que u - v es la velocidad de dm relativa a m, simbolizada como v$_{rel}$, esta ecuación final se puede reescribir como:
    $\mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m{\mathrm {d} \mathbf {v} \over \mathrm {d} t}$

### Caso ablación/eyección de masa

En un sistema en el que la masa está saliendo de un cuerpo, la deducción varía ligeramente. A un tiempo t, una masa m viaja a una velocidad v, por lo tanto el momento inicial del sistema viene dado por:
    $\mathbf {p} _{\mathrm {1} }=m\mathbf {v}$
Como el cuerpo principal va a estar perdiendo masa, dm será negativo, esto significa que para un tiempo t + dt el momento de sistema será:
    $\mathbf {p} _{\mathrm {2} }=(m+\mathrm {d} m)(\mathbf {v} +\mathrm {d} \mathbf {v} )+\mathbf {u} (-\mathrm {d} m)=m\mathbf {v} +m\mathrm {d} \mathbf {v} +\mathbf {v} \mathrm {d} m+\mathrm {d} m\mathrm {d} \mathbf {v} -\mathbf {u} \mathrm {d} m$
donde u es la velocidad de la masa eyectada. Por tanto, en dt el momento del sistema es:
    $\mathrm {d} \mathbf {p} =\mathbf {p} _{\mathrm {2} }-\mathbf {p} _{\mathrm {1} }=(m\mathbf {v} +m\mathrm {d} \mathbf {v} +\mathbf {v} \mathrm {d} m-\mathbf {u} \mathrm {d} m)-(m\mathbf {v} )=m\mathrm {d} \mathbf {v} -(\mathbf {u} -\mathbf {v} )\mathrm {d} m$
Ya verificamos que el dp es el mismo en ambos casos; es decir, la conclusión es la misma y la fórmula se mantiene.

## Formas de expresión

Por definición, la aceleración está dada por a = dv/dt, entonces la ecuación de movimiento de un sistema de masa variable se puede escribir:
    $\mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m\mathbf {a}$
Para cuerpos que no son tratados como masas puntuales, a debe ser reemplazado por a$_{cm}$, la aceleración del centro de masa del sistema, es decir:
    $\mathbf {F} _{\mathrm {ext} }+\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}=m\mathbf {a} _{\mathrm {cm} }$
La fuerza debida a un empuje se suele definir como $\mathbf {F} _{\mathrm {empuje} }=\mathbf {v} _{\mathrm {rel} }{\frac {\mathrm {d} m}{\mathrm {d} t}}$ de manera que:
    $\mathbf {F} _{\mathrm {ext} }+\mathbf {F} _{\mathrm {empuje} }=m\mathbf {a} _{\mathrm {cm} }$
Esta forma nos muestra que un cuerpo puede tener aceleración debida a un empuje, incluso si no actúan fuerzas en él. (F$_{ext}$ = 0). Nótese finalmente que si tomamos F$_{net}$ como la suma de F$_{ext}$ y F$_{empuje}$ entonces la ecuación retoma la forma usual de la segunda ley de Newton
    $\mathbf {F} _{\mathrm {net} }=m\mathbf {a} _{\mathrm {cm} }$