About StarSimu |
 |
StarSimu is a file that runs on Microsoft® Excel, one of the most
widely used electronic spreadsheet programs
for IBM PC and compatible computers.
This file allows you to calculate the physical variables in a stellar
interior, which in themselves determine the structure of a star. Moreover
you are
absolutely free to specify all the relevant parameters of the simulation.
Hystorical background
StarSimu was born as a proposal for a short dissertation
in stellar astrophysics at the University of Pisa, Italy.
It was prompted by the book "Physics of the stars"
by the Russian astrophysicist S.A. Kaplan. The calculations
concerning the polytropic models mentioned in this book have been
incorporated into StarSimu.
Useful insights have been provided by the book "Astrophysics with
a PC" by Paul Hellings, an excellent guide
for amateur astronomers who want to explore astrophysics with a home
computer. Among the topics that the text covers are the shape of comet
tails, meteor dynamics, evolution of binary-star orbits,
homogeneous stellar models, stellar atmospheres, star formation in the
galaxy, and cosmological models for the Universe.
Assumptions made
-
The static configuration of a gaseous sphere held together by
self-gravitation must satisfy the so-called condition of hydrostatic
equilibrium, a local equation applicable at local points in
the star. The essence of the equation of hydrostatic equilibrium is
that the variation of the pressure P throughout the star is sufficient
to just balance the gravitational field of the star.
Furthermore, the pressure is determined by an equation of state
applicable to the local conditions in the stellar interior.
The above considerations do not in themselves determine the structure
of a star. At each point in a stellar interior any specified increase
of pressure that may be required to
balance the gravitational force is obtainable from an unlimited number of
combinations of density (rho) and temperature (T) values. So one or more
additional conditions relating the physical
variables in the stellar interior are needed.
An approximation introduced at the beginning of this century
is that of eliminating T by assuming that the pressure can be expressed
as some power of the
density only. This "polytropic" pressure-density relationship
is written as P = K*rho^Y, where the exponent Y
is a free parameter which enables us to
consider model stars (called polytropes) of various types, and K is
constant which depends upon the nature itself of the polytrope.
-
The polytrope relation was selected more on mathematical
than physical grounds since it cannot
be deduced from physical arguments. But certain
idealized physical circumstances for a star lead naturally to a polytropic
pressure-density relationship. Here follow
some classical polytropes that were important in the past:
- Y = 5/3: model of a fully convective star (i.e. a star in
which the energy is transported by convection) where radiation
pressure is not important; also used for a
white dwarf supported by the pressure of a non-relativistic degenerate
electron gas.
- Y = 4/3: the standard model of Eddington, used to describe a
star in which the gas pressure is a constant fraction of the total
pressure throughout the star (this particular polytrope corresponds
more closely to stars in radiative equilibrium, i.e. stars for which
the energy is transported by radiative transfer rather than by convection);
also used to describe a white dwarf with relativistic degeneration.
- Y = 1: an approximation for the distribution of stars in a
globular cluster (it represents a self-gravitating isothermal gas, and
a globular cluster can be seen as an isothermal gaseous sphere
in which the particles are not molecules, but stars).
- A polytropic model is only space dipendent, as it gives the variation of the
pressure (P), the density (rho), and the mass inside r (Mr) as a function
of r, the distance to
the center of the model. No time dependent processes - such as nuclear
reactions, chemical mixing, energy transport, or mass loss by stellar
wind - are included. Moreover the model does not contain any information
on the temperature throughout the star.
The star is divided into successive concentric thin shells, or layers,
from the center to the surface. The relationships to be satisfied by the
physical variables of the stellar model are then imposed at each layer.
The distance between two layers - or the
thickness of a single layer - is given by dr, which will be taken as
constant throughout the star.
The thickness dr of a layer should be small compared to the radius
of the star. For a central density of 10^10, 10^11, or 10^12 g/cm³
our suggested values for dr are 1.6, 0.8, or 0.4 km, respectively.
Adopting them the pressure vanishes after
about 800-900 layers,
meaning that we have reached the surface of the star.
The values of r and Mr at that point are therefore the total
radius and the total mass of the star.
Main features
- Input parameters. The input data of the model are the polytropic
exponent Y, the constant K, the central density rho_c (in g/cm³),
and the thickness dr (in km) of a layer.
- Output parameters. The output data are M, the total
mass of the star expressed in solar units (1 solar mass = 2*10^33 g),
and the stellar radius R (in km).
- Physical quantities. For each layer, StarSimu
gives the density rho
(in g/cm³), the pressure P (in dyne/cm²), the
mass Mr (in grams), as functions of the radius r (in cm). It also gives
the same quantities referred to their maximum value, i.e. the
values of rho/rho_c, P/Pc (the subscript "c" means "central"),
and Mr/M as functions of r/R.
- Number of layers. Our file is setted for
calculating up to 1000 layers,
but the spreadsheet's built-in copying feature will allow you
to add as many layers as you want.
- Pictures and graphs. They can be easily produced thanks to the
sophisticated graphing capabilities of all the current spreadsheet
programs.
Back to StarSimu homepage
Copyright © 1998 Torre Luciana Astronomical
Observatory. All rights reserved.