Calculation of a symmetric Gnomon that approximately corrects the Equation of Time, 30.3. 2009 1
The Gnomon of a sundial can be shaped such that the Equation of Time is automatically corrected and the Sundial shows the Civil Time. Because the sun has the same declination twice a year, but the Equation of Time is different at these dates, the gnomon has to be changed twice a year, like it is e.g. the case for the dial of Martin Bernhardt: http://www.praezisions-sonnenuhr.de/ In order to avoid changing the gnomon twice a year I developed a semi transparent double gnomon that incorporates both: http://riegler.home.cern.ch/riegler/sundial/mainpage.htm The absolute value of the equation of time at a given declination is not very different, so by approximating the analemma with a symmetric curve one can arrive at a single Gnomon that doesn t have to be changed. The error introduced by this approximation is at maximum 1.7 minutes. The calculation of this Gnomon is described in the next pages. 2
Equation of Time (Degrees) Equation of Time (Minutes) Solar Declination (Degrees) Solar Declination Days after December 21 st Equation of Time (EoT): Civil Time = Sundial Time + EoT Days after December 21 st Equation of Time in Degrees =EOT (Minutes)/4 Days after December 21 st 3
Equation of Time Error (Minutes) Equation of Time (Degrees) Equation of Time (Degrees) EoT vs. Declination Dec. 21 to Jun. 21 Jun. 21 to Dec. 21 Solar Declination (Degrees) EoT vs. Declination Averaged to make it symmetric Solar Declination (Degrees) Error from Symmetrization Maximum 1.7 Min. around June 21 Solar Declination (Degrees) 4
Symmetrized Equation of Time τ (Degrees) Using the symmetric Approximation the Equation of Time τ is now a unique function of the Solar Declination δ. We call the function absolute value of the approximated equation of Time τ (δ). Solar Declination δ (Degrees) 5
Calculation of the Gnomon: At 12:00 + EoT the Sun is exactly in southern direction A Ray connected to Point P and pointing South with a declination of δ must touch the Gnomon. This defines a ray for every declination δ y South R x P x=-r*sin(τ (δ)) y=-r*cos(τ (δ)) 12:00 EoT (degrees) 6
Parametric Representation of the sunrays: South x=- R*Sin(τ (δ)) + k*0 y=- R*Cos(τ (δ)) + k*1 z= 0 + k*tan(δ) δ =-23.45 to 23.45 k= 0 to 1500 mm Magnifications: 7
Sunrays: x (k, δ)=- R*Sin(τ (δ)) + k*0 y (k, δ) =- R*Cos(τ (δ)) + k*1 z (k, δ) = 0 +k*tan(δ) δ =-23.45 to 23.45 k= to 1500 mm Gnomon: The Gnomon must be the rotationally symmetric body that touches the surface of sun rays. How can we calculate this body? To find the radius of the Gnomon at height z0 we cut the sunrays with a plane at z=z0 which gives a 1-dimensional curve. Example z0=150mm Then we find the circle centered at zero that touches this curve. The radius r0 of this circle is the radius of the Gnomon at height z0. The parametrization of this curve at z=z0 is: z0= k*tan(δ) k=z0/ Tan(δ) X(δ)=- R*Sin(τ (δ)) y (δ) =- R*Cos(τ (δ)) + z0/ Tan(δ) δ =-23.45 to 23.45 To find the circle, centered at zero, that touches this curve, we simply have to find the minimum distance between the point (0,0) and the curve, i.e. r(δ) 2 = x(δ) 2 +y(δ) 2 = = R 2 *Sin(τ (δ)) 2 +[R*Cos(τ (δ)) -z0/tan(δ)] 2 Minimum r0 8
On the following pages the touching circles at different heights z0 above the equatorial plane are shown. The radius R of the dial is assumed to be 500mm. 9
Z0=1mm r(δ) δ 10
Z0=50mm r(δ) δ 11
Z0=100mm r(δ) δ 12
Z0=150mm r(δ) δ 13
Z0=200mm r(δ) δ 14
Z0=216mm r(δ) δ 15
Z0=-50mm r(δ) δ 16
Z0=-100mm r(δ) δ 17
Z0=-150mm r(δ) δ 18
Z0=-200mm r(δ) δ 19
Z0=-213mm r(δ) δ 20
Above and below a certain limiting declination there is no solution to the minimization problem. The reason is that the correction for one day is shadowing the correction of the next day, so the correcting Gnomon can only exist in a certain range of declinations. The solar declination ranges from -23.45 to 23.45 degrees. The rage where a solution for the correcting Gnomon can be found is -23.30 to 23.42 degrees. This means that only ±10 days around Decmeber 21 st and ±2 days around June 21 st the indicator isn t precise. The optimum shape to minimize the error during these days is to simply cut the gnomon at this height. 21
Limit at positive declination: 23.42 degrees, z0=216mm r(δ) Z0=-213mm δ Limit at negative declination: 23.3 degrees, z0=-213mm r(δ) δ The total height of the gnomon is therefore 213+216 = 429mm 22
Final Gnomon 23
Final Gnomon 24
Final Gnomon z 0 (mm) r 0 (mm) z 0 (mm) r 0 (mm) z 0 (mm) r 0 (mm) z 0 (mm) r 0 (mm) -213. 8.49-212. 9.98-211. 11.31-210. 12.46-209. 13.48-208. 14.44-207. 15.32-206. 16.14-205. 16.9-204. 17.6-203. 18.27-202. 18.91-201. 19.52-200. 20.1-199. 20.65-198. 21.17-197. 21.66-196. 22.14-195. 22.6-194. 23.06-193. 23.49-192. 23.91-191. 24.31-190. 24.7-189. 25.08-188. 25.44-187. 25.78-186. 26.11-185. 26.43-184. 26.73-183. 27.03-182. 27.32-181. 27.6-180. 27.87-179. 28.14-178. 28.38-177. 28.62-176. 28.85-175. 29.08-174. 29.3-173. 29.52-172. 29.73-171. 29.92-170. 30.12-169. 30.29-168. 30.47-167. 30.64-166. 30.8-165. 30.97-164. 31.11-163. 31.26-162. 31.40-161. 31.54-160. 31.67-159. 31.8-158. 31.92-157. 32.03-156. 32.13-155. 32.24-154. 32.34-153. 32.43-152. 32.52-151. 32.61-150. 32.69-149. 32.77-148. 32.84-147. 32.91-146. 32.97-145. 33.02-144. 33.08-143. 33.12-142. 33.17-141. 33.21-140. 33.25-139. 33.28-138. 33.31-137. 33.34-136. 33.36-135. 33.38-134. 33.4-133. 33.41-132. 33.42-131. 33.43-130. 33.43-129. 33.43-128. 33.42-127. 33.42-126. 33.41-125. 33.39-124. 33.38-123. 33.36-122. 33.34-121. 33.31-120. 33.29-119. 33.25-118. 33.22-117. 33.19-116. 33.15-115. 33.1-114. 33.06-113. 33.01-112. 32.96-111. 32.9-110. 32.84-109. 32.79-108. 32.73-107. 32.66-106. 32.6-105. 32.53-104. 32.47-103. 32.39-102. 32.32-101. 32.25-100. 32.17-99. 32.09-98. 32. -97. 31.91-96. 31.82-95. 31.73-94. 31.64-93. 31.55-92. 31.45-91. 31.35-90. 31.26-89. 31.16-88. 31.05-87. 30.94-86. 30.83-85. 30.72-84. 30.61-83. 30.49-82. 30.38-81. 30.26-80. 30.14-79. 30.02-78. 29.900-77. 29.78-76. 29.65-75. 29.52-74. 29.39-73. 29.26-72. 29.13-71. 28.99-70. 28.86-69. 28.72-68. 28.58-67. 28.44-66. 28.29-65. 28.15-64. 28.01-63. 27.86-62. 27.71-61. 27.56-60. 27.40-59. 27.25-58. 27.09-57. 26.93-56. 26.78-55. 26.62-54. 26.46-53. 26.31-52. 26.15-51. 25.99-50. 25.82-49. 25.65-48. 25.48-47. 25.31-46. 25.14-45. 24.97-44. 24.79-43. 24.62-42. 24.44-41. 24.27-40. 24.09-39. 23.91-38. 23.73-37. 23.54-36. 23.36-35. 23.18-34. 23. -33. 22.81-32. 22.63-31. 22.45-30. 22.27-29. 22.08-28. 21.90-27. 21.71-26. 21.52-25. 21.32-24. 21.12-23. 20.93-22. 20.73-21. 20.54-20. 20.34-19. 20.15-18. 19.95-17. 19.75-16. 19.56-15. 19.36-14. 19.16-13. 18.96-12. 18.76-11. 18.56-10. 18.36-9. 18.16-8. 17.95-7. 17.75-6. 17.54-5. 17.34-4. 17.13-3. 16.93-2. 16.72-1. 16.52 25
Final Gnomon z 0 (mm) r 0 (mm) z 0 (mm) r 0 (mm) z 0 (mm) r 0 (mm) z 0 (mm) r 0 (mm) 1. 16.11 2. 15.9 3. 15.69 4. 15.48 5. 15.27 6. 15.06 7. 14.85 8. 14.65 9. 14.44 10. 14.24 11. 14.03 12. 13.82 13. 13.61 14. 13.4 15. 13.19 16. 12.98 17. 12.77 18. 12.55 19. 12.34 20. 12.13 21. 11.92 22. 11.71 23. 11.5 24. 11.29 25. 11.08 26. 10.87 27. 10.66 28. 10.45 29. 10.24 30. 10.03 31. 9.82 32. 9.6 33. 9.39 34. 9.18 35. 8.97 36. 8.76 37. 8.55 38. 8.34 39. 8.12 40. 7.91 41. 7.7 42. 7.49 43. 7.28 44. 7.06 45. 6.86 46. 6.65 47. 6.44 48. 6.23 49. 6.02 50. 5.81 51. 5.61 52. 5.4 53. 5.19 54. 4.98 55. 4.78 56. 4.57 57. 4.37 58. 4.16 59. 3.96 60. 3.75 61. 3.55 62. 3.35 63. 3.15 64. 2.95 65. 2.75 66. 2.54 67. 2.34 68. 2.14 69. 1.94 70. 1.75 71. 1.55 72. 1.35 73. 1.15 74. 0.96 75. 0.76 76. 0.57 77. 0.38 78. 0.19 79. 0 80. 0.19 81. 0.37 82. 0.56 83. 0.74 84. 0.92 85. 1.1 86. 1.29 87. 1.48 88. 1.66 89. 1.86 90. 2.05 91. 2.23 92. 2.41 93. 2.58 94. 2.76 95. 2.92 96. 3.09 97. 3.26 98. 3.44 99. 3.61 100. 3.78 101. 3.95 102. 4.12 103. 4.28 104. 4.44 105. 4.60 106. 4.76 107. 4.92 108. 5.08 109. 5.24 110. 5.39 111. 5.54 112. 5.7 113. 5.84 114. 5.99 115. 6.14 116. 6.28 117. 6.43 118. 6.57 119. 6.71 120. 6.85 121. 6.98 122. 7.12 123. 7.25 124. 7.39 125. 7.52 126. 7.64 127. 7.77 128. 7.9 129. 8.02 130. 8.14 131. 8.26 132. 8.38 133. 8.5 134. 8.62 135. 8.74 136. 8.85 137. 8.95 138. 9.06 139. 9.16 140. 9.26 141. 9.36 142. 9.45 143. 9.55 144. 9.64 145. 9.73 146. 9.82 147. 9.9 148. 9.99 149. 10.07 150. 10.15 151. 10.22 152. 10.29 153. 10.36 154. 10.42 155. 10.49 156. 10.55 157. 10.6 158. 10.66 159. 10.71 160. 10.76 161. 10.8 162. 10.85 163. 10.89 164. 10.93 165. 10.96 166. 10.98 167. 11. 168. 11.02 169. 11.04 170. 11.05 171. 11.07 172. 11.07 173. 11.07 174. 11.07 175. 11.06 176. 11.05 177. 11.04 178. 11.02 179. 10.99 180. 10.96 181. 10.92 182. 10.87 183. 10.83 184. 10.77 185. 10.72 186. 10.65 187. 10.58 188. 10.5 189. 10.41 190. 10.32 191. 10.22 192. 10.11 193. 9.99 194. 9.87 195. 9.73 196. 9.6 197. 9.44 198. 9.27 199. 9.1 200. 8.92 201. 8.73 202. 8.53 203. 8.3 204. 8.05 205. 7.79 206. 7.51 207. 7.2 208. 6.87 209. 6.52 210. 6.12 211. 5.67 212. 5.18 213. 4.60 214. 3.94 215. 3.09 216. 1.86