层流—紊流共存流场中岩溶裂隙网络演化过程的数值模拟方法

研究岩溶水系统的演化过程对许多资源与环境问题有着重要意义。岩溶地区水资源、油气资源的预测与开发、水土流失成因及防治等课题都与此密切相关。岩溶管道是地下水对裂隙的逐渐溶蚀扩展形成的,岩溶演化初期,所有裂隙宽度不大,流场整体成层流状态。随着溶蚀的进行,一部分裂隙优先扩大使其中的水流进入紊流状态,而另一部分裂隙中水流仍呈层流状态。文章提出一种数值方法,能够模拟层流-紊流共存流场的岩溶演化过程。利用蒙特卡洛(Monte-Carlo)方法模拟初始裂隙网络,通过非连续介质方法模拟裂隙网络中的渗流。裂隙扩展的速度通过岩石表面溶蚀速度经验公式计算,使用迭代方法求解层流-紊流共存条件下的水头非线性方程。构建了能够利用解析法求解的模型,把数值解和解析解的结果进行对比,验证了本研究的数值法及软件的可行性。

湍流边界层近壁区多个相干结构的数值模拟

从流动稳定性理论中的一般共振三波概念出发 ,提出一种湍流边界层近壁区多个相干结构的理论模型 ,采用高精度差分格式和Fourier谱展开相结合的方法 ,求解三维不可压Navier Stokes方程 ,直接数值模拟近壁湍流多个相干结构的演化问题 .

有向外力场作用下机器人路径规划方法

针对有向外力场作用下基于格点的机器人路径规划方法不能得到最优解的问题,采用了基于水平集方法的路径规划方法.将机器人的粒子跟踪转变为曲线的数值演化,通过解哈密尔顿-雅克比方程得到最优时间路径,改进后向路径追踪的数值计算精度保证了规划路径的可行性.仿真结果表明,算法可以有效应对复杂的有向外力场,并且在强外力场中仍保持路径的可行性.采用的连续路径规划方法可以突破传统机器人路径规划算法基于格点搜索的限制,并有效利用空间中存在的外力场.

张裂型盆地地热参数的垂向变化与琼东南盆地热流分布特征

海洋热流数据是开展海洋地球动力学研究和油气资源评价的基础数据.为深入认识琼东南盆地的地热特征,本文首先利用耦合沉积作用与岩石圈张裂过程的数值模型分析了张裂型盆地主要地热参数的垂向变化特征;并通过钻孔资料的详细分析,获得了琼东南盆地44口钻孔的热流数据;结合海底地热探针获取的热流数据,对琼东南盆地地热特征及其主要影响因素进行了简要分析.结果表明：沉积作用的热披覆效应对表层热流有较明显的抑制作用,由于沉积物生热效应与披覆效应的共同作用,同一钻孔处海底表层热流与钻孔深度3000~4000m处热流或与海底间的平均热流差异很小,可以一起用于分析琼东南盆地的热流分布特征;莺歌海组、乐东组热导率随深度变化小于黄流组及其下地层热导率的变化,钻孔沉积层平均热导率约为1.7 W·（m·K）-1,钻孔地层生热率一般低于2.5μW·m-3,平均生热率为1.34μW·m-3,平均地温梯度主要介于30~45℃/km,热流介于50~99mW·m-2,陆架区热流主要集中于60~70mW·m-2,深水区钻孔具有较高的地温梯度和热流值;从北部陆架与上陆坡区往中央坳陷带,热流值从50~70mW·m-2,增高为65~85mW·m-2,并且往东有升高趋势,在盆地东部宝岛凹陷、长昌凹陷与西沙海槽北部斜坡带构成一条热流值高于85mW·m-2的高热流带.进一步分析认为,琼东南盆地现今热流分布特征是深部热异常、强烈减薄岩石圈的裂后冷却作用、晚期岩浆热事件、地壳与沉积层的生热贡献以及沉积作用的热披覆效应等多种主要因素综合作用的结果.

