© Charles Chandler

 

The first step is to pre-process models that only specify abundances, to determine the volumes of each element.

1. Assume that each element has been subjected to enough pressure to compress it down to its liquid density, and no further (i.e., the Coulomb barrier is not being crossed).
2. Divide the liquid density by the atomic mass to find the number of atoms in a given volume of an element.
3. Multiply the atoms/volume by the abundance factor, to get the volume.

Knowing the volume of each element, we can determine the average density.

 

Assuming that the elements are not mixed, and that the densest elements have settled to the bottom, we can then estimate the radii of each elementary layer.

 

Abundances

Elements

Old Code

 

Knowing which elements are present at which depths, the FEA engine can then be invoked, to calculate the force of gravity, given the density of each parcel. From that force, the pressure can be determined, and the model temperature can be checked.

 

Code