Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. © 2013 AIP Publishing LLC.

Trace theory

The concept of traces has been introduced for describing non-sequential behaviour of concurrent systems via its sequential observations. Traces represent concurrent processes in the same way as strings represent sequential ones. The theory of traces can be used as a tool for reasoning about nets and it is hoped that applying this theory one can get a calculus of the concurrent processes anologous to that available for sequential systems. The following topics will be discussed: algebraic properties of traces, trace models of some concurrency phenomena, fixed-point calculus for finding the behaviour of nets, modularity, and some applications of the presented theory.

谓词μ演算和模态图的语义一致性

模态图是谓词μ演算的一种有效的图形表示形式。证明了谓词μ演算和模态图的语义一致性，详细讨论了谓词μ演算公式、嵌套谓词等式系和模态图之间的关系，并给出了一种优化的从线性公式到嵌套谓词等式系的转换算法。

牦牛粪便对青藏高原东部高寒草甸土壤养分系统的影响

川西北高寒草甸位于青藏高原东部地区，是我国四大牧区之一，也是长江和黄河等江河源区的重要水源涵养功能区。近几十年来，大量牦牛粪便被牧民作为生物能源、肥料或者食用菌产业的原料而利用，为草原生态系统的养分平衡增加了生态风险。鉴于在青藏高原地区针对牦牛粪便的相关研究尚未见报道，本文从粪便和土壤养分相互变化趋势的角度出发，研究了粪便在夏秋两季的分解状况和粪便其下及其周围土壤的养分变化。通过此研究，将有助于了解牦牛粪便在自然状态下的分解速率以及粪便对土壤养分及土壤微生物生物量的影响，为粪便对草地土壤生态系统的贡献提供概念性认识，同时也为高寒草甸草地这种脆弱生态系统的可持续管理提供理论依据。 针对牦牛粪便自身养分释放及其对土壤在时间和空间尺度上的影响，本文通过模拟牛粪堆积，在不同时间和固定区域内对牛粪和土壤进行了夏秋两季的采样测定，分析了牛粪及土壤NO3--N、NH4+-N、速效K、无机P、有机C、全N 和全P 含量随时间的变化趋势。得到如下结果： 1. 在研究区域内，牛粪对草地生态系统具有较强的养分（N、P）贡献能力。据初步统计，其估值大致为氮素699~932 kg ha-1，P 素为110~147 kg ha-1。牛粪（湿重、干重）在夏秋季节的分解速率具有较大差异，夏季显著快于秋季。夏季牛粪湿重、干重在2 个月左右之后分别降为初期的14%和24%，而秋季同期分别降为35%和52%。牛粪养分（NO3--N、NH4+-N、速效K、无机P、有机C、全N 和全P）的分解速率在夏季也要显著快于秋季。秋季经过2 个月左右的分解之后，牛粪以上含量分别降至初始态的32%、60%、36%、64%、58%、63%和43%，远高于夏季的同期水平。 2. 在不同季节，牛粪周围的土壤养分（NO3--N、NH4+-N、速效K 和无机P）含量变化随时间改变呈现相同的趋势。但是，牛粪周围不同远近下土壤养分随时间的动态变化幅度略有不同。粪下土壤养分含量随时间变化波动较大，距粪便越远，随时间变化的变幅越为平缓。总体来说，在夏季由于粪便分解较快，对土壤养分的持续作用时间不及秋季，秋季粪便分解变慢，表现出养分的缓释特征。其次，牛粪对粪下土壤影响的持续时间也长于对周围土壤的作用时间。 3. 粪便对土壤养分（NO3--N、NH4+-N、速效K、无机P）影响的范围在不同季节具有差异。夏季要高于秋季，但对周围土壤养分影响的持续时间低于秋季。在夏季短期内，牛粪对土壤NO3--N 和速效K 含量的影响范围能够超过30cm，而对NH4+-N 和无机P 的影响范围则介于10cm 和30cm 之间。在2 个月左右之后，牛粪对周围土壤养分的影响能力基本消失。在秋季，牛粪对周围土壤养分影响范围难以达到30cm 处。粪便在夏季对其下土壤和周围土壤的有机C、全N 和全P 含量并无显著影响，但在秋季能显著增加其下土壤有机C、全N 和全P 含量。 4. 牛粪在秋季对土壤SMB-C、SMB-N 和SMB-P 含量的影响能够持续2 个月以上，由于秋冬季节牛粪分解缓慢，因此推断这种效应持续时间至少能够1 年左右。另一方面，牛粪在秋季对土壤SMB-C、SMB-N 和SMB-P 含量的影响范围主要集中在其下土壤，而周围的影响效应并不明显。 The grassland on the eastern fringe of Qinghai-Tibetan Plateau was one of the four greatest pasture zones in our country and the main water conservation function zones in the hesastream of Yangtse River and Huanghe river. Rencent years, lots of dung in this area was used as biological energy, fertilizer or material of fungoid growing, leading to high risk of nutrient banlance in grassland ecosystem. In view of the researches on the impact of yak dung in this area are relatively rare, the present study focused on the relationship of dung and soil nutrient transformation in sunmer and autumn, which could profoundly illuminate the mechanism of dung decomposition and the effect of dung on soil chemical properties and soil microbe biomass. The present study also contributed to a basic understand and provided scientific management in the high-frigid ecosystem. Decomposition of yak dung and its effect on soil chemical properties in eastern grassland of Qinghai-Tibetan Plateau were determined. The study simulated the real dung pats, took dung and soil samples at different time and fixed-point in summer and autumn. The samples were analysed for NO3--N, NH4+-N, available K, inorganic P, total organic C (TOC), total N (TN), total P (TP). It was concluded that: 1. In study area, the yak dung supplied to ecosystem substantial nutrient. It is estimated that the N contribution of dung was approximately 699~932 kg ha-1, P contribution was approximately 110~147 kg ha-1. The rate of yak decomposition was more rapid in summer than autumn, the wet and dry weight of yak decreased to 14% and 24% respectively after 2 months when dung excreted in summer, with 35% and 52% in autumn. The content of NO3--N、NH4+-N、available K、inorganic P、TOC、TN and TP in dung decomposed more rapid in summer too. After 2 months when dung excreted in autumn, the content of above nutrient decreased to 32%、60%、36%、64%、58%、63% and 43% respectively，which were significantly higher than summer. 2. The content of NO3--N、NH4+-N、available K and inorganic P in soil around dung had the same transformation trend in each season, whereas it was distinguishing at different gradient of distance from dung, the nutrient in soil below dung had the most significant change while the more far from dung, the less change in soil. It was concluded that the yak dung had prolong impacts on soil in autumn compared with summer, besides, it aslo showed that the yak dung had protract effect on soil below dung compared with soil around dung. 3. The yak dung had expansive impact on soil around dung in summer whereas had relatively short effect compared with autumn. In short-term at summer, there was a significant increase about the content of NO3--N and available K around 30cm radius from dung pat while the content of NH4+-N and inorganic P between 10cm to 30cm. After 2 months, the impact almost disappeared. In autumn, the effect was hard to reach 30cm. The yak dung had no significant effects on the content of TOC、TN and TP in soil below or around dung in summer whereas there was a obvious increase in soil below dung pat in autumn. 4. The duration of effect of yak dung on soil microbial biomass(SMB) C、N and P was at least 2 months, maybe even more than 1 year. On the other hand, the impact of dung on SMB-C、SMB-N and SMB-P mainly acted on soil below dung while no obious effect on soil around dung.

