用切贝谢夫配点法消除地球自振常微分方程组在地心处的奇异性

用切贝谢夫配点法求解地球自振常微分方程组，无需进一步改化即可消除这组方程在地心处的奇异性，并能获得高精度的结果     

地面总辐射的异常减少对气候影响的数值试验

在与实际台风场近似一致的环境流场条件下，应用柱会标中的正压无辐散模式，讨论了台风中扰动的稳定性问题。研究表明，台风暖心结构的台风外围适当强度的冷空气活动有利限台风中扰动发展，且台风愈强，台风中扰动愈易发展。     

不同方案求解非线性化学动力学方程组的比较

在一个简化的气相化学模式的基础上，比较了Ｈｙｂｒｉｄ、ＱＳＳＡ、Ｓｋｌａｒｅｗ三种计算方案求解非线性化学动力学方程组的差异，并引入误差分配和线性组合两种质量守恒技术，分析它们对模拟结果的影响．研究结果表明：不同方案对有机烃类浓度的计算几乎没有影响，但对其它气态物浓度有一定的影响．用ＱＳＳＡ和Ｓｌａｒｅｗ方案预测的结果更为相近，而采用Ｈｙｂｒｉｄ方案得到的气态物浓度的平衡时间略为提前．从计算效率来看，ＱＳＳＡ方案优于其它两种方案，两种质量守恒技术均能改善模拟结果．     

A SIMPLIFIED APPROACH TO MODELING 3D SEDIMENT-LADEN TURBULENT FLOWS

ＡＳＩＭＰＬＩＦＩＥＤＡＰＰＲＯＡＣＨＴＯＭＯＤＥＬＩＮＧ３ＤＳＥＤＩＭＥＮＴ－ＬＡＤＥＮＴＵＲＢＵＬＥＮＴＦＬＯＷＳＤｏｎｇｈｕｏＺＨＯＵａｎｄＳａｍＳ．Ｙ．ＷＡＮＧＡｂｓｔｒａｃｔ：Ａ３－ｄｎｕｍｅｒｉｃａｌｍｏｄｅｌｆｏｒｓｉｍｕｌａｔｉｎｇｓ...     

Stochastic perturbation analysis of groundwater flow. Spatially variable soils,semi-infinite domains and large fluctuations

As is well known, a complete stochastic solution of the stochastic differential equation governing saturated groundwater flow leads to an infinite hierarchy of equations in terms of higher-order moments. Perturbation techniques are commonly used to close this hierarchy, using power-series expansions. These methods are applied by truncating the series after a finite number of terms, and products of random gradients of conductivity and head potential are neglected. Uncertainty regarding the number or terms required to yield a sufficiently accurate result is a significant drawback with the application of power series-based perturbation methods for such problems. Low-order series truncation may be incapable of representing fundamental characteristics of flow and can lead to physically unreasonable and inaccurate solutions of the stochastic flow equation. To support this argument, one-dimensional, steady-state, saturated groundwater flow is examined, for the case of a spatially distributed hydraulic conductivity field. An ordinary power-series perturbation method is used to approximate the mean head, using second-order statistics to characterize the conductivity field. Then an interactive perturbation approach is introduced, which yields improved results compared to low-order, power-series perturbation methods for situations where strong interactions exist between terms in such approximations. The interactive perturbation concept is further developed using Feynman-type diagrams and graph theory, which reduce the original stochastic flow problem to a closed set of equations for the mean and the covariance functions. Both theoretical and practical advantages of diagrammatic solutions are discussed; these include the study of bounded domains and large fluctuations.     

反演二维弹性波动方程组多参数的渐近展开法

采用摄动法对二维弹性波动方程组的系数反问题进行渐近近似,给出了反演系数的具体计算方法。用摄动法解微分方程组的系数反问题还是首次。此外,还采用了地球物理分层模型,这种分层模型便于更精细地反映地球不同圈层的物理特征,更接近于实际情况。同时对进一步提高反演分辨率、增加反演参数、对病态方程组如何克服不适定性等问题,作了具体阐述,具有重要的理论意义与应用价值。     

