A comprehensive study on the parameters setting in smoothed particle hydrodynamics (SPH) method applied to hydrodynamics problems

Smoothed particle hydrodynamics (SPH) is a meshfree, Lagrangian particle method which has advantages in handling solids with extremely large deformation. Like any other numerical methods, cares must be taken to ensure its desirable accuracy and stability through considering several correction techniques in calculation. The selection of values for parameters in those correction approaches is a key step in SPH simulation, which is always difficult for new beginners to deal well with effectively. This paper examines the common inconsistency and instability problems in SPH method and studies its computational efficiency when applied to hydrodynamics problems with material strength like soil column collapse. We analyzed in detail how the correction techniques mitigate these inconsistency and instability problems. Also, the numerical testing results associate with different values for the parameters used in the correction techniques are provided for better understanding the influence of these parameters and for finding out the desirable values. It is found that (1) the SPH method is easily subjected to an inconsistency problem in the boundary area due to the boundary deficiency, and it can be treated well by adopting “virtual particles” contributing to the particle summations. (2) The numerical oscillation in SPH simulation can be mitigated effectively by artificial viscosity with the suggested parameter values. (3) The tension cracking treatment, artificial viscosity and artificial stress work well in removing the tensile instability problem in SPH method. In addition, the nearest neighboring particle searching (NNPS) algorithm, spacing ratio, smoothing length and time step influence the efficiency and accuracy of SPH method significantly. It is shown that SPH method with suggested parameters values can produce a very good result compared with the experimental result.     

Radial pulsations of helium stars with masses from 10 to 50<Emphasis Type="Italic">M</Emphasis><Subscript>⊙</Subscript>

We have performed hydrodynamic calculations of the radial pulsations of helium stars with masses 10M_⊙ ≤ M ≤ 50M_⊙, luminosity-to-mass ratios 5 × 10³L_⊙/M_⊙ ≤ L/M ≤ 2.5 × 10⁴L_⊙/M_⊙, and effective temperatures 2 × 10⁴ K ≤ T_eff ≤ 10⁵ K for helium and heavy-element mass fractions of Y=0.98 and Z=0.02, respectively. We show that the high-temperature boundary of the instability region for radial pulsations at L/M ? 10⁴L_⊙/M_⊙ extends to T_eff≈10⁵ K. The amplitude of the velocity variations for outer layers is several hundred km s^?1, while the brightness variations in the B band of the UBV photometric system are within the range from several hundredths to half a magnitude. At constant luminosity-to-mass ratio, the radial pulsation period is determined only by the effective temperature of the star. In the ranges of luminosity-to-mass ratios 10⁴L_⊙/M_⊙ ≤ L/M ≤ 2 × 10⁴L_⊙/M_⊙ and effective temperatures 5 × 10⁴ K ≤ T_eff ≤ 9 × 10⁴ K, the periods of the radial modes are within 6 min ?Π?10³ min.     

利用空间测地技术结果估计地球平均下地幔粘性(英)

基于IEC－4G冰期后地壳反弹模型,和地球上Laurentia,Fennoscandia,Antarctica,andGreenland四大冰盖最近21000年以来的冰融参数,计算了对地球最大主转动惯量的影响△I33,并进而由现代空间测地技术观测资料分析得到的地球自转非潮汐加速项为约束,估计了地球平均下地幔（670km以下）粘性vLM为（0．9～2．5）·1022Pas,这个结果表明了vLM应具有1022Pas量级．     

Bayesian analysis of stage–fall–discharge models for gauging stations affected by variable backwater

Parsimonious stage–fall–discharge rating curve models for gauging stations subject to backwater complications are developed from simple hydraulic theory. The rating curve models are compounded in order to allow for possible shifts in the hydraulics when variable backwater becomes effective. The models provide a prior scientific understanding through the relationship between the rating curve parameters and the hydraulic properties of the channel section under study. This characteristic enables prior distributions for the rating curve parameters to be easily elicited according to site‐specific information and the magnitude of well‐known hydraulic quantities. Posterior results from three Norwegian and one American twin‐gauge stations affected by variable backwater are obtained using Markov chain Monte Carlo simulation techniques. The case studies demonstrate that the proposed Bayesian rating curve assessment is appropriate for developing rating procedures for gauging stations that are subject to variable backwater. Copyright © 2009 John Wiley & Sons, Ltd.     

