同位网格有限体积法在变密度浆液流场计算中的应用

为研究自膨胀浆液的扩散机理,基于现代计算流体动力学理论,建立了一种变密度浆液二维流场计算方法.该方法采用结构化同位网格剖分计算区域,利用有限体积法离散浆液流动控制方程,借助动量插值技术构建压力修正方程;采用压力耦合方程组的半隐式算法(SIMPLE)迭代求解动量离散方程和压力、速度修正方程.通过变密度浆液在二维矩形狭槽中自由膨胀算例对该算法进行检验.结果表明,该方法求解正确,并具有较高的计算精度,为建立膨胀性高聚物注浆材料在二维裂隙中流动扩散仿真分析方法奠定了基础.     

用HH-SIMPLE方法求解明渠中三维浅水底孔泄流问题

本文针对明渠中三维底孔喷射流问题,提出了一个新的求解 N-S 方程的有限差分方法——HH-SIMPLE 方法.用该方法求解三维水流问题,可以保持 SIMPLE 方法的优点,即可将动量方程中的速度与压力解锅;求解离散化动量方程时可用隐式迭代格式,而无需作“刚盖近似”或“静水压力假设”,因而可很好地处理自由表面.本文中用该法求解的算例在国际计算流体力学领域中目前尚为空白,其三维计算结果可精确地描述在上游有来流情况下的浅水底孔喷流的流场、自由水面和压力场问题.计算结果和试验资料吻合较好.     

摄动迭代及其在 Burgers 方程数值解中的应用

本文推广了文[1]中的摄动迭代格式,并将其应用于二维耦合 Burgers 方程组的数值求解。结果表明,这种计算格式对求解一类隐式非线性差分方程十分有效.     

风生环流的准三维数值模拟

本文在水平流速垂向分布表达式采用二次多项式的基础上,对二维浅水方程组中的动量方程进行修正,建立了风生环流的准三维数学模型.应用破开算子法求解考虑动量修正的深度平均方程,对破开后的两个分步分别采用隐式格式和预估-校正格式离散,得到了准三维数值解.验证计算表明该数值模式能用于模拟河口,水库和港湾内的风生环流,算法的稳定性好,收敛快.     

求解渗流自由面的有限元混合不动点法

基于有自由面渗流分析的高斯点,建立了求解渗流问题的非光滑非线性方程组模型和求解此类问题的有限元混合不动点算法,此方法属于固定网格法,只需划分一次网格,不需要对数据做任何近似处理,完全利用程序迭代计算渗流自由面.讨论了非光滑方程组解的存在性和该不动点算法的收敛性,通过节点压强插值绘制出渗流自由面.算例结果表明,该方法简单且收敛速度快.对不动点法的收敛性分析为迭代法的收敛提供了理论依据.     

模转换问题的二维标量有限差分波束传输法

为了使用标量有限差分波束传输法求解耦合的波动方程,使用两套网格分别代表两种耦合场变量(TE模或TM模),建立无条件稳定的隐式差分方程组,精确求得完全模转换的转换长度,针对二维LiNbO3TETM转换器给出算例。     

齿顶及啮合间隙对双螺杆压缩机流场特性影响

为了研究转子啮合间隙和齿顶间隙对双螺杆压缩机流场特性的影响,应用计算流体力学技术对压缩机内流场进行数值求解。通过有限体积法求解压缩机流场模型,采用RNG/K-ε湍流模型,流场迭代方法选用压力隐式算子分割法,耗散率方程和湍动能方程采用一阶迎风格式,能量和动量方程采用二阶迎风格式,基于改进的弹簧近似方法进行网格滑移和重构,对不同啮合间隙和齿顶间隙的二维压缩机内部流场参数进行数值求解。研究结果表明,压缩机啮合间隙、齿顶间隙越小,其出口压力越大,内部温度以及内部气体流速越小。     

2维Ginzburg-Landau方程的分裂LOD高阶紧致格式

采用分裂技巧研究了2维的Ginzburg-Landau方程构造高效的数值格式.把2维Ginzburg-Landau方程变成线性和非线性问题以避免求解耦合的非线性方程组.为减少存储量和计算量,对线性问题进一步运用局部1维方法,把它分解为2个1维问题求解.所得到的数值格式具有高效、高精度等数值特征.最后,用数值算例模拟了2维Ginzburg-Landau方程所描述的物理现象,新方法具有较大的优越性.     

剪切变稀自由表面流动问题的不可压SPH模拟

