基于Fluent的空气射流切削式反循环钻头参数优化

针对黄土地层,提出利用空气代替钻井液的反循环气体喷射技术。但现有技术对反循环气体喷射钻头切削土体的能力没有具体优化讨论。利用fluent软件,优化反循环喷射钻头结构参数,对不同底喷孔直径、底喷孔数量、底喷孔扩大段直径、底喷孔扩大段数量、喷射孔数量5个结构参数进行模拟分析。分析结果表明,在进风流量一定的情况下,增加底喷孔的数量或增大底喷孔的直径会削弱钻头的空气射流切削能力和反循环能力。底喷孔直径为3 mm,底喷孔个数为2个,底喷孔扩大段直径为8 mm,底喷孔扩大段长度为5 mm,内喷孔数量为5个时,反循环喷射钻头喷射能力和反循环能力最优。     

金融领域的随机建模与基于软件R的Monte Carlo模拟(2):Cox-Ross-Rubinstein模型

主要讨论了 Cox-Ross-Rubinstein (CRR) 模型和广义的CRR模型,并研究了如何基于CRR模型和广义的CRR模型利用Monte Carlo模拟计算资产价格以及期权价值.     

金融领域的随机建模与基于软件R的Monte Carlo模拟(5):Black-Scholes的世界

主要研究了Black-Scholes模型.与Black-Scholes期权定价公式相比,将再次强调和证实利用R软件Monte Carlo模拟的强大作用.     

Stability of a planet in the HD 41004 binary system

The Hill stability criterion is applied to analyse the stability of a planet in the binary star system of HD 41004 AB, with the primary and secondary separated by 22 AU, and masses of 0.7 M_⊙ and 0.4 M_⊙, respectively. The primary hosts one planet in an S‐type orbit, and the secondary hosts a brown dwarf (18.64 M_J) on a relatively close orbit, 0.0177 AU, thereby forming another binary pair within this binary system. This star‐brown dwarf pair (HD 41004 B+Bb) is considered a single body during our numerical calculations, while the dynamics of the planet around the primary, HD 41004 Ab, is studied in different phase‐spaces. HD 41004 Ab is a 2.6 M_J planet orbiting at the distance of 1.7 AU with orbital eccentricity 0.39. For the purpose of this study, the system is reduced to a three‐body problem and is solved numerically as the elliptic restricted three‐body problem (ERTBP). The Hill stability function is used as a chaos indicator to configure and analyse the orbital stability of the planet, HD 41004 Ab. The indicator has been effective in measuring the planet's orbital perturbation due to the secondary star during its periastron passage. The calculated Hill stability time series of the planet for the coplanar case shows the stable and quasi‐periodic orbits for at least ten million years. For the reduced ERTBP the stability of the system is also studied for different values of planet's orbital inclination with the binary plane. Also, by recording the planet's ejection time from the system or collision time with a star during the integration period, stability of the system is analysed in a bigger phase‐space of the planet's orbital inclination, ≤ 90°, and its semimajor axis, 1.65–1.75 AU. Based on our analysis it is found that the system can maintain a stable configuration for the planet's orbital inclination as high as 65° relative to the binary plane. The results from the Hill stability criterion and the planet's dynamical lifetime map are found to be consistent with each other. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

洪湖流域自然农耕条件下营养盐沉积输移演化模拟研究

选择洪湖流域为研究区域,应用分布式流域模型SWAT,从长期演化的角度出发,对流域营养盐输移演化进行模拟。模拟的边界条件设置为自然农耕背景,探讨百年尺度洪湖流域营养盐输移演变规律。对TN、TP的模拟分汛期与非汛期两个时间段,分别进行200年的模拟运算。模拟结果显示模型可以很好的体现流域营养盐沉积输移演变的基本特征。TN、TP浓度变化有明显的季节特征,随时间变化特征与沉积记录相一致。模拟结果的时间序列分析,TN、TP显示出不同的变化趋势。流域模型为百年时间尺度上流域营养盐沉积输移演化的自然趋势研究,提供了很好评价依据。     

A comparison of seven methods for the inverse modelling of groundwater flow. Application to the characterisation of well catchments

Inverse modelling is a key step in groundwater-related hydrological studies. Several inversion techniques were developed during the last decades, but hardly any comparison between them was presented. We compare seven modern inverse methods for groundwater flow: the Regularised Pilot Points Method (both the estimation, RPPM-CE, and the Monte Carlo (MC) simulation variants, RPPM-CS), the MC variant of the Representer Method (RM), the Sequential Self-Calibration Method (SSC), the Moment Equations Method (MEM), the Zonation Method (ZM) and a non-iterative Semi-Analytical Method (SAM). These methods are applied to a two-dimensional synthetic example, depicting steady-state groundwater flow around a pumping well. Their relative performance is assessed in terms of their ability to characterise the log-transmissivity and hydraulic head fields and to predict the extent of the well catchment, both for a mildly and a strongly heterogeneous transmissivity field. The main conclusions drawn from the comparison are: (1) MC-based methods (RPPM-CS, SSC and RM) yield very similar results, regardless the degree of heterogeneity and despite they use different parameterisation schemes and objective functions; (2) statistical moments of the target quantities provided by MEM and RPPM-CE are similar to those of MC-based methods; (3) ZM and SAM are negatively affected by strong heterogeneity; and (4) in general, observed differences between the performances of all methods are not very large. MC-based inverse methods need considerably more CPU time than the other tested approaches. An advantage of MC-based methods is that they allow computing the posterior probability distribution of the target quantities, which can be directly fed to probabilistic risk-assessment procedures.     