基于循环注蒸气和烟气激励的海上稠油开采技术

通过实验和油藏数值模拟,研究循环注蒸气和烟气激励的加热降黏、气体溶解降黏、气体扩大加热腔与减少热损失、蒸气与气体协同增产效应等开采机理,优化注入蒸气和烟气温度、烟气蒸气体积比、周期蒸气和气体注入量等参数,建议注入蒸气与烟气的温度高于300℃,在大气压下烟气蒸气体积比为0.4～0.6.在渤海稠油油田开展4口井循环注蒸气和烟气激励开采的现场试验.结果表明,循环注蒸气和烟气激励开采后稠油井平均日产油量达到冷采的3倍以上,最高日产油量超过100m3,周期产油量达到1.0×104 m3左右.基于循环注蒸气和烟气激励的海上开采稠油技术可以提高海上稠油开采速度和采收率,有潜力成为未来海上稠油油田的主要热采开发方式.     

改进的HASM-AD算法及在空间变量模拟的应用分析

高精度曲面建模(HASM)可以显著提高空间曲面模拟的精度,但是计算速度低限制了该模型的进一步应用。为了提高HASM模型的计算效率,本文对HASM-AD算法作了改进,通过在计算过程中为采样点添加索引,避免了计算过程中对采样点信息的重复查找操作;同时,在遍历独立计算单元时实时计算第一类、第二类基本量及克式符号,避免了全局存储上述变量所需要的额外内存消耗。数值试验表明,由于将全局线性方程组求解问题转化为局部独立计算单元(5×5栅格)内的方程组求解,改进的HASM-AD算法显著提高了计算效率,同时降低了模型运行过程中的内存消耗。最后,本文以全国陆地降水空间分布模拟分析作为实例,验证改进的HASM-AD算法模拟精度及计算效率,模拟结果表明,改进的HASM-AD算法模拟结果精度优于其他HASM算法(以HASM-PCG为例),并且计算效率优势更为明显,实现全国10km分辨率的降水分布模拟耗时仅为4s。表明改进的HASM-AD算法提高了计算速度,并且适于大尺度的空间变量模拟应用。     

中国省域煤矿事故时空分异特征及影响因素研究

煤矿事故给国家和人民带来巨大的生命财产损失,科学把握煤矿事故发生的宏观规律及区域煤矿安全的影响因素有助于促进煤矿安全生产。本文尝试采用GIS空间分析方法开展区域煤矿事故的研究,基于空间自相关性、时空演化规律方法研究了2006-2015年中国省域煤矿事故的时空分异特点,并分析了省域煤矿安全的影响因素及相对风险度。主要结论：①中国各省区之间煤矿事故存在明显的空间自相关性,但随着各省区煤矿安全生产形势改善情况的不同,其自相关性在不断下降;②以空间重心转移曲线分析中国煤矿事故发生的时空演化规律,10年间煤矿事故空间重心向东北方向转移。③重特大事故频发的省区由于监管措施的有力执行往往具有较低的煤矿百万吨死亡率;平均煤矿生产能力低的省区具有较高的煤矿百万吨死亡率;④建立基于解释变量的贝叶斯空间模型,以煤炭百万吨死亡率为指标评估中国各省区煤矿事故相对风险度,中国煤矿安全生产相对危险度存在明显的空间分布规律且各省区差异较大,相对危险度较高的省区主要集中在南方地区,包括湖北、福建等,危险度较低的省区主要集中在中东部地区,包括山西、内蒙古等。     

CFD-based numerical analysis of a variable cross-section cylinder

Using ANSYS-CFX, a general purpose fluid dynamics program, the vortex-induced vibration (VIV) of a variable cross-section cylinder is simulated under uniform current with high Reynolds numbers. Large eddy simulation (LES) is conducted for studying the fluid-structure interaction. The vortex shedding in the wake, the motion trajectories of a cylinder, the variation of drag and lift forces on the cylinder are analyzed. The results show that the vortices of variable cross-section cylinder are chaotic and are varying along the cylinder. In places where cross-sections are changing significantly, the vortices are more irregular. The motion trail of the cylinder is almost the same but irregular. The drag and lift coefficients of the cylinder are varying with the changes of diameters.