Sieve mesh size and quantitative chironomid paleoclimatology

The minimum sieve mesh size for fossil chironomid analysis is usually set at 100 μm, to ensure adequate recovery of small species. Yet taking into account the labor intensity of sorting and identifying fossil chironomid remains, the large numbers of samples requiring processing in paleoclimate studies with high temporal resolution, and the increasing need to engage non-specialist analysts in this work, it seems appealing to sieve samples through a larger mesh size that would mainly retain easier-to-see, easier-to-extract and easier-to-identify fossil specimens. In this study we evaluated the influence of sieve mesh size on chironomid-based quantitative paleoenvironmental reconstructions in African lakes. We developed two chironomid-salinity calibration data sets based on either the >100 μm or >150 μm size fractions of surface-sediment fossil assemblages, and compared the performance of salinity-inference models derived from them. We find that, despite ∼35% additional fossil loss, restriction to >150 μm data did not appreciably affect individual taxon optima or tolerances, and resulted in only a modest reduction of mean taxon richness per lake. Parameters of statistical model performance were as good, or better than those of models based on the >100 μm data, albeit only after excluding two atypically dilute (<100 μS/cm) lakes from the calibration. Application of inference models based on >150 μm or >100 μm data to a 200-year fossil record from Lake Abiyata, Ethiopia, produced very similar trends, amplitudes and uncertainty ranges of inferred past salinity change. Restriction to >150 μm data reduced the mean fossil yield of core samples by ∼20% on average, i.e. fossil loss was markedly lower than in surface-sediment samples. Using the larger mesh size reduced sample processing time by up to 50%, partly by removing a significant proportion of visually obstructive organic debris, and partly by allowing a greater fraction of chironomid fossils to be identified directly in the sorting tray. The fraction of 1st instar group taxa in surface-sediment samples was reduced from 13% to 3%, increasing the mean taxonomic resolution of fossil assemblages, and thus their ecological specificity.     

三维稳定渗流的 p 型自适应有限元分析

建立了三维稳定渗流的p型自适应有限元分析的误差估计器。充分利用p型有限元前处理少、计算精度高和收敛速度快的优点，分析了p型有限元法求解渗流问题的技术特点，给出了求解三维稳定渗流的自适应升阶策略。通过具体算例，研究了单元阶次及网格尺度对计算结果的影响，验证了三维p型有限元法在求解渗流问题时的可行性。     

Efficient discrete modelling of composite structures for rockfall protection

This paper presents a discrete framework for the modelling of composite structures for rockfall protection. The model is applied to analyse the dynamic response of a cylindrical damping module upon impact of a boulder. The damping module consists of a cylindrical wire mesh, two steel rings, a boundary rope, a geotextile lining and a granular filling material. The chain-link wire mesh, the steel rings and the boundary rope are represented with deformable cylinder elements. The geotextile lining is incorporated into the openings of the wire mesh by using deformable facets. The filling material is represented using spherical particles.     

A novel interface element with asymmetric nodes and its application on concrete-faced rockfill dam

Based on the Goodman element, the Guyan reduction method is introduced to develop the interface element with asymmetric nodes, which accomplishes the coarse-fine mesh transition between soil and structure in finite element models by providing a different number of nodes on the two sliding surfaces. The number of Gauss points is greater than that of the traditional Goodman element, which ensures its accuracy. The developed interface element is employed in the dynamic elasto-plastic analysis of CFRDs. The results indicate that the use of the constructed interface element can significantly decrease the number of elements with little influence on the accuracy.     

An auto-adaptive moving mesh method for the numerical simulation of piping erosion

A numerical approach for simulating piping erosion is proposed based on an auto-adaptive moving mesh. This approach adapts a flownet mesh (FM) to discretize the domain. It reduces the simulation of two-dimensional transport of moving particles to one-dimensional problem by taking advantages of flow tube concept. Furthermore, it employs a stepwise procedure to decouple the solid-liquid interaction equations for modeling piping erosion. Specifically, the FM is updated and evolved automatically in accordance to simulated erosion since the flownet is refined in high-velocity regions, where the piping erosion progresses. Therefore, the FM is auto-adaptive to the piping erosion paths.     

Bathymetry Inversion with the Gravity-Geologic Method: A Study of Long-Wavelength Gravity Modeling Based on Adaptive Mesh