山体隆升历史与地貌演化过程的数值模拟约束——以青藏高原东北缘河西走廊中段的周边年轻上升山地为例

地貌演化数值模型是分析构造-气候相互作用,再现地表形态演化过程的重要工具.本文利用地貌演化数值模型模拟了青藏高原东北缘河西走廊中段周边山体的演化过程.河西走廊气候干旱,外动力过程相似,但走廊中段南北两侧各山体间的形态特征存在明显差异(北侧山地:金塔南山东段、金塔南山西段、合黎山;南侧山地:榆木山).首先,各山体内选取典型区域分析了真实地形的特征指标(准周期性谷间距、地形高差与出水口数量);其次利用Landlab平台构建了典型区域的地貌演化数值模型,并通过数值实验正演了各山体地形的演化过程;最后将真实地形与不同时间节点的模拟地形进行比对,以准周期性谷间距等地形特征指标为判识标准,推断了河西走廊中段山体开始强烈隆升的年代和抬升速率.模拟结果与区域已发表的地质地貌证据在年轻上升山地(未达到稳定状态的幼年期地形)地区呈现较好的一致性,表明地貌演化数值模型结合地形特征指标是约束年轻上升山地强烈隆升年代和速率的可靠工具.本文提出的全新研究思路与方法有望在年轻上升山地推广并成为构造地貌学研究的常用手段.

基于水平集方法的时变流场中弱驱动多机器人队形保持方法

针对弱驱动多机器人群体编队控制和最优路径在环境流场影响下不能兼顾的问题,提出了将最优时间路径规划与领航者-跟随者编队控制相结合的队形保持方法。水平集方法得到领航者的最优时间路径,队形保持以此路径为基准。根据某一时刻领航者的位置和跟随者的可到达集合,通过求解队形优化函数来确定跟随者的位置,通过迭代求解该队形优化问题得到完整的路径。采用蒙特卡洛仿真了4种组合实验,结果显示基于水平集的路径规划方法和本研究提出的队形保持在路径的可行性、机器人完成任务的平均速度和队形的形状保持方面都有明显的优势。

Theoretical and numerical study on reinforcing effect of rock-bolt through composite soft rock-mass

Anchoring mechanism and failure characteristics of composite soft rock with weak interface usually exhibit remarkable difference from those in single rock mass.In order to fully understand the reinforcement mechanism of composite soft roof in western mining area of China,a mechanical model of composite soft rock with weak interface and rock bolt which considering the transverse shear sliding between different rock layers was established firstly.The anchoring effect was quantified by a factor defined as anchoring effect coefficient and its evolution equation was further deduced based on the deformation relationship and homogenized distribution assumption of stress acting on composite structure.Meanwhile,the numerical simulation model of composite soft rock with shear joint was prompted by finite element method.Then detailed analysis were carried out for the deformation features,stress distribution and failure behavior of rock mass and rock bolt near the joint under transverse load.The theoretical result indicates that the anchoring effect of rock-bolt through weak joint changes with the working status of rock mass and closely relates with the physical and geometric parameters of rock mass and rock bolt.From the numerical results,the bending deformation of rock bolt accurately characterized by Doseresp model is mainly concentrated between two plastic hinges near the shear joint.The maximum tensile and compression stresses distribute in the plastic hinge.However,the maximum shear stress appears at the positions of joint surface.The failure zones of composite rock are produced firstly at the joint surface due to the reaction of rock bolt.The above results laid a theoretical and computational foundation for further study of anchorage failure in composite soft rock.

Numerical simulation of wheel wear evolution for heavy haul railway

