土石坝溃决模拟及水流计算研究进展

对土石坝溃决模拟技术的最新研究成果进行了总结和评价,在此基础上对该研究领域今后的研究工作提出了若干建议,包括:应研究提出不同坝型溃决可能性分析评价方法;针对不同坝型开展溃坝离心模型试验,揭示在不同致灾因子作用下,土石坝溃决机理和溃口发展过程;开展土石坝初始管涌形成以后发展过程的试验研究,揭示孔流转变为堰流的控制条件;针对溃坝水流的流线曲率较大、溃坝水流大多是非恒定超临界流以及筑坝材料粒径级配范围变化大等特点,开展大型溃坝水工水力学模型试验,揭示不同坝型的溃口流量过程、泥沙输移及下游河道洪水演进规律。     

计算机模拟引擎及其在铀成矿过程模拟中的应用研究

本文介绍了计算机模拟引擎的设计与开发 ,铀成矿过程模拟的实施步骤 ,并对其在伊犁盆地南缘的应用作了探讨。     

Catchment scale geostatistical simulation and uncertainty of soil erodibility using sequential Gaussian simulation

Soil erodibility (K) affects sediment delivery to streams and needs to be appropriately quantified and interpolated as a fundamental geographic variable for implementing suitable catchment management and conservation practices. The spatial distribution of K for erosion modelling at non-sampling grid locations has traditionally been estimated using interpolation algorithms such as kriging which do not adequately represent the uncertainty of estimates. These methods cause smoothing effects through overestimating the low values and underestimating the large values. In this study observed values were used to implement a sequential Gaussian simulation (SGS) procedure to evaluate the certainty of modelled data. Soil erodibility values were computed using 41 soil samples taken from the top 10 cm soil layer regularly distributed across four catchments, 367–770 ha in area, within Kangaroo River State forest, New South Wales (NSW). One hundred realisations were applied in the simulation process to provide spatial uncertainty and error estimates of soil erodibility. The results indicated that values simulated by the SGS algorithm produced similar K values for the neighbouring cells. At the pixel level, the SGS approach generated a reliable estimation of soil erodibility in most areas. Spatial variation of the K factor in this study was strongly related to soil landscape differences across the catchments; within catchments slope gradient did not have a substantial impact on the numerical values of the K factor using pixel-by-pixel comparisons of raster grid maps.     

淹没水跃的数值模拟

利用标准双方程紊流模型数值模拟了进口Fr₁为819、淹没度为024、闸下出流形成的淹没水跃。详细介绍了数值离散格式特别是用VOF方法跟踪非恒定自由表面的数值方法。将淹没水跃的数值模拟结果包括自由表面位置、平均流速、雷诺应力分量的分布与物模试验结果进行了对比,表明该数学模型及数值方法能够较为准确地模拟淹没水跃的大尺度紊流结构及其微观结构。     

Computer simulation of crystal structures

It is impossible to calculate from first principles the details of crystal structures. However, it is possible to rely on information which is certain, namely on interatomic distances, and calculate the details of a crystal structure based on these. Individual distances between atoms can be predicted accurately if one considers the coordination numbers of anions and cations, the extended electrostatic valence rule, and the effects of shared edges between different coordination polyhedra. The predicted interatomic distances are used as observations in a distance least squares refinement in which the positional parameters and the cell edges are adjusted until the calculated interatomic distances correspond as closely as possible to the predicted distances. The topology of the simulated structure has to be known or assumed. Appropriate weighting can facilitate proper modelling. Applications include: a) solution of pseudosymmetric structures; b) preliminary refinement of trial structures; c) geometric refinement without X-ray intensities; d) comparison of hypothetical structures with observed polymorphs of the same compositions; e) simulation of temperature dependence of structures; f) simulation of pressure dependence; g) calculation of structures isomorphous to known structures; h) calculation of thermal ellipsoids; i) calculation of local environments deviating from the overall symmetry of a structure; j) testing of hypotheses about the behavior of structures at varying conditions.     

Spatial simulation of geologic variables

Analytical formulation of a general simulation model for generation of anisotropic as well as isotropic synthetic patterns in one, two, or three dimensions is proposed. It has significance for the purpose of modeling geologic properties such as ore grades, reservoir porosity, mineral distribution, fracture spacings, aperture, orientations, etc. General procedures for such a simulation by the autoregressive process are given for model parameters estimation and synthetic pattern generation. The model works on the square net basis and generates sequential patterns first along any desired direction for unidimensional simulation and then two-dimensional patterns are constructed with reference to two orthogonal unidimensional sequences. Applications to synthetic two-dimensional pattern are shown for isotropic cases with different model parameters. The extension of the model to three-dimensional space is readily available.     

