Numerical Analysis of Wellbore Instability in Gas Hydrate Formation During Deep-Water Drilling

Gas hydrate formation may be encountered during deep-water drilling because of the large amount and wide distribution of gas hydrates under the shallow seabed of the South China Sea. Hydrates are extremely sensitive to temperature and pressure changes, and drilling through gas hydrate formation may cause dissociation of hydrates, accompanied by changes in wellbore temperatures, pore pressures, and stress states, thereby leading to wellbore plastic yield and wellbore instability. Considering the coupling effect of seepage of drilling fluid into gas hydrate formation, heat conduction between drilling fluid and formation, hydrate dissociation, and transformation of the formation framework, this study established a multi-field coupling mathematical model of the wellbore in the hydrate formation. Furthermore, the influences of drilling fluid temperatures, densities, and soaking time on the instability of hydrate formation were calculated and analyzed. Results show that the greater the temperature difference between the drilling fluid and hydrate formation is, the faster the hydrate dissociates, the wider the plastic dissociation range is, and the greater the failure width becomes. When the temperature difference is greater than 7℃, the maximum rate of plastic deformation around the wellbore is more than 10%, which is along the direction of the minimum horizontal in-situ stress and associated with instability and damage on the surrounding rock. The hydrate dissociation is insensitive to the variation of drilling fluid density, thereby implying that the change of the density of drilling fluids has a minimal effect on the hydrate dissociation. Drilling fluids that are absorbed into the hydrate formation result in fast dissociation at the initial stage. As time elapses, the hydrate dissociation slows down, but the risk of wellbore instability is aggravated due to the prolonged submersion in drilling fluids. For the sake of the stability of the wellbore in deep-water drilling through hydrate formation, the drilling fluid with low temperatures should be given priority. The drilling process should be kept under balanced pressures, and the drilling time should be shortened.

Radio-frequency measurements of energetic-electron-driven emissions using high-frequency magnetic probe on XuanLong-50 spherical torus

A high-frequency magnetic probe is designed and developed on the XuanL ong-50(EXL-50)spherical torus to measure high-frequency magnetic field fluctuation.The magnetic loop,radio filters,radio-frequency limiter,and data acquisition system of the probe are comprehensively examined.The fluctuation data from the EXL-50 plasma are analyzed in the time–frequency domain using fast Fourier transforms.Moreover,distinct high-frequency instabilities are detected using this diagnostic system.In particular,significant frequency chirping is observed,which is consistent with the bumpon-tail drive instability predicted using the Berk–Breizman model.

Effects of different numerical algorithms on simulation of chemical dissolution-front instability in fluid-saturated porous rocks

Many scientific and engineering problems need to use numerical methods and algorithms to obtain computational simulation results because analytical solutions are seldom available for them.The chemical dissolution-front instability problem in fluid-saturated porous rocks is no exception.Since this kind of instability problem has both the conventional(i.e.trivial)and the unconventional(i.e.nontrivial)solutions,it is necessary to examine the effects of different numerical algorithms,which are used to solve chemical dissolution-front instability problems in fluid-saturated porous rocks.Toward this goal,two different numerical algorithms associated with the commonly-used finite element method are considered in this paper.In the first numerical algorithm,the porosity,pore-fluid pressure and acid/solute concentration are selected as basic variables,while in the second numerical algorithm,the porosity,velocity of pore-fluid flow and acid/solute concentration are selected as basic variables.The particular attention is paid to the effects of these two numerical algorithms on the computational simulation results of unstable chemical dissolution-front propagation in fluid-saturated porous rocks.The related computational simulation results have demonstrated that:1)the first numerical algorithm associated with the porosity-pressure-concentration approach can realistically simulate the evolution processes of unstable chemical dissolution-front propagation in chemical dissolution systems.2)The second numerical algorithm associated with the porosity-velocity-concentration approach fails to simulate the evolution processes of unstable chemical dissolution-front propagation.3)The extra differential operation is the main source to result in the failure of the second numerical algorithm.

Non-Linear Numerical Simulation of Coal and Gas Outburst as a Coupling Instability Problem of Solid-Fluid Biphase Media

