PHP Science Labs Wiki

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
This arbitrary-precision PHP function computes the time in seconds that it takes to accelerate from rest to reach a given distance in meters.  It also takes into account the time-dilating affect of general relativity.

<code>
/*
   ===========================================================================
   This function returns the statistics pertaining to accelerating from rest
   until reaching distance (D meters) at a specified rate (m/s^2), taking
   into account the effect of relativistic time dilation in the case of
   extreme values, such as interstellar or intergalactic distances.

   Given SI arguments:
   DistM = Distance accelerated  (meters)
   Accel = Rate of acceleration  (meters/s^2)

   -----------------------------------------------------------------
   The formulas below were derived from basic relativistic geometry.

   Returns string of space-delimited values: 'm a D T A b t'

   Where:
   D = Distance in light seconds equivalent to (DistM).
     = DistM / 299792458

   a = Rate of acceleration (m/s^2) experienced by traveler.

   g = Standard gravitational acceleration rate on Earth
     = 9.80665 m/s^2

   A = Relativistic acceleration potential.
     = a * T / c
     = a * T / 299792458
       If there was no natural speed limit, then this would be the
       the relativistic speed we would achieve after accelerating
       at rate (a) for (T) seconds, which can reach values greater
       than the speed of light, which is 'impossible', hence the
       use of the term 'potential'.

   T = Earth observer (static frame) time interval to reach (D).
     = sqrt(D*D + 2*D/A)
       The Earth is taken to be static relative to the moving traveler.

   b = Relativistic speed (v/c) reached at distance (D).
     = T / (D + 1/A)

   t = Traveler's proper (moving frame) time interval to reach (D).
     = log((1 + b) / sqrt(1 - b*b)) / A
       Where log(x) refers to the natural (Napierian) logarithm.

      (T - t) will always equate to >= 0 due to relativistic time dilation.

       Relativistic time dilation has been experimentally confirmed.

   ERRORS:
   Error message returned on invalid argument.
   ===========================================================================
*/

   function bc_Time_to_Accel_to_Dist ($DistMetersStr, $AccelStr='9.80665')
{
   $Q = 64;          // Set internal working decimals limit.
   $c = '299792458'; // Speed of light m/s

// Read distance in meters and acceleration rate arguments.
   $m = trim($DistMetersStr); // m
   $a = trim($AccelStr);      // m/s^2
   $S = '1';

// Error if invalid argument.
   if (!Is_Numeric($m) or !Is_Numeric($a))
      {return "ERROR:\n\nBoth arguments must be numeric.";}

// Compute light seconds distance from meters.
   $D = bcDiv($m, bcMul($c, $S, $Q),$Q);

// Compute relativistic acceleration potential.
   $A = bcDiv(bcMul($a, $S, $Q), $c, $Q);

   $_D2 = bcMul($D,  $D, $Q);
   $_2D = bcMul('2', $D, $Q);
   $_2D_over_A = bcDiv($_2D, $A, $Q);

// Compute Earth time interval to accelerate
// to distance D light-time units.
   $T = bcSqRt(bcAdd($_D2, $_2D_over_A, $Q),$Q);

// Relativistic speed beta (v/c) achieved at distance (D).
   $b = bcDiv($T, bcAdd($D, bcDiv('1', $A, $Q), $Q),$Q); // beta

   $b2 = bcMul($b,$b, $Q);
   $_1_minus_b2 = bcSub('1', bcMul($b,$b, $Q),$Q);

   $k = bcDiv(bcAdd('1', $b, $Q), bcsqrt($_1_minus_b2, $Q),$Q);

// Compute traveler's proper time to
// reach final distance (D).  This time
// value is limited to double-precision.
   $t  = log($k) / $A;


// Do some output formatting.
   $Q  = 32; // Output decimals limit.
   $m  = RTrim(RTrim(bcAdd('0', $m,  $Q), '0'), '.');
   $a  = RTrim(RTrim(bcAdd('0', $a,  $Q), '0'), '.');
   $D  = RTrim(RTrim(bcAdd('0', $D,  $Q), '0'), '.');
   $T  = RTrim(RTrim(bcAdd('0', $T,  $Q), '0'), '.');
   $A  = RTrim(RTrim(bcAdd('0', $A,  $Q), '0'), '.');
   $b  = RTrim(RTrim(bcAdd('0', $b,  $Q), '0'), '.');
   $t  = RTrim(RTrim(bcAdd('0', $t,  $Q), '0'), '.');

   return "$m $a $D $T $A $b $t";

} // end of  bc_Time_to_Accel_to_Dist (...)


</code>