The prediction of the wheel wear is a fundamental problem in heavy haul railway. A numerical methodology is introduced to simulate the wheel wear evolution of heavy haul freight car. The methodology includes the spatial coupling dynamics of vehicle and track, the three-dimensional rolling contact analysis of wheel-rail, the Specht's material wear model, and the strategy for reproducing the actual operation conditions of railway. The freight vehicle is treated as a full 3D rigid multi-body model. Every component is built detailedly and various contact interactions between parts are accurately simulated, taking into account the real clearances. The wheel-rail rolling contact calculation is carried out based on Hertz's theory and Kalker's FASTSIM algorithm. The track model is built based on field measurements. The material loss due to wear is evaluated according to the Specht's model in which the wear coefficient varies with the wear intensity. In order to exactly reproduce the actual operating conditions of railway,dynamic simulations are performed separately for all possible track conditions and running velocities in each iterative step.Dimensionless weight coefficients are introduced that determine the ratios of different cases and are obtained through site survey. For the wheel profile updating, an adaptive step strategy based on the wear depth is introduced, which can effectively improve the reliability and stability of numerical calculation. At last, the wear evolution laws are studied by the numerical model for different wheels of heavy haul freight vehicle running in curves. The results show that the wear of the front wheelset is more serious than that of the rear wheelset for one bogie, and the difference is more obvious for the outer wheels. The wear of the outer wheels is severer than that of the inner wheels. The wear of outer wheels mainly distributes near the flange and the root; while the wear of inner wheels mainly distributes around the nominal rolling circle. For the outer wheel of front wheelset of each bogie, the development of wear is gradually concentrated on the flange and the developing speed increases continually with the increase of traveled distance.

Numerical simulation study of the failure evolution process and failure mode of surrounding rock in deep soft rock roadways

Based on the safety coefficient method,which assigns rock failure criteria to calculate the rock mass unit,the safety coefficient contour of surrounding rock is plotted to judge the distribution form of the fractured zone in the roadway.This will provide the basis numerical simulation to calculate the surrounding rock fractured zone in a roadway.Using the single factor and multi-factor orthogonal test method,the evolution law of roadway surrounding rock displacements,plastic zone and stress distribution under different conditions is studied.It reveals the roadway surrounding rock burst evolution process,and obtains five kinds of failure modes in deep soft rock roadway.Using the fuzzy mathematics clustering analysis method,the deep soft surrounding rock failure model in Zhujixi mine can be classified and patterns recognized.Compared to the identification results and the results detected by geological radar of surrounding rock loose circle,the reliability of the results of the pattern recognition is verified and lays the foundations for the support design of deep soft rock roadways.

Registration of image feature points using differential evolution

This paper introduces a robust global nonlinear optimizer—differential evolution(DE),which is a simple evolution algorithm to search for an optimal transformation that makes the best alignment of two sets of feature points.To map the problem of matching into the framework of DE,the objective function is proportional to the registration error which is measured by Hausdorff distance,while the parameters of transformation are encoded in floating-point as the functional variables.Three termination criteria are proposed for DE.A simulation of 2-dimensional point sets and a similarity transformation are presented to compare the robustness and convergence properties of DE with genetic algorithm’s (GA).And the registration of an object and its contour model have been demonstrated by using of DE to natural images.

Exact Soliton Solutions to a Generalized Nonlinear Schrdinger Equation

The （1＋1）-dimensional F-expansion technique and the homogeneous nonlinear balance principle have been generalized and applied for solving exact solutions to a general （3＋1）-dimensional nonlinear Schr6dinger equation （NLSE） with varying coefficients and a harmonica potential. We found that there exist two kinds of soliton solutions. The evolution features of exact solutions have been numerically studied. The （3＋1）D soliton solutions may help us to understand the nonlinear wave propagation in the nonlinear media such as classical optical waves and the matter waves of the Bose-Einstein condensates.

The Impacts of Initial Perturbations on the Computational Stability of Nonlinear Evolution Equations

The impacts of initial perturbations on the computational stability of nonlinear evolution equations for non-conservative difference schemes and non-periodic boundary conditions are studied through theoretical analysis and numerical experiments for the case of onedimensional equations.The sensitivity of the difference scheme to initial values is further analyzed.The results show that the computational stability primarily depends on the form of the initial values if the difference scheme and boundary conditions are determined.Thus,the computational stability is sensitive to the initial perturbations.

