恒星线性非绝热非径向脉动问题的退耦化处理方法

本文发展了一种解恒星线性非绝热非径向脉动问题的退耦化方法。这个方法把非绝热非径向脉动问题的六阶线性微分方程,分解为由一个代数方程联系起来的一个四阶线性微分方程和另一个二阶线性微分方程进行数值求解。这样的一个退耦处理,有利于克服以前在数值解这类问题时常常遇到的收敛域小和收敛速度慢等困难,并且为数值解方程时所采用的Henyey方法提供了一个自然和方便的初始猜测解。     

星际气体激波与恒星线性密度波的总合自洽解

本文从二维气体动力学和泊松方程出发,研究星际气体激波和线性恒星密度波自洽的密度波理论,求出准稳、紧卷、双臂的线性恒星密度波和星际气体激波的总合自洽解。所得结果表明,在采用不完全的流体力学线性方程组及全局的引力方程的条件下,得到的是收敛的自洽解,该解给出了与局部渐近解相同的色散关系、非谐波形式的总密度分布及总引力分布,所得结果还表明,若采用完全的流体力学线性化方程组,将得出不稳定的线性恒星密度波。     

粒子在标量张量理论的二阶后闵可夫斯基近似弱引力场中的引力偏转

标量张量理论是目前为止最成功的相对论性引力理论之一,经受住了已有实验和观测的检验.随着实验和观测精度的提高,理论将得到进一步的检验.为此将该理论的可观测效应计算到二阶后牛顿近似非常必要.标量张量理论的二阶后牛顿近似结果已经由Xie等人给出.为得到有质量粒子在弱引力场中的偏转情况,利用二阶后闵可夫斯基近似下的度规解得到了粒子的测地线方程,并通过迭代的方法求得在非束缚条件下,粒子在弱引力场中的轨迹的解析解.进一步利用此结果计算了粒子在弱引力场中的偏转角(与光线的偏转角不同)并与前人的结果进行了比较.     

白天赤道电激流区梯度漂移不稳定性的非局域分析(英文)

八十年代,人们在赤道电激流区观测到数公里尺度的强非线性波动现象。对这一现象目前还没有满意的理论解释。 根据一些作者对速度剪切对不稳定性的影响的讨论,我们认为速度剪切对长波的影响不是决定性的,描述梯度漂移不稳定性的基本方程可以作相当大的简化。对一个简单的等离子体密度分布模型,该微分方程可以解析地求解。 本文结果表明,大尺度波的相速度和增长率不仅依赖于波数,也依赖于模数—非局域解在不稳定区域中的节点数。与局域理论相比,非局域效应大大地减小了增长率。结果还表明,如果考虑复合耗散,增长率的长波截断与观测结果是吻合的,另外,波的相速度与观测结果也相合。 线性分析的结果还不能完全解释大尺度波现象,但是本文的解析解对今后的非线性分析很有帮助。 类似的分析方法将应用到夜间赤道电激流区和中纬Es层。     

太阳光谱中的磁场敏感谱线

本文通过斯托克斯参数转移方程的数值解,在比较普遍的情况下,计算了FeIλλ6302.499、6173.341、5250.216三条磁场敏感谱线的斯托克斯参数轮廓,并研究了它们的磁场增强、温度敏感等一系列性质。海野和三郎求解斯托克斯参数转移方程时,其中有两个假设为:(1)能源函数B_2随光学深度r_2呈线性分布;(2)大气结构为密尔恩-爱丁顿模型。这两个假设与实际情况相差太远,本文放弃了它们。与海野的结果比较起来,我们的数值解与观测结果符合较好。     

数值历表对动力系统参数的偏导数

本文推导了天体运动方程的数值解对积分初始条件、天体质量等动力系统参数的偏导数所满足的微分方程和初始条件。     

Chandler摆动的非线性分岔耗散

