22.08.2015 Views

Arquivo do trabalho - IAG - USP

Arquivo do trabalho - IAG - USP

Arquivo do trabalho - IAG - USP

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Formação estelar em nuvens turbulentas e o transporte de Fluxo MagnéticoSantos-Lima et al. (2010)), usan<strong>do</strong> um resolve<strong>do</strong>r HLL e Runge-Kutta de segunda ordempara a integração temporal. A equação de Poisson é resolvida por um algorítimo basea<strong>do</strong>em multigrid (Press, Teukolsky & Flannery, 1978). Veja uma descrição mais detalha <strong>do</strong>código numérico no Ap. B.5.2 Condições IniciaisNa Tabela 5.1 temos as condições iniciais usadas para to<strong>do</strong>s os modelos simula<strong>do</strong>s come sem auto-gravidade. Os modelos sem auto-gravidade são denomina<strong>do</strong>s ”R”, os autogravitantescomeçan<strong>do</strong> fora <strong>do</strong> equilíbrio são chama<strong>do</strong>s ”N”e o modelo começan<strong>do</strong> emequilíbrio magneto-hidrostático é denomina<strong>do</strong> modelo ”E”. Os modelos ”R”, embora sejamindependentes da escala são também descritos aqui por suas condições iniciais emunidades físicas para comparação com os outros modelos, exceto pelo parâmetro adimensionalA que mede a força <strong>do</strong> potencial gravitacional A = (G M ⋆ )/(L c 2 s ), onde G é aconstante gravitacional de Newton, M ⋆ é a massa <strong>do</strong> potencial estelar, L é o tamanho <strong>do</strong><strong>do</strong>mínio computacional (que também dá a unidade de distância <strong>do</strong> código), e c s é a velocidade<strong>do</strong> som (que é também a unidade de velocidade no código). Os outros parâmetros naTabela são a massa total M tot , isto é, a soma das massas das estrelas e <strong>do</strong> gás, a densidadenumérica inicial gasosa da nuvem n 0 , a qual para o modelo em equilíbrio refere-se à densidadecentral inicial <strong>do</strong> gás da nuvem, o campo magnético inicial B 0 admiti<strong>do</strong> uniforme,a razão pressão térmica-pressão magnética β, a velocidade turbulenta v turb , a velocidadede Alfvén v A , a massa de Jeans inicial da nuvem em presença de campo magnético m J,B ,e a razão entre as densidades iniciais de energia turbulenta e magnética E turb /E mag parao núcleo (r c ≤ 0.3pc). As velocidades turbulenta e de Alfvén são dadas em unidades davelocidade <strong>do</strong> som isotérmica c s = 3×10 4 cm/s.A Tabela 5.2 apresenta o tempo dinâmico, t dyn = L/c s . Este pode ser compara<strong>do</strong> como tempo de queda-livre <strong>do</strong> gás da nuvem, da<strong>do</strong> por,√3πt ff =32 G ρe ao tempo de difusão <strong>do</strong> campo magnético, t diff ,(5.6)89

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!