Mostrar código-fonte de Equilíbrio hidrostático

{{Mecânica do contínuo}}
Em [[mecânica de fluidos]], diz-se que um [[fluido]] está em '''equilíbrio hidrostático''' ou '''mecânico''' quando ele estiver em repouso, isto é, quando a velocidade de [[escoamento]] não mudar ao longo do tempo em todos os pontos do fluido. Isso significa que um fluido em equilíbrio (perfeito) é um referencial inercial, já que nele a força resultante que atua em cada ponto é nula.

O equilíbrio hidrostático ocorre em [[estrela]]s quando há o balanço entre a atração do [[campo gravitacional]] das camadas superiores pelas inferiores, e entre a diferença de [[pressão]] entre essas camadas, que provoca uma força cujo sentido é contrário ao centro da estrela. Se a força de pressão e a força da gravidade não forem equivalentes, a estrela sofrerá mudanças radicais em sua estrutura.<ref name= if-ufrgs>{{citar web|url=http://astro.if.ufrgs.br/evol/node13.htm|título=Equilíbrio Hidrostático|data=|website=IF-UFRGS|acessodata=30 de outubro de 2018|publicado=30 de setembro de 2008|ultimo=S.O.|primeiro=Kepler}}</ref> Estrelas com menos massa permanecem em equilíbrio por mais tempo, pois consomem seu estoque de [[hidrogênio]] mais lentamente e, consequentemente, produzindo por mais tempo a força externa necessária para manter sua estrutura. O [[Sol]], por exemplo, durará mais cerca de 4,5 bilhões de anos. Além de estrelas, o equilíbrio hidrostático é importante para todos os corpos celestes.

== Equilíbrio hidrostático para fluidos na Terra ==
{{Artigo principal|Lei de Stevin}}
A equação do equilíbrio hidrostático para fluidos em um campo gravitacional constante, como na [[superfície terrestre]], é:
{{{indent|:}}}{|cellpadding="{{{cellpadding|5}}}" style="border:{{{border|2}}}px solid {{{border colour|#0073CF}}};background: {{{background colour|#F5FFFA}}}; text-align: center;"
|{{{title|'''Equação do Equilíbrio Hidrostático''' ''(gravidade uniforme)''}}}
{{{equation|<math>\frac{dP}{dz}=-\rho(z)g</math>}}}
|}

Em que <math>z</math> é a altura com relação a uma origem (geralmente o solo), <math>P=P(z)</math> é a pressão em função da altura e <math>\rho(z)</math> é a densidade em um altura particular.

:{| class="toccolours collapsible collapsed" width="50%" style="text-align:left"
!Demonstração<ref name= uff-termo>{{citar web|url=https://www.youtube.com/watch?v=tag-pS1KKt0|título=Aula 2.4 - Exercícios do capitulo 1 (II): atmosfera exponencial|data=|website=Youtube|acessodata=11 de novembro de 2018|publicado=30 de julho de 2015|ultimo=Martins|primeiro=Jorge Sá}}</ref>
|-
|
Para haver equilíbrio hidrostático, todas as camadas de um fluido devem estar sujeitas a uma força nula. Considerando uma camada delimitada pelas altura <math>z</math> e <math> z + \Delta z</math>, com uma área transversal <math>A</math> da camada. Há duas força de pressão atuando sobre este elemento: uma causada pela pressão fluido acima de <math>z+\Delta z</math> (vertical e para baixo); e outra causada pela pressão do fluido abaixo da altura <math>z</math> (vertical e para cima). Essas forças podem ser obtidas pelo produto da pressão na camada pela área transversal:

:<math> F_{superior} = - P(z+\Delta z).A</math>
:<math> F_{inferior} = P(z).A</math>

Já o peso da camada de ar é igual ao produto de sua massa <math>M</math> pela aceleração <math>g</math> da gravidade. Lembrando que <math>M =\rho V = \rho A \Delta z</math>, em que <math>V</math> é o volume da camada e <math>\rho</math> é a densidade do fluido, o peso é expresso da seguinte forma:

:<math> F_{g} = - Mg = - \rho g A \Delta z</math>

Para que haja equilíbrio, a soma de todas essas forças tem que ser nula:

:<math> F_{inferior} + F_{superior} + F_{g} = 0 </math>
:<math> - P(z+\Delta z).A + P(z).A - \rho g A \Delta z = 0 </math>
:<math> - \left( P(z+\Delta z) - P(z) \right)A = \left(\rho g \right)A \Delta z</math>

Dividindo ambos os lados da equação por <math>-A \Delta z</math> e identificando <math>P(z+\Delta z) - P(z) = \Delta P</math>, obtém-se:

:<math> \frac{\Delta P}{\Delta z}= -\rho g</math>

Diminuindo a espessura da camada até torná-la infinitesimal (<math>\Delta z \to 0</math>), encontra-se, finalmente:

{{{indent|:}}}{|cellpadding="{{{cellpadding|5}}}" style="border:{{{border|2}}}px solid {{{border colour|#000000}}};background: {{{background colour|#f5f5f5}}}; text-align: center;"
|{{{title|}}}
{{{equation|<math>\frac{dP}{dz}=-\rho(z)g</math>}}}
|}
|}

Caso a densidade seja constante, essa equação reduz-se ao [[Lei de Stevin|princípio de Stevin]].

== Equilíbrio hidrostático em estrelas ==
[[File:Equilíbrio Hidrostático.png|thumb|300px|right|Ilustração do balanço entre a pressão interna do gás e a auto-gravidade de uma estrela. Essa descrição aplica-se para o nosso Sol, que não se expande nem se contrai, porque está em equilíbrio.]]
{{Artigo principal|Formação estelar}}

Estrelas são formadas dentro de concentrações relativamente densas de gás e [[poeira interestelar]], denominadas [[Nuvem molecular|nuvens moleculares]]. Isso ocorre devido à força gravitacional que colapsa um fragmento da nuvem, que então passa a ser considerado um objeto individual com campo gravitacional próprio, independente do restante da nuvem, em primeira aproximação. Conforme esse fragmento se contrai, o gás que compõem sua região central aquece e aumenta sua densidade, formando a [[protoestrela]]. O disco externo, em rotação devido à conservação do [[momento angular]] do fragmento, evolui mais lentamente e formará, possivelmente, um sistema planetário. A matéria continua a cair em direção à protoestrela, aumentando rapidamente seu tamanho.<ref name= pdf-ufrgs>{{citar web|url=http://www.if.ufrgs.br/~tiberio/disciplinas/fis02010/aula21.pdf |titulo=Evolução estelar e equilíbrio hidrostático|website=IF-UFRGS|acessodata=9 de novembro de 2018|publicado=|ultimo=Vale|primeiro=Tibério B.}}</ref>

A queda de matéria em direção ao centro é então interrompida quando se iniciam as reações de [[fusão nuclear]] no centro da protoestrela. Tais reações produzem um radiação que causa um gradiente de pressão entre camadas da atmosfera. Esse gradiente sempre aponta para longe do centro, opondo-se a força gravitacional e impedindo o colapso do objeto. Diz-se a partir desse momento que foi atingido o equilíbrio hidrostático e a protoestrela passa a ser uma [[estrela]] de tamanho bem definido, já que o objeto não mais se contrai.<ref name= pdf-ufrgs/>

=== Equação ===
A equação fundamental do equilíbrio hidrostático estelar é:<ref name= if-ufrgs/>
{{{indent|:}}}{|cellpadding="{{{cellpadding|5}}}" style="border:{{{border|2}}}px solid {{{border colour|#FF0000}}};background: {{{background colour|#F5FFFA}}}; text-align: center;"
|{{{title|'''Equação do Equilíbrio Hidrostático''' ''(não-relativística)''}}}
{{{equation|<math>\frac{dP}{dr}=-\frac{G M(r)\rho(r)}{r^2}</math>}}}
|}

Em que <math>r</math> é a distância ao centro da estrela, <math>P=P(r)</math> é a pressão em função dessa distância, <math>\rho(r)</math> é a densidade naquela particular distância, e <math>M(r)</math> é a massa da estrela do centro até aquela distância, calculada através da seguinte integral:<ref name= if-ufrgs/>

:<math>M(r)=\int_0^r 4\pi r'^2\rho(r')dr'</math>

:{| class="toccolours collapsible collapsed" width="80%" style="text-align:left"
!Demonstração<ref name= if-ufrgs/>
|-
|
Para haver equilíbrio hidrostático, todas as forças que atuam em qualquer elemento de volume dentro da estrela têm que ser compensadas exatamente, já que uma força resultante não nula implicaria movimentos e, portanto, mudanças na estrutura. As únicas forças que precisamos considerar são a força gravitacional, para dentro, e a força de pressão, para fora.

Considerando um elemento de volume cilíndrico, a uma distância <math>r</math> do centro da estrela, com seu eixo na direção do centro, com uma seção transversal <math>dS</math> e um comprimento <math>dr</math>. A força de pressão atuando sobre este elemento, causada pela diferença de pressão nas paredes interna e externa, é expressa pelo produto da área infinitesimal do cilindro por essa diferença de pressão:
:<math> dF = - dP.dS</math>
Pela [[regra da cadeia]], <math> dP = \frac{dP}{dr}.dr</math>. O sinal negativo da equação acima significa então que a pressão aumenta conforme o raio <math>r</math> diminui.

Já a força gravitacional atuante nesse volume infinitesimal é dada pelo produto do elemento de massa <math>dm</math> pela aceleração gravitacional naquele ponto <math>g=\frac{GM}{r^2}</math>, em que <math>M = M(r)</math> é a massa total abaixo do elemento de volume que, portanto, o atrai. Essa massa das camadas inferiores pode ser calculada se conhecida a função densidade <math>\rho(r)</math> da estrela, já que <math>\rho(r)=\frac{dM}{dV}=\frac{dm}{dV}</math>, em que <math>dV</math> é o elemento de volume. Assumindo a simetria esférica <math>(V(r')=\frac{4\pi r'^3}{3}</math>), portanto, a massa das camadas inferiores é calculada pela seguinte integral:
:<math> M(r) = {\int_{0}^{r}} dM = {\int_{0}^{r}} {\frac{dM}{dV}.\frac{dV}{dr'}} dr' = {\int_{0}^{r}} {\rho(r).4\pi r'^2} dr'</math>

Finalmente, a expressão do elemento de força gravitacional atuando no particular volume deve ser (fazendo uso da regra da cadeia e de que <math>dV=drdS</math>):
:<math>dF_{g}= dm.g =\frac{dm}{dV}\frac{GM(r)}{r^2} dV = \rho(r)\frac{GM(r)}{r^2}drdS</math>

Para que haja equilíbrio segue que <math> dF = dF_{g}</math>:

:<math> - \frac{dP}{dr}.drdS = \rho(r)\frac{GM(r)}{r^2}drdS</math>

Isso implica que:

{{{indent|:}}}{|cellpadding="{{{cellpadding|5}}}" style="border:{{{border|2}}}px solid {{{border colour|#000000}}};background: {{{background colour|#f5f5f5}}}; text-align: center;"
|{{{title|}}}
{{{equation|<math>\frac{dP}{dr}=-\frac{G M(r)\rho(r)}{r^2}</math>}}}
|}
|}

A equação de equilíbrio hidrostático acima é embasada pela [[mecânica Newtoniana]]. Em regimes relativísticos, em que tal abordagem não é válida, é necessário usar a [[Equação Tolman-Oppenheimer-Volkoff|equação de Tolman-Oppenheimer-Volkoff]].

==== Densidade constante ====
Como a equação do equilíbrio hidrostático é uma [[Equação integral|equação integro]][[Equação diferencial|-diferencial]], é necessário conhecer a função densidade <math>\rho(r)</math> para resolvê-la. A solução mais simples e não-trivial (que seria a solução para <math>\rho(r)=0</math>) é a solução para densidade <math>\rho(r) = constante</math>.

Usando os métodos de resolução de equações diferenciais, descobre-se que a pressão <math>P</math> em função da distância <math>r</math> ao centro da estrela é dado pela seguinte expressão:

:<math>P(r)=\frac{2 \pi G}{3} (\rho R)^2 \left [1-\left (\frac{r}{R} \right)^2 \right] </math>

Em que <math>R</math> é o raio da estrela.

:{| class="toccolours collapsible collapsed" width="80%" style="text-align:left"
!Demonstração
|-
|
Para essa demonstração, basta usar a condição de <math>\rho(r) = constante</math> para obter a massa <math>M(r)</math>:

:<math>M(r) = \int_0^r 4\pi r'^2\rho(r')dr' = 4\pi\rho\int_0^r  r'^2 dr' = \frac{4\pi\rho}{3} r'^3 \Biggr|_{0}^{r} = \frac{4\pi\rho}{3} r^3</math>

Conhecida a massa total em função de <math>r</math>, usa-se então a equação principal:

:<math>\frac{dP}{dr}=-\frac{G M(r)\rho}{r^2} = -\frac{G\rho}{r^2}\frac{4\pi\rho}{3} r^3 = -\frac{4\pi G}{3}\rho^2 r </math>

Mudando de variável temporariamente (<math>r \to r'</math>) para integrar e colocar <math>r</math> como uma variável no extremo de integração inferior, obtemos:

:<math>\int_r^R \frac{dP}{dr'} dr' = -\int_r^R \frac{4\pi G}{3}\rho^2 r' dr'</math>

:<math> P(r') \Biggr|_{r}^{R} = -\frac{4\pi G}{3}\rho^2 \frac{r'^2}{2} \Biggr|_{r}^{R}</math>

:<math> P(R) - P(r)  = -\frac{2\pi G}{3}\rho^2 (R^2 - r^2)</math>

Em que <math>R</math> é o raio da estrela. Como à distância do raio da estrela a atmosfera deve sumir, a pressão deve ser nula, ou seja, <math>P(R)=0</math>. Desenvolvendo o que sobrou na equação, finalmente obtemos:

:<math> - P(r)  = -\frac{2\pi G}{3}\rho^2 [R^2] \left[1 - \left(\frac{r^2}{R^2} \right) \right]</math>

{{{indent|:}}}{|cellpadding="{{{cellpadding|5}}}" style="border:{{{border|2}}}px solid {{{border colour|#000000}}};background: {{{background colour|#f5f5f5}}}; text-align: center;"
|{{{title|}}}
{{{equation|<math>P(r)=\frac{2 \pi G}{3} (\rho R)^2 \left [1-\left (\frac{r}{R} \right)^2 \right]</math>}}}
|}
|}

{{Referências}}
{{controle de autoridade}}
{{Portal3|Física|Astronomia}}

[[Categoria:Mecânica de fluidos]]
[[Categoria:Conceitos de astronomia]]
[[Categoria:Definição de planeta]]