Impact Dynamics on Granular Plate

In this paper, a theoretical and numerical study on the impact of a rubber solid on the free surface of a granular plate is presented, showing a simulation of an aircraft wheel on impact with a flexible landing surface. This physical action, when we use a theological approach, becomes a fundamental parameter to investigate wear and tear, and consequently strength to micro and macro pavements failure. The study has developed initially from a microscopic point of view and subsequently on macroscale. The effects are strictly linked with material degradation associated with damage evolution. The problem is developed by energetic approach on an elastic-plastic element using the functional energy containing two contributions, bulk and surface. The model simulates the behaviour of flexible runway pavements during the landing phase.

EVOLUTION OF MOIST POTENTIAL VORTICITY DURING A WARM-ZONE HEAVY RAINFALL EVENT IN THE PEARL RIVER DELTA

First,based on routine meteorological data,the synoptic characteristics of a heavy warm-sector rainfall that occurred on June 13,2008 in the Pearl River Delta were analyzed.Second,a mesoscale numerical model,Weather Research and Forecasting(WRFV2.2),was used to simulate the heavy rainfall. Diagnostic analyses were done of moist potential vorticity(MPV)for its horizontal components(MPV2) and vertical components(MPV1)based on the simulation results of WRFV2.2 to identify the mechanism of the rainfall development.The results showed that the heavy rainfall occurred when there were high MPV1 in the upper levels and low MPV1 and high MPV2 in the lower levels.Disturbances of high MPV1 in the upper levels came from the southwest or northwest,those of low MPV1 in the lower levels came from the southwest,and those of high MPV2 came from the south.Disturbances of low MPV1 at low levels were the direct cause of convective instability.Enhanced vertical shear of meridional wind led to increased MPV2 at lower levels,strengthened baroclinicity,and active warm and wet flows.These distributions of MPV helped to trigger the release of unstable energy and produce warm-sector heavy rainfall.As it integrates the evolution of dynamic and thermal fields,MPV is able to reveal the development of this heavy rainfall effectively.

Numerical study on the evolution of the shock-accelerated SF_6 interface:Influence of the interface shape

The early phases of the shock interaction process on two-dimensional interfaces with different shapes are numerically investigated in this study,which are closely related to the shock refraction and reflection,vorticity production and transport.The numerical method employs an adaptive unstructured quadrilateral mesh,which can capture the wave pattern and interface evolution very well.Simulations are carried out under the conditions of an incident shock Mach number of 1.2 and the light/heavy (air/SF 6) interface.Five different shapes are considered in the simulations:rectangle,ellipse,diamond and two kinds of triangle.The results show that the interfacial shapes can influence the wave patterns particularly on the shape and evolution of refracted shock waves.The generation and the distribution of vorticity on the interfaces with five different shapes also have dissimilarities.The circulation deposition on five interfaces is quantitatively investigated and compared with theoretical model.A good agreement is found between the numerical results and the predictions by the theoretical model.Some characteristic scales of the interface are tracked.Under the influence of nonlinear-acoustic effect and vorticity effect,the interfaces present different evolution modes.

Instabilities and pattern evolution in a vertically heated annulus

The convection in an annular container with heated bottom,cooled top and insulated side walls are studied by both linear instability analysis and direct numerical simulation.The onset of convection is investigated by linear stability analysis and corresponding pattern selection mechanisms are discussed.The nonlinear evolution of different flow patterns and the convective heat transfer are simulated.The transition to oscillatory flow is also given by stability analysis where the base flow is a steady three dimensional flow.The stability predictions are in good agreement with the numerical simulations,including both the growth rate and the dimensionless frequency.

Massive Scalar Field Evolution in the Dyadosphere Spacetime of Charged Black Hole

Scalar field quasinormal modes in the dyadosphere spacetime of charged black hole are studied by using the third-order WKB approximation. From numerical results obtained, we find that the scalar field mass u plays an important role in studying the quasinormal frequencies. With the scalar field mass increases, the real parts increase and the magnitudes of the imaginary parts decrease. Partieulary, these change are almost linearly.