Mercury(II) ion sensing with PVC-matrix based membrane sensors

Thiosemicarbazone derivatives have been used as ion carriers for the preparation of PVC-matrix based mercury(II)-selective membrane sensors. The electrodes give near-Nernstian responses in the linear concentration range of 1.0×10-1-5.0×10-6 M with detection limits of the order of 10-6 M. The stable potentiometric signals are obtained within a short time period of 20-25s. The effect of different plasticizers has been studied and dioctylsebacate (DOS) found to give a better response in comparison to other plasticizers. Selectivity coefficient values (log KPotHg,M) have been evaluated using fixed interference method. Better selectivity for mercury(II) ions is observed over many of the monovalent (Na+, K+ and NH4+) and divalent ions (Mg2+, Ca2+, Zn2+, Pb2+, Ni2+, Co2+, etc.). The sensors have also been used as indicator electrodes in potentiometric titration of mercury(II) ions with EDTA and its determination in synthetic water samples.

WAVELENGTH POSITIONING ERRORS AND TRUE DETECTION LIMIT IN RELATION TO MULTICOMPONENT ANALYSIS TECHNIQUES FOR THE CORRECTION OF LINE INTERFERENCES IN INDUCTIVELY COUPLED PLASMA ATOMIC EMISSION-SPECTROMETRY

