高水头船闸第二分流口自分流体型水力特性研究

高水头船闸输水水流能量巨大,水力指标高,对输水系统布置和型式提出了更高要求。大藤峡水利枢纽工程船闸最大设计水头40.25 m,输水系统第二分流口拟采用空腔自分流型式,在国内外尚无先例。通过采用三维数值模拟的技术手段对该分流口型式及输水系统整体分流性能进行系统研究,详细分析在各典型充、泄水时刻第二分流口部位的水动力特性,论证在工程实际应用中的可行性。研究结果表明,在合理布置各分支廊道的出水支孔尺寸后,该分流口型式在充、泄水过程中分流效果较好,但在分流口内部存在典型流速、压力分布区域,特别是在泄水工况下分流口大墩头附近存在较大范围三角低流速紊动区域,可能会对分流口结构造成不利影响。研究成果可为全面掌握该分流口型式的水力特性以及进一步改进优化提供技术参考。     

入海水道二期对洪泽湖地区洪水风险的影响研究

为探究淮河入海水道二期工程对洪泽湖的泄洪能力以及周边滞洪区的影响,在对研究区河网概化及地形处理的基础上,建立了研究区河网一维、洪泽湖湖区与周边滞洪区二维耦合的水动力数值模型,并采用1991,2003,2006,2007四年实测历史洪水资料对模型进行率定及验证。基于现状工况、规划工况以及不同洪水量级,设计了5个对比方案对研究区进行洪水演进数值模拟,对不同洪水量级进行横向对比,对不同工况进行纵向对比,并对不同方案进行淹没面积分析和洪水影响分析。结果表明,入海水道二期工程的启用可以有效降低洪泽湖水位,提高洪泽湖的防洪标准和周边滞洪区的启用标准,减少同等洪水量级下的受影响人口数和区域GDP,为洪泽湖的防洪安全建设以及周边蓄滞洪区的稳定发展创造了条件。     

高水头船闸一字闸门水动力特性数值模拟

针对某山区高水头船闸30 m淹没水深下一字闸门的水动力特性进行研究,借助FLUENT软件提供的Realizable k-ε紊流模型对门体启闭过程非恒定流场进行动态模拟,通过建立由二维平面混合网格纵向拉伸成的三维棱柱体网格,结合2.5D重构和弹性光顺动网格法,在降低网格数量、提高计算效率的同时,实现一字闸门启闭全过程三维流场模拟。引入VOF法对自由水面迭代求解,得到一字闸门启闭过程门前后水位差及动水阻力矩变化规律,将计算结果与模型试验进行比较,两者吻合较好。模拟计算发现:门底间隙对闸门运行整个阶段动水阻力矩均有影响,动水阻力矩峰值随门底间隙增大而变小。     

基于分形理论的陶瓷材料显微组织模拟

利用VB开发工具,开发出基于二维和三维系统单相、复相陶瓷材料显微组织及其演变的仿真模拟软件,其主要包括三个模块:二维系统仿真模块、三维系统仿真模块和显微结构定量分析模块;可实现二维和三维的显微组织仿真图像输出,并且可方便地观察空间任意角度截面显微组织及其演变;应用分形理论、拓扑学等方法,对显微结构进行了定量分析,得到晶粒数及平均体积等参数与仿真时间MCS的变化关系。     

Particle-in-Cell algorithms for emerging computer architectures

We have designed Particle-in-Cell algorithms for emerging architectures. These algorithms share a common approach, using fine-grained tiles, but different implementations depending on the architecture. On the GPU, there were two different implementations, one with atomic operations and one with no data collisions, using CUDA C and Fortran. Speedups up to about 50 compared to a single core of the Intel i7 processor have been achieved. There was also an implementation for traditional multi-core processors using OpenMP which achieved high parallel efficiency. We believe that this approach should work for other emerging designs such as Intel Phi coprocessor from the Intel MIC architecture.     

CUDA programs for the GPU computing of the Swendsen–Wang multi-cluster spin flip algorithm: 2D and 3D Ising,Potts, and XY models

We present sample CUDA programs for the GPU computing of the Swendsen–Wang multi-cluster spin flip algorithm. We deal with the classical spin models; the Ising model, the q$q$-state Potts model, and the classical XY model. As for the lattice, both the 2D (square) lattice and the 3D (simple cubic) lattice are treated. We already reported the idea of the GPU implementation for 2D models (Komura and Okabe, 2012). We here explain the details of sample programs, and discuss the performance of the present GPU implementation for the 3D Ising and XY models. We also show the calculated results of the moment ratio for these models, and discuss phase transitions.     

Motion-based method for estimating time required to attach self-adhesive insulators

A self-adhesive insulator is a component of a home appliance that is used to suppress vibration or prevent humidity affecting the internal parts of the appliance. There is a wide range of types and designs available, allowing them to be applied to areas having different shapes. At the design stage, once an insulator design has been developed sufficiently to identify its dimensions and features, the attaching time and baseline cost must be estimated with reasonable accuracy to enable a comparison of vendor quotes. However, the current estimation method is not sufficiently accurate in terms of the baseline cost. This paper presents a motion-based time-estimating scheme with which the time required for the attachment of such insulators can be calculated more accurately. The scheme has been developed by analyzing the motions needed to attach 350 insulators and then designating representative motions and their time values. For this purpose, a modular arrangement of predetermined time standards (MODAPTS) is adopted. Motion-based time-estimation method is useful in terms of simplicity and accuracy. It enables design engineers to estimate the time required for the attachment based only on a drawing of the insulator and a few MODAPTS rules. Estimates made with this method should vary from the actual value by no more than 9.5%.     

苯并呋喃酮在抑制聚丙烯降解过程中的作用

采用多次挤出的方法对三种不同结构苯并呋喃酮与受阻酚及亚磷酸酯的三元复配稳定体系在等规聚丙烯 (PP)中的稳定化作用进行了研究。结果表明:含5,7-二叔丁基-3-(3,4-二甲基苯基)-3氢-苯并呋喃-2-酮 (OXBF1)的复配体系3#配方5次挤出之后的熔体流动速率相对于其同分异构体5,7-二叔丁基-3-(2,5-二甲基苯基)-3氢-苯并呋喃-2-酮(PXBF1)、5,7-二叔丁基-3-(2,4-二甲基苯基)-3氢-苯并呋喃-2-酮 (MXBF1)复配体系的5#配方和6#配方分别下降了24．8％和26％;并且,其抑制所稳定PP加工发黄的作用优于其他两种结构,使得3#配方5次挤出之后的黄度指数(YI)相对于5#配方和6#配方分别下降了9．3％和11．9％;此外,亚磷酸酯对氢过氧化物的还原作用同样提高了苯并呋喃酮高温含氧条件下在PP中的热氧稳定化作用。     

GPELab,a Matlab toolbox to solve Gross–Pitaevskii equations I: Computation of stationary solutions

This paper presents GPELab (Gross–Pitaevskii Equation Laboratory), an advanced easy-to-use and flexible Matlab toolbox for numerically simulating many complex physics situations related to Bose–Einstein condensation. The model equation that GPELab solves is the Gross–Pitaevskii equation. The aim of this first part is to present the physical problems and the robust and accurate numerical schemes that are implemented for computing stationary solutions, to show a few computational examples and to explain how the basic GPELab functions work. Problems that can be solved include: 1d, 2d and 3d situations, general potentials, large classes of local and nonlocal nonlinearities, multi-components problems, and fast rotating gases. The toolbox is developed in such a way that other physics applications that require the numerical solution of general Schrödinger-type equations can be considered.