三维复杂域非结构化六面体网格全自动生成方法

使用自适应有限单元法在三维复杂域内求解水工及岩土工程问题,讨论了应用行波法在多介质三维复杂块体内六面体单元网格的自动生成过程。其特点是对不规则区域具有较强的适应性,并考虑了节理单元的自动生成。     

A high order Godunov scheme with constrained transport and adaptive mesh refinement for astrophysical and geophysical MHD

We propose to extend the well-known MUSCL-Hancock scheme for Euler equations to the induction equation modeling the magnetic field evolution in kinematic dynamo problems. The scheme is based on an integral form of the underlying conservation law which, in our formulation, results in a “finite-surface” scheme for the induction equation. This naturally leads to the well-known “constrained transport” method, with additional continuity requirement on the magnetic field representation. The second ingredient in the MUSCL scheme is the predictor step that ensures second order accuracy both in space and time. We explore specific constraints that the mathematical properties of the induction equations place on this predictor step, showing that three possible variants can be considered. We show that the most aggressive formulations reach the same level of accuracy than the other ones, at a lower computational cost. More interestingly, these schemes are compatible with the Adaptive Mesh Refinement (AMR) framework. It has been implemented in the AMR code RAMSES. It offers a novel and efficient implementation of a second order scheme for the induction equation. The scheme is then adaptated to solve for the full MHD equations using the same methodology. Through a series of test problems, we illustrate the performances of this new code using two different MHD Riemann solvers (Lax–Friedrich and Roe) and the need of the Adaptive Mesh Refinement capabilities in some cases. Finally, we show its versatility by applying it to the ABC dynamo problem and to the collapse of a magnetized cloud core.     

一种基于视觉特征的地形模型数据提取与快速显示方法

DEM模型是一种重要的地形模型，由于其数据量极其巨大，如何实现快速显示是一个难题。针对DEM地形模型，提出基于视点与观察分辨率的数据抽取及基于法矢量的模型简化的快速绘制方法算法，最后进行实验验证。结果表明，该方法数据压缩量大，绘制速度快，且逼真度无明显变化。     

专题地图的生产制作

介绍使用计算机制作专题地图的工艺流程,及计算机制图在印刷前的生产制作的工作步骤和方法,指出了在制作中需要注意的事项,对出版高质量的地图在各环节如何加以控制提出了建议。     

球面四元三角网的复杂拓扑关系计算

针对已有球面四元三角网拓扑关系计算方法不能处理带空洞复杂目标的情况,提出了基于公共三角格网来计算球面四元三角网复杂拓扑关系,该方法能准确判断出带空洞复杂目标的拓扑关系,为全球空间分析和应用提供了新的思路。     

论网格与网格地理信息系统

文中论述和分析了网格、空间信息网格和地图上的格网等概念。在分析目前WebGIS的主要关键技术及实际功能和效果的基础上,从软件和数据两个方面分析了目前WebGIS存在的问题及其原因。针对目前WebGIS对于网格技术的不适应性,提出了关于基于多智能代理的网格地理信息系统的思考,论述了代理、移动代理和智能代理及其优势,给出了系统的体系结构,并对系统客户端、服务器及其智能代理技术进行了说明,研究了系统中移动代理的设计、迁移机制、巡回路线规划和通信机制等。最后给出了初步的结论和讨论意见。     

Stream tube integration method for analysing subsurface fluid flow

A stream tube integration method is introduced to solve transient subsurface fluid flow problems. The method combines a geometry-embedded form of Darcy's Law and the notion of location of average. Two types of problems, transient radial flow to a well of finite radius in an areally infinite aquifer and in a double porosity system, are solved by the stream tube integration method and the integral finite difference method. Results of the solutions show that the stream tube integration method, with fixed coarse mesh, are more accurate and better behaved than the integral finite difference method, with fine mesh. The fixed mesh stream tube integration method is readily extended to the moving mesh method. With much coarse mesh, the moving mesh technique can obtain the same accurate results as the fixed mesh stream tube integration method. It is suggested that the stream tube integration method is a viable way to state, solve, interpret and verify numerical solutions. The method provides efficient computation and improved accuracy for analysing subsurface fluid flow. © 1998 John Wiley & Sons, Ltd.     

低速冲击下高强钢绞线网-高性能砂浆加固钢筋混凝土梁试验研究

研究了高强钢绞线网-高性能砂浆加固钢筋混凝土梁的抗低速冲击性能,利用落锤冲击试验法对4根简支钢筋混凝土梁进行对比试验研究,其中1根为普通梁,3根为加固梁,通过分析加固前后梁的破坏形态,结合挠度变化、钢筋的应变以及加速度的时程曲线,对比得到了钢筋混凝土梁加固后的抗冲击性能。研究结果表明：加固后钢筋混凝土梁的抗冲击能力显著提高,钢绞线和砂浆作为外加层,不仅增加了梁的抗弯剪能力和刚度,提高了结构的耗能能力,而且限制了裂缝的发展,使梁的整体性和延性更好;在相同冲击作用下,钢绞线直径的增加能适当提高结构的抗冲击性能。     

Multi‐method virtual electromagnetic experiments for developing suitable monitoring designs: A fictitious CO2 sequestration scenario in Northern Germany

We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO₂ sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO₂ injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO₂ dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO₂, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO₂ plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.     

高密度聚乙烯/纳米辉钼矿复合材料性能研究

以径厚为30~45 nm的辉钼矿作为导电填料,采用熔融共混法制备高密度聚乙烯/纳米辉钼矿复合材料,应用高阻计ZC-36、SEM、XRD、TG等研究复合材料的电性能、晶体结构、界面微观形貌、热性能等。研究结果表明:当纳米辉钼矿含量为25%,复合材料的体积电阻率为5.1×109Ωcm,纳米辉钼矿高定向排列,呈链状集合体;复合材料tonset分解温度提高了35℃;辉钼矿的晶体结构发生变化,(002)晶面2θ右移了1.12°,面间距缩小0.257 nm。