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

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

Impulsive time-free transfers between halo orbits

A methodology is developed to design optimal time-free impulsive transfers between three-dimensional halo orbits in the vicinity of the interior L ₁ libration point of the Sun-Earth/Moon barycenter system. The transfer trajectories are optimal in the sense that the total characteristic velocity required to implement the transfer exhibits a local minimum. Criteria are established whereby the implementation of a coast in the initial orbit, a coast in the final orbit, or dual coasts accomplishes a reduction in fuel expenditure. The optimality of a reference two-impulse transfer can be determined by examining the slope at the endpoints of a plot of the magnitude of the primer vector on the reference trajectory. If the initial and final slopes of the primer magnitude are zero, the transfer trajectory is optimal; otherwise, the execution of coasts is warranted. The optimal time of flight on the time-free transfer, and consequently, the departure and arrival locations on the halo orbits are determined by the unconstrained minimization of a function of two variables using a multivariable search technique. Results indicate that the cost can be substantially diminished by the allowance for coasts in the initial and final libration-point orbits.An earlier version was presented as Paper AIAA 92-4580 at the AIAA/AAS Astrodynamics Conference, Hilton Head Island, SC, U.S.A., August 10–12, 1992.     

Adiabatic invariants for the nonconservative Kepler's problem

This paper considers adiabatic invariants for the classical Kepler problem with resisting forces. The analysis is based on the theory of integrating factors and theory of adiabatic invariants in the Krylov-Bogoliubov-Mitropolski variables. The adiabatic invariants are series with respect to a small parameter. Also, for every particular case of nonconservative forces, it is shown that, with a complete set of adiabatic invariants, an approximate solution of the problem can be obtained. Four problems are analyzed in detail where approximate solutions are compared with numerical.     

On the global solution of the N-body problem

In connection with the publication (Wang Qiu-Dong, 1991) the Poincaré type methods of obtaining the maximal solution of differential equations are discussed. In particular, it is shown that the Wang Qiu-Dong'sglobal solution of the N-body problem has been obtained in Babadzanjanz (1979). First the more general results on differential equations have been published in Babadzanjanz (1978).     

Conditioning mean steady state flow on hydraulic head and conductivity through geostatistical inversion

Nonlocal moment equations allow one to render deterministically optimum predictions of flow in randomly heterogeneous media and to assess predictive uncertainty conditional on measured values of medium properties. We present a geostatistical inverse algorithm for steady-state flow that makes it possible to further condition such predictions and assessments on measured values of hydraulic head (and/or flux). Our algorithm is based on recursive finite-element approximations of exact first and second conditional moment equations. Hydraulic conductivity is parameterized via universal kriging based on unknown values at pilot points and (optionally) measured values at other discrete locations. Optimum unbiased inverse estimates of natural log hydraulic conductivity, head and flux are obtained by minimizing a residual criterion using the Levenberg-Marquardt algorithm. We illustrate the method for superimposed mean uniform and convergent flows in a bounded two-dimensional domain. Our examples illustrate how conductivity and head data act separately or jointly to reduce parameter estimation errors and model predictive uncertainty.This work is supported in part by NSF/ITR Grant EAR-0110289. The first author was additionally supported by scholarships from CONACYT and Instituto de Investigaciones Electricas of Mexico. Additional support was provided by the European Commission under Contract EVK1-CT-1999-00041 (W-SAHaRA-Stochastic Analysis of Well Head Protection and Risk Assessment).     

A BME solution of the inverse problem for saturated groundwater flow

In most real-world hydrogeologic situations, natural heterogeneity and measurement errors introduce major sources of uncertainty in the solution of the inverse problem. The Bayesian Maximum Entropy (BME) method of modern geostatistics offers an efficient solution to the inverse problem by first assimilating various physical knowledge bases (hydrologic laws, water table elevation data, uncertain hydraulic resistivity measurements, etc.) and then producing robust estimates of the subsurface variables across space. We present specific methods for implementing the BME conceptual framework to solve an inverse problem involving Darcys law for subsurface flow. We illustrate one of these methods in the case of a synthetic one-dimensional case study concerned with the estimation of hydraulic resistivity conditioned on soft data and hydraulic head measurements. The BME framework processes the physical knowledge contained in Darcys law and generates accurate estimates of hydraulic resistivity across space. The optimal distribution of hard and soft data needed to minimize the associated estimation error at a specified sampling cost is determined. This work was supported by grants from the National Institute of Environmental Health Sciences (Grant no. 5 P42 ES05948 and P30ES10126), the National Aeronautics and Space Administration (Grant no. 60-00RFQ041), the Army Research Office (Grant no. DAAG55-98-1-0289), and the National Science Foundation under Agreement No. DMS-0112069.     

Special Solutions in a Simplified Restricted Three-Body Problem with Gravitational Radiation Taken into Account

The distinctive feature of the relativistic restricted three-body problem within the c ^–5 order of accuracy (2 post-Newtonian approximation) is the presence of the gravitational radiation. To simplify the problem the motion of the massive binary components is assumed to be quasi-circular. In terms of time these orbits have linearly changing radii and quadratically changing phase angles. By substituting this motion into the Newtonian-like equations of motion one gets the quasi-Newtonian restricted quasi-circular three-body problem sufficient to take into account the main indirect perturbations caused by the binary radiation terms. Such problem admits the Lagrange-like quasi-libration solutions and rather simple quasi-circular orbits lying at large distance from the binary.     

GIS中线元的误差熵带研究

基于现有的线元位置不确定性模型大多与置信水平的选取有关,而置信水平的选取带有一定程度的主观性,因而不能惟一确定,引入信息熵理论,提出了线元的误差熵带模型,并将它与“E－带”进行了比较,计算了落入其内的概率。该模型根据联合熵惟一确定,与置信水平的选取无关。     

The Non-Linear Stability of L4 in the Restricted Three-Body Problem when the Bigger Primary is a Triaxial Rigid Body

The non-linear stability of L ₄ in the restricted three-body problem has been studied when the bigger primary is a triaxial rigid body with its equatorial plane coincident with the plane of motion. It is found that L ₄ is stable in the range of linear stability except for three mass ratios: