Numerical Stabilization of Orbital Motion

Mainly, the author focuses on Baumgarte's method and its applications in satellite, asteroid, stellar and planetary problems. In the paper arguments are given for the use of energy relations for stabilization in the elliptical two-body problem. Stabilizing properties of Baumgarte's equations and others are discussed. A simple approach is proposed for stabilizing the equations of almost circular motion. By using Baumgarte's technique, the author derives stabilized equations of perturbed restricted three-body problem. It is shown experimentally that stabilization in the problems mentioned above can raise the accuracy of numerical integration by several orders.     

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

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

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Ground motion estimation during the Kashmir earthquake of 8th October 2005

In this article, analytical methods have been used to estimate ground motion during the 8 October 2005, Kashmir earthquake. Peak ground acceleration (PGA) values at several stations in the epicentral region have been estimated by empirical analytical source mechanism models. As an alternate analysis, PGA estimates have also been obtained using the stochastic finite fault seismological model. The estimated PGAs are compared with that obtained from damage values. A PGA contour map in the near-source region is provided. It is found that very near to the epicenter, PGA would have reached more than 1 g. It is demonstrated that empirical analytical models can be effectively used to estimate ground motion due to rupture of active faults.     

Numerical study of the motions and drift force of a floating OWC device

The motion and the drift force of a floating OWC (oscillating water column) wave energy device in regular waves are studied taking account of the oscillating surface-pressure due to the pressure drop across the duct of the air chamber. The potential problem inside the chamber is formulated by making use of the Green integral equation associated with the Rankine-type Green function while the outer problem with the Kelvin-type Green function. The added mass, wave damping and excitation coefficients as well as the motion and drift force of the OWC device are calculated for various values of parameter related to the pressure drop.     

船舶受冲击力作用时的水动力试验研究

根据冲量相似的要求进行船舶在冲击作用下的模型试验，并将模型试验结果与理论计算及实际试验进行了对比和分析。     

Effect of wave-directionality on second-order loads induced by the set-down

At first-order of approximation a sea-state may be considered as an infinite sum of Airy components with angular frequencies ω_i and wave-number vectors k_i. A second-order analysis shows the co-existence of long waves appearing at the difference frequencies ω_i—ω_j with wave-number vectors k_i—k_j. In shallow water they become appreciable in amplitudes and may induce slow-drift motion of moored structures.For small values of ω_i—ω_j,k_i—k_j may take all kinds of directions for an angular-spread wave system. Then it may be questioned how the in-line and transverse second-order accelerations compare to those obtained for a mono-directional wave-system.This analysis is carried out here by relating the spectra of the second-order horizontal accelerations to the directional wave-spectrum. Numerical applications are first performed for deep water. They show that at low frequencies, even for very narrowly spread wave systems, the transverse component is larger than the in-line component. In shallow water both components are dratically reduced as compared to the mono-directional case. As a consequence one may question the validity of model-testings or numerical models which take no account of the directionality of the wave-system.     

Rheological Properties and Incipient Motion of Cohesive Sediment in the Haihe Estuary of China

The Haihe cohesive sediment, which is typical in China, is studied systematically for its basic physical and incipientmotion properties. Following the requirements of dredging works in the Haihe Estuary, cohesive sediment samples weretaken from three locations. Laboratory experiments were conducted to determine the rheological properties of these sam-ples and to examine the incipient motion of the cohesive sediment. It is found that the cohesive sediment has an obviousyield stress τb, which increases with the mud density in a manner of an exponential function, and so does the viscosityparameter η. The cohesive sediment behaves like a Bingham fluid when its density is below 1.38 ～ 1.40 g/cm3, andwhen denser than these values, it may become a power-law fluid. The incipient motion experiment also revealed that theincipient velocity of the cohesive sediment increases with the density in an exponential manner. Therefore, the incipientmotion is primarily related to the density, which is different from the case for non-cohesive sediment in which the incipi-ent motion is correlated with the diameter of sand particles instead. The incipient motion occurs in two different ways de-pending on the concentration of mud in the bottom. For sufficiently fine particles and a concentration lower than1.20 g/cm3, the cohesive sediment appears as fluidized mud, and the incipient motion is in the form of instability of aninternal wave. For a higher concentration, the cohesive sediment appears as general quasi-solid-mud, and the incipientmotion can be described by a series of extended Shields curves each with a different porosity for newly deposited alluvial mud.     

Effect of Berm Width and Elevation on Irregular Wave Run-Up

This paper discusses the effect of berm width and elevation of composite slope on irregular wave run-up. Based on the data obtained from model tests, the formula and distribution of irregular wave run-up on composite slope are derived. The changing of wind speed, width and elevation of the berm are considered comprehensively. The wave run-up with various exceedance probability can be es-timated utilizing the distribution curves of irregular wave run-up.     

Research on Characteristics of Density Current Under the Action of Waves

In this paper,the characteristics of density current under the action of waves are describedwith the help of flume experiment and theoretical analysis.The study shows that turbid water under the ac-tion of the waves can present three types of motion,i.e.significant stratification,fragile stratification andstrong mixing.The motion of turbid water presents significant stratification when(H/D)/△ρ/ρ~(1/2)≤4.5,generally this state is known as density current.The formulas of motionvelocity,thickness,and discharge of density current moving on horizontal bottom are derived by use of ba-sic equations such as momemtum equation,equation of energy conservation and continuity equation offluid.The time-average velocity and the thickness of density current under the action of waves have a rela-tionship with such parameters as relative density(△ρ/ρ),wave height(H),and water depth(D).Whenthese parameters are determined,the time-average thickness and motion velocity of density current are al-so determined.The relat