Third Order Effect of Rotation on Stellar Oscillations of a B Star

We aim at investigating the effect of rotation up to the third order in the angular velocity of a star on the p and g modes, based on the formalism developed by Soufi et al. Our ultimate goal is the study of oscillations of β Cephei stars which are often rapidly rotating stars. Our results show that the third-order perturbation formalism presented by Soufi et al. should be corrected for some missing terms and some misprints in the equations. As a first step in our study of β Cephei stars, we quantify by numerical calculations the effect of rotation on the oscillation frequencies of a uniformly rotating zero-age main-sequence star with 12 M<,??>. For an equatorial velocity of 100km s-1, it is found that the second-and third-order corrections for (l, m)=(2, 2), for instance, are of the order of 0.01% of the frequency for radial order n=6 and reaches up to 0.5% for n=14.     

Estimates of size for the large octopus Macroctopus maorum from measures of beaks in prey remains

Macroctopus maorum, also known as Octopus maorum and Pinnoctopus cordiformis, grows to 12 kg and is the largest octopod in Australasia. In New Zealand, this species features in the diets of seals and albatrosses. Beaks of depredated octopus resist digestion and accumulate in the stomach until regurgitated or defecated. Beaks of M. maorum were differentiated reliably from those of other octopods in prey remains for octopus > c.100 g. Previous studies estimated the size of depredated M. maorum from measures of lower beaks using equations derived either from a small sample of M. maorum or from related species. Here, allometric equations to estimate ventral mantle length (4.6–27.5 cm), dorsal mantle length (5.1–29.0 cm), total length (20–198 cm), and body mass (36–12 000 g) were derived from measures of both upper and lower beaks of 90 specimens of M. maorum. These regressions were imprecise and reflected innate variability in octopus morphometrics, a characteristic reflected in the results of analyses of published data for M. maorum from northern New Zealand and for two other large octopods, Octopus vulgaris and Enteroctopus magnificus.     

南黄海中部海底沉积物原位声速与物理性质相关关系

基于在南黄海中部获得的海底沉积物原位声速数据,分析讨论了原位声速与密度、含水量、孔隙比、孔隙度等沉积物物理性质参数的相关关系。通过回归分析建立了海底沉积物原位声速预测方程。结果表明,原位声速与上述物理参数之间具有良好的相关性,相关系数r均大于0.93。对比分析了原位声速预测方程与甲板声速预测方程的差异,两者最小相差15.8 m/s,最大相差31.3 m/s,平均差值约为22.6 m/s。将原位声速预测方程与其他研究者建立的声速预测方程进行了对比,结果表明不同声速预测方程存在明显差异,初步分析了存在差异的原因。     

Elastic wave modelling method based on the displacement–velocity fields: an improving nearly analytic discrete approximation

In this paper, we present an improvement to our previously published nearly analytic discrete method (NADM) to solve acoustic and elastic wave equations. We compare the numerical errors of the improved NADM with the original NADM and also the fourth-order Lax–Wendroff correction and present examples of three-component wave fields in 2D transversely isotropic media with strong velocity contrasts. Comparing with the original NADM, we find that the improved method requires significantly less storage space and can increase the time accuracy from second order of the original NADM to fourth order, while the space accuracy remains the same as that of the original one. Theoretical analyses and numerical results suggest that our improved NADM is suitable for large-scale numerical modelling as it can effectively suppress numerical dispersion and source-generated noises caused by discretising wave equations when too-coarse grids are used.     

Comparing predicted and observed ground motions from subduction earthquakes in the Lesser Antilles

This brief article presents a quantitative analysis of the ability of eight published empirical ground-motion prediction equations (GMPEs) for subduction earthquakes (interface and intraslab) to estimate observed earthquake ground motions on the islands of the Lesser Antilles (specifically Guadeloupe, Martinique, Trinidad, and Dominica). In total, over 300 records from 22 earthquakes from various seismic networks are used within the analysis. It is found that most of the GMPEs tested perform poorly, which is mainly due to a larger variability in the observed ground motions than predicted by the GMPEs, although two recent GMPEs derived using Japanese strong-motion data provide reasonably good predictions. Analyzing separately the interface and intraslab events does not significant modify the results. Therefore, it is concluded that seismic hazard assessments for this region should use a variety of GMPEs in order to capture this large epistemic uncertainty in earthquake ground-motion prediction for the Lesser Antilles.     

甘肃玛曲大水金矿成矿过程中资源量生长的数量规律

