三维圆柱体-弹性板非线性流固耦合数值计算

采用彼得罗夫-盖尔金有限元法离散不可压缩流体,并对大变形弹性结构体使用连续盖尔金有限元法进行离散.使用任意的拉格朗日-欧拉(ALE)方法处理流固网格的大幅变形,并且采用全解耦的隐式分区方法分别求解流固两相,进而将其应用于流固耦合问题的数值计算.根据一系列的网格划分方案,通过改流体和固体网格之间的网格数量系统地进行了网格收敛性测试,同时将开发的求解方法应用于三维标准的圆柱体-柔性板问题,并与文献中标准解进行对比.结果表明:开发的求解器在流固耦合问题中具有较好的准确性和可靠性.     

移动边界下塑料熔体充填模拟的稳定性方法

针对具有移动边界的塑料复合成形过程模拟中因网格变形而导致的数值稳定性问题,提出了塑料熔体充填模拟的稳定性方法．该方法基于任意拉格朗日欧拉法建立移动边界下充填模拟的控制方程,采用四面体单元中引入中心速度稳定的速度压力耦合计算方案,提升了移动网格下塑料熔体充填模拟的稳定性和精度;进一步针对计算过程中网格变形导致四面体单元反转错误的问题,采用引入反转单元罚函数的线弹性网格优化方法,提升了移动边界下网格变形的稳定性;数值稳定性实验表明该方法在网格极端变形和不同单元纵横比下均能稳定求解．最后采用具有柔性移动边界特征的注射压缩-热压复合成形工艺对该方法进行验证,发现移动边界下塑料熔体流动前沿的模拟误差为3.3%,实现了对此类成形工艺的精确模拟．     

翼型颤振分析的时域方法

使用全耦合法预测了翼型NLR7301的颤振边界.在流体域的动网格上,通过求解非定常的雷诺平均N-S方程,计算出翼型表面的非定常气动力.在耦合了翼型的结构动力学方程后,得到翼型表面的变形位移,通过变形网格技术对流体域的网格进行更新.在此基础上利用时间推进的方法模拟了翼型的自激振动过程.     

侧向挤压过程的耦合无网格-有限元法数值模拟

应用刚塑性耦合无网格-有限元法,对平面侧向挤压过程进行数值模拟.采用修正的罚函数法约束体积不可压缩条件,避免体积闭锁现象.采用边界奇异核方法,直接、准确地施加本质边界条件.针对侧向挤压工艺的特点,将无网格节点布置在变形剧烈的区域,充分发挥其处理大变形问题的能力;在变形较小的区域采用有限元单元,利用其较高的计算效率.数值模拟结果表明,该耦合方法对于侧向挤压这类局部成形问题的求解具有良好的效果.     

镦挤过程自适应耦合无网格-有限元法分析

提出自适应耦合无网格-有限元法,并对镦挤成形过程进行数值分析.基于刚塑性流动理论推导轴对称问题的耦合无网格-有限元法求解公式,避免了有限元法因网格畸变而需重新划分网格的缺点,相对纯无网格法明显提高了计算效率,并提出实现自适应耦合的具体算法,包括动态区域划分算法与无网格节点加密算法.数值算例结果表明,自适应耦合方法能够将模拟一次进行,并得到与有限元法接近的结果,在处理局部大变形问题时具有显著效果.     

自由运动船体迎浪状态下上浪问题的数值模拟

基于FLUENT软件进行二次开发,得到了一个用于求解甲板上浪问题的完整的数值模拟方案.该方案融合了动网格功能、数值波浪水池的消波和推板边界造波、船体运动的自动求解、捕捉自由液面大变形的流体体积分数法、高效的网格划分方案和稳定精确的时间步长方案,最终实现了波浪和船体运动的耦合,成功模拟了甲板上浪现象,并将该方案应用到自由运动船体在顶浪状态下的上浪问题研究中,模拟得到了物体在波浪上的运动及甲板上浪的过程,同时还得到了上浪产生的冲击压力.     

自由界面问题的拉格朗日粒子和流体体积耦合算法

针对自由界面问题,构建了拉格朗日粒子和流体体积(volume of fluid, VOF)耦合算法.拉格朗日粒子方法可以准确追踪运动界面,但是一般很难保证流体的质量守恒性.VOF方法可以保证很好的质量守恒性,但是不容易计算界面的几何信息.因此,本文构造了一种耦合算法,吸收两种方法的优点.耦合算法中还引入了四叉树自适应网格技术,可以在大变形区域提高界面的分辨率,并能减少计算量.利用耦合算法模拟经典的Zalesak旋转盘问题和单涡剪切流动,数值结果和文献已有结果吻合较好,验证了耦合算法的稳定性、有效性和准确性.     