采空场覆岩变形数值模拟与相似模拟比较研究

用有限元数值计算尝试了离层模拟,并分别采用数值模拟与相似模拟的方法研究了木架山矿区143剖面倾斜采空场处理时的地表移动与覆岩破坏规律。研究发现,随着矿柱的崩落,覆岩跨度不断增加,引起的地表移动范围也不断增大,地表垂直沉降范围大致为跨度的1.5～1.7倍,垂直沉降曲线始终对称于最大沉降点;水平移动的范围大于垂直沉降的范围,最大水平移动量大约为最大沉降量的40 %,且水平移动不对称。另外,当覆岩发生整体下沉时离层主要存在于跨度两端,离层由层理剪切破坏引起。     

Numerical simulation of pile installation

During installation of a displacement pile, the soil around the pile is heavily distorted. The resulting changes in soil density and stress state around the pile determine the ultimate pile capacity. In most finite element models, the installation phase is not explicitly modelled. In this paper the full installation phase is modelled in two different ways, in a numerical code capable of large deformations.     

Two-dimensional simulation by turning bands

Journel (1974) developed the turning-bands method which allows a three-dimensional data set with specified covariance to be obtained by the simulation of several one-dimensional realizations which have an intermediate covariance. The relationship between the threedimensional and one-dimensional covariance is straightforward and allows the one-dimensional covariance to be obtained immediately. In theory a dense uniform distribution of lines in three-dimensional space is required along which the one-dimensional realizations are generated; in practice most workers have been content to use the fifteen axes of the regular icosahedron. Many mining problems may be treated in two dimensions, and in this paper a turning-bands approach is developed to generate two-dimensional data sets with a specified covariance. By working in two dimensions, the area on which the data is simulated may be divided as finely as desired by the lines on which the one-dimensional realizations are first generated. The relationship between the two-dimensional and one-dimensional covariance is derived as a nontrivial integral equation. This is solved analytically for the onedimensional covariance. The method is applied to the generation of a two-dimensional data set with spherical covariance.     

Atomistic simulation of oxide surfaces

This paper is concerned with the atomistic simulation of the surfaces of MgO, CaO, SrO, Li₂O, c-ZrO₂, SrTiO₃, -Al₂O₃ and -Fe₂O₃. Details are reported of the structure and energy of pure, non-defective surfaces and the constant electrostatic potential in the bulk that results from surface relaxation and electronic polarisation. Cation impurity segregation is discussed in some detail with an emphasis on calculated heats of segregation, equilibrium surface coverages, the lattice structure of segregated surfaces and the energy of impure surfaces.Paper from Conference on Quantum Theory and Experiment, July 1986     

区域地下水动态数值模拟

在系统综述区域地下水数值模拟的基础上,运用有限差分和有限体积法,结合济南泉域裂隙岩溶地下水含水介质系统,采用二维非稳定流数学模型进行了数值模拟,对泉域地下水动态变化进行了评价,为泉域地下水资源评价与管理提供了一种有效的分析手段,为区域地下水资源的科学管理提供了有意义的参考。     

Hierarchical simulation of aquifer heterogeneity: implications of different simulation settings on solute-transport modeling

The fine-scale heterogeneity of porous media affects the large-scale transport of solutes and contaminants in groundwater and it can be reproduced by means of several geostatistical simulation tools. However, including the available geological information in these tools is often cumbersome. A hierarchical simulation procedure based on a binary tree is proposed and tested on two real-world blocks of alluvial sediments, of a few cubic meters volume, that represent small-scale aquifer analogs. The procedure is implemented using the sequential indicator simulation, but it is so general that it can be adapted to various geostatistical simulation tools, improving their capability to incorporate geological information, i.e., the sedimentological and architectural characterization of heterogeneity. When compared with a standard sequential indicator approach on bi-dimensional simulations, in terms of proportions and connectivity indicators, the proposed procedure yields reliable results, closer to the reference observations. Different ensembles of three-dimensional simulations based on different hierarchical sequences are used to perform numerical experiments of conservative solute transport and to obtain ensembles of equivalent pore velocity and dispersion coefficient at the scale length of the blocks (meter). Their statistics are used to estimate the impact of the variability of the transport properties of the simulated blocks on contaminant transport modeled on bigger domains (hectometer). This is investigated with a one-dimensional transport modeling based on the Kolmogorov-Dmitriev theory of branching stochastic processes. Applying the proposed approach with diverse binary trees and different simulation settings provides a great flexibility, which is revealed by the differences in the breakthrough curves.     

