VB2019 Simultion of our Galagy, including sagitarius and 10 mini Black holes
Introduction and Purpose
The purpose of the program "VB Sagittarius" is to simulate the black hole Sagittarius A* and 10 mini BH's:
This program is part of what is called the "Sagitarrius Project". To read more select this link
4.4 Example of type 2 experiment: The Star s2 around the BH Sagittarius A*
Temporary important is this: Carpe Diem and Sagittarius.
Please study that document in relation with Display 2, shown below, and observe that the w parameter changes when Display 2 is selected or not selected. The same is true for inclination. The omega parameter is constant, which it should not.
The program consists of 4 Forms: Control, Display, Menu and Figure
- The Control Form is used to select a planet simulation and to perform a simulation.
- The Display Form shows the results of a simultion and the result of a simulation saved in a Data Bae
- The Menu Form is to select additional parameters used by a planet simulation
- The Figure Form is to select a simulation saved in a Data Base
Control Form - Operation
The control form is used to select and perform a simulation.
The control form is also use to create a certain number of predefined figures.
- First select a particular simulation to be performed.
In the example this is simulation 7, which demonstrates "Sagittarius and 10 Mini BH's"
- Next Select "Parameter Selection Menu" in order to change cetain parameters relevant for this simulation. This is an option.
- Next Select Start
- Next the simulation will show the Display Form as shown in Display 2.
An additinal display shows the evolution of the simulation.See Display 3.
This option is currently not available.
- First select a particular figure to be created.
Select Menu Selection Form. Select SubTest and return back to control form
- Next Select Start
- Next the simulation will show the Display Form as shown in Figure 2
Display Form - Operation
The Display Form is divided in 5 parts
Now we go to something more specific data.
- A top header section which shows general information.
- A top box left showing information about 10 mini BH's
- A top box right showing information about 10 mini BH's
- A bottum box left showing information about 10 mini BH's
- A bottum box right showing specific information about mini BH 1.
This infomation is related to the first line of each of the other 3 Boxes.
- The Column marked Date shows 10 dates. This shows the information the last time when the BH was either at aphelion (farthest distance) or at perihelion (closest) distance. For line 1, BH 1, this in the year 3247.
When you select Display 2 you get a different image. Now line 1, BH 1, this changes in the year 2254.
- In that same column there is a section marked # of revs . This shows the total number of revolutions of the BH since the start of the simulation. The revolutions are counted based on the number of times when the mini BH reaches aphelion. For BH 1 this are 75 revolutions. When Display 2 is selected this changes to 14.
- The bottom box right shows a column marked Year. Specific of interest are the 4 values in the row starting with "t 0",
"t -1", "t -2" and "t max". This is the angle when the BH is at aphelion. Specific this is the longitude of the pericenter, W, which is the sum of omega and w.
In this particular case for BH 1 omega = 228 degrees and w = 293 degrees (See bottum left box) which gives a total of 521 degrees or 161 degrees, which is a good match.
The reason for using w in stead of w is because the original program was used to simulate the planets around the Sun.
- The column "time in sec" shows the time when the BH was aphelion. This was in cycle 193597952, 193597951 and 193597950. To calculate these numbers you have to divide the "time in sec" value by 200 which is delta time dt between each calculation cycle. This value is shown in the header behind dtime.
- Using the first three values at "t 0","t -1" and "t -2" it is possible to calulate the maximum distance which is shown in the row "t max". Unfortunate the actual distances shown in column dmax don't demonstrate this.
This reflects work in progress
Both Display 3a and 3b show the same situation. The situation enlarged with a factor 2.
The most remarqueble chance is in the star S62, in green. Also interesting is the star S2. What is observed is the mutual influence between the 10 stars or 10 mini Black Holes.
When Display 3a is not selected the display shows the situation in the year 2213.
When Display 3a is selected the display shows the situation in the year 16050.
When Display 3b is not selected the display shows the situation in the year 2218.
When Display 3b is selected the display shows the situation in the year 16171.
The simulation is solely done using Newton's law.
Menu Form - Operation
Figure Form - Operation
VB2019 Sagittarius Software Description
Source and destination
The program supports 1 BH in the center and 10 sun's or mini BH's. Together 11 objects.
The names of the of the objects are stored in the array BHname0$(0-10)
This array stores the names in a standard order from low (S1) to high (S62)
The numbers of BH's are stored in the array source(0-10). The standard order of these numbers are from 0 to 10.
The VB2019 application supports 3 databases or files. The main reason is to save the results of a specific simulation and to allow to temporarely stop a simulation and to continue at a later date.
What is following is more specific information about the three files used.
- A Command data base. This file contains only one record. The file contains the parameters of a simulation. The two most important parameters are Subtest and test
The two subroutines used are READCOM and WRITECOM. The filenames used end with CMn. The letter n identifies the subtest.
- A DATA file. This file contains only one record. The file contains all the parameters for all the objects, like position and velocity to continue with a simulation exactly in the same state where it was temporarily halted. The two subroutines used are READDATA and WRITEDATA
- A DataBase file. This file contains n records and is updated at regular intervals. The two subroutines used are READDB and WRITEDB
- READCOM en WRITECOM use a BINARY file structure.
At the start of each subroutine the file is opened with: FileOpen(1,filenmcm$,OpenMode.Binary.
To read a record READCOM uses the statement FileGet(1,Comrecord
To write a record WRITECOM uses the statement FilePUT(1,Comrecord
At the end the file is closed with FileClose(1)
- READDATA en WRITEDATA use a Random file structure.
At the start of each subroutine the file is opened with: FileOpen(2,filenmou$,OpenMode.Random,,,ReclenDATA).
To read a record READDATA uses the statement FileGet(2,Datarecord, recordnumberou)
To write a record WRITEDATA uses the statement FilePUT(2,Datarecord, recornumberou)
At the end of each subroutine the file is not closed
- READDB en WRITEDB use a BINARY file structure.
At the start of each subroutine the file is opened with: FileOpen(3,filenmdb$,OpenMode.Binary.
To read a record READDB uses the statement FileGet(3,Bufrecord
To write a record WRITEDB uses the statement FilePUT(3,Bufrecord
At the end the file is closed with FileClose(3)
Angular orbital elements
The background is the book Explanatory supplement to the astronomical almanac paragraph 1.412 Orbital Parameters, page 27, specific Figure 1.412.1 Angular orbital elements.
In this section we use instead of mini BH the word planet.
In Figure 1.412.1 the path of the planet is an ellipse and drawn counter clock wise. This defines the the orbital plane. The Blackhole is at point M. The path results in 4 states:
- State 1 is the point N. This is point when the planet crosses the reference plane. There are two of these points (See state 3). The major characteristic of point N is that the distance towards the BH at point M (object 0) is decreasing.
Point N defines the parameter omega which has a value between 0 and 360 degrees.
- State 2 defines point P when the planet is at Perihelium or shortest distance from the BH
Point P defines the parameter w or the argument of Pericenter, which has a value between 0 and 360 degrees.
Point P also allows the calcultion of the inclination which has a value between -180 and 180
Point P is also used to calculate the minimum distance from the BH.
- State 3 is the point N'. This is point when the planet crosses the reference plane. There are two of these points (See state 1). The major characteristic of point N is that the distance towards the BH at point M (object 0) is increasing.
Point N' is not used to define omega.
- State 4 defines the point when the planet reaches Aphelium or the furtest distance from the BH.
Point P is also used to calculate the maximum distance from the BH.
When both the maximum distance and the minimum distance are calculated the eccentricity can be updated.
id1 a e i (°) omega w (°) Tp (yr) P (yr) Kmag q (AU) v (%c) dv m0
S1 0.5950 0.5560 119.14 342.04 122.30 2001,800 166.0 14.70 2160.7 0.55 0.03 12.40
S2 0.1251 0.8843 133.91 228.07 66.25 2018,379 16.1 13.95 118.4 2.56 0.00 13.6
S4 0.3570 0.3905 80.33 258.84 290.80 1957,400 77.0 14.40 1779.7 0.57 0.01 12.2
S6 0.6574 0.8400 87.24 85.07 116.23 2108,610 192.0 15.40 860.3 0.94 0.00 9.2
S8 0.4047 0.8031 74.37 315.43 346.70 1983,640 92.9 14.50 651.7 1.07 0.0 13.2
S9 0.2724 0.6440 82.41 156.60 150.60 1976,710 51.3 15.10 793.2 0.93 0.02 8.2
S12 0.2987 0.8883 33.56 230.10 317.90 1995,590 58.9 15.50 272.9 1.69 0.01 7.6
S13 0.2641 0.4250 24.70 74.50 245.20 2004,860 49.0 15.80 1242.0 0.69 0.01 10.
S14 0.2863 0.9761 100.59 226.38 334.59 2000,120 55.3 15.70 56.0 3.83 0.06 10.
S62 0.0905 0.9760 72.76 122.61 42.62 2003,330 9.9 16.10 16.4 7.03 0.04 10.
What the result of the simulations show is that the configuration is very stable.
Created: 8 August 2020
Updated: 15 November 2020
Back to my home page: Contents of This Document