HydroStatic Star

Purpose

This module aims to reproduce a situation of hydrostatic equilibrium (HSE), a condition where a volume of a fluid is at rest or at constant velocity. This occurs when compression due to gravity is balanced by a pressure gradient force.

we can translate the image above into the condition:

frac{dP}{dh} + \rho \cdot g = 0\\

Implementation

Two fundamental objects are used by this module:

• Ambients
• Profiles
• PointGravity

For this module, the first step to recreate the solution is to determine a 1D profile in hydrostatic equilibrium, this profile can then be applied to a spherically symmetric 3D model using interpolation.

With a little bit of imagination we can think of a 1D column of gas as an array of values:

COLUMN TOP

h5	ρ5	P5
h4	ρ4	P4
h3	ρ3	P3
h2	ρ2	P2
h1	ρ1	P1
h0	ρ0	P0

COLUMN BASE

We can relate each level of the column to an position in an array which contains height (h), density(ρ) and pressure(P).