Based on the theory of continuum mechanics of multi-pbase media, a mathematical model and non-linear FEM equation of the coupling instability problem of solid-fluid biphase media for coal-methane outburst under finite deformation are established. The critical conditions of the surface instability are presented as the singularity of the total stiffness matrices of the coal body for coal-methaue outburst. That means the deformtion or the coal body emerges bifurcatiou phenomena. The numerical simulation of a typical outburst is made.

A note for the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions

In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensional discrete model of elastic wave motion, the module of reflection factor will be greater than 1 in high frequency band when artificial wave velocity is greater than 1.5 times the ratio of discrete space step to discrete time step. Based on the proof, the frequency band in which instability occurs is discussed in detail, showing such high-frequency waves are meaningless for the numerical simulation of wave motion.

A Numerical Study on Piezoelectric Energy Harvesting by Combining Transverse Galloping and Parametric Instability Phenomena

This paper aims to numerically investigate the effects of parametric instability on piezoelectric energy harvesting from the transverse galloping of a square prism. A two degrees-of-freedom reduced-order model for this problem is proposed and numerically integrated. A usual quasi-steady galloping model is applied, where the transverse force coefficient is adopted as a cubic polynomial function with respect to the angle of attack. Time-histories of nondimensional prism displacement, electric voltage and power dissipated at both the dashpot and the electrical resistance are obtained as functions of the reduced velocity. Both, oscillation amplitude and electric voltage, increased with the reduced velocity for all parametric excitation conditions tested. For low values of reduced velocity, 2:1 parametric excitation enhances the electric voltage. On the other hand, for higher reduced velocities, a 1:1 parametric excitation(i.e., the same as the natural frequency) enhances both oscillation amplitude and electric voltage. It has been also found that, depending on the parametric excitation frequency, the harvested electrical power can be amplified in 70% when compared to the case under no parametric excitation.

Numerical studies of the Kelvin-Hemholtz instability in a coronal jet

Kelvin-Hemholtz（K-H）instability in a coronal EUV jet is studied via 2.5D MHD numerical simulations.The jet results from magnetic reconnection due to the interaction of the newly emerging magnetic field and the pre-existing magnetic field in the corona.Our results show that the Alfv e′n Mach number along the jet is about 5–14 just before the instability occurs,and it is even higher than 14 at some local areas.During the K-H instability process,several vortex-like plasma blobs with high temperature and high density appear along the jet,and magnetic fields have also been rolled up and the magnetic configuration including anti-parallel magnetic fields forms,which leads to magnetic reconnection at many X-points and current sheet fragments inside the vortex-like blob.After magnetic islands appear inside the main current sheet,the total kinetic energy of the reconnection outflows decreases,and cannot support the formation of the vortex-like blob along the jet any longer,then the K-H instability eventually disappears.We also present the results about how the guide field and flux emerging speed affect the K-H instability.We find that a strong guide field inhibits shock formation in the reconnecting upward outflow regions but helps secondary magnetic islands appear earlier in the main current sheet,and then apparently suppresses the K-H instability.As the speed of the emerging magnetic field decreases,the K-H instability appears later,the highest temperature inside the vortex blob gets lower and the vortex structure gets smaller.

Inherent Numerical Instability in Computing Invariant Measures of Markov Chains

Invariant measures of Markov chains in discrete or continuous time with a countable set of states are characterized by its steady state recurrence relations. Exemplarily, we consider transition matrices and Q-matrices with upper bandwidth n and lower bandwidth 1 where the invariant measures satisfy an (n + 1)-order linear difference equation. Markov chains of this type arise from applications to queueing problems and population dynamics. It is the purpose of this paper to point out that the forward use of this difference equation is subject to some hitherto unobserved aspects. By means of the concept of generalized continued fractions (GCFs), we prove that each invariant measure is a dominated solution of the difference equation such that forward computation becomes numerically unstable. Furthermore, the GCF-based approach provides a decoupled recursion in which the phenomenon of numerical instability does not appear. The procedure results in an iteration scheme for successively computing approximants of the desired invariant measure depending on some truncation level N. Increasing N leads to the desired solution. A comparison study of forward computation and the GCF-based approach is given for Q-matrices with upper bandwidth 1 and 2.

Numerical simulation and mechanism analysis of freak waves

