Explaining Formation of Astronomical Jets Using Dynamic Universe Model

Astronomical jets are observed from the centres of many Galaxies including our own Milkyway. The formation of such jet is explained using SITA simulations of dynamic universe model. For this purpose the path traced by a test neutron is calculated and depicted using a set up of one densemass of the mass equivalent to mass of Galaxy center, 90 stars with similar masses of stars near Galaxy center, mass equivalents of 23 Globular Cluster groups, 16 Milkyway parts, Andromeda and Triangulum Galaxies at appropriate distances. Five different kinds of theoretical simulations gave positive results. The path travelled by this test neutron was found to be an astronomical jet emerging from Galaxy center.


Introduction
Dynamic Universe model is a singularity free tensor based math model. The tensors used are linear without using any differential or integral equations, there will not be any divided by zero errors, multiple real solutions or any other imaginary solutions for a single set of input data. Only one calculated output set of values exists. Data means properties of each mass like its three dimensional coordinates, velocities, accelerations and it's mass. In this paper, the set of Galaxy center of the Milkyway along with ninety numbers Wolf-Rayet stars is simulated. This whole set is under the continuous and dynamical influence of twenty three globular cluster groups, sixteen Milkyway parts, Andromeda and Triangulum Galaxies. It is a total of 133 masses in this work. This galactic center can be of any another Galaxy also but the outer parts will vary. The dynamical gravitation effect of all this set of masses on a test neutron is calculated and depicted as graphs. Various cases of simulations and their output graphs are discussed in output section. For any N-body problem calculations, the more accurate our input data the better will the calculated results; one should take extreme care, while collecting the input data.
Basic constituent parts of this paper are input data collection, procedure for calculations, various simulations and their output graphs, and lastly discussions/conclusions. In the 'Procedure for calculations' section a discussion can be seen about the total time taken for calculations and number of iterations in each type of simulations.

Input data collection
For conducting these simulations/calculations using a N-body problem solution called dynamic universe model is used. The required real observational data is collected from various sources. That included many research papers and web-pages. All these are referenced. Even Wikipedia also was checked many times for better understanding of some of the practices involved.

Reference to coordinate system
Sun is the reference as usual with it's mass as 1.99E+30 kg, and the Helio Centric Ecliptic Xyz coordinate values of Sun as on 01.01.2000@00.00:00 hrs are Xecliptic=0.0, yecliptic=0.0 and Zecliptic=0.0 Meters. Even though SUN is the reference it will not be used in the calculations as Sun is very far off from Galaxy center.
The stars S1 to S100 will be staying in nearby locations to Galaxy (MILKY Way) center and they will be used in calculations. There is no Coordinate system which is centered on Milkyway center, which would have been a better choice for doing these calculations.

Galactic center
The distance of Galactic Center is of the order of 26000 light-years from solar system. Its direction is having a Right Ascension of 17 h 45 m 40.045 s and Declination of (-) 29°0'27.9". Its mass can be estimated as 4.1 million M ☉ or about 8.2 × 10 36 kg. The orbits of WR star S14 [1], specify that the radius is no more than 6.25 light-hours. In another paper its mass is of the order of 4.31 million M ☉ [2]. Lower value of both the estimates is taken here. The author needs to say one more thing at this point, about the super massive Blackhole [3] at the center of Milkyway according to BBC news. He could not find any further published paper about it on internet. A densemass is selected instead of Black hole. This paper requires only mass and coordinates of this densemass, size is not important for these calculations.

Wolf-Rayet stars (WR stars or Galaxy Center Stars)
Masses of WR stars: Masses of WR stars near Galaxy core are named like S1to S100 etc. Their masses are simulated using the Excel formula This formula generates random numbers between 24 and 16. This generated random number is multiplied with solar mass of 1.99e30 kg to give a star mass which is minimum 16 times to 24 times solar mass [4]. these WR-Stars are simulated using the following Excel formula.
-For x ecliptic the formula applied=$K$31 × (R AND () × (0.00000478802)+1) -For y ecliptic the formula applied=$L$31 × (R AND () × (0.00000478802)+1) -For z ecliptic the formula applied=$M$31 × (R AND () × (0.00000478802)+1) (4) Distances: Distances of the above simulated stars are calculated using above simulated coordinates: -For the Distances the Excel formula used =((K32)^2+(L32)^2+(M32)^2)^0.5 (5) Test Particle… Neutron Initial position and mass of the test particle Neutron: It was chosen as it is electrically neutral. The mass of this was taken from Wikipedia. Position of Neutron was taken near the first Galaxy center star. A value of 1000000000 meters for x, 2000000000 for y and for z 3000000000 meters were added to the coordinates of first star and these values were taken as initial x, y, z positions for the Neutron.
Initial velocity for the Neutron: An initial velocity of 1% of velocity of light was taken. That is V=30000 km/s=30000000 m /sec.