大水金矿床是著名的超大型矿床.在借鉴前人矿床地质研究成果的基础上, 通过对矿体规模数据的多重分形特征及吨-品位关系研究, 揭示了成矿过程中资源量生长的数量规律.矿体长度服从简单分形分布, 矿床资源量的吨-品位关系能用指数函数或直线方程更好地拟合, 矿床与矿体资源量双对数图、吨-品位的双对数图的多重分形特征表明: 在分形构造成矿空间中, 特定成矿阶段的矿床资源量的生长方式是自组织非线性增殖过程, 与自组织有关的随机性使吨-品位关系趋于指数关系; 在整个成矿过程中, 矿床资源量的生长方式是各阶段矿体资源量的线性叠加.进一步根据熵最大化原理及效用理论导出了与统计结果吻合的矿床吨-品位方程, 并得到柴家庄金矿矿体地质数据的验证.      

Theoretical analysis and simulation for avoiding antenna oscillations without distributed resistance

The so-called Wu and King antenna pattern is widely used in GPR because of its simple design and construction features. The main disadvantage is its limited efficiency due to transmitter energy losses which occur through lumped resistors. Based on the analysis of the electromagnetic field behaviour across the antenna, it is possible to replace the effect of the resistors by either storing the energy of the electric pulses, or damp them by means of one matched resistor, which will theoretically improve the efficiency of the antenna. In this paper we provide a theoretical analysis using a modified transmission line model together with simulation based on delayed potentials among other electromagnetic software, and measurement results using an impulse transmitter with fast MOSFET switches and a matched resistor that support this idea.     

Inverse analysis of stochastic moment equations for transient flow in randomly heterogeneous media

We present a nonlinear stochastic inverse algorithm that allows conditioning estimates of transient hydraulic heads, fluxes and their associated uncertainty on information about hydraulic conductivity (K) and hydraulic head (h  ) data collected in a randomly heterogeneous confined aquifer. Our algorithm is based on Laplace-transformed recursive finite-element approximations of exact nonlocal first and second conditional stochastic moment equations of transient flow. It makes it possible to estimate jointly spatial variations in natural log-conductivity (Y=lnK)$(Y=\ln K)$, the parameters of its underlying variogram, and the variance–covariance of these estimates. Log-conductivity is parameterized geostatistically based on measured values at discrete locations and unknown values at discrete “pilot points”. Whereas prior values of Y at pilot point are obtained by generalized kriging, posterior estimates at pilot points are obtained through a maximum likelihood fit of computed and measured transient heads. These posterior estimates are then projected onto the computational grid by kriging. Optionally, the maximum likelihood function may include a regularization term reflecting prior information about Y. The relative weight assigned to this term is evaluated separately from other model parameters to avoid bias and instability. We illustrate and explore our algorithm by means of a synthetic example involving a pumping well. We find that whereas Y and h can be reproduced quite well with parameters estimated on the basis of zero-order mean flow equations, all model quality criteria identify the second-order results as being superior to zero-order results. Identifying the weight of the regularization term and variogram parameters can be done with much lesser ambiguity based on second- than on zero-order results. A second-order model is required to compute predictive error variances of hydraulic head (and flux) a posteriori. Conditioning the inversion jointly on conductivity and hydraulic head data results in lesser predictive uncertainty than conditioning on conductivity or head data alone.     

Multilayer shallow water flow using lattice Boltzmann method with high performance computing

A multilayer lattice Boltzmann (LB) model is introduced to solve three-dimensional wind-driven shallow water flow problems. The multilayer LB model avoids the expensive Navier–Stokes equations and obtains stratified horizontal flow velocities as vertical velocities are relatively small and the flow is still within the shallow water regime. A single relaxation time BGK method is used to solve each layer coupled by the vertical viscosity forcing term. To increase solution stability, an implicit step is suggested to obtain flow velocities. The main advantage of using the LBM is that after selecting appropriate equilibrium distribution functions, the LB algorithm is only slightly modified for each layer and retains all the simplicities of the LBM within the high performance computing (HPC) environment. The performance of the parallel LB model for the multilayer shallow water equations is investigated on CPU-based HPC environments using OpenMP. We found that the explicit loop control with cache optimization in LBM gives better performance on execution time, speedup and efficiency than the implicit loop control as the number of processors increases. Numerical examples are presented to verify the multilayer LB model against analytical solutions. We demonstrate the model’s capability of calculating lateral and vertical distributions of velocities for wind-driven circulation over non-uniform bathymetry.