为解决传统弱可压光滑粒子动力学(SPH)方法中,密度的微小变化能引起压力的非物理振荡,发展了一种新的不可压SPH方法,并将其推广到剪切变稀自由表面流动问题的数值模拟中。该方法采用压力Poisson方程求解流体压力,流体控制方程则通过分数步长法进行求解。首先,利用不可压SPH方法模拟了牛顿液滴在水平固壁上的铺展现象,并通过与文献结果的比较,验证了不可压SPH方法模拟自由表面流的有效性;随后,对牛顿液滴和剪切变稀液滴在倾斜固壁上的铺展现象进行了不可压SPH模拟,并讨论了两种液滴铺展现象的不同流动特征。数值结果表明,不可压SPH方法不仅能够有效克服弱可压SPH方法所产生的压力振荡问题,并且能够成功推广到剪切变稀自由表面流的数值模拟中。     

对几种SPH方法的对比研究

用光滑粒子流体动力学(Smoothed Particle Hydrodynamics,简称SPH)方法对函数及其一阶和二阶空间导数核近似进行了详细研究.讨论了几种SPH方法的计算思路,给出了一维、二维和三维情况下分别用CSPH(Corrective SPH)、MSPH(Modified SPH)、SSPH(symmetric SPH)方法,对函数及其一阶和二阶空间导数进行核近似的计算方法,通过一维和二维数值算例,对四种不同的SPH方法进行了比较和误差分析,结果表明MSPH和SSPH方法极大地提高了边界处的精度而且SSPH方法的误差最小.     

Characterization of Ordovician carbonate reservoirs, southeastern Saskatchewan, Canada

The discovery of the prolific Ordovician Red River reservoirs in 1995 in southeastern Saskatchewan was the catalyst for extensive exploration activity which resulted in the discovery of more than 15 new Red River pools. The best yields of Red River production to date have been from dolomite reservoirs. Understanding the processes of dolomitization is, therefore, crucial for the prediction of the connectivity, spatial distribution and heterogeneity of dolomite reservoirs.The Red River reservoirs in the Midale area consist of 3～4 thin dolomitized zones, with a total thickness of about 20 m, which occur at the top of the Yeoman Formation. Two types of replacement dolomite were recognized in the Red River reservoir: dolomitized burrow infills and dolomitized host matrix. The spatial distribution of dolomite suggests that burrowing organisms played an important role in facilitating the fluid flow in the backfilled sediments. This resulted in penecontemporaneous dolomitization of burrow infills by normal seawater. The dolomite in the host matrix is interpreted as having occurred at shallow burial by evaporitic seawater during precipitation of Lake Almar anhydrite that immediately overlies the Yeoman Formation. However, the low δ18O values of dolomited burrow infills (-5.9‰～ -7.8‰, PDB) and matrix dolomites (-6.6‰～ -8.1‰, avg. -7.4‰ PDB) compared to the estimated values for the late Ordovician marine dolomite could be attributed to modification and alteration of dolomite at higher temperatures during deeper burial, which could also be responsible for its 87Sr/86Sr ratios (0.7084～0.7088) that are higher than suggested for the late Ordovician seawaters (0.7078～0.7080). The trace amounts of saddle dolomite cement in the Red River carbonates are probably related to "cannibalization" of earlier replacement dolomite during the chemical compaction.     

An Index mechanism for multi-scale view in spatial databases

There are numerous geometric objects stored in the spatial databases. An importance function in a spatial database is that users can browse the geometric objects as a map efficiently. Thus the spatial database should display the geometric objects users concern about swiftly onto the display window. This process includes two operations:retrieve data from database and then draw them onto screen. Accordingly, to improve the efficiency, we should try to reduce time of both retrieving object and displaying them. The former can be achieved with the aid of spatial index such as R-tree, the latter require to simplify the objects. Simplification means that objects are shown with sufficient but not with unnecessary detail which depend on the scale of browse. So the major problem is how to retrieve data at different detail level efficiently. This paper introduces the implementation of a multi-scale index in the spatial database SISP (Spatial Information Shared Platform) which is generalized from R-tree. The difference between the generalization and the R-tree lies on two facets: One is that every node and geometric object in the generalization is assigned with a importance value which denote the importance of them, and every vertex in the objects are assigned with a importance value,too. The importance value can be use to decide which data should be retrieve from disk in a query. The other difference is that geometric objects in the generalization are divided into one or more sub-blocks, and vertexes are total ordered by their importance value. With the help of the generalized R-tree, one can easily retrieve data at different detail levels.Some experiments are performed on real-life data to evaluate the performance of solutions that separately use normal spatial index and multi-scale spatial index. The results show that the solution using multi-scale index in SISP is satisfying.     