Direction:
Direction of test Neutron is taken towards Milky Way center. The direction cosines are calculated using the formula for equation of line through point 0 0 0 , , x y z and 1 1 1 , , x y z :

Time-step
Time-step is the elapsed time between iterations. Sometimes a time-step of 1 second is given initially to stabilize the system for two iterations. A time step of 100 years was taken in some of these calculations.

Initial conditions
The initial conditions like Mass and Initial distance from Sun in x direction, are shown in Figures 1 and 2. The spike at Mass no.132 indicates Andromeda Galaxy. In Figure 3 the mass distribution taken here for the WR stars is shown.The graphs are provided for visual understanding of the initial conditions for doing these simulations.

Starting conditions used in various files
Starting conditions of Various files used in this paper are given below. Neutron mass and all the other distances were not changed, they were maintained the same in all these theoretical experiments as discussed above.

Mathematical Background
The mathematics of Dynamic Universe is available in many published papers. The following Equation (7) is the basis for many calculations. This basic equation is sufficient for almost all practical purposes ( ) 1 1 This concept can be extended to still higher levels in a similar way.

SITA (Simulation of Inter-Intra-Galaxy Tautness and Attraction forces)
SITA is a totally non-general relativistic algorithm. Here in NO way GR effects are taken into consideration. No space-time continuum. No λ factor to introduce repulsion between Galaxies at any distance. In this SITA Simulation Universe is assumed to be dynamically moving and rotating. This is not a static model as assumed by Newton. Additionally on SITA, a inhomogeneous and anisotropic lumpy universe was assumed.

Procedure
The quest is simple. We need to find out on which input conditions the path traced by neutron will bend perpendicularly. Disappointingly high different numbers of varieties of theoritical experiments with different input conditions were conducted. Using the above defined initial values the SITA calculations procedure gave the data outputs. A large amount of output data was produced. Here different input conditions means different neutron velocities ranging from full light velocity to 1% of velocity of light in different directions. This output data was analysed using graphs and checked results. Mistakes were corrected. Finally we could get results which were thought probably never possible. All the experiments and outputs were logged. Only final cases are discussed in the next section.
Another hurdle faced is the laptop getting heated up. A new computer has to be purchased with a forced cooling fan. For a simple retired person it is a difficult job to spare some money. But any way, each experiment lasted for a duration of 2 to 6 hours of rigorous calculations by the new computer. Power failures combined with UPS (uninterrupted power supply unit) failures caused hevock many times. Everything needed to be restarted

Results from Simulations Abbreviations
In most of the graphs in this document, we will use GC for Galaxy center or Milkyway center, NS for near star or a Star moving close to GC, PM: point mass or mass, the equivalent mass of the body situated at its gravitational center. Pos means the ecliptic Cartesian xyz coordinate position of the point mass or Neutron, and sx, sy and sz denote positions of xyz axes. Unit of all distances is meters. Vel means the ecliptic xyz velocities associated with the point mass at that instant of time and position and vx, vy and vz denote velocities of point mass in xyz axes. Unit of all velocities is meters/second. Another word rev is used for revised, bend is for bending, itr is used for iteration number and perp is used for perpendicular.

Results consolidation:
The details of various final values for the different theoretical experiments are consolidated in the file: … "vvtc Eureka perp results consolidation.xls".
Positive bending results: Using the input data and starting conditions as discussed in the previous section, many theoretical experiments were conducted. Finally there is a ray of hope. The first moon is visible. The galaxy Densemass astronomical jets emerging from the Milkyway centre is coming true. There was a perpendicular movement observed in a xy pos movement graph. (sy vs sx position).