A numerical wave model based on the modified fourth-order nonlinear Schroe dinger equation （mNLSE） in deep water was developed to simulate the formation of freak waves and a standard split - step, pseudo-spectral method was used to solve the equation. The validation of the model is firstly verified, then the simulation of freak waves was performed by changing sideband conditions, and the variation of wave energy was also analyzed in the evolution. The results indicate that Benjamin - Feir instability （ sideband instability） is an important mechanism for freak wave formation.

Universality of Periodic Oscillation Induced in Side Branch of a T-Junction in Numerical Simulation

The flow instability through the side branch of a T-junction is analyzed in a numerical simulation. In a previous experimental study, the authors clarified the mechanism of fluid-induced vibration in the side branch of the T-junction in laminar steady flow through the trunk. However, in that approach there were restrictions with respect to extracting details of flow behavior such as the flow instability and the distribution of wall shear stress along the wall. Here the spatial growth of the velocity perturbation at the upstream boundary of the side branch is investigated. The simulation result indicates that a periodic velocity fluctuation introduced at the upstream boundary is amplified downstream, in good agreement with experimental result. The fluctuation in wall shear stress because of the flow instability shows local extrema in both the near and distal walls. From the numerical simulation, the downstream fluid oscillation under a typical condition has a Strouhal number of 1.05, which approximately agrees with the value obtained in experiments. Therefore, this periodic oscillation motion is a universal phenomenon in the side branch of a T-junction.

High-Frequency Trading:Deception and Consequences

This commentary is based on the work of Cooper,Davis,and Van Vliet(2016)and the commentary focuses on what problem high-frequency trading poses.It lists key literature on high-frequency trading that is missing and points out that the poker analogy to defend deception in financial markets is weak and misleading.The article elaborates on the negative impact created by spoofing and quote stuffing,the two typical deceptive practices used by high-frequency traders.The recent regulations regarding high-frequency trading,in response to the“Flash Crash”of 2010,are preventive,computerized and more effective.They reflect ethical requirements to maintain fair and stable financial markets.

Investigations of high-frequency induction hardening process for piston rod of shock absorber

The microhardness of piston rods treated with different induction hardening processes was tested. The experimental results reveal that the depth of the hardened zone is proportional to the ratio of the moving speed of the piston rod to the output power of the induction generator. This result is proved correct through the Finite Element Method (FEM) simulation of the thermal field of induction heating. From tensile and impact tests, an optimized high frequency induction hardening process for piston rods has been obtained, where the output power was 82%×80 kW and the moving speed of workpiece was 5364 mm/min. The piston rods, treated by the optimized high frequency induction hardening process, show the best comprehensive mechanical performance.

Numerical investigation on evolution of cylindrical cellular detonation

Cylindrical cellular detonation is numerically investigated by solving two- dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.

NUMERICAL STUDY OF BUOUNCY- AND THERMOCAPILLARY-DRIVEN FLOWS IN A CAVITY

Thermocapillary-and buoyancy-driven convection in open cavities with differentially heated endwalls is investigated by numerical solutions of the two- dimensional Navier-Stokes equations coupled with the energy equation. We studied the thermocapillary and buoyancy convection in the cavities, filled with low-Prandtl- number fluids, with two aspect-ratios A=1 and 4, Grashof number up to 10~5 and Reynolds number |Re|≤10~4. Our results show that thermocapillary can have a quite significant effect on the stability of a primarily buoyancy-driven flow, as well as on the flow structures and dynamic behavior for both additive effect (i.e., positive Re) and opposing effect (i.e., negative Re).

Improving the numerical stability of the MAGIC model

MAGIC is one of the most widely used models for forecasting long-term acidification. The model's code, however, has been experiencing numerical instability, though this might not be widely known to the public users. The major instability comes from the analytical solution to two cubic equations for calculating SO42- concentration and the exchangeable fraction of Al on the soils. The mathematical algorithm for calculating the concentration of SO42- from a quadratic equation is also found instable. This paper is aimed at improving the instability above through proved numerical algorithms.

Numerical and Experimental Investigation of the Aero-Hydrodynamic Effect on the Behavior of a High-Speed Catamaran in Calm Water

