一种从经纬度网格资料获得正方形网格资料的插值方案

本文提出了一种把经纬度网格点上的气象要素值转化到正方形网格点上的多项式插值方案。对插值方案作了稳定性讨论和定量及定性分析。结果表明,这种插值方案方便可行,而且有广泛适用性。     

Simulation of adaptive optics for the Vacuum Tower Telescope

A simulation model of the adaptive optics of the German Vacuum Tower Telescope (VTT), Observatorio del Teide, Tenerife, is presented. The model uses modules from the integrated model of the Euro50 extremely large telescope, and includes submodels of a Shack-Hartmann wavefront sensor, a de-formable mirror, a tip-tilt mirror, high-voltage amplifier low-pass filters, a reconstructor and a controller. We investigate the impact on the closed loop bandwidth of changes in controller configuration and certain system parameters, such as low pass filter bandwidth and camera integration and readout time. Control strategies were tested on simple models before implementation on the full VTT model. Using the models, different control strategies are compared.     

Numerical calculations in the orbital determination of an artificial satellite for a long arc

Numerical methods are usually used for the computation of ephemerides with perturbations for the precise orbital determination of an artificial satellite. But their numerical stability will be encountered in a long arc. In this case the use the improved Encke special perturbation methods has been suggested. The results of this paper show that Encke's method does indeed have a certain effectiveness, but cannot yet completely resolve the numerical stability, and the more efficient method is to use the energy integral or its variational relation to control the growth of the along-track error in general numerical calculations so that the aim of stabilization can be achieved.     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.     

Obtaining basic simulation data for a heterogeneous field with stratified marine soils

Thinly stratified sedimentary deposits in a heterogeneous field were investigated to obtain basic physical data for the simulation of water flow. A procedure is described which translates a thinly stratified soil profile into a number of functional layers using functional hydrological properties. A functional layer is defined as a combination of one or more soil horizons and should (i) be recognizable during a soil survey using an auger and (ii) show significantly different functional hydrological properties when compared with another functional layer. This procedure gave three easily recognizable functional layers. Sets of hydrological characteristics of these three functional layers were obtained by physical measurements of the soil and by estimation, using textural data for classification into a standard Dutch series. The performance of several combinations of these sets was tested by comparing simulated and measured soil matric potentials for seven plots during one year. The best simulation results were obtained if measured soil hydraulic characteristics were used for relatively homogeneous functional layers and if the soil hydraulic characteristics were estimated at each location for the most heterogeneous layer.     

Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

Three-dimensional discontinuous deformation analysis (3-D DDA) using a new contact resolution algorithm

This paper presents a new contact calculating algorithm for contacts between two polyhedra with planar boundaries in the three-dimensional discontinuous deformation analysis (3-D DDA). In this algorithm, all six type contacts in 3-D (vertex-to-face, vertex-to-edge, vertex-to-vertex, face-to-face, edge-to-edge, and edge-to-face) are simply transformed into the form of point-to-face contacts. The presented algorithm is a simple and efficient method and it can be easily coded into a computer program. In this paper, formulations of normal contact, shear contact and frictional force submatrices based on the new method are derived and the algorithm has been programmed in VC⁺⁺. Examples are provided to demonstrate the new contact rule between two blocks.     

Nonlinear dynamic finite element simulation of Alibey earth dam

In this study, dynamic behavior and earthquake resistance of Alibey earth dam was investigated. The dam was modeled with four node plane-strain finite elements (FE) and displacement-pore pressure coupled FE analyses were performed. Nonlinear material models such as pressure dependent and independent multi yield materials were implemented during the analyses. Transient dynamic FE analyses were performed with Newmark method. The Newton-Raphson solution scheme was adopted to solve the equations. Liquefaction and/or cyclic mobility effects were considered during the analysis. For the FE analyses, OpenSees (Open System for Earthquake Engineering Simulation) framework was adopted.     

On CO<Subscript>2</Subscript> fluid flow and heat transfer behavior in the subsurface,following leakage from a geologic storage reservoir

Geologic storage of CO₂ is expected to produce plumes of large areal extent, and some leakage may occur along fractures, fault zones, or improperly plugged pre-existing wellbores. A review of physical and chemical processes accompanying leakage suggests a potential for self-enhancement. The numerical simulations presented here confirm this expectation, but reveal self-limiting features as well. It seems unlikely that CO₂ leakage could trigger a high-energy run-away discharge, a so-called “pneumatic eruption,” but present understanding is insufficient to rule out this possibility. The most promising avenue for increasing understanding of CO₂ leakage behavior is the study of natural analogues.