The present paper deals with the evaluation of the relative error (DELTA(A)) in estimated analyte concentrations originating from the wavelength positioning error in a sample scan when multicomponent analysis (MCA) techniques are used for correcting line interferences in inductively coupled plasma atomic emission spectrometry. In the theoretical part, a quantitative relation of DELTA(A) with the extent of line overlap, bandwidth and the magnitude of the positioning error is developed under the assumption of Gaussian line profiles. The measurements of eleven samples covering various typical line interferences showed that the calculated DELTA(A) generally agrees well with the experimental one. An expression of the true detection limit associated with MCA techniques was thus formulated. With MCA techniques, the determination of the analyte and interferent concentrations depend on each other while with conventional correction techniques, such as the three-point method, the estimate of interfering signals is independent of the analyte signals. Therefore. a given positioning error results in a larger DELTA(A) and hence a higher true detection limit in the case of MCA techniques than that in the case of conventional correction methods. although the latter could be a reasonable approximation of the former when the peak distance expressed in the effective width of the interfering line is larger than 0.4. In the light of the effect of wavelength positioning errors, MCA techniques have no advantages over conventional correction methods unless the former can bring an essential reduction ot the positioning error.

工业机器人的小脑模型控制

本文为动力学控制工业机器人提出了一种综合学习算法,这种学习算法可将以前所学的信息用于新的控制输入.这种控制方法不需要事先知道机器人动力学,它易于应用于特殊的控制问题或修改以适应实际系统中的变化,控制方法在时间上是有效的,且很适合于定点实现.学习控制算法的有效性通过4自由度的直接驱动机器人前两个关节在重复运动中的计算机仿真实验得到了验证.

工业机器人机械手的小脑模型自适应控制

本文为动力学控制工业机器人机械手提出一种综合控制算法。该控制算法,利用小脑模型算术计算机模块模拟机器人机械手的动力学方程并计算实现期望运动所需力矩作为前馈力矩控制项;利用自适应控制器实现反馈控制,以消除由输入扰动和参数变化而引起的机器人机械手运动误差。这种控制方法在时间上是有效的,且很适合于定点实现。控制方法的有效性通过四自由度的直接驱动机器人前两个关节的计算机仿真实验得到验证。

膨胀土化学改性实验研究