In this paper,the effect of water and air fluids on the behavior of a planing catamaran in calm water was studied separately in calm water by using experimental and numerical methods.Experiments were conducted in a towing tank over the Froude number range of 0.49–2.9 with two degrees of freedom.The model vessel displacement of 5.3 kg was implemented in experimental tests.Craft behavior was evaluated at the displacements of 5.3,4.6,and 4 kg by using the numerical method.The numerical simulation results for the hull’s resistance force were validated with similar experimental data.The fluid volume model was applied to simulate two-phase flow.The SST k-ωturbulence model was used to investigate the effect of turbulence on the catamaran.The results showed that in the planing mode,the contribution of air to pressure resistance increased by 55%,40%,and 60%at the mentioned displacements,whereas the contribution of air to friction resistance was less than 15%on average.The contribution of the air to the total lift force at the abovementioned displacements exceeded 70%,60%,and 50%in the planing mode but was less than 10%in the displacement mode.At the displacements of 5.3 and 4 kg,the area under the effect of maximum pressure moved around the center of gravity and caused porpoising longitudinal instability at the Froude numbers of 2.9 and 2.4,respectively.However,at the displacement of 4.6 kg,this effect did not occur,and the vessel maintained its stability.

库水升降作用下灰岩岸坡稳定性数值研究

基于现场调研与采样,采用室内力学试验测定灰岩力学参数与干湿循环下强度劣化规律;基于UDEC建立库岸边坡数值模型,研究了三峡库区反倾岸坡、顺层岸坡、近水平层状岸坡及含危岩体岸坡在库水升降作用下的失稳演化过程及消落带劣化深度与岩层厚度对岸坡稳定性的影响规律。结果表明:库水升降作用下岸坡的失稳过程可以分为裂隙萌生发育、裂隙发展贯通和坡体失稳破坏阶段;随着消落带劣化深度的增加,4种岸坡稳定性系数降低;随着岩层厚度的增加,4种岸坡的稳定性系数升高。近水平层状岸坡稳定性系数受以上2种因素影响最大,而随着干湿循环次数增加发生的变化最小;当劣化深度大于10 m以及岩层厚度大于5 m时,含危岩体岸坡稳定性系数随影响因素水平变化而产生的变化量较小。长远来看,受库水影响的4种典型岩质岸坡稳定性从大到小依次为反倾岸坡、近水平层状岸坡、顺层岸坡、含危岩体岸坡。

复合板起皱失稳数值模拟模型的建立

复合板被广泛应用于汽车工业、航空航天、国防军工等各大领域,起皱失稳是复合板塑性成形过程中影响零件质量的主要问题之一,而目前能够准确复现复合板起皱形貌的数值模拟模型的建立方式尚不明晰。因此本文以方板对角拉伸试验(YBT)作为对比验证对象,在ABAQUS软件中采用特征值屈曲分析分别与动力显式和静力通用两种常规算法相结合的buckle-dynamic算法以及buckle-static算法,对选用不同建模方式和单元类型建立的模型,以起皱高度、厚向位移和应变云图等为评估依据,进行了起皱形貌模拟对比分析。结果表明,采用buckle-static算法对复合层建模方式建立的常规壳模型分析所得结果的综合性能最佳。若模拟工况复杂或计算不易收敛,分析方法也可选用buckle-dynamic算法。

压气机内部旋转不稳定的研究综述

旋转不稳定是压气机工作在高负荷近失速工况时的一种常见现象。研究旋转不稳定在降低压气机工作噪声、减小流致振动以及保障航空发动机稳定工作等方面具有重要意义。首先对旋转不稳定现象进行了回顾,详细讨论了旋转不稳定的特征。其次,重点调研了旋转不稳定的起源和机理,将旋转不稳定产生的原因归纳为叶尖泄漏流、涡脱落以及流动剪切等类别。此外,回顾了模拟旋转不稳定的数值方法,讨论了多种流动控制手段对旋转不稳定的作用效果。最后,对旋转不稳定的研究现状进行了总结,对未来的研究趋势进行了展望。

运煤巷道富水泥灰岩破碎带注浆技术及数值模拟研究

运煤巷道受富水泥灰岩破碎带的影响,易发生失稳变形,本文以某运煤巷道为工程背景,结合大量隧道注浆加固工程经验,制定出针对富水泥灰岩破碎带巷道的施工注浆方案,并以实际工程参数为依据,结合FLAC3D软件进行数值模拟,对巷道注浆加固效果进行验证,为类似工程问题的设计与施工提供一定的参考依据。