重力插值方法研究

为了根据离散观测数据构制连续空间重力变化图像，分析和讨论了３种数值插值方法，计算结果表明多面函数方法插值精度最高。由于逐步回归分析筛选核函数中心点的计算繁琐，文中提出根据分形理论和Ｓｈａｎｎｏｎ取样定理来确定核函数中心点。对滇西试验场进行模拟试算，插值精度可达到４～５（１０－８ｍｓ－２）。     

On water ice formation in interstellar clouds

A model is proposed for the formation of water ice mantles on grains in interstellar clouds. This occurs by direct accretion of monomers from the gas, be they formed by gas or surface reactions. The formation of the first monolayer requires a minimum extinction of interstellar radiation, sufficient to lower the grain temperature to the point where thermal evaporation of monomers is just offset by monomer accretion from the gas. This threshold is mainly determined by the adsorption energy of water molecules on the grain material; for hydrocarbon material, chemical simulation places this energy between 0.5 and 2 kcal mol⁻¹, which sets the (true) visible extinction threshold at a few magnitudes. However, realistic distributions of matter in a cloud will usually add to this an unrelated amount of cloud core extinction, which can explain the large dispersion of observed (apparent) thresholds. Once the threshold is crossed, all available water molecules in the gas are quickly adsorbed, because the grain cools down and the adsorption energy on ice is higher than on bare grain. The relative thickness of the mantle, and, hence, the slope of  τ₃( A _v)  depend only on the available water vapour, which is a small fraction of the oxygen abundance. Chemical simulation was also used to determine the adsorption sites and energies of O and OH on hydrocarbons and study the dynamics of formation of water molecules by surface reactions with gaseous H atoms, as well as their chances to stick in situ.     

The Yarkovsky effect as a heat engine

We show how the Yarkovsky effect can be understood as a heat engine. The output of the engine, manifested in the rate of change in semimajor axis of the body, has a maximum at an intermediate heat capacity, depending on the rotation rate of the body. This maximum arises because the work output depends on the product of the solar heat absorbed by the body and transported from its morning to evening side (this am-pm heat flux increases with heat capacity) and the Carnot efficiency (which declines with heat capacity).     

Microwave observations of Saturn's rings: anisotropy in directly transmitted and scattered saturnian thermal emission

We present a new Very Large Array (VLA) image of Saturn, made from data taken in October 1998 at a wavelength of λ3.6 cm. The moderate ring opening angle (B≈15°) allows us to explore direct transmission of microwave photons through the A and C rings. We find a strong asymmetry of photons transmitted through the A ring, but not in the C ring, a new diagnostic of wake structure in the ring particles. We also find a weak asymmetry between east and west for the far side of the ansae. To facilitate quantitative comparison between dynamic models of the A ring and radio observations, we extend our Monte Carlo radiative transfer code (described in Dunn et al., 2002, Icarus 160, 132-160) to include idealized wakes. We show the idealized model can reproduce the properties of dynamic simulations in directly transmitted light. We examine the model behavior in directly transmitted and scattered light over a range of physical and geometric wake parameters. Finally, we present a wake model with a plausible set of physical parameters that quantitatively reproduces the observed intensity and asymmetry of the A ring both across the planet and in the ansae.     

Evolution of Polytropic Reimann Disks in Halos

The equilibrium of elliptical Riemann disks with a polytropic equation of state and their evolution under the influence of viscosity and gravitational radiation inside spheroidal halos with a relative surface mass density k. The evolutionary trajectory of a disk inside a halo with k<0.5, which is analogous to the evolution of an isolated disk, differs from that of a disk inside a denser halo.     

On the relation between orbital structure and observed bar morphology

We investigate the morphological relation between the orbits of the central family of periodic orbits ( x ₁ family) and the bar itself using models of test particles moving in a barred potential. We show that different bar morphologies may have as a backbone the same set of x ₁ periodic orbits. We point out that by populating initially axisymmetric stellar discs exponentially with test particles in circular, or almost circular motion, we may end up with a response bar which reveals a shape different in crucial details from that of the individual stable x ₁ orbits. For example, a bar model in which the x ₁ orbits are pure ellipses may have a much more complicated response morphology. This depends on the particular invariant curves around x ₁, which are populated in each model.     

Study on numerical simulation to water lowering in an open coal mine digging

In view of the situation of excavation of open coal mine for the underground water disaster,we should carry out simulation studies for the numerical value of the water lowering project and improve the accuracy and the level of the water lowering project.On the basis of the hydrological geological conditions of certain open mine digging,a more reasonable seepage numerical model was built according to MODFLOW.It was simulated in advance that the process of the confined water level descending with the time,and combining with the actual observations to test the correctness of the model.The calculation showed that the results coincided well with the results of actual measurement.Based on this,different water lowering numerical simulations were built for the open coal mine digging.It could be simulated and forecast that the changes of the groundwater level in drainage process within and outside the mine pit,and it was quantitatively assessed that the possible water lowering result of the opencast water drainage process,which provide an important basis for the actual water lowering project and the possible project disposal.