在已有的Ｃｈａｎｄｌｅｒ摆动非线性共振激发模型的基础上,引入微小的非线性项,得到不可忽略的非线性分岔解。实际上,Ｃｈａｎｄｌｅｒ摆动方程是可用线性阻尼非线性平面单摆来逼近的,由Ｄｕｆｆｉｎｇ方程可以求得方程的解。     

地震位移场对极移的激发

本文讨论了地震位移场对极移的激发问题.与已有的文献中采用的液核静力方程不同,我们采用文[1]中提出改进的液核静力方程.这样核幔间的连续性条件都可以得到满足.本文所用的地震参数取自文[2]. 在微分方程求解过程中,把整个球体分成两部分:内部是球心附近的一个小球,外部为一球壳.在内部我们应用幂级数渐近解,外部采用Runge-Kutta数值解.这样可以避免数值解过程中球心处的奇异性问题.计算结果表明,虽然用同样的地震参数,由于液核方程的不同,我们的结果比文[2]的结果大三倍左右.这说明,不同的液核方程的影响是比较大的.     

窄吸积环动力学不稳定性的数值解

本文在细环近似和共转半么位移不为零的条件下，考虑吸积环中非轴对称动力学不稳定性的线性扰动过程，采用数值计算方法求得不稳定性的线性增长率和共转半径位移随波数的变化关系，发现线性增长率受共转半径不为零的影响较小，而共转半径位移项随波数的色散关系与线性Ｋｄｖ方程的色散关系相同，说明窄吸收环在动力学不稳定性数值模拟中出现的“行星状解”很可能是类似于Ｋｄｖ方程中的孤子解。     

Vinti型中间轨道根数的摄动方程

本文给出了Vinti型中间轨道根数a,e,s,M_s,ψ_s,Ω_s的摄动方程,结果准到二阶,并消除了方程中的泊松项。     

VLBI数据处理中的剩余钟行为与剩余大气效应建模

抽样选取了 2 7次上海天文台佘山站参与的天测与测地VLBI实验 ,分别进行了单次解算。通过分析解算参数随剩余钟行为和剩余大气效应分段拟合长度的变化 ,得到以下初步结论 :(1)选取不同分段拟合长度时 ,站坐标解算结果和时延残差加权均方根存在差异 ,最大分别至厘米级和数十皮秒 ,因而分段拟合长度不能随意选取。 (2 )分段拟合时段长度存在某一合理取值范围 ,它不宜过长 ,否则钟和大气的剩余效应短周期变化不能很好地模型化。为了保证待估参数解算时有足够的自由度 ,拟合时段不宜过短 ,否则将导致法方程近于或出现奇异 ,达不到较好控制噪声的效果。 (3)由于各次实验、同一次实验中的不同台站相应的钟和大气条件存在差异 ,有必要对每次实验以及每次实验中的各观测台站分别分析 ,寻找合适的分段拟合长度。这在实际操作中显然相当烦琐。 (4)一般而言 ,在剩余钟行为拟合长度缺省值 6 0min情况下 ,剩余大气效应分段拟合长度以介于 10min至 4 0min为宜 ;在剩余大气效应分段拟合长度缺省值 2 0min情况下 ,剩余钟行为分段拟合长度以介于 2 0min至 10 0min为宜     

Ekman layer damping of r modes revisited

We investigate the damping of neutron star r modes due to the presence of a viscous boundary (Ekman) layer at the interface between the crust and the core. Our study is motivated by the possibility that the gravitational wave driven instability of the inertial r modes may become active in rapidly spinning neutron stars, for example, in low-mass X-ray binaries, and the fact that a viscous Ekman layer at the core–crust interface provides an efficient damping mechanism for these oscillations. We review various approaches to the problem and carry out an analytic calculation of the effects due to the Ekman layer for a rigid crust. Our analytic estimates support previous numerical results, and provide further insight into the intricacies of the problem. We add to previous work by discussing the effect that compressibility and composition stratification have on the boundary-layer damping. We show that, while stratification is unimportant for the r-mode problem, composition suppresses the damping rate by about a factor of 2 (depending on the detailed equation of state).     

