|
Accuracy of representation of astronomical periods or ratios by gear trains comprises several components. These include:
Accuracy versus precision"Accuracy" is a measure of how well a stated value of a quantity corresponds to the true value of that quantity. "Precision" is usually a measure of the of the repeatability of the values of a measured quantity. It is quite possible for a measurement to be quite precise (that is, successive determinations of a quantity can produce the same number to a large number of significant figures) without it being particularly accurate due to bias(es) in the measurements. These latter can be due to a poorly calibrated measuring device (a meter stick short by a few millimeters, for example) or to wrong or missing corrections to raw measurements (using the wrong temperature or expansion coefficient in accounting for the change in length in a meter stick due to a temperature differing from that for which its length was calibrated). In current standard astronomical practice, the values given for each of the coefficients in the various quantities or ratios which GearTrains attempts to match are given to the number of significant figures in which the last figure is uncertain by an amount less than ±10. "Uncertain", in this instance, is intended to represent the accuracy of the quantity represented, usually based on an elaborate statistical analysis of the measurements made, including estimates of both the repeatability of the measurement device(s), based on an analysis of the measured values, and the accuracy of the measurement device(s), based on an analysis of the design, construction, and conditions of use of these device(s). Accuracy of mathematical expressionsAs given in Seidelmann [1], the current recognized authority and reference for descriptions of motions and aspects of astronomical objects in the solar system, these expressions are usually truncated series expansions in powers of time elapsed after a reference epoch, a Taylor series truncated to a few terms. The time argument, T, is usually given in Julian centuries (see Synodic month for an explanation of this unit of time). The Taylor series were derived by fitting coefficients of their terms to a series of (weighted) observations taken over some period of time. The Taylor series expansions are then used to predict the value of the period or ratio of interest into the future. Astronomical phenomena of interest in horology include
Synodic Month The expression for the length of the synodic month (in days) is:
See Synodic month for the definition of T and further discussion of this expression and its source, Seidelmann [1]. Tropical Year The expression for the length of the tropical year (in days) is:
The tropical year is the time between two successive passages of the Sun through the Vernal Equinox. The Vernal Equinox is a (nearly) fixed point in the sky corresponding to that intersection between the celestial equator (the apparent great circle of intersection between the plane of the Earth’s equator and the apparent celestial sphere) and the ecliptic (the apparent great circle of the Sun‘s annual path through the sky) at which the Sun passes from the southern side of the celestial equator to the northern; the Sun does this each year at the beginning of Spring, on or about 21 March. Again, see Synodic month for the definition of T and further discussion of this expression‘s source, Seidelmann [1]. Ratio of lengths of Mean Solar to Mean Sidereal Day The mean solar day can be thought of as the average over the year of the time interval between two successive passages of the Sun across the visible meridian of a fixed observer on the Earth. The visible portion of an observer‘s meridian is an arc which extends from the north point on his horizon through the point directly overhead (the zenith) to the south point on his horizon. The mean sidereal day is the average time interval of one full rotation of the Earth with respect to the Vernal Equinox. Due to the Earth‘s motion in its orbit about the Sun, the Sun appears to move eastward across the sky with respect to the (nearly) fixed background stars, making one complete turn in a tropical year from the Vernal Equinox back to the Vernal Equinox. The Earth rotates eastward on its axis, so when it has made a full turn with respect to the Vernal Equinox, a mean sidereal day, it will have to rotate just a bit farther to catch the (mean) Sun, a mean solar day. Thus, the mean solar day is slightly longer than the mean sidereal day, about 3m 56s longer. The ratio of lengths of the mean solar to the mean sidereal day (or minutes or seconds) is:
Again, see Synodic month for the definition of T and further discussion of this expression‘s source, Seidelmann [1]. Validity intervals of mathematical expressionsThe bounds of these intervals for expressions (1) through (3) in Accuracy of mathematical expressions are not stated in the current authoritative reference, Seidelmann [1], but may be available in his references for the expressions he provides. We have not as yet examined his references. In the previous edition of the Explanatory Supplement [2] the quantities are represented by expressions (1) through (3):
at the standard epoch of J1900 (1900 January 0.5 = JD 241 5020.0 [JD is the abbreviation for Julian Date. See Synodic month for a discussion of this time scale.]). T19 is the number of Julian centuries after this epoch, which is calculated as:
T19 is analogous to T as discussed in Synodic month, but applies to times with respect to JD 2415020.0, the standard epoch near the beginning of the 20th century, rather than to JD 2451545.0 (2000 January 1.5, J2000), the current standard epoch (see Synodic month for further discussion of this topic). It is interesting to propagate the expressions in Reference [1] to J1900 to compare current values of these quantities with those given in the Reference [2]. Synodic month The value given in expression (4) is the result of direct calculation from the theory of the motion of the Moon current in the mid 1950s. Reference [2] notes that the change in the mean length of the synodic month "does not exceed a few hundredths of a second per century, and depend(s) partly upon the variations in the rate of rotation of the Earth". Probably due to the lack of high stability atomic clocks at that time and the difficulty of obtaining accurate (and precise) measurements of the position of the Moon because it is an extended body (not a point), no time-varying terms are provided in Reference [2] for this value. The appropriate value of T in expressions (1), (2), and (3) for propagating from J2000 to J1900 is T = -1.0. Substituting this value in expression (1) yields the following value:
This propagated value agrees with that given in expression (4) to the number of significant figures given in the latter expression. It is interesting to note how many more significant figures the length of the synodic month can now be determined compared with a half-century ago. Tropical year The theory of the apparent motion of the Sun is far less complex than that for the Moon, and in 1960, at the time of the publication of Reference [2], this theory was substantially unchanged since about 1900. The length of the tropical year at J1900 will be the constant term in expression (5), 365.24219879 days. The value for that epoch calculated from expression (2), again using T = -1.0, is:
This value differs by about 300 in the eighth (last) decimal place in expression (5), which is a surprise. As noted under Accuracy versus precision, it is usual to publish measured values of quantities with the number of significant figures such that the uncertainty in the value of the least significant figure published, here the eighth decimal place, is of the order of ±10. This difference is thirty times that. It is unlikely that improved observational and timing techniques over the past 50 years would have brought such a change. A possible explanation would be that there is a misprint in one or the other of References [1] or [2], but a not-exhaustive investigation failed to reveal any misprint. A mystery. Ratio of the lengths of the mean solar to the mean sidereal day This ratio is also determined from the theory of the apparent motion of the Sun. The value of this ratio at J1900 as given in Reference [2] will be the constant term in expression (6), 1.002737909265. The value for that epoch calculated from expression (3), again using T = -1.0, is:
The constant term in expression (6) has 13 significant figures, while that in expression (3) has 16. Since this propagation was made using a pocket calculator which displays only 12 significant figures, the last digit in expression (10) is uncertain by ±0.5 or so. The agreement in the 12th significant figure is within ±3, which is reasonable. So, this analysis suggests that expressions (1) through (3) are likely to be valid for at least ±100 years from J2000, their common epoch. Another way of estimating the interval of validity of a given series would be to assume that it is unlikely to extend substantially longer into the future (or the past) than the time at which the value of a term of a given power of T has the same magnitude as the corresponding value of the next term in the series, of some higher power of T. On this basis, for expression (1), the magnitude of the terms in T and T2 will be equal in a bit less than 25 Julian centuries, or about 2,500 years. Because the change in the length of the synodic month is substantially linear, it might be safe to assume that expression (1) would be a reasonable representation of the true length of the synodic month for at least 10% of this period, or about 250 years. For expression (2), the magnitude of the terms in T2 and T3 will be equal for T = ±2.76 Julian centuries, or ±276 years from J2000. Because the coefficients of these terms are so nearly equal in magnitude, they are likely not strongly determined from the available data. However, the term in T dominates the higher order terms by a factor of about at least 8,400 in this expression and differs relatively little from that in expression (5) from reference [2]. Thus, the terms in T and T2 will be equal in about 8,400 Julian centuries, or about 840,000 years. All this suggests that expression (2) would likely be useful for a few centuries or so into the future or past. For expression (3), the magnitude of the terms in T and T2 will be equal in almost exactly 10,000 Julian centuries, or about 1 million years. From the discussion above, this would suggest that this ratio is well-determined, and would also likely be useful for a few centuries or so into the future or past. In any case, it should be recognized that expressions (1) through (3) referred to above should be used with caution in trying to determine the value of the period or ratio they represent for times substantially before or after their reference epoch. Also, in any calculation using these expressions to determine the relevant length or ratio, the resulting value should be rounded to the number of significant figures given in the constant term of each expression. Accuracy of calculationsBecause most computer languages and representations of integers and floating point numbers in computer hardware have limits to the size and/or the number of significant figures of representation, it is important to keep this in mind in order to maximize the significance of representation (accuracy) of numerical results. In the case of gear trains, there may be twelve (or even more) gears in a given train, so order of computation of the components of the reduction ratio is critical. The best approach is to use integer multiplication of the tooth counts separately for all of the gears in the numerator and all of the gears in the denominator of the ratio, then to do one floating-point division to determine the gear set ratio (GSR). For a given target ratio (TR) to be matched by a gear train, the fractional Error of this train would be
Note that if GSR and TR are very nearly equal, the number of significant figures in Error will be reduced substantially. The Java language and many computer hardware architectures use the IEEE 754 standard for integers and floating-point numbers. The longest IEEE 754 integer ("long" or "double") is 63 bits plus a sign bit, or
This would seem to be adequate for even a 16-gear train with a maximum tooth count of as much as 225, since
Floating-point division is a bit less accurate. The IEEE 754 floating point number fraction (mantissa) is a 52-bit integer:
Thus, if the accuracy of representation of a gear train is about 1/4.5 x 1015, or about a second in 142 million years, the Error will be lost in roundoff and cannot be calculated. With 12-digit scientific calculators, the Error will be lost in roundoff at about 1 second in 31,000 years, a factor of 4,500 smaller as would be expected. The software used by GearTrains has been carefully designed to minimize roundoff issues by judicious use of Java BigInteger and BigDecimal numbers in its calculations. BigInteger and BigDecimal, properly used, can perform calculations to any number of significant figures desired. References[1] Seidelmann, P. Kenneth "Explanatory Supplement to the Astronomical Almanac", The Nautical Almanac Office [now: Astronomical Applications Department ], U.S. Naval Observatory ; University Science Books , Sausalito, CA 1992 [2] ____, "Explanatory Supplement to the Astronomical Ephemeris and the American Ephemeris and Nautical Almanac", Her Majesty’s Stationery Office, London, 1961 | |||||||||||||||||||||||||||||||||||||||||