Discussion
Introduction to dynamic universe model Dynamic universe model of cosmology uses tensors without differential and integral equations, and gives unique solutions. SITA is the name of one method of programming used for calculation of tensor evaluations in dynamic universe model. Tensors are generally tough to understand interpret and appreciate. This is mainly because it is the number of equations that each tensor that will be subdivided into. The overall concept is difficult to comprehend. Now let us see what dynamic universe model of cosmology is… It uses tensor mathematics based on Newtonian physics. This mathematics used here is simple and straightforward. All the mathematics and the Excel based software details are explained in the three books published by the author [5][6][7]. In the first book, the solution to N-body problem-called dynamic universe model (SITA) is presented; which is singularity-free, inter-body collision free and      dynamically stable. This is the basic theory of dynamic universe model published in 2010 [5]. The second book in the series describes the equations and SITA software in EXCEL emphasizing the singularity free portions. It explains more than 21,000 different equations [6]. The third book describes the SITA software in EXCEL in the accompanying CD/DVD emphasizing mainly HANDS ON usage of a simplified version in an easy way. The third book contains explanation for 3000 equations instead of earlier 21,000 [7]. With this same SITA setup, many physical problems were solved, which are otherwise not possible. For using this SITA, we have to give the initial values of masses and cartesian X Y Z "coordinates of positions, velocities, and accelerations". Feeding accelerations is not compulsory. Velocities are also not very important, after little iteration of calculations, all the three dimensional Velocities and accelerations will be formed automatically. A point to be noted here is that the dynamic universe model never reduces to general relativity on any condition. The fourth book [8,9] in the series on dynamic universe model: SITA, gave simulations that predicted the existence of the large number of blue-shifted galaxies in 2004, ie., more than about 35~40 blue-shifted galaxies known at the time of Astronomer Edwin Hubble in 1930s. The far greater numbers of blueshifted galaxies was confirmed by the Hubble space telescope (HST) observations in the year 2009. This dynamic universe model can be used for solving general N-body problem. This method solved many unsolved problems earlier like galaxy disk formations, missing mass in galaxies, pioneer anomaly, non-collapsing large scale mass structures, new horizons trajectory predictions etc. We will start with Equation 25, which is a tensor equation, which will be devided into 21000 linear equations on which SITA calculations work. This is the main equation giving many results, that are not possible otherwise today. Now it attempts to solve new problems like a. Variable mass rocket trajectory problem.
c. Astronomical jets observed from Milkyway Center.
d. Dynamic universe model's prediction "no dark matter" in the universe came true! [11].
To support dynamic universe model the we can find the following supporting observations. a. Mathematical prediction of existence of blue shifted galaxies came true.
c. The first redshifted Quasar 3C273 is blue shifted (supporting existence of blue shifted galaxies). d. The most distant Quasar eso1122 found to have a blue Shift of 0.110473 (supporting existence of blue shifted galaxies).

Solving new problems
Variable mass rocket trajectory problem: When the rocket is moving due to thrust of the fuel in it, its mass will continuously reduce. The positions and gravitation forces of Moon, Sun, Earth and other planets will dynamically change. Dynamic Universe Model predicts such trajectory continuously. This will be very useful for optimizing the overall mission efficiency for reducing the thruster rocket sizes and increasing payload capacity of rocket.

Explaining large variation in the gravitational bending results of VLBI:
In this, the effect of universal gravitational force is calculated on a radio photon by using a singularity free and collision free N-body problem solution called dynamic universe model. Here the capabilities of this dynamic universe model are extended into micro world i.e. to light photons and radio wavelength photons and neutrinos etc. By doing so a real world very long baseline interferometry (VLBI) observations   are explained [10]. The VLBI techniques give gravitational bending results in a wide range of values. Now dynamic universe model explains reason for such variation.The basic difference is that where present day physics considers gravitation effects of only sun or the main gravitating body only, the dynamic universe model considers gravitational effect of sun, planets, globular clusters, milky-way, local systems etc., and finds the universal gravitational force vector at that instant of time for that configuration of the Universe.
Can the gravitational effect of sun and moon be neglected on me when i am standing on earth? No. For example, tide caused by sun and moon in oceans--We observe high tide and low tide in the mornings and evenings, or on full-moon-day and no-moon-day. These tides are caused by gravitation of sun and moon only. So we can not neglect gravitation effect of sun and moon on earth. For better accuracies we have to consider planets also….
There is a large variation in the gravitational bending results of very long baseline interferometry (VLBI) , in the field of radio astronomical observations of quasars, galaxies etc. This variation is clearly visible when the solar gravitational bending/deflection angle is plotted against solar elongation angle [11][12][13][14].
Astronomical jets observed from center of disk galaxies like milkyway: Many galaxies show the formation of astronomical jets from its center including our milkyway. These set of N-body simulations show the root cause of these formations.