Groebner-basis solution of the three-dimensional resection problem (P4P)

The three-dimensional (3-D) resection problem is usually solved by first obtaining the distances connecting the unknown point P{X,Y,Z} to the known points P_i{X_i,Y_i,Z_i}i=1,2,3 through the solution of the three nonlinear Grunert equations and then using the obtained distances to determine the position {X,Y,Z} and the 3-D orientation parameters {,, }. Starting from the work of the German J. A. Grunert (1841), the Grunert equations have been solved in several substitutional steps and the desire as evidenced by several publications has been to reduce these number of steps. Similarly, the 3-D ranging step for position determination which follows the distance determination step involves the solution of three nonlinear ranging (`Bogenschnitt') equations solved in several substitution steps. It is illustrated how the algebraic technique of Groebner basis solves explicitly the nonlinear Grunert distance equations and the nonlinear 3-D ranging (`Bogenschnitt') equations in a single step once the equations have been converted into algebraic (polynomial) form. In particular, the algebraic tool of the Groebner basis provides symbolic solutions to the problem of 3-D resection. The various forward and backward substitution steps inherent in the classical closed-form solutions of the problem are avoided. Similar to the Gauss elimination technique in linear systems of equations, the Groebner basis eliminates several variables in a multivariate system of nonlinear equations in such a manner that the end product normally consists of a univariate polynomial whose roots can be determined by existing programs e.g. by using the roots command in Matlab.Acknowledgments.The first author wishes to acknowledge the support of JSPS (Japan Society of Promotion of Science) for the financial support that enabled the completion of the write-up of the paper at Kyoto University, Japan. The author is further grateful for the warm welcome and the good working atmosphere provided by his hosts Professors S. Takemoto and Y. Fukuda of the Department of Geophysics, Graduate School of Science, Kyoto University, Japan.     

一类双曲-抛物型方程的广义解

讨论了一类奇摄动双曲-抛物型方程广义初边值问题,在适当的条件下,用Galerkin方法研究了广义解的存在性、唯一性,同时得到了解的渐近估计式。     

Multiple Parameter Regularization: Numerical Solutions and Applications to the Determination of Geopotential from Precise Satellite Orbits

Kaula’s rule of thumb has been used in producing geopotential models from space geodetic measurements, including the most recent models from satellite gravity missions CHAMP. Although Xu and Rummel (Manuscr Geod 20 8–20, 1994b) suggested an alternative regularization method by introducing a number of regularization parameters, no numerical tests have ever been conducted. We have compared four methods of regularization for the determination of geopotential from precise orbits of COSMIC satellites through simulations, which include Kaula’s rule of thumb, one parameter regularization and its iterative version, and multiple parameter regularization. The simulation results show that the four methods can indeed produce good gravitational models from the precise orbits of centimetre level. The three regularization methods perform much better than Kaula’s rule of thumb by a factor of 6.4 on average beyond spherical harmonic degree 5 and by a factor of 10.2 for the spherical harmonic degrees from 8 to 14 in terms of degree variations of root mean squared errors. The maximum componentwise improvement in the root mean squared error can be up to a factor of 60. The simplest version of regularization by multiplying a positive scalar with a unit matrix is sufficient to better determine the geopotential model. Although multiple parameter regularization is theoretically attractive and can indeed eliminate unnecessary regularization for some of the harmonic coefficients, we found that it only improved its one parameter version marginally in this COSMIC example in terms of the mean squared error.     

顾及地球引力摄动的卫星和卫星星下点位置计算

本文叙述了顾及地球引力摄动的卫星化置计算方法；利用天球坐标系与地球坐标系之间的坐标变换关系，给出了由卫星位置计算卫星星下点位置的严密计算公式。