A computer generator for randomly layered structures

AcomputergeneratorforrandomlylayeredstructuresYUJia shun1,2,HEZhen hua2(1.TheInstituteofGeologicalandNuclearSciences,NewZealand;2.StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,China)Abstract:Analgorithmisintrod…     

海南岛地体及其毗邻陆缘晚中生代—新生代古地磁研究和构造演化

本文叙述了对海南岛及其毗邻大陆边缘白垩纪到第四纪地层岩石进行古地磁研究的全部工作过程。通过分析岩石中剩余磁矢量的磁偏角及磁倾角的变化,提出海南岛白垩纪以来经历的构造演化模式如下:早期伴随顺时针旋转而向南迁移,后期伴随逆时针转动并向北运移。联系该地区及邻区的地质、地球物理资料,对海南岛上述的构造地体运动提出以下认识:北部湾内早期有一拉张作用,主要是该作用使湾内地壳显著伸长减薄,形成北部湾盆地。从而导致了海南岛的早期构造运动,而海南岛后期的构造运动则主要是受南海海底扩张的影响。海南地体运动规律的阐明对于了解北部湾油气盆地的形成演化有重要的理论和实际意义。     

Exciton Seebeck in Molecular Systems

Various applications relevant to the exciton dynamics,such as the organic solar cell,the large-area organic light-emitting diodes and the thermoelectricity,are operating under temperature gradient.The potential abnormal behavior of the exicton dynamics driven by the temperature difference may affect the efficiency and performance of the corresponding devices.In the above situations,the exciton dynamics under temperature difference is mixed with     

Achieving Unusual States of Matter Under High Pressure

The elongation method,originally proposed by Imamura was further developed for many years in our group.As a method towards O(N)with high efficiency and high accuracy for any dimensional systems.This treatment designed for one-dimensional(ID)polymers is now available for three-dimensional(3D)systems,but geometry optimization is now possible only for 1D-systems.As an approach toward post-Hartree-Fock,it was also extended to     

《Chinese Science Bulletin》SUBJECT INDEX TO VOLUME 59(2014)

正~~     

An on-line scaling method for improving scalability of a database cluster

The explosive growth of the Internet and database applications has driven database to be more scalable and available, and able to support on-line scaling without interrupting service. To support more client's queries without downtime and degrading the response time, more nodes have to be scaled up while the database is running. This paper presents the overview of scalable and available database that satisfies the above characteristics. And we propose a novel on-line scaling method. Our method improves the existing on-line scaling method for fast response time and higher throughputs. Our proposed method reduces unnecessary network use, i.e. , we decrease the number of data copy by reusing the backup data. Also, our on-line scaling operation can be processed parallel by selecting adequate nodes as new node. Our performance study shows that our method results in significant reduction in data copy time.     

An improved R-tree based on childnode''s probability

R-Tree is a good structure for spatial searching. But in this indexing structure,either the sequence of nodes in the same level or sequence of traveling these nodes when queries are made is random. Since the possibility that the object appears in different MBR which have the same parents node is different, if we make the subnode who has the most possibility be traveled first, the time cost will be decreased in most of the cases. In some case, the possibility of a point belong to a rectangle will shows direct proportion with the size of the rectangle. But this conclusion is based on an assumption that the objects are symmetrically distributing in the area and this assumption is not always coming into existence. Now we found a more direct parameter to scale the possibility and made a little change on the structure of R-tree, to increase the possibility of founding the satisfying answer in the front sub trees. We names this structure probability based arranged R-tree (PBAR-tree).     

Integrating GIS Web services based on mediating architecture

The geographic information service is enabled by the advancements in general Web service technology and the focused efforts of the OGC in defining XML-based Web GIS service. Based on these models, this paper addresses the issue of services chaining,the process of combining or pipelining results from several interoperable GIS Web Services to create a customized solution. This paper presents a mediated chaining architecture in which a specific service takes responsibility for performing the process that describes a service chain. We designed the Spatial Information Process Language (SIPL) for dynamic modeling and describing the service chain, also a prototype of the Spatial Information Process Execution Engine (SIPEE) is implemented for executing processes written in SIPL. Discussion of measures to improve the functionality and performance of such system will be included.