移动交界面流固耦合传热的数值稳定性分析

研究了交界面移动情况下流固耦合稳态传热的数值稳定性问题.考虑Dirichlet-Robin组合边界条件,用速度表征交界面的移动情况,流体域和固体域分别采用有限体积法和有限单元法进行离散及数值求解,利用Goudonov-Ryabenkii理论正则模态分析方法重点研究了交界面移动时数值方法的稳定性,最终获得了一条由耦合系数和移动速度组成的最优曲线,并且证明了当耦合系数和移动速度在这条曲线上取值时,离散的求解域能够达到最快的收敛速度及绝对的稳定性特征.为设计人员进行数值仿真时选取合理的参数提供了参考.     

考虑动态流体网格的颗粒-流体耦合算法

传统颗粒-流体耦合算法难以考虑流体动态边界问题,在模拟大变形案例时因流固边界不匹配易造成计算误差,影响结果的准确性.针对该问题,引入流体网格动态更新方法,推导动网格下的达西渗流方程和颗粒-流体相互作用方程,在离散元商业软件PFC~(2D)基础上,开发了考虑流体动网格的颗粒-流体耦合算法.将该算法用于模拟饱和土不排水剪切双轴试验,通过与常体积法的结果对比,验证了该方法的有效性.最后,采用该算法模拟了不同围压下的不排水双轴试验,计算结果的规律与室内试验具有较好的一致性.该算法考虑了流体动网格的问题,在模拟大变形下的三轴压缩试验、一维固结试验等案例时可获得相适应的流固边界,有助于提高模拟精度,因此可为类似研究提供参考.     

WC-CMPS法对溃坝溢流问题的数值模拟

无网格粒子法(拉格朗日法)是计算流体力学领域最新一代的计算方法,由于不需要生成网格,无网格粒子法能够模拟任意形状的界面及其变形和碎裂,尤其适用于计算如溃坝溢流及波浪破碎等具有复杂自由表面运动的流体力学问题。其中,光滑粒子流体动力学方法(SPH)和移动粒子半隐式法(MPS)是无网格计算方法的典型代表。本文采用基于MPS法的修正弱可压缩移动粒子半隐式(Weakly Compressible Corrected Moving Particle Semi-implicit Method)法构建了用于求解自由表面流问题的数值计算模型,并用于模拟干床面与湿床面上的溃坝溢流问题。计算结果与实验结果的对比分析表明,采用基于WC-CMPS无网格法建立的数值计算模型对干、湿床面上溃坝溢流的溢流形态、水柱高度以及波前速度等均给出了良好的计算结果。     

Markedness and its Interpretation in Language Use

Language markedness is a common phenomenon in languages, and is reflected from hearing, vision and sense, i.e. the variation in the three aspects such as phonology, morphology and semantics. This paper focuses on the interpretation of markedness in language use following the three perspectives, i.e. pragmatic interpretation, psychological interpretation and cognitive interpretation, with an aim to define the function of markedness.     

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.     

A computer generator for randomly layered structures

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

低渗透非均质砂岩油藏启动压力梯度研究

理论推导与室内实验相结合,建立了低渗透非均质砂岩油藏启动压力梯度确定方法。首先借助油藏流场与电场相似的原理,推导了非均质砂岩油藏启动压力梯度计算公式。其次基于稳定流实验方法,建立了非均质砂岩油藏启动压力梯度测试方法。结果表明:低渗透非均质砂岩油藏的启动压力梯度确定遵循两个等效原则。平面非均质油藏的启动压力梯度等于各级渗透率段的启动压力梯度关于长度的加权平均;纵向非均质油藏的启动压力梯度等于各渗透率层的启动压力梯度关于渗透率与渗流面积乘积的加权平均。研究成果可用于有效指导低渗透非均质砂岩油藏的合理井距确定,促进该类油藏的高效开发。     

On Tragic Consciousness In The Works By Ernest Hemingway

As an American modern novelist who were famous in the literary world, Hemingway was not a person who always followed the trend but a sharp observer. At the same time, he was a tragedy maestro, he paid great attention on existence, fate and end-result. The dramatis personae's tragedy of his works was an extreme limit by all means tragedy on the meaning of fearless challenge that failed. The beauty of tragedy was not produced on the destruction of life, but now this kind of value was in the impact activity. They performed for the reader about the tragedy on challenging for the limit and the death.     

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

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

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.     

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

正~~     

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     

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