Modeling of long-wavelength gravity anomaly is crucial for bathymetry inversion with a gravity-geologic method. We propose a new method, named as iGGM, to approximate the long-wavelength gravity anomalies by using a finite element method based on an adaptive triangular mesh which is constructed by uneven control points. The mesh size is suitably controlled to ensure that there are several control points in each grid. By using iGGM, the bathymetry in the South China Sea (Test Area #1: 112°E–119°E, 12°N–20°N) and East China Sea (Test Area #2: 125°E–130°E, 25°N–30°N) is estimated. The performance of the method was evaluated by comparing the predictions with Earth topographical database 1 (ETOPO1) model and shipborne depths in the test points. Results show that the depths derived by iGGM have a strong correlation with the shipborne depths. In the test points, the mean values of their differences are smaller than 10 m. The standard deviations of their differences are smaller than 120 m and their correlation is stronger than 0.98. Meanwhile, the results provided by the iGGM model are comparable with that obtained by the ETOPO1 model, e.g., the differences between iGGM and ETOPO1 models in test points for Test Areas 1 and 2 are 116 and 70 m in standard deviation, respectively.     

Coupled computations for an elastic‐perfectly plastic soil using adaptive mesh refinement

An adaptive mesh refinement algorithm has been developed for non‐linear computations in geomechanics, based on a smoothed stress–strain finite element formulation. This uses estimates of error in the incremental shear strain invariant to guide the regeneration of unstructured meshes at regular intervals during loading. Following each mesh‐update, no re‐analysis of previous increments with the new mesh is necessary. Algorithm performance has been investigated by analysing a passive earth pressure problem using a linear elastic‐perfectly plastic Mohr–Coulomb soil model. Perfectly drained behaviour has been considered, as have partially drained situations using hydromechanical coupling, while undrained behaviour has been approximated using time steps close to zero. In all cases, mesh adaptivity has been successful in capturing regions of high strain gradient. The results have been compared with analytical solutions. Accurate computations of limit load and shear band orientation have been obtained for a wide range of material dilation angles. Copyright © 2000 John Wiley & Sons, Ltd.     

金刚石钻头参数与钢筋砼的适应性研究

根据钢筋砼的特性,从理论分析结合实践经验合理选择了金刚石钻头的主要性能参数,ＧＷ２号配方对钢筋砼的适用性较好,胎体性能稳定,并选用了粒度为４０－５０目、浓度为４５％的较高品级金刚石与该胎体相匹配,在生产使用中取得了较好的钻进效果 。     

Unstructured mesh generation and landcover-based resistance for hydrodynamic modeling of urban flooding

Urban flood inundation modeling with a hydrodynamic flow solver is addressed in this paper, focusing on strategies to effectively integrate geospatial data for unstructured mesh generation, building representation and flow resistance parameterization. Data considered include Light Detection and Ranging (LiDAR) terrain height surveys, aerial imagery and vector datasets such as building footprint polygons. First, a unstructured mesh-generation technique we term the building-hole method (BH) is developed whereby building footprint data define interior domain boundaries or mesh holes. A wall boundary condition depicts the impact of buildings on flood hydrodynamics. BH provides an alternative to the more commonly used method of raising terrain heights where buildings coincide with the mesh. We term this the building-block method (BB). Application of BH and BB to a flooding site in Glasgow, Scotland identifies a number of tradeoffs to consider at resolutions ranging from 1 to 5 m. At fine resolution, BH is shown to be similarly accurate but execute faster than BB. And at coarse resolution, BH is shown to preserve the geometry of buildings and maintain better accuracy than BB, but requires a longer run time. Meshes that ignore buildings completely (no-building method or NB) also support surprisingly good flood inundation predictions at coarse resolution compared to BH and BB. NB also supports faster execution times than BH at coarse resolution because the latter uses localized refinements that mandate a greater number of computational cells. However, with mesh refinement, NB converges to a different (and presumably less-accurate) solution compared to BH and BB. Using the same test conditions, Hunter et al. [Hunter NM, Bates PD, Neelz S, Pender G, Villanueva I, Wright NG, Liang D, et al. Benchmarking 2D hydraulic models for urban flood simulations. ICE J Water Manage 2008;161(1):13–30] compared the performance of dynamic-wave and diffusive-wave models and reported that diffusive-wave models under-predicted the longitudinal penetration of the flood zone due to important inertial effects. Here, we find that a relatively coarse-mesh implementation of a dynamic-wave model suffers from the same drawback because of numerical diffusion. This shows that whether diffusion is achieved through the mathematics or numerics, the effect on flood extent is similar. Finally, several methods of distributing resistance parameters (e.g., Manning n) across the Glasgow site were evaluated including methods that utilize aerial imagery-based landcover classification data, MasterMap^® landcover classification data and LiDAR-based feature height data (e.g., height of shrubs or hedges). Results show that landcover data is more important than feature height data in this urban site, that shadows in aerial imagery can cause errors in landcover classification which degrade flood predictions, and that aerial imagery offers a more detailed mapping of trees and bushes than MasterMap^® which can locally impact depth predictions but has little impact on flood extent.     

基于非结构化风格的三维大地电磁矢量有限无模拟

Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.