An efficient and robust approach to predict ship self-propulsion coefficients

Achieving a reliable and accurate numerical prediction of the self-propulsion performance of a ship is still an open problem that poses some relevant issues. Several CFD methods, ranging from boundary element methods (BEM) to higher-fidelity viscous Reynolds averaged Navier–Stokes (RANS) based solvers, can be used to accurately analyze the separate problems, i.e. the open water propeller and the hull calm water resistance. However, when the fully-coupled self-propulsion problem is considered, i.e. the hull advancing at uniform speed propelled by its own propulsion system, several complexities rise up. Typical flow simplifications adopted to speed-up the simulations of the single analysis (hull and propeller separately) lose their validity requiring a more complex solver to tackle the fully-coupled problem. The complexity rises up further when considering a maneuver condition. This aspect increases the computational burden and, consequently, the required time which becomes prohibitive in a preliminary ship design stage.The majority of the simplified methods proposed in literature to include propeller effects, without directly solve the propeller flow, in a high-fidelity viscous solver are not able to provide all the commonly required self-propulsion coefficients. In this work, a new method to enrich the results from a body force based approach is proposed and investigated, with the aim to reduce as much as possible the computational burden without losing any useful result. This procedure is tested for validation on the KCS hull form in self-propulsion and maneuver conditions.     

椭球谐和球谐系数之间一个简单的转换关系

本文推导的椭球谐系数和球谐系数相互之间转换关系的核心思想是在ε~2量级下利用Legendre函数的正交性,从球谐系数求解的积分表示出发,将积分中的椭球坐标变量与球坐标变量相互转换,从而得出椭球谐系数与球谐系数之间的转换关系。本文导出的转换关系有以下优点:①对于第二类Legendre函数的计算采用Laurent级数表示,使计算第二类Legendre函数更为简单;②保留了ε~2量级下,导出的转换关系相比文献[2]的形式更简单,满足物理大地测量边值问题线性化的要求;③顾及了余纬和归化余纬的区别。     

平面无结节网衣水阻力系数的试验研究

通过物理模型试验,统计分析得到了无结节网衣的水阻力系数在不同状态下随雷诺数、倾角和布置状态的变化规律,并通过网衣沉子的组合试验验证了所得经验公式的正确性和适用性,为深水网箱等渔具的水动力计算提供了依据.     

Analytical method of radiation by a wave energy device with dual rectangular floating bodies

The system with one floating rectangular body on the free surface and one submerged rectangular body has been applied to a wave energy conversion device in water of finite depth. The radiation problem by this device on a plane incident wave is solved by the use of an eigenfunction expansion method, and a new analytical expression for the radiation velocity potential is obtained. The wave excitation force is calculated via the known incident wave potential and the radiation potential with a theorem of Haskind employed. To verify the correctness of this method, an example is computed respectively through the bound element method and analytical method. Results show that two numerical methods. are in good agreement, which shows that the present method is applicable. In addition, the trends of hydrodynamic coefficients and wave force are analyzed under different conditions by use of the present analytical method.     

巴丹吉林沙漠北缘拐子湖流沙下垫面近地层湍流强度和陆面过程特征

利用巴丹吉林沙漠北缘拐子湖流沙下垫面2013年7、10月和2014年1、4月的湍流通量资料,计算并分析了研究区近地层湍流强度,同时针对风速分量、温度、水汽和CO₂归一化标准差随稳定度的变化关系和总体输送系数等陆面过程特征进行分析。结果表明：（1）风速各分量的湍流强度均随风速的增加逐渐减小,风速处于2 m·s^-1以下时湍流发展最为旺盛。湍流强度主要由水平方向风速分量决定,垂直方向风速的作用较小,且近中性和不稳定层结利于湍流的发展。与其他地区相比,平坦且没有建筑物的沙漠地区,机械湍流较弱,湍流强度相应较小。（2）风速各分量的归一化标准差与稳定度（z/L）均满足1/3次方函数规律,其中垂直方向风速分量的拟合曲线方程较好。（3）动量输送系数C_d具有明显的夏季高、冬季低的变化状态且各月的日变化形态均呈夜间低、日间高的循环形态。热量输送系数C_h的不同月份日变化间并没有明显的排列次序,且日出日落前后具有明显的波动。不稳定层结时,C_d和C_h均随风速的增加逐渐减小;稳定层结时,C_d和C_h均随着风速的增加逐渐上升。     

Compass系统与GPS系统L1频段信号干扰仿真分析

全球导航卫星系统间的频率兼容性研究已经成为一个热点问题。首先介绍了国际电联建议的卫星导航系统间干扰分析方法;然后根据该方法对Compass系统和GPS系统在L1频段的信号进行了干扰仿真分析,并比较了系统间互干扰和系统自干扰的影响。仿真计算结果表明,除了互用信号L1C和B1C外,其余信号间的自干扰在总干扰中占主导,互干扰的影响相对较小。     

VTI介质qP波方程高精度有限差分算子

波动方程有限差分法是一种使用广泛的地震波数值模拟方法.但是有限差分法本身固有存在着数值频散问题,会降低地震波场模拟的精度与分辨率.为了克服常规有限差分算子的数值频散,本文针对VTI介质地震波数值模拟问题,构造了频率-空间域qP波波动方程高精度有限差分优化算子,根据最优化理论中高斯-牛顿法确定了高精度有限差分算子的优化系数.利用常规差分算子和高精度优化差分算子对归一化相速度的频散关系精度进行了对比分析,并对均匀各向同性介质和均匀VTI介质中的qP波地震波场进行了有限差分数值模拟,通过频散关系精度分析和波场数值模拟结果表明：有限差分优化算子具有较高的波场数值模拟精度,有效压制了传统有限差分算子数值模拟中的数值频散现象,提高了有限差分算子精度,为VTI介质频率-空间域qP波正演模拟奠定了基础.     

