This evening's study is returning back to the Schwarzchild solution and thinking about what spatial warping will occur with a whole galaxy full of stars.
The Einstein equations are not linear, so we can't just superimpose solutions and expect it all to work out nicely. I don't know if a transform has ever been worked out to do this; a machine which inputs two relativistic gravitational fields, and outputs their sum.
What I notice about the Schwarzchild star is that the stretching is one dimensional towards the mass, and any significant space-time curvature is highly localised to a few million kms around the star. Beyond that distance, the curvature is very low, so solutions to the Einstein equations may be treated as approximately linear in this region. Meaning that we may superimpose space-warping solutions without incurring much error for modelling the large scale mechanics of stellar clusters, and principally where I'm heading here is the relativistic dynamics of our own galaxy.
The second feature of the Schwarzchild solution is that the time warping apparently mirrors the space warping. I'll maintain this relationship. This figure is a scalar quantity, ie is independent of the direction of the spatial warping. So having superimposed our space warping solutions, we can calculate the total spatial divergence at a point, and invert this figure for the degree of time warping for that region.
At this stage I'm unsure how to combine these three scalar components of space warping: as a root-sum-square, or as a new/old volume ratio. These are both coordinate system independent combinations. I'll think about that later.
My Method:
I'll treat the galaxy as a 2-dimensional disc of mass M, radius R which is composed of a uniform spread of an infinite number of infinitessimal masses dM, each creating a weak relativistic distortion of space and time. For any point along the radius of that disc (distance r from the core), we integrate the sum of these effects to find the total spatial divergence. This will be done numerically for various radii r from 0 upto 3R, and profile graphs created.
I'll move to AI units now.
ds/dr= (1-2Gm/rc^2)^(-1/2)
Let A = 2Gm/c^2. Although compiled with universal constants, this represents the mass of the galaxy.
ds/dr = (1-A/r)^(-1/2)
Far from the star, where A/r is small, a binomial expansion of the derivative will be approximate.
ds/dr ~ 1 + (-1/2).(-A/r) + (-1/2)(-3/2).(-A/r)^2 / 2! + ...
ds/dr ~ 1 + A/2r + O[1/r^2]
This error term is insignificant when r is large, so may be omitted.
So we take ds/dr = 1 + A/2r
The expansion is simply (ds-dr)/dr = A/2r
This is linear in m, so the integration of over the disc area is natural as dM tends to zero.
How on Qo'Nos do I do this integral?

Hmmm..
Right, what is the polar form of an off centre circle? Cosine rule...
H = r.cos(x) + [ R^2 - ( r.sin(x) )^2 ]^(1/2)
INTEGRAL (x=0 to 2.PI) INTEGRAL (h=0 to H) [ A/2h .h] dh.dx
NB: The second h in the integrand represents the arc expansion of polar coordinates
= INTEGRAL (x=0 to 2.PI) [ AH/2 ] dx
For the two component directions, we now introduce component terms into the integrand. These are absolute values, since the stretching is the same whether the star is infront or behind us. I'm assuming that adjacent gravitational fields amplify rather than cancel-out their space warping effects in the null gravity point inbetween them.
Galactic Radial component = A/2 . INTEGRAL (x=0 to 2.PI) [ H.|cos(x)| ] .dx
Galactic Angular component = A/2 . INTEGRAL (x=0 to 2.PI) [ H.|sin(x)| ] .dx
These will obviously reduce to two half integrals.
Let sin(x)=y.R/r
Galactic Radial component = A/2 . INTEGRAL (x=0 to 2.PI) [ H.|cos(x)| ] .dx
= A/2 . {0 + 2.R^2/r .LOOP_INTEGRAL (x=-PI/2 to PI/2) (1-y^2)^(1/2) .dy
= A.R^2/r . [r/R .sin(x). (1-r^2/R^2. sin(x)^2 )^(1/2) + (1/2).INV_SIN(r.sin(x)/R) ] (with x=PI/2)
Letting a=r/R, the proportion of distance to the rim of the galactic disc.
So, the radial component of space-warping at distance r=aR from galactic core is...
= A.R.[ (1-a^2)^(1/2) + INV_SIN(a) / (2a) ]
----------------------------------------------
How nice.

I don't think that R should be there as a product though

maybe A/R but not AR.
Although to be fair, I got through about 4 sheets of paper before finding the right substitution. It's been a while.

This is probably all nonsense anyway.
Now for the angular component:
This time let cos(x)=y.(R^2-r^2)^(1/2)/r
Galactic Angular component = A/2 . INTEGRAL (x=0 to 2.PI) [ H.|sin(x)| ] .dx
= A/2 . {0 + 2.(R^2-r^2)/r .LOOP_INTEGRAL (x=0 to PI) (1+y^2)^(1/2) .dy
= A.(R^2-r^2)/r . [(y/2).(1+y^2)^(1/2) + (1/2).LOG[y + (1+y^2)^(1/2)]] (limits x=0 , x= PI/2)
= A.(R^2-r^2)/r . { a/(1-a^2) + (1/2).LOG [(1+a)/(1-a)] }
So, the angular component of space-warping at distance r=aR from galactic core is...
= A/R . { 1 + (1-a^2)/(2a).LOG [(1+a)/(1-a)] }
----------------------------------------------
That one looks better.
Pictures to follow.
This calculus has been quite fun actually. I'm surprised I still remember it all, so it must be like riding a bike.