The simulation of ore-dressing plants

An efficient and automated ore-dressing plant simulator has been developed.In this simulator, stream variables can be used to describe a large number of unique solid particles that have a wide range of properties. Consequently, a large number of simultaneous equations need to be solved.This paper demonstrates how particles can be categorized and the mass balance for complex recycle systems rapidly achieved. The simulator that has been developed uses unit modules that are linked together by an executive computer programme. The executive, by treating each module in turn, obtains a set of simultaneous equations that are solved iteratively. The sequence in which the modules are treated is determined by partitioning and tearing algorithms that select tear streams, thus rendering the system acyclic. The algorithms are designed to choose the tear set (as there are usually several) that will allow convergence in the fewest possible iterations. This choice is based on a criterion that can be extended to apply to all the convergence methods considered. A simultaneous convergence technique has been developed that is applicable specifically to linear ore-dressing plants, but can also be applied to non-linear systems. This technique, which is called the reduced Newton method, is shown to exhibit local convergence under reasonable conditions and is generally superior to direct substitution and to the bounded Wegstein methods.     

波生沿岸流数值模拟

为了更好地研究近岸海域波生沿岸流,建立了基于高阶Boussinesq水波方程的波生沿岸流时域数值模型。控制方程在中等水深范围内具有较好的色散性和变浅作用性能,同时具有二阶完全非线性特征,适合描述近岸区域波浪强非线性运动。通过采用松弛造波方法实现了非线性波浪的无反射入射,采用周期性侧边界条件模拟开敞边界。通过数值试验,讨论了模型中主要参数对数值结果的影响。利用率定后的参数模拟了均匀坡度海岸上产生的沿岸流,通过和实验数据的对比验证了模型的准确性和适用性。利用模型数值模拟了不同波浪入射条件(包括周期、波高和波浪入射角度)对波生沿岸流的影响。     

散粒体系统动力学仿真

针对多面体粒子组成的散粒体系统,基于离散单元法,建立了每个散粒体的邻居目录,在邻居目录内以射线交叉法判断散粒体间的接触,给出了散粒体系统动力学的仿真方法,有效地实现了对散粒体系统接触碰撞动力学的数值仿真。该方法具有重要的理论价值,可推广应用于众多工程领域的散粒体系统堆积、碰撞动力学仿真     

Numerical simulation of percussive drilling

This paper presents a novel dynamical model to analyze the long‐term response of a percussive drilling system. This departs from existing approaches that usually consider a single activation and bit/rock interaction cycle for the analysis of the process performance. The proposed model integrates the axial dynamics of an elastic piston and an elastic drill bit, a motion‐dependent pressure law to drive the piston, and a generalized bit/rock interaction law representative of the dynamic indentation taking place at the bit/rock interface. It applies to down‐the‐hole percussive drilling as well as top‐hole, with minor modifications. The model does not account for the angular motion or the hole cleaning, however. The model is first formulated mathematically; then, a finite‐dimensional approximation is proposed for computations. Numerical analyses of the model response, for a low‐size down‐the‐hole percussive system, follow. The period‐1 stationary response for the reference configuration is studied in detail, and parametric analyses assessing the influence on the rate of penetration of the bit/rock interaction parameters, the feed force, and the percussive activation parameters are conducted. These analyses reveal that the multiscale nature of the process is well captured by the model and recover expected trends for the influence of the parameters. They also suggest that a significant increase of the penetration rate can be achieved by increasing the percussive frequency. Copyright © 2015 John Wiley & Sons, Ltd.     

Computer simulation of reef growth

Light is one of the major controls on reef growth and carbonate production. The growth of present reef builders depends largely upon the amount of light available for photosynthesis. As light decreases with water depth, so does reef growth. The computer model presented extends this principle by combining two functions, one for photosynthesis and the other for the extinction of light in water. The model is used to simulate the growth of Alacran Reef, Mexico, two reefs of the Great Barrier Reef and the reefs of the windward platform of St Croix, US Virgin Islands. The model also gives an accurate simulation of the growth of fore-reef walls in Belize, in agreement with the accretion hypothesis developed for this feature.     

Numerical simulation of ground anchors

A procedure of finite-element modeling and beam-column modeling of ground anchors was proposed in this study to investigate the load transfer mechanism in ground anchors. The procedure included the modeling of soil, grout, and strand tendon and the interface modeling of soil–grout and grout–strand in ground anchors. A series of finite element analyses and beam-column analyses were performed using the proposed models on ground anchors. The numerical predictions were compared with observed measurements in a field load test. The results indicated that the numerical simulation of load transfer mechanism on ground anchors can provide reasonable predictions.