地基大口径望远镜圆顶设计

地基大口径望远镜圆顶设计 Dome design of large ground-based telescope 1 first-author 刘祥意（1987-），男，吉林敦化人，博士，副研究员，2012年于哈尔滨工业大学获得硕士学位，2015年于中国科学院长春光机所获得博士学位，主要从事地基大口径望远镜结构设计及热控方面的研究。E-mail： liuxiangyi107@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921195&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921199&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921197&type= 22.01332855 32.00400543 刘祥意 （1987-），男，吉林敦化人，博士，副研究员，2012年于哈尔滨工业大学获得硕士学位，2015年于中国科学院长春光机所获得博士学位，主要从事地基大口径望远镜结构设计及热控方面的研究。E-mail： liuxiangyi107@163.com liuxiangyi107@163.com 1 2 corresp E-mail： zcwang911@163.com E-mail： zcwang911@163.com 王志臣（1980-），男，黑龙江大兴安岭人，副研究员，2003年、2007年于哈尔滨工业大学分别获得学士和硕士学位，主要从事地基大口径望远镜光机结构方面的研究。E-mail： zcwang911@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921200&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921202&type= 22.01332855 32.00400543 王志臣 （1980-），男，黑龙江大兴安岭人，副研究员，2003年、2007年于哈尔滨工业大学分别获得学士和硕士学位，主要从事地基大口径望远镜光机结构方面的研究。E-mail： zcwang911@163.com zcwang911@163.com 1 corresp 中国科学院 长春光学精密机械与物理研究所，吉林 长春130033 中国科学院 长春光学精密机械与物理研究所，吉林 长春130033 Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China 中国科学院大学，北京 100049 中国科学院大学，北京 100049 University of Chinese Academy of Sciences， Beijing 100049， China University of Chinese Academy of Sciences， Beijing 100049， China 圆顶视宁度是影响地基大口径望远镜成像质量重要的影响因素之一。视宁度产生的原因主要是由于空气折射率的不均匀变化导致，首先通过理论分析确定影响空气折射率的主要因素为空气的温度、速度和压力，然后借助CFD流体仿真软件，得到望远镜圆顶周围附近空气的温度场、速度场和压力场的变化，进而推算出望远镜成像光路中圆顶及望远镜附近周围空气折射率的变化情况，以评估不同结构形式圆顶所产生的圆顶视宁度对望远镜成像的影响。通过仿真分析可知，升降式圆顶更利于望远镜周围空气的温度、速度和压力的快速变化及稳定，能够减少圆顶视宁度的影响，为4 m级及更大口径望远镜的圆顶方案选择及设计提供指导。 The dome seeing is an important influencing factor for imaging quality of ground-based large-aperture telescopes. Seeing is mainly due to the uneven change in the refractive index of air. First， theoretical analysis shows that the main factors affecting the refractive index of air are the temperature， velocity， and pressure of the air. Then， with the help of CFD fluid simulation software， the temperature field， velocity field， and pressure field changes of the air near the telescope dome are obtained， and the change of the refractive index of the air around the dome and telescope is calculated. Thus， the effect of dome seeing produced by domes of different structural forms on telescopic imaging can be evaluated. The simulation analysis shows that the lifting dome is more conducive to the rapid change and stability of the temperature， speed， and pressure of the air around the telescope to reduce the impact of dome seeing. From the perspective of dome seeing， this paper provides a guide for the scheme selection and design of a dome at the 4-meter level and larger telescopes. 望远镜 地基大口径 视宁度 圆顶 telescope large aperture seeing dome 1 引 言 1 地基大口径望远镜作为造价不菲精密的天文观测与空间探测仪器，其安装及工作环境的要求很高。首先体现在对望远镜站址选择上，站址一般选择在全年气候干燥少雨的高海拔地区，这样既能增加望远镜的全年工作时间，也能减少地表大气对流层对望远镜观测的影响，获得较好的大气相干长度。一般而言，优良望远镜站址的大气相干长度为10~30 cm。其次体现在望远镜对站址内塔台等人为建筑的建设上，不仅要求具有较高的抗地质灾害、抗暴风、暴雨和暴雪等气象灾害的能力，而且随着人们对视宁度认识的提高，还要求尽量较少望远镜的塔台所导致视宁度对望远镜观测的影响。塔台圆顶作为直接储存望远镜的建筑部分，不仅需要在强风、强雨雪等恶劣气象条件下对望远镜起到保护作用，而且需要在望远镜工作状态下尽量不影响观测 ［ 1 1 - 2 2 1-2 ］ 。 国外的大口径望远镜在设计建设初期，除了要对望远镜的备选观测站址的气象条件、大气视宁度等进行长时间的统计研究，还要综合考虑望远镜工作条件、资金和自然环境等因素，因此对望远镜塔台及圆顶的结构形式要进行深入的分析。例如，口径为3.5 m的GREGOR望远镜的圆顶为半球对开式升降圆顶 ［ 3 3 ］ ，口径为3.67 m的AEOS望远镜的圆顶为多段升降旋开式圆顶 ［ 4 4 ］ ，口径为3.5 m的WIYN望远镜的圆顶为多边形敞开式随动圆顶 ［ 5 5 ］ ，为口径为4.2 m的SOAR望远镜的圆顶为5/8球形结构的敞开式随动圆顶 ［ 6 6 ］ 。 我国的大口径望远镜研究起步较晚，技术相对落后，目前国内已建成的单镜主镜口径最大口径为4 m级 ［ 7 7 - 9 9 7-9 ］ ，相对于国外的单镜口径8 m级望远镜 ［ 10 10 ］ 还有一定的技术差距。国内1 m级望远镜的圆顶形式主要有两种，兴隆基地1 m望远镜采用的是随动圆顶，云南天文台NVST望远镜采用敞开式圆顶，而4 m望远镜的圆顶形式还需要进一步的探索研究。本文针对4 m口径的大口径望远镜圆顶的设计，分别分析了升降式圆顶与敞开式随动圆顶两种结构形式下视宁度的影响，为4 m望远镜的圆顶设计提供指导。 2 视宁度的影响 1 视宁度对望远镜成像质量的影响主要是光程上空气折射率的波动引起的 ［ 11 11 - 12 12 11-12 ］ 。根据Gladstone-Dale公式可知 ［ 13 13 - 14 14 13-14 ］ ，空气的折射率可以由空气密度 <math id="M1"><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921207&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921205&type= 2.11666679 3.64066648 表示，即： <math display="block" id="M2"><mi>N</mi><mo>=</mo><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">G</mi><mi mathvariant="normal">D</mi></mrow></msub><mi>ρ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921211&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921208&type= 22.77533340 4.57200003 ， 其中 <math id="M3"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">G</mi><mi mathvariant="normal">D</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921215&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921213&type= 6.01133299 3.72533321 为Gladstone-Dale常数，其大小取决于望远镜的观测波长 <math id="M4"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921217&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921216&type= 1.77800000 2.79399991 （μm），表达式如下： <math display="block" id="M5"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">G</mi><mi mathvariant="normal">D</mi></mrow></msub><mo>=</mo><mn mathvariant="normal">2.33</mn><mo>×</mo><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">4</mn></mrow></msup><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><mn mathvariant="normal">7.52</mn><mo>×</mo><msup><mrow><mn mathvariant="normal">10</mn></mrow><mrow><mo>-</mo><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><msup><mrow><mi>λ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921219&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921218&type= 59.60533142 9.56733322 . 通常而言，望远镜的观测波段为0.1~4 μm，则 <math id="M6"><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">G</mi><mi mathvariant="normal">D</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921221&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921220&type= 6.01133299 3.72533321 约为2.23×10 - 4 。 一般而言，望远镜及其圆顶周围的风速不会超过20 m/s，小于0.3 Ma，空气可以视为不可压缩的气体，所以空气的压强为动压与静压之和： <math display="block" id="M7"><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mi>p</mi><mo>+</mo><mn mathvariant="normal">0.5</mn><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921223&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921222&type= 28.02466965 4.82600021 ， <math display="block" id="M8"><mi>V</mi><mo>=</mo><mroot><mrow><msup><mrow><mi>u</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mi>w</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921225&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921224&type= 30.90333366 4.40266657 ， 其中： <math id="M9"><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921228&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921227&type= 3.21733332 4.57200003 为标准状态下的空气密度， u ， v ， w 分为风速三个正交矢量方向的速度。 根据静态方程 <math id="M10"><mi>p</mi><mo>≅</mo><mi>ρ</mi><mi>R</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921232&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921230&type= 15.32466698 4.57200003 和 式（3） 式（3） ，空气的密度可以表示为： <math display="block" id="M11"><mi>ρ</mi><mo>≅</mo><mfrac><mrow><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>R</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>R</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><msup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921237&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921235&type= 32.51200104 8.38199997 ， 其中： <math id="M12"><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921239&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921238&type= 3.72533321 3.72533321 为空气的温度， R 为空气常数。 根据式（1）~ 式（5） 式（5） ，空气折射率的计算公式为： <math display="block" id="M13"><mi>N</mi><mo>≅</mo><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">G</mi><mi mathvariant="normal">D</mi></mrow></msub><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>R</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>R</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><msup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mo>=</mo><mn mathvariant="normal">1</mn><mo>+</mo><mi mathvariant="normal">Δ</mi><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921241&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921240&type= 66.71733856 9.22866726 . 由 式（6） 式（6） 可知，相对于真空条件下的折射率，光在空气中的折射率取决于 <math id="M14"><mi mathvariant="normal">Δ</mi><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921245&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921244&type= 5.67266655 2.87866688 ，而 <math id="M15"><mi mathvariant="normal">Δ</mi><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921245&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921244&type= 5.67266655 2.87866688 主要与空气的速度、温度和压力有关。因此，光程上空气折射率的变化率的一致性由空气的速度、温度和压力的变化率决定，即： <math display="block" id="M16"><mi mathvariant="normal">d</mi><mi>n</mi><mo>=</mo><mfrac><mrow><msub><mrow><mi>K</mi></mrow><mrow><mi mathvariant="normal">G</mi><mi mathvariant="normal">D</mi></mrow></msub></mrow><mrow><mi>R</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo stretchy="false">[</mo><mi mathvariant="normal">d</mi><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>V</mi><mi mathvariant="normal">d</mi><mi>V</mi><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mfrac><mrow><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">]</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921247&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921246&type= 74.76066589 8.38199997 . 随着CFD流体仿真技术的发展，借助常用的CFD流体仿真软件CFX，计算得到望远镜圆顶周围附近空气的温度、速度和压力变化，再通过 式（7） 式（7） 计算望远镜在圆顶附近周围空气折射率的变化情况，以评估圆顶视宁度对望远镜成像的影响。 3 CFD仿真分析 1 3.1　有限元模型的建立 2 本文的塔台圆顶设计用于即将建设的国内大口径望远镜，其主镜口径为4 m，次镜口径为0.36 m，望远镜的总高约为11.5 m。塔台圆顶的设计参考了国外大口径望远镜的塔台圆顶设计。国外4 m级望远镜的塔台高度如 表1 表1 所示，综合考虑望远镜拟定站址的当地气象条件，拟定望远镜的塔台高度为15 m，圆顶的内径为 Φ 23 m，圆顶内中心高为15 m。根据参数分别建立敞开式随动圆顶和升降式圆顶的三维模型，如 图1 图1 和 图2 图2 所示。 望远镜 国家 位置 口径/m 塔台高度/m DCT 美国 Arizona 4.2 15 Soar 智利 Chile 4.1 15 ARIES 印度 India 3.6 13 EST 西班牙 Canary Islands 4.0 40 ATST 西班牙 La Palma 4.0 28 为了研究望远镜及圆顶周围空气在有风条件下的空气参数，需要在望远镜及圆顶周围建立空气流域场。为了减少边界条件不确定性对CFD仿真计算所造成的误差，根据望远镜及塔台的外形轮廓尺寸长 L 、宽 W 和高 H ，以塔台的底面为中心，分别向前后、左右和高度方向建立长度分别为5 L ，5 W 和5 H 的空气流域场，建立的模型如 图3 图3 和 图4 图4 所示。 3.2　网格划分 2 综合考虑仿真计算的收敛稳定性，计算的精度和准确性，网格划分采用非结构网格划分。为了在不降低运算精度的情况下尽可能减少网格的数量，提高计算效率，对结构复杂尺寸相对较小的望远镜及塔台采用较小的网格尺寸进行划分，而对于外围的较大的空气流域场则采用较大的网格尺寸进行划分。故先在望远镜周围建一个较小的空气流域，然后在较小的空气流域外建一个大的空气流域，以便对网格尺寸进行控制，如 图3 图3 和 图4 图4 所示。望远镜及塔台周围的平均网格尺寸为0.1 m，外围较大空气流域的平均网格尺寸为0.5 m，两种圆顶的有限元网格分别如 图5 图5 和 图6 图6 所示。其中，敞开式随动圆顶的有限元网格总数为2 754 484个，升降式圆顶的有限元网格总数为2 746 021个。 3.3　边界条件设置 2 望远镜的拟定站址为青海省海拔为3 000 m左右的某山，大气压强为67.2 kPa，空气密度为0.902 kg/m 3 ，平均风速最大为10 m/s，所以设置入口风速为10 m/s，空气流域的上边界和左右边界设置为freely slip walls，下边界设置为no-slip walls。由于望远镜自身的热惯性远远大于空气的热惯性，在工作过程中，当外部环境温度下降时，望远镜的温度变化很难跟上环境温度的变化。为了模拟这一现象，以外部环境温度变化速率±1 ℃/h（20%~90%的温度变化情况）为参考，假设望远镜、塔台、圆顶以及地面的初始温度高于环境温度1 ℃，环境的初始参考温度为20 ℃。 流体的实际运动状态大多为湍流。常用的湍流数值模拟方法主要3种，分别为直接数值模拟方法、大涡流模拟和雷诺时均模拟法，其中雷诺时均模拟法在实际工程中应用最为广泛。综合考虑各方面因素，如实际条件、计算消耗和计算精度等问题，根据常用的3种雷诺时均模拟法湍流模型的特点（见 表2 表2 ），最终确定选择SST模型，湍流强度设置为3.5%。 模 型 特 点 标准 <math id="M17"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921277&type= 1.43933344 2.37066650 - <math id="M18"><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921279&type= 2.11666679 2.37066650 模型 考虑低雷诺数、可压缩性和剪切流，在层流、平板和圆柱绕流较为理想 标准 <math id="M19"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921277&type= 1.43933344 2.37066650 - <math id="M20"><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921279&type= 2.11666679 2.37066650 模型 适用范围广，精度合理，属于半经验公式，但遇到强流线弯曲、漩涡、旋转流时，精度较低 SST模型 使 <math id="M21"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921277&type= 1.43933344 2.37066650 - <math id="M22"><mi>ω</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921279&type= 2.11666679 2.37066650 模型在广泛的流动领域中有了更高的精度和可行度 3.4　CFD仿真结果 2 经过瞬态分析得到敞开式随动圆顶和升降式圆顶望远镜中心剖面的流体仿真结果，如 图7 图7 和 图8 图8 所示。对比 图7 图7 和 图8 图8 可知，敞开式随动圆顶由于圆顶的作用，望远镜周围空气的速度、温度和相对压力的变化更大，且更加复杂，根据 式（7） 式（7） 可知，由此导致的空气折射率变化的一致性更差，圆顶视宁度的影响更大。而升降式圆顶则可以让周围的空气顺畅流通，望远镜周围空气的温度、速度和压力的变化相对较快，从而在一定程度上降低了空气折射率的不均匀变化，降低了圆顶视宁度的影响。 1;2;3; 1;2;3; 为了进一步分析两种圆顶结构方案对望远镜观测区域内空气折射率变化的影响，提取 图9 图9 所示4条光线上仿真结果的温度、速度和相对压力的数据。每条光线的长度为4.5 m，利用matlab绘制相应的变化曲线，如 图10 图10 和 图11 图11 所示。 由 图10 图10 和 图11 图11 可知，敞开式随动圆顶的变化范围更大，相应的速度、温度和相对压力变化为0~10 m/s，19.8~20.2 ℃和 - 28~51 Pa，而升降式圆顶的变化为2~12 m/s，19.96~20.07 ℃和17~25 Pa。根据所得速度、温度和相对压力数据，计算得到空气折射率相对于真空的折射率 <math id="M23"><mi mathvariant="normal">Δ</mi><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921317&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921316&type= 5.67266655 2.87866688 的变化曲线如 图12 图12 和 图13 图13 所示。由图可知，敞开式随动圆顶折射率变化 <math id="M24"><mi mathvariant="normal">Δ</mi><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921317&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=37921316&type= 5.67266655 2.87866688 的最大值与最小值之差为4.6×10 - 7 ，而升降式圆顶的折射率变化为1.2×10 - 7 。根据仿真结构可知，敞开式随动圆顶由于圆顶的作用，望远镜周围的空气会产生不稳定的涡流，导致空气的温度、速度和压力急剧变化。敞开式随动圆顶内光路中空气折射率的变化约为升降式圆顶内光路中空气中折射率变化的4倍。 4 结 论 1 本文以4 m级望远镜为研究对象，主要从圆顶视宁度的角度对望远镜的塔台圆顶结构进行了分析。视宁度产生的原因主要是由于空气折射率的不均匀变化，通过理论分析确定影响空气折射率的主要因素为空气的温度、速度和压力，然后建立敞开式随动圆顶和升降式圆顶的有限元分析模型，利用CFD流体仿真分析软件CFX进行分析，根据仿真结构可知，升降式圆顶内光路中空气中折射率变化约为敞开式随动圆顶内光路中空气折射率的变化1/4，可以减少圆顶视宁度的影响，更适合作为4 m级及更大口径望远镜的圆顶结构形式。 参考文献 张景旭 . 地基大口径望远镜系统结构技术综述 ［J］. 中国光学 ， 2012 ， 5 （ 4 ）： 327 - 336 . doi: 10.3788/CO.20120504.0327 http://dx.doi.org/10.3788/CO.20120504.0327 地基大口径望远镜系统结构技术综述 ZHANG J X . Overview of structure technologies of large aperture ground-based telescopes ［J］. Chinese Optics ， 2012 ， 5 （ 4 ）： 327 - 336 . （in Chinese） . doi: 10.3788/CO.20120504.0327 http://dx.doi.org/10.3788/CO.20120504.0327 Overview of structure technologies of large aperture ground-based telescopes 程景全 . 天文望远镜原理和设计射电、红外、光学、X射线和 γ 射线望远镜 ［M］. 北京 ： 中国科学技术出版社 ， 2003 . 天文望远镜原理和设计射电、红外、光学、X射线和γ射线望远镜 CHENG J Q . Principles of Astronomical Telescope Design ［M］. Beijing ： China Science and Technology Press ， 2003 . （in Chinese） Principles of Astronomical Telescope Design DENKER C ， LAGG A ， PUSCHMANN K G ， et al . The GREGOR Solar Telescope ［J］. IAU Special Session ， 2012 ， 6 ： E2 . 03. The GREGOR Solar Telescope Jr ROBERTS L C ， FIGGIS P D . Thermal conditioning of the AEOS telescope ［C］. Astronomical Telescopes and Instrumentation. Proc SPIE 4837 ， Large Ground-Based Telescopes ， Waikoloa ， Hawai'i， USA . 2003 ， 4837 ： 264 - 271 . Thermal conditioning of the AEOS telescope BLANCO D R ， JOHNS M W . Thermal design of the WIYN 3 .5-meter telescope enclosure［C］. 1994 Symposium on Astronomical Telescopes and Instrumentation for the 21st Century. Proc SPIE 2199 ， Advanced Technology Optical Telescopes V ， Kailua ， Kona， HI ， USA . 1994， 2199 ： 743 - 753 . Thermal design of the WIYN TERAN J ， PORTER D S ， HILEMAN E A ， et al . Unique dome design for the SOAR telescope project ［C］. Astronomical Telescopes and Instrumentation. Proc SPIE 4004 ， StructuresTelescope， Enclosures， Controls， Assembly/Integration/Validation， and Commissioning， Munich， Germany. 2000 ， 4004 ： 155 - 163 . Unique dome design for the SOAR telescope project 胡海飞 ， 赵宏伟 ， 刘振宇 ， 等 . 4 m口径SiC反射镜原位检测用静压支撑系统 ［J］. 光学 精密工程 ， 2017 ， 25 （ 10 ）： 2607 - 2613 . doi: 10.3788/OPE.20172510.2607 http://dx.doi.org/10.3788/OPE.20172510.2607 4 m口径SiC反射镜原位检测用静压支撑系统 HU H F ， ZHAO H W ， LIU ZH Y ， et al . Hydrostatic support system for in situ optical testing of a 4 m aperture SiC mirror ［J］. Opt. Precision Eng. ， 2017 ， 25 （ 10 ）： 2607 - 2613 . （in Chinese） . doi: 10.3788/OPE.20172510.2607 http://dx.doi.org/10.3788/OPE.20172510.2607 Hydrostatic support system for in situ optical testing of a 4 m aperture SiC mirror 吴小霞 ， 李剑锋 ， 宋淑梅 ， 等 . 4 m SiC轻量化主镜的主动支撑系统 ［J］. 光学 精密工程 ， 2014 ， 22 （ 9 ）： 2451 - 2457 . doi: 10.3788/OPE.20142209.2451 http://dx.doi.org/10.3788/OPE.20142209.2451 4 m SiC轻量化主镜的主动支撑系统 WU X X ， LI J F ， SONG S M ， et al . Active support system for 4 m SiC lightweight primary mirror ［J］. Optics and Precision Engineering ， 2014 ， 22 （ 9 ）： 2451 - 2457 . （in Chinese） . doi: 10.3788/OPE.20142209.2451 http://dx.doi.org/10.3788/OPE.20142209.2451 Active support system for 4 m SiC lightweight primary mirror 吴小霞 ， 王鸣浩 ， 明名 ， 等 . 大口径SiC轻量化主镜热变形的定标 ［J］. 光学 精密工程 ， 2012 ， 20 （ 6 ）： 1243 - 1249 . doi: 10.3788/OPE.20122006.1243 http://dx.doi.org/10.3788/OPE.20122006.1243 大口径SiC轻量化主镜热变形的定标 WU X X ， WANG M H ， MING M ， et al . Calibration of thermal distortion for large aperture SiC lightweight mirror ［J］. Optics and Precision Engineering ， 2012 ， 20 （ 6 ）： 1243 - 1249 . （in Chinese） . doi: 10.3788/OPE.20122006.1243 http://dx.doi.org/10.3788/OPE.20122006.1243 Calibration of thermal distortion for large aperture SiC lightweight mirror BOCCAS M . Protected-silver coatings for the 8-m Gemini telescope mirrors ［J］. Thin Solid Films ， 2006 ， 502 （ 1/2 ）： 275 - 280 . Protected-silver coatings for the 8-m Gemini telescope mirrors SEBAG J ， VOGIATZIS K . Estimating dome seeing for LSST ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 9150， Modeling， Systems Engineering， and Project Management for Astronomy VI ， Montréal， Quebec， Canada . 2014 ， 9150 ： 268 - 277 . Estimating dome seeing for LSST LAI O ， WITHINGTON J K ， LAUGIER R ， et al . Direct measure of dome seeing with a localized optical turbulence sensor ［J］. Monthly Notices of the Royal Astronomical Society ， 2019 ， 484 （ 4 ）： 5568 - 5577 . Direct measure of dome seeing with a localized optical turbulence sensor ZAGO L . Engineering handbook for local and dome seeing ［C］. Proc SPIE 2871， Optical Telescopes of Today and Tomorrow ， 1997 ， 2871 ： 726 - 736 . Engineering handbook for local and dome seeing WYNGAARD J C ， IZUMI Y ， COLLINS S A . Behavior of the refractive-index-structure parameter near the ground ［J］. JOSA ， 1971 ， 61 （ 12 ）： 1646 - 1650 . Behavior of the refractive-index-structure parameter near the ground BIDA T A ， DUNHAM E W ， BRIGHT L P ， et al . Site testing for the discovery channel telescope ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 5489， Ground-Based Telescopes， Glasgow， United Kingdom . 2004 ， 5489 ： 196 - 206 . Site testing for the discovery channel telescope KRABBENDAM V L ， RUTHVEN G P ， BENNETT V P ， et al . Active optical system design for the 4 .2-m SOAR telescope［C］. Astronomical Telescopes and Instrumentation. Proc SPIE 4003， Optical Design， Materials， Fabrication， and Maintenance ， Munich， Germany . 2000， 4003 ： 122 - 135 . Active optical system design for the SAGAR R ， NAJA M ， MAHESWAR G ， et al . Science at high-altitude sites of ARIES - astrophysics and atmospheric sciences ［J］. Proceedings of the Indian National Science Academy ， 2014 ， 80 （ 4 ）： 759 . Science at high-altitude sites of ARIES - astrophysics and atmospheric sciences COZAR-CASTELLANO J ， CAGIGAL M N ， VEGA N ， et al . Overview of telescope structure， enclosure， and pier preliminary design of the European Solar Telescope ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 12187， Modeling， Systems Engineering， and Project Management for Astronomy X ， Montréal， Québec， Canada . 2022 ， 12187 ： 693 - 703 . Overview of telescope structure， enclosure， and pier preliminary design of the European Solar Telescope DALRYMPLE N E ， Jr OSCHMANN J M ， HUBBARD R P . ATST enclosure： seeing performance， thermal modeling， and error budgets ［C］. SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 5497， Modeling and Systems Engineering for Astronomy， Glasgow， United Kingdom . 2004 ， 5497 ： 497 - 507 . ATST enclosure： seeing performance， thermal modeling， and error budgets 10.37188/OPE.20223023.3031 TH743 A 1004-924X（2022）23-3031-08 刘祥意,王志臣,王志.地基大口径望远镜圆顶设计[J].光学精密工程,2022,30(23):3031-3038. 刘祥意,王志臣,王志.地基大口径望远镜圆顶设计[J].光学精密工程,2022,30(23):3031-3038. LIU Xiangyi,WANG Zhichen,WANG Zhi.Dome design of large ground-based telescope[J].Optics and Precision Engineering,2022,30(23):3031-3038. LIU Xiangyi,WANG Zhichen,WANG Zhi.Dome design of large ground-based telescope[J].Optics and Precision Engineering,2022,30(23):3031-3038. 国家自然科学基金资助项目（No.U2031126） 2022-11-08 2022-11-18 2022-12-10 2022-12-20 光学精密工程 23 30