En dinámica de gases , la ecuación de Chaplygin , llamada así en honor a Sergei Alekseevich Chaplygin (1902), es una ecuación diferencial parcial útil en el estudio del flujo transónico . [1] Es
![{\displaystyle {\frac {\partial ^{2}\Phi }{\partial \theta ^{2}}}+{\frac {v^{2}}{1-v^{2}/c^{ 2}}}{\frac {\partial ^{2}\Phi }{\partial v^{2}}}+v{\frac {\partial \Phi }{\partial v}}=0.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Aquí, está la velocidad del sonido , determinada por la ecuación de estado del fluido y la conservación de la energía. Para gases politrópicos, tenemos , donde es la relación de calor específico y es la entalpía de estancamiento, en cuyo caso la ecuación de Chaplygin se reduce a![{\displaystyle c=c(v)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle c^{2}/(\gamma -1)=h_{0}-v^{2}/2}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \gamma}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle h_{0}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {\partial ^{2}\Phi }{\partial \theta ^{2}}}+v^{2}{\frac {2h_{0}-v^{2}}{2h_ {0}-(\gamma +1)v^{2}/(\gamma -1)}}{\frac {\partial ^{2}\Phi }{\partial v^{2}}}+v{ \frac {\partial \Phi }{\partial v}}=0.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La ecuación de Bernoulli (ver la derivación a continuación) establece que la velocidad máxima ocurre cuando la entalpía específica está en el valor más pequeño posible; se puede tomar la entalpía específica como cero correspondiente a la temperatura del cero absoluto como valor de referencia, en cuyo caso es la velocidad máxima alcanzable. Las integrales particulares de la ecuación anterior se pueden expresar en términos de funciones hipergeométricas . [2] [3]![{\displaystyle 2h_{0}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Derivación
Para el flujo potencial bidimensional, la ecuación de continuidad y las ecuaciones de Euler (de hecho, la ecuación de Bernoulli compresible debido a la irrotacionalidad) en coordenadas cartesianas que involucran las variables velocidad del fluido , entalpía específica y densidad son![{\displaystyle (x,y)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle h}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle\rho}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{alineado}{\frac {\partial }{\partial x}}(\rho v_{x})+{\frac {\partial }{\partial y}}(\rho v_{y })&=0,\\h+{\frac {1}{2}}v^{2}&=h_{o}.\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
con la ecuación de estado actuando como tercera ecuación. Aquí está la entalpía de estancamiento, es la magnitud del vector velocidad y es la entropía. Para el flujo isentrópico , la densidad se puede expresar solo como función de la entalpía , que a su vez, utilizando la ecuación de Bernoulli, se puede escribir como .![{\displaystyle \rho =\rho (s,h)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle h_ {o}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v^{2}=v_{x}^{2}+v_{y}^{2}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle s}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \rho =\rho (h)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \rho =\rho (v)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Como el flujo es irrotacional, existe un potencial de velocidad y su diferencial es simplemente . En lugar de tratar y como variables dependientes, utilizamos una transformación de coordenadas de modo que y se conviertan en nuevas variables dependientes. De manera similar, el potencial de velocidad se reemplaza por una nueva función ( transformación de Legendre ) [4]![{\displaystyle \phi}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle d\phi =v_{x}dx+v_{y}dy}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v_{x}=v_{x}(x,y)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v_{y}=v_{y}(x,y)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x=x(v_{x},v_{y})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle y = y (v_ {x}, v_ {y})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \Phi =xv_{x}+yv_{y}-\phi }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
tal entonces su diferencial es , por lo tanto![{\displaystyle d\Phi =xdv_{x}+ydv_{y}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle x={\frac {\partial \Phi }{\partial v_{x}}},\quad y={\frac {\partial \Phi }{\partial v_{y}}}.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Introduciendo otra transformación de coordenadas para las variables independientes de a según la relación y , donde es la magnitud del vector velocidad y es el ángulo que forma el vector velocidad con el eje -, las variables dependientes se convierten en![{\displaystyle (v_{x},v_{y})}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle (v,\theta)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v_{x}=v\cos \theta }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v_{y}=v\sin \theta }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle v}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle \theta}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\ Displaystyle v_ {x}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\begin{aligned}x&=\cos \theta {\frac {\partial \Phi }{\partial v}}-{\frac {\sin \theta }{v}}{\frac {\partial \Phi }{\partial \theta }},\\y&=\sin \theta {\frac {\partial \Phi }{\partial v}}+{\frac {\cos \theta }{v}}{\ frac {\partial \Phi }{\partial \theta }},\\\phi &=-\Phi +v{\frac {\partial \Phi }{\partial v}}.\end{aligned}}}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
La ecuación de continuidad en las nuevas coordenadas se convierte en
![{\displaystyle {\frac {d(\rho v)}{dv}}\left({\frac {\partial \Phi }{\partial v}}+{\frac {1}{v}}{\frac {\partial ^{2}\Phi }{\partial \theta ^{2}}}\right)+\rho v{\frac {\partial ^{2}\Phi }{\partial v^{2}} }=0.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Para flujo isentrópico, donde es la velocidad del sonido. Usando la ecuación de Bernoulli encontramos![{\displaystyle dh=\rho ^{-1}c^{2}d\rho }](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle c}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {d(\rho v)}{dv}}=\rho \left(1-{\frac {v^{2}}{c^{2}}}\right)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
dónde . Por lo tanto, tenemos![{\displaystyle c=c(v)}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
![{\displaystyle {\frac {\partial ^{2}\Phi }{\partial \theta ^{2}}}+{\frac {v^{2}}{1-{\frac {v^{2} }{c^{2}}}}}{\frac {\partial ^{2}\Phi }{\partial v^{2}}}+v{\frac {\partial \Phi }{\partial v} }=0.}](data:image/gif;base64,R0lGODlhAQABAIAAAAAAAP///yH5BAEAAAAALAAAAAABAAEAAAIBRAA7)
Ver también
Referencias
Wikimedia Commons tiene medios relacionados con la ecuación de Chaplygin .
- ^ Chaplygin, SA (1902). Sobre corrientes de gas. Colección completa de obras.(Ruso) Izd. Akád. Nauk SSSR, 2.
- ^ Sedov, LI, (1965). Problemas bidimensionales en hidrodinámica y aerodinámica. Capítulo X
- ^ Von Mises, R., Geiringer, H. y Ludford, GSS (2004). Teoría matemática del flujo de fluidos compresibles. Corporación de mensajería.
- ^ Landau, LD ; Lifshitz, EM (1982). Mecánica de fluidos (2 ed.). Prensa de Pérgamo. pag. 432.