A new numerical method for solving radiation driven winds from hot stars

We present a general method for solving the non‐linear differential equation of monotonically increasing steady‐state radiation driven winds. We graphically identify all the singular points before transforming the momentum equation to a system of differential equations with all the gradients explicitly given. This permits a topological classification of all singular points and to calculate the maximum and minimum mass‐loss of the wind. We use our method to analyse for the first time the topology of the non‐rotating frozen‐in ionisation m‐CAK wind, with the inclusion of the finite disk correction factor, and find up to 4 singular points, three of the x‐type and one attractor‐type. The only singular point (and solution passing through) that satisfies the boundary condition at the stellar surface is the standard m‐CAK singular point. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

地球大气边界层结构和天文视宁度

依据声雷达资料和其他边界层资料讨论了大气边界层结构、地形等对视宁度的影响和有关选址问题．     

CCD联合平差方法及其初步实例检验

在常规照相天体测量工作中 ,对暗天体的定位是采用逐级定标星过渡的方法实现的。由于星表系统差和局部差的存在 ,最终定位结果可能包含难以把握的误差积累。在CCD小视场观测情况下时常难以找到足够数量的定标星。鉴于此 ,推导了CCD观测联合平差方法的严格矢量表达式 ,以期通过共同星的过渡 ,实现大天区统一平差。原理上此方法可实现相对于河外天体的定位 ,以削弱星表系统差和局部差对定位结果的影响。讨论了CCD联合平差方法实现中的问题 ,包括底片常数模型选择、切点位置改正、共同星的较差改正和法方程解算方法等 ,并分析了具体的处理措施。用实例初步检验了CCD联合平差方法的效果 ,表明通过对多张CCD观测的联合平差 ,可以削弱单张CCD观测参考星数目过少、单张CCD观测参考星分布不均、个别参考星存在位置偏差和局部区域参考星存在系统性位置偏差等不利因素对归算结果的影响 ,进而达到扩大视场和提高归算精度的目的     

The boundary-layer flow of a radiating gas by the integral equation method

The flow in the boundary layer of a very hot two-component plasma is analysed when the radiative heat flux is given by the exact integral equation expression. The basic nonlinear integro-differential equation is solved by perturbing it about the differential approximation for radiation. In this way some light is shed on the order of accuracy of the differential approximation of radiation. In fact an error of about 4% may be incurred by invoking the differential approximation.     

Determination of the potential in a synodic system

Determination of the potential field in a fixed (inertial) system may be accomplished by the solution of a homogeneous linear partial differential equation when a family of orbits of a body moving in the field is given. This partial differential equation was presented and thoroughly analyzed earlier. The present paper discusses the same problem in a rotating system where the centrifugal and Coriolis effects render the pertinent partial differential equation in general non-homogeneous and non-linear. A linear, though non-homogeneous, partial differential equation for the determination of the synodic potential is obtained only in the special case of iso-energetic families of orbits.     

疏散星团NGC6530天区的恒星自行和成员概率

利用上海天文台的照相底片资料,确定了疏散星团NGC6530天区364颗恒星的自行和成员概率,并对有关自行测定的方法、结果和精度等问题作了较为详细的介绍和讨论。使用的底片历元差为87年,全部恒星自行中误差的均方根值为1．09mas/a。     

ADS辅助实现线极化到圆极化的转换

给出C波段双圆极化低噪声接收机中圆极化输出的实现方案,研究了利用90°电桥进行线- 圆极化转换时两个线极化通道增益及相位的平衡度对圆极化隔离的影响。分析了如何通过调节一对同轴电缆的长度差来补偿相位不平衡,并借助ADS软件计算了该长度差,实现了线-圆极化转换,取得接收机两个圆极化之间的隔离度好于20dB的效果。