Expansive soil is a kind of typical unsaturated soil with characteristics of high swelling-shrinking deformation, cracks and over consolidation. It is very harmful to civil engineering, As a new processing method deal with expansive soil, Chemistry treatment has widespread applied in developed countries such as Europe and America, and also gained remarkable result. Based on the embankment filling soil improving testing projects in Meng-Xin freeway, this paper proposed a new processing method to expansive soil embankment wrapped with PAS-treated soil, experimental study of expansive soil chemical improved by PAS is been carried out. The water content change is the external factor which causes expansive soil to have swelling-shrinkage deformation. this reflected that the soil body swelling-shrinkage characteristic mainly depends on its mineral ingredient and the soil-water mutual function. This paper takes expansive soil as one kind of ordinary high plastic clay from angle of clay-water mutual function explained the expansive soil swelling-shrinkage deformation mechanism on microscopic. And take this swelling-shrinkage mechanism as the master line, Cooperates with the China Academy of Chemistry, we developed the new method PAS treatment, trough ionic exchange, joint, package and flocculation, the stronger static electricity function weakened the level through adsorption and the stronger static electricity function, PAS can weakened the negative charge repulsion between levels, causes the electric potential to reduce, diffusion layer thickness to be thinner, and improves the water affinity performance of expansive soil effectively. Moreover the space network architecture compromised with PAS and soil enhanced the joint strength between the clay particles , enable the soil body to have comparatively high strength and the distortion rate. pointed proposed the PAS modified principle. Combine with the construction of experimented road, this paper sums up and presents the construction craft and technology requirement of PAS treatment to expansive soil embankment. Through many experimental studied the basic physical property, the intensity characteristic and water stability changes of expansive soil and PAS-treated soil. The results of study indicate that adding lime into the expansive soil can reduce the content of clay gain obviously, reduce the plasticity notably, increase the strength greatly, control the property of swelling and shrinking effectively, and can meliorate the stability of sucking water clearly. Simultaneity PAS don’t change the cultivate capacity of the soil, the modified slope of the embankment can adopt plant fixed slope method as ecology protection. Finally the processing effect of use different treatment has analyzed through numerical simulation, summarized the PAS chemical wrapping treatment process in the actual project application, and appraised its processing effect and the project efficiency. The research indicated that PAS chemical treatment is one effective method to improve expansive soil. Compare with long-distance replacement, especially in the high plastic expansive soil massive distribution area, PAS treatment has the very greatly economical superiority to be promoted. The study in the paper not only afforded technique method to Meng-Xin expressway construction but also important for improvement of the expressway construction theory in swelling soil areas. Key words: PAS; expansive soil; swelling-shrinkage deformation mechanism; wrapping embankment; chemical modified treatment.

基于波动方程的层间多次波预测方法研究

In exploration geophysics，velocity analysis and migration methods except reverse time migration are based on ray theory or one-way wave-equation. So multiples are regarded as noise and required to be attenuated. It is very important to attenuate multiples for structure imaging, amplitude preserving migration. So it is an interesting research in theory and application about how to predict and attenuate internal multiples effectively. There are two methods based on wave-equation to predict internal multiples for pre-stack data. One is common focus point method. Another is inverse scattering series method. After comparison of the two methods, we found that there are four problems in common focus point method: 1. dependence of velocity model; 2. only internal multiples related to a layer can be predicted every time; 3. computing procedure is complex; 4. it is difficult to apply it in complex media. In order to overcome these problems, we adopt inverse scattering series method. However, inverse scattering series method also has some problems: 1. computing cost is high; 2. it is difficult to predict internal multiples in the far offset; 3. it is not able to predict internal multiples in complex media. Among those problems, high computing cost is the biggest barrier in field seismic processing. So I present 1D and 1.5D improved algorithms for reducing computing time. In addition, I proposed a new algorithm to solve the problem which exists in subtraction, especially for surface related to multiples. The creative results of my research are following: 1. derived an improved inverse scattering series prediction algorithm for 1D. The algorithm has very high computing efficiency. It is faster than old algorithm about twelve times in theory and faster about eighty times for lower spatial complexity in practice; 2. derived an improved inverse scattering series prediction algorithm for 1.5D. The new algorithm changes the computing domain from pseudo-depth wavenumber domain to TX domain for predicting multiples. The improved algorithm demonstrated that the approach has some merits such as higher computing efficiency, feasibility to many kinds of geometries, lower predictive noise and independence to wavelet； 3. proposed a new subtraction algorithm. The new subtraction algorithm is not used to overcome nonorthogonality, but utilize the nonorthogonality's distribution in TX domain to estimate the true wavelet with filtering method. The method has excellent effectiveness in model testing. Improved 1D and 1.5D inverse scattering series algorithms can predict internal multiples. After filtering and subtracting among seismic traces in a window time, internal multiples can be attenuated in some degree. The proposed 1D and 1.5D algorithms have demonstrated that they are effective to the numerical and field data. In addition, the new subtraction algorithm is effective to the complex theoretic models.

各向同性与各向异性介质中地震波旅行时正演模拟与反演成像