Investigation of highly efficient algorithms for solving linear equations in the discontinuous deformation analysis method

Large‐scale engineering computing using the discontinuous deformation analysis (DDA) method is time‐consuming, which hinders the application of the DDA method. The simulation result of a typical numerical example indicates that the linear equation solver is a key factor that affects the efficiency of the DDA method. In this paper, highly efficient algorithms for solving linear equations are investigated, and two modifications of the DDA programme are presented. The first modification is a linear equation solver with high efficiency. The block Jacobi (BJ) iterative method and the block conjugate gradient with Jacobi pre‐processing (Jacobi‐PCG) iterative method are introduced, and the key operations are detailed, including the matrix‐vector product and the diagonal matrix inversion. Another modification consists of a parallel linear equation solver, which is separately constructed based on the multi‐thread and CPU‐GPU heterogeneous platforms with OpenMP and CUDA, respectively. The simulation results from several numerical examples using the modified DDA programme demonstrate that the Jacobi‐PCG is a better iterative method for large‐scale engineering computing and that adoptive parallel strategies can greatly enhance computational efficiency. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatial and multivariate analysis of geochemical data from metavolcanic rocks in the Ben Nevis area, Ontario

A study of the lithogeochemistry of metavolcanics in the Ben Nevis area of Ontario, Canada has shown that factor analysis methods can distinguish lithogeochemical trends related to different geological processes, most notably, the principal compositional variation related to the volcanic stratigraphy and zones of carbonate alteration associated with the presence of sulphides and gold. Auto- and cross-correlation functions have been estimated for the two-dimensional distribution of various elements in the area. These functions allow computation of spatial factors in which patterns of multivariate relationships are dependent upon the spatial auto- and cross-correlation of the components. Because of the anisotropy of primary compositions of the volcanics, some spatial factor patterns are difficult to interpret. Isotropically distributed variables such as CO ₂ are delineated clearly in spatial factor maps. For anisotropically distributed variables (SiO ₂ ), as the neighborhood becomes smaller, the spacial factor maps becomes better. Interpretation of spatial factors requires computation of the corresponding amplitude vectors from the eigenvalue solution. This vector reflects relative amplitudes by which the variables follow the spatial factors. Instability of some eigenvalue solutions requires that caution be used in interpreting the resulting factor patterns. A measure of the predictive power of the spatial factors can be determined from autocorrelation coefficients and squared multiple correlation coefficients that indicate which variables are significant in any given factor. The spatial factor approach utilizes spatial relationships of variables in conjunction with systematic variation of variables representing geological processes. This approach can yield potential exploration targets based on the spatial continuity of alteration haloes that reflect mineralization.List of symbols c _i Scalar factor that minimizes the discrepancy between andU _i - D Radius of circular neighborhood used for estimating auto- and cross-correlation coefficients - d Distance for which transition matrixU is estimated - d _ij Distance between observed valuesi andj - E Expected value - E _i Row vector of residuals in the standardized model - F(d _ij) Quadratic function of distanced _ij F(d _ij)=a+bd _ij+cd _ij ² - L Diagonal matrix of the eigenvalues ofU - _i Eigenvalue of the matrixU;ith diagonal element ofL - N Number of observations - p Number of variables - Q Total predictive power ofU - R Correlation matrix of the variables - R _0j Variance-covariance signal matrix of the standardized variables at origin;j is the index related tod andD (e.g.,j=1 ford=500 m,D=1000 m) - R _1j Matrix of auto- and cross-correlation coefficients evaluated at a given distance within the neighborhood - R _m ² Multiple correlation coefficient squared for themth variable - S _i Column vectori of the signal values - s _k ² Residual variance for variablek - T _i Amplitude vector corresponding toV _i;ith row ofT=V ^–1 - T Total variation in the system - U Nonsymmetric transition matrix formed by post-multiplyingR ₀₁ ^–1 byR _ij - U _i Componenti of the matrixU, corresponding to theith eigenvectorV _i;U _i= _iV_iT_i - U* _i ComponentU _i multiplied byc _i - U _ij Sum of componentsU _i+U _j - V _i Eigenvector of the matrixU;ith column ofV withUV=VL - w Weighting factor; equal to the ratio of two eigenvalues - X _i Random variable at pointi - x _i Value of random variable at pointi - y _i Residual ofx _i - Z _i Row vectori for the standardized variables - z _i Standardized value of variable     

Ultraviolet absorption cross-sections of chloro and chlorofluoro-methanes at stratospheric temperatures

Absorption cross-sections of nine halomethanes (CCl₄, CHCl₃, CH₂Cl₂, CH₃Cl, CFCl₃, CF₂Cl₂, CF₃Cl, CHFCl₂, and CHF₂Cl), measured between 174 and 250 nm for temperatures ranging from 225 to 295 K, are presented with uncertainties ranging from 2 to 4% and compared with previous determinations made for comparable temperature ranges.The largest temperature effect which takes place near the absorption threshold, decreases the absorption cross-section up to 50% for highly chlorinated methanes, but is negligible for molecules highly stabilized by hydrogen and/or fluorine. Extrapolated values for temperatures of aeronomical interest are presented, as well as parametrical formulas which give absorption cross-section values for given wavelength and temperature ranges.