The effect of intensity of turbulence in umbrella cloud on tephra dispersion during explosive volcanic eruptions: Experimental and numerical approaches

During an explosive volcanic eruption, tephra fall out from the umbrella region of the eruption cloud to the ground surface. We investigated the effect of the intensity of turbulence in the umbrella cloud on dispersion and sedimentation of tephra by performing a series of laboratory experiments and three dimensional (3-D) numerical simulations. In the laboratory experiments, spherical glass-bead particles are mixed in stirred water with various intensities of turbulence, and the spatial distribution and the temporal evolution of the particle concentration are measured. The experimental results show that, when the root-mean-square of velocity fluctuation in the fluid (W_rms) is much greater than the particle terminal velocity (v_t), the particles are homogeneously distributed in the fluid, and settle at their terminal velocities at the base of the fluid where turbulence diminishes. On the other hand, when W_rms is as small as or smaller than v_t, the particle concentration increases toward the base of the fluid during settling, which substantially increases the rate of particle settling. The results of the 3-D simulations of eruption cloud indicate that W_rms is up to 40 m/s in most of the umbrella cloud even during a large scale plinian eruption with a magma discharge rate of 10⁹ kg/s. These results suggest that relatively coarse pyroclasts (more than a few mm in diameter) tend to concentrate around the base of the umbrella cloud, whereas fine pyroclasts (less than 1/8 mm in diameter) may be distributed homogeneously throughout the umbrella cloud during tephra dispersion. The effect of the gradient of particle concentration in the umbrella cloud explains the granulometric data of the Pinatubo 1991 plinian deposits.     

Physical and numerical simulations of the presence of a forsterite transitional zone at high temperature and pressure: insight for scCO2 geological storage

Geological sequestration is one of the most effective ways to reduce greenhouse gases, such as carbon dioxide (CO₂). The deep oceanic crust dominated by ultrabasic rock could store CO₂ permanently. However, the storage mechanism has not been thoroughly understood because of the limited amount of research and experiments conducted. This study explored the reactive mechanisms of water–rock–gas in an ultrabasic system under different conditions. Forsterite, the most dominant mineral found in ultrabasic reservoirs, was used to conduct laboratory physical simulation experiments. Two experimental systems were designed including an scCO₂–forsterite–water system and an N₂–forsterite–water system. All experiments were performed for 1000 h at an experimental temperature of 150°C and a pressure of 150 bar, respectively, to mimic the geological conditions. The liquid products from the experiments were analysed by inductively coupled plasma-optical emission spectrometry, whereas the solid samples were analysed by scanning electron microscopy with energy disperse spectroscopy. Results showed that: (1) in the early stage during scCO₂/N₂–forsterite–water interaction, forsterite was dissolved with a reactive transitional zone forming on the surface, which caused H⁺ to enter into the silicate framework and accelerated the reaction; (2) in the N₂ system, the dissolution of forsterite was inhibited by the Mg²⁺ concentration after reaching its saturation in the late stage; and (3) in the scCO₂ system, magnesite was precipitated as a secondary mineral during the late stage, which promoted the dissolution of forsterite. As a result, the degree of dissolution of forsterite in the scCO₂ system was far higher than in the N₂ system. The experimental results are consistent with the numerical simulation using TOUGHREACT, a geochemical simulation procedure, which showed that CO₂ promotes the dissolution of forsterite greatly at high temperature and pressure.     

组合动力锚水动力特性数值模拟研究

作为海上浮式结构物的一种新型锚固基础,动力锚具有自安装、不受水深影响、适用范围广的特点。在动力锚的基础上研发的组合动力锚结合了动力锚自安装和板锚法向受荷的特点,具有安装快速、适用多种类型海床、承载效率高等性质。组合动力锚在水中自由下落时的水动力学特性（下落速度、方向稳定性等）会受到锚链、尾翼宽度和助推器质量等因素影响。若下落速度过小或方向稳定性过差,则会影响锚的安装成功率。采用计算流体动力学方法模拟流体对锚的冲击和锚在水中自由下落过程,以优化组合动力锚的尾翼尺寸;其次研究锚链作用力和助推器质量对组合动力锚下落速度和偏角的影响规律。计算结果表明：组合动力锚的拖曳阻力系数为0.45左右,尾翼宽度最优尺寸为翼板宽度的1.25倍。连接在锚眼处的锚链会减小组合动力锚的下落速度并加剧锚的偏转,需综合锚的下落速度和偏角来确定锚在水中下落高度。