By seismic tomography, interesting results have been achieved not only in the research of the geosphere with a large scale but also in the exploration of resources and projects with a small scale since 80'. Compared with traditional inversion methods, seismic tomography can offer more and detailed information about subsurface and has been being paid attention by more and more geophysicists. Since inversion based on forward modeling, we have studied and improved the methods to calculate seismic traveltimes and raypaths in isotropic and anisotropic media, and applied the improved forward methods to traveltime tomography. There are three main kinds of methods to calculate seismic traveltime field and its ray path distribution, which are ray-tracing theory, eikonal equation by the finite-difference and minimum traveltime tree algorithm. In ray tracing, five methods are introduced in the paper, including analytic ray tracing, ray shooting, ray bending, grid ray tracing and rectangle grid ray perturbation with three points. Finite-difference solution of eikonal equation is very efficient in calculation of seismic first-break, but is awkward in calculation of reflection traveltimes. We have put forward a idea to calculate traveltimes of reflected waves using a combining way of eikonal equation method and other one in order to improve its capability of dealing with reflection waves. The minimum traveltime tree algorithm has been studied with emphases. Three improved algorithms are put forward on the basis of basic algorithm of the minimum traveltime tree. The first improved algorithm is called raypath tracing backward minimum traveltime algorithm, in which not only wavelets from the current source but also wavelets from upper source points are all calculated. The algorithm can obviously improve the speed of calculating traveltimes and raypaths in layered or blocked homogeneous media and keep good accuracy. The second improved algorithm is raypath key point minimum traveltime algorithm in which traveltimes and raypaths are calculated with a view of key points of raypaths (key points of raypths mean the pivotal points which determine raypaths). The raypath key point method is developed on the basis of the first improved algorithm, and has better applicability. For example, it is very efficient even for inhomogeneous media. Another improved algorithm, double grid minimum traveltime tree algorithm, bases upon raypath key point scheme, in which a model is divided with two kinds of grids so that the unnecessary calculation can be left out. Violent undulation of curved interface often results in the phenomenon that there are no reflection points on some parts of interfaces where there should be. One efficacious scheme that curved interfaces are divided into segments, and these segments are treated respectively is presented to solve the problem. In addition, the approximation to interfaces with discrete grids leads to large errors in calculation of traveltimes and raypaths. Noting the point, we have thought a new method to remove the negative effect of mesh and to improve calculation accuracy by correcting the traveltimes with a little of additional calculation, and obtained better results.

Foundations of Actor Semantics

The actor message-passing model of concurrent computation has inspired new ideas in the areas of knowledge-based systems, programming languages and their semantics, and computer systems architecture. The model itself grew out of computer languages such as Planner, Smalltalk, and Simula, and out of the use of continuations to interpret imperative constructs within A-calculus. The mathematical content of the model has been developed by Carl Hewitt, Irene Greif, Henry Baker, and Giuseppe Attardi. This thesis extends and unifies their work through the following observations. The ordering laws postulated by Hewitt and Baker can be proved using a notion of global time. The most general ordering laws are in fact equivalent to an axiom of realizability in global time. Independence results suggest that some notion of global time is essential to any model of concurrent computation. Since nondeterministic concurrency is more fundamental than deterministic sequential computation, there may be no need to take fixed points in the underlying domain of a power domain. Power domains built from incomplete domains can solve the problem of providing a fixed point semantics for a class of nondeterministic programming languages in which a fair merge can be written. The event diagrams of Greif's behavioral semantics, augmented by Baker's pending events, form an incomplete domain. Its power domain is the semantic domain in which programs written in actor-based languages are assigned meanings. This denotational semantics is compatible with behavioral semantics. The locality laws postulated by Hewitt and Baker may be proved for the semantics of an actor-based language. Altering the semantics slightly can falsify the locality laws. The locality laws thus constrain what counts as an actor semantics.

Birkhoff normal forms for Fourier integral operators II

Iantchenko, A.; Sj?strand, J., (2001) 'Birkhoff normal forms for Fourier integral operators II', American Journal of Mathematics 124(4) pp.817-850 RAE2008

Asymptotic behaviour of the elastic solution near the tip of a crack situated at a nonideal interface

Mishuris, G; Kuhn, G., (2001) 'Asymptotic behaviour of the elastic solution near the tip of a crack situated at a nonideal interface', Zeitschrift f?r Angewandte Mathematik und Mechanik 81(12) pp.811-826 RAE2008

Twierdzenia o punktach stałych pewnych klas operatorów nieliniowych wraz z ich zastosowaniami

Wydział Matematyki i Informatyki