Characteristics of nonlinear evolution of wavepackets in boundary layers

The nonlinear evolution of a finite-amplitude disturbance in a 3-D supersonic boundary layer over a cone was investigated recently by Liu et al. using direct numerical simulation （DNS）. It was found that certain small-scale 3-D disturbances amplified rapidly. These disturbances exhibit the characteristics of second modes, and the most amplified components have a well- defined spanwise wavelength, indicating a clear selectivity of the amplification. In the case of a cone, the three-dimensionality of the base flow and the disturbances themselves may be responsible for the rapid amplification. In order to ascertain which of these two effects are essential, in this study we carried out DNS of the nonlinear evolution of a spanwise localized disturbance （wavepacket） in a flat-plate boundary layer. A similar amplification of small-scale disturbances was observed, suggesting that the direct reason for the rapid amplification is the three-dimensionality of the disturbances rather than the three-dimensional nature of the base flow, even though the latter does alter the spanwise distribution of the disturbance. The rapid growth of 3-D waves may be attributed to the secondary instability mechanism. Further simulations were performed for a wavepacket of first modes in a supersonic boundary layer and of Tollmien-Schlichting （T-S） waves in an incompressible boundary layer. The re- suits show that the amplifying components are in the band centered at zero spanwise wavenumber rather than at a finite spanwise wavenumber. It is therefore concluded that the rapid growth of 3-D disturbances in a band centered at a preferred large spanwise wavenumber is the main characteristic of nonlinear evolution of second mode disturbances in supersonic boundary layers.

Numerical analysis on shock-cylinder interaction using immersed boundary method

The problem of shock interaction with a rigid circular cylinder has been investigated using a compressible immersed boundary method coupled with high-order weighted-essentially non-oscillatory(WENO) scheme.First,the accuracy of the developed code is validated.Then,influences of the incident shock Mach number on the flow-field structure and dynamic drag coefficient,as well as time evolution of the flow field are studied.For different shock Mach number,the flow structure shows very different features.At a given dimensionless time,both the normalized shock detachment distance and the normalized vertical distance from the highest point of the primary reflected shock to the centerline of the cylinder decreases with increasing shock Mach number.However,location of the upper triple point varies non-monotonically with shock Mach number.For a case with given shock Mach number,the trajectory of the upper triple point and the time evolution of the normalized vertical distance from the highest point of the primary reflected shock to the centerline of the cylinder can both be predicted by linear correlation.Nevertheless,the time evolution of the normalized shock detachment distance is biased to be non-linear.Meanwhile,time evolution of force exerted on the cylinder is quite unsteady for a case with given shock Mach number and given cylinder diameter.For small shock Mach number,there exists a negative valley,and it disappears when the incident shock Mach number increases to a large value,e.g.,1.7.Furthermore,correlations to predict the occurrence of the peak drag and its value under different shock Mach numbers have been proposed.

Interface Characteristics and Evolution Mechanism of W/CuCrZr in Hot Melt Explosion Welding

Explosion welding was carried out on the basis of vacuum hot melt W/CuCrZr composite plate.Metallurgical microscope,scanning electron microscope and energy dispersive X-ray spectroscope were used to observe the microscopic morphology of the bonding interface.At the same time,combined with finite element calculations,the evolution mechanism of the interface of the hot melt explosion welded W/CuCrZr composite plate was explored.The results show that the interface bonding of the hot melt explosion welded W/CuCrZr composite plate is good and there is a cross-melting zone with 3–8μm in thickness,but cracks are developed on the W side.The numerical simulation reproduces the changes of pressure,stress,strain and internal energy at the bonding interface in the process of hot melt explosion welding.The location of the crack generated in the experiment coincides with the high stress position calculated by numerical simulation.The high pressure and high temperature near the hot melt explosion welding interface further promote the bonding of the interface.