Observations that support dynamic universe model
Mathematical prediction of existence of blue shifted galaxies: The forecasts of dynamic universe model in 2004 came true observationally in 2009 [8]. Dynamic universe model calculated to anticipate the existence of a large number of blue shifted galaxies which are much more higher than the 40 numbers known from the time of Astronomer Hubble in 1930s. It was confirmed by Hubble space telescope (HST) observations in the year 2009. Today the HST observed more than 10000 directly. In addition one can safely assume and prove various other types of galaxies such as Quasars, UV galaxies, X-ray, γ-Ray sources and other blue galaxies etc., are also blue shifted galaxies. That means about 31.7% of galaxy count are blue shifted. One should not neglect such large number of blue shifted galaxies. It appears it is a Godly devotion to Bigbang cosmologies! The author submitted a paper to PRD in 2004.

SN1987A-Neutrino emission:
Many flavours of Neutrinos are generated from sun and stars. One will find other flavours are generated from SN1987A. This covers the whole spectrum of Neutrinos. There are no unkown additional flavours generated from Bigbang.

Other cosmologies-comparison
Newton's static universe model requires fine balancing of bodies in all directions, so that all bodies stay in static equilibrium of attraction forces. This was described as such equilibrium as though a set of needles is finely balancing on their noses, any small disturbance will cause all to fall. Here in our dynamic universe model, gravitational attraction forces are balanced, by centrifugal forces not by balancing attraction forces. SITA proves that bodies will not collapse but revolve about each other. Dynamic universe model will not have Big-bang singularity, as we are proposing a nonexpanding anisotropic and heterogeneous universe model without considering the general relativity. This is a dynamic universe model without space-time continuum. No Big bang singularity. Hence singularity theorem is not applicable here. Hawking and Penrose in their singularity theorem said that ' Isotropic and homogeneous expanding universe, there must be a Big bang singularity some time in the past according to general theory of relativity. PCP was not considered true here as in steady state universe we need not assume any homogeneity and isotropy here at any point of time. This is a non expanding universe and matter need not be created to keep the density constant. The steady state cosmological model was presented by Hoyle The perfect cosmological principle (PCP) stated by Hoyle is that, Isotropy and homogeneity and other statistical properties of the universe are time independent. Universe has no beginning. No starting point for time scale. Matter is required to be created to keep the density constant in the expanding universe. { In a recent paper Aguirre and Gratton time like geodesic are not complete in Hoyle's Steady-state model. They proposed a geodesically complete Steady-state model, in which two universes are simultaneously present. In one of them, the universe is expanding and time is moving forwards, and in the other, it is contracting and time is moving backwards}. Friedmann-Robertson-Walker models are popular. These are standard Bigbang models. Naturally all the problems inherent in the Bigbang models are present here also. In the absence of other working cosmological models, many workers choose these next. Missing mass, lesser age of the universe, anisotropy of cosmic microwave background, Bigbang singularity etc., are some of the problems present in these models. Bowen and Ferreira said, In models by de Sitter or any other matter filled models, there will be mass loss by scalar charges in these types of expanding universe models. That means a point like particle carrying charge q, acts like a source for mass less scalar field. It looses its mass in time. There is one more popular line of thought, which is being seen now a days. They are cyclic universe models presented by many workers. We will see some the recent work done by Steinhardt and Turok, in which the universe starts from Bigbang to end up in Big-crunch only to start again in Bigbang to start the cycle. They attempted to say a little about, what happened before Bigbang. Hawking and Penrose for detailed work Hawking and Ellis in their singularity theorem, showed that Bigcrunch heads towards a cosmic singularity, where General relativity fails. After big crunch what happens, nobody knows. There is a basic problem in all these models, including String theory and M-theory; the matter density is significantly low, which makes these models impractical. In these models the universe is flat but not closed. So the question comes what happens to all these radiation? Steinhardt and Turok, presented another model of cyclic universe, to overcome the problem of failure of general relativity after Big-Crunch. They pushed the Big-crunch singularity into 5 th dimension, so that other three spatial and one time coordinates will be intact. It may be Steady state model or cyclic universe model; one thing is there in common. Both types of models ask for the creation of matter from vacuum. Earlier on this point the Bigbang people were criticizing the Steady state people. Now let's see about rotation models presented various authors from Gödel to Korotky and Obukhov.There were many authors. Gödel metric described the solution of General relativity with homogeneous spacetime and with casualty condition violated. All these people gave mainly a line element as a solution to Einstein's General relativity and tested that solution. Nobody talked about revolution. Mainly they argued about the rotation of universe, saying "when everything rotating, why not universe also?" But they have not considered the revolution of parts of the universe. Another difficulty faced by Korotky and Obukhov, is that it is impossible to combine pure rotation with expansion of universe in a solution of general relativity for a pure simple source. There were many authors who faced problems like closed time like curves (CTCs) Obukhov Yu N and Saulo Carneiro. The problems like non linearity of coordinate axes and interdependency between coordinate axes is still present inherently in all these models.
There is a fundamental difference between galaxie/systems of galaxies and systems that normally use statistical mechanics, such as molecules in a box. The molecules repel each other. But in gravitation we have not yet experienced any repulsive forces. Binny and Tremaine only attraction forces were seen. Einstein introduced cosmological constant λ to introduce repulsive forces at large scales like inter galactic distances in his General relativity based cosmological considerations in for expanding universe 1917. This was not liked by many, and created turbulence in the scientific world. One of the reasons for his cosmological constant λ is that he disliked the picture at infinity given by Newtonian gravitation. Though his ideas about infinity were good, the cosmological constant λ and repulsive forces created havoc in the scientific community for at least last hundred years! Almost every worker/scientist in this field faced problems either conceptually or mathematically. Singularities were big hurdles for many of us.

Dynamic Universe Model: Evidences
Presence of blue shifted galaxies in the universe, is the main evidence. Hubble deep space houses thousands of blue shifted galaxies which is one of the greatest mysteries for expanding universe models could not explain.
Our galaxy the Milky way is moving with a speed 454 ± 125 km/sec towards l=63 ± 15° and b=-11 ± 14°relative to distant part of samples and 474 ± 164 km/sec towards l=167 ± 20° and b=5 ± 20° relative to nearer part of samples. The local group comprising of Milky way, NGC6822, Andromida galaxy and other dwarf elliptical galaxies, Magellanic clouds rotate about their centers and revolve around a common center. SM Faber and David Burstain in their paper " Motions of galaxies in the neighborhood of Local group "{presented in a symposium,' Large scale motions of universe' Princeton described the streeming motions towards the Great Attractor (located at l=309 and b=+18) by the local group, Virgo cluster, Ursa major, Centaurus, Camelopardalis, Perseus-Pisces etc ., clusters with speeds ranging up to 1000km/sec. Please note the difference in directions of movement as well as speeds. All these clusters form a super cluster which also rotate and revolve about each other. Groups of super clusters form Filament structures and to grate walls and so on. This is how our universe is Lumpy and anisotropic even at large scale.
Another piece of supporting evidence for the dynamic universe model was there. There is a considerable discussion was as to whether GA: the Great attractor exists at all. For example Mathewson, Ford and Buckhorn have measured the peculiar velocities 1355 spiral galaxies. They find no backside in fall into GA region, rather a bulk flow of about 400 km/sec on the scales of 100 ho-1 MPC. Thus there is a considerable doubt about the existing of an attracting mass there. Both the parties find streaming motions or bulk flow. If there is no attracting mass, then why they are moving? this super cluster must be in revolution motion. Birch has discovered the asymmetric distribution of the angles of rotation of polarization vectors of 132 radio sources and tried to explain this via the global rotation. We think that the asymmetric distribution of the angles of rotation of polarization vectors is due to the galaxies or parts of clusters revolving in different directions.

Conclusions
Astronomical jets emerging from galaxy center of milkyway or any other galaxy are simulated here using dynamic universe model in this paper. This N-body model consists of 133 bodies of point masses, viz., 90 Wolf-Rayet stars, Galaxy Center, 23 Globular Cluster groups, 16 Milkyway parts, Andromeda and Triangulum Galaxies, and a Test neutron the path which was traced.
Lots of theoretical experiments were conducted with many types of assumed initial conditions, which did not give expected bending results and were discarded. That means only in certain conditions astronomical jets are formed as shown in the results section.
In all these simulations presented in this paper the time step is 1 sec, sharp bending of neutron path either after 1030 or 1650 iterations, is visible in x and y axes graphs only, and z axis don't show any bending. Another important observation is either Galaxy core mass was increased or near star masses were reduced slightly. The directions of the neutron initial velocities are slightly different from the first simulation to the others.
Probably this may be too early to conclude from these N-body simulation experiments. So our initial proposition is correct. In some conditions as discussed above, these bodies produce such gravitational effects on the neutrons travelling in between these bodies, these neutrons will take such will take some paths which are perpendicular to the central plane of Galaxy near the Galaxy center. Though it is heavier for my purse, purchase of new computer became necessary to take the continual computational load [21,22].