基于通量非结构网格有限体积法的Level Set方程求解

在非结构网格上,采用基于通量的格心有限体积法求解Level Set方程,并给出了其数值离散格式.对Zalesak圆盘、旋转收缩圆及圆面剪切三个经典界面运动问题进行了数值模拟,其质量盈亏量与收敛阶的计算结果表明,该方法克服了传统数值方法造成的质量不守恒的缺点,能准确有效地求解拓扑结构发生变化的复杂运动界面问题,且计算量较格点格式大为减少.     

粘弹流动的间断有限元模拟

针对传统有限元法求解Oldroyd-B本构方程时需加入稳定化方案的缺点,本文基于非结构网格给出了统一间断有限元求解框架.该框架包含采用IPDG(interior penalty discontinuous Galerkin)求解质量方程和动量方程,与采用RKDG(RungeKutta discontinuous Galerkin)求解本构方程这两个核心.数值结果表明:该方法在求解Oldroyd-B本构方程时无需加入稳定化方案,实施比有限元法简便,且具有较高的计算精度,可有效地模拟包含应力奇异点的复杂粘弹流动问题,进而揭示非牛顿粘弹流动的基本特征.     

非定常不可压Navier-Stokes方程基于Crank-Nicolson格式的两水平变分多尺度方法

不可压缩粘性流是密度不发生变化的流体运动．它们被用来描述许多重要的物理现象,例如：天气、洋流、绕翼型流动和动脉内的血液流动．Navier-Stokes方程是不可压缩粘性流的基本方程．因此,求解Navier-Stokes方程的数值方法在近几十年得到了广泛的关注．本文主要给出非定常不可压Navier-Stokes方程基于Crank-Nicolson格式的两水平变分多尺度方法．该方法分为两步：第一步,在粗网格上求解稳定的非线性Navier-Stokes系统;第二步,在细网格上求解稳定的线性问题去校正粗网格上的解．通过该方法推导的速度的误差估计关于时间是二阶收敛的．数值实验验证了在粗细网格匹配合理的情形下,本文的方法与直接在细网格上使用单网格的变分多尺度方法相比,可以节约大量的计算时间．     

饱和两相介质近场波动问题的一种时域全显式数值计算方法

基于 u -p形式的饱和两相介质弹性波动方程,开展了饱和两相介质近场波动问题时域显式数值计算方法的研究。通过对波动方程中的质量矩阵和孔隙流体压缩矩阵进行对角化处理,消除了方程中的动力耦联,实现了波动方程的解耦。分别应用中心差分法和Newmark常平均加速度法求解固相位移和速度,基于向后差分法求解孔隙流体压力,推导得到了饱和两相介质动力响应的时域显式逐步积分的计算列式,建立了饱和两相介质近场波动问题的一种新的时域全显式数值计算方法。进行了该文方法中矩阵对角化合理性的验证。将该方法的数值解与相应的解析解进行对比,二者符合良好,验证了该方法的正确性。将该文建立的时域数值计算方法与透射人工边界方法相结合,应用于饱和两相介质的近场波动问题,进行了饱和土场地地震响应的计算研究,计算结果符合弹性波动理论的基本规律,表明该方法对于饱和两相介质近场波动问题时域计算求解的适用性。基于该方法中时域递推计算格式的传递矩阵,进行了该方法稳定性特性的研究。该文建立的数值计算方法具有时域全显式算法的基本特征。方法中对动力响应的全部分量均采用递推和迭代的模式进行求解,避免了求解耦联的动力方程组。该方法具有较高的计算效率,是进行饱和两相介质近场波动问题时域计算求解的一种有效的算法。     

饱和多孔介质不同动力耦合形式数值分析EI北大核心CSCD

Biot饱和多孔介质波动行为的数值模拟在众多工程领域中具有重要的意义和作用,由于固相与液相耦合方程难以解耦,使该问题的数值模拟难度较大。针对饱和多孔介质中部分耦合u-p及全耦合u-p-U方程形式的特征,推导了相应动力耦合控制方程的有限元弱形式,并引入不同耦合形式的饱和多孔介质时域黏性边界,综合利用Comsol Multiphysics提供的偏微分方程应用模式进行二次开发求解,通过一维饱和多孔介质动力响应的解析解和数值解验证了模型求解技术的合理性和可行性,基于u-p-U耦合形式探讨了冲击荷载作用下干砂饱和砂地基动力固结中应力波传播特性。计算结果表明慢纵波对动力固结的影响比较显著,合理的冲击荷载持续时间有利于固结效果的改善。     

边界层流场时空模式动力稳定性的数值计算及相互转换

数值模拟了边界层流场的时间稳定性和空间稳定性问题,分析了时空两种模式之间的转换关系.在数值模拟中,采用变步长的Runge-Kutta方法,求解基于Chebyshev插值点的边界层Blasius方程,然后对边界条件进行坐标变换,使Chebyshev谱方法适用于边界层稳定性方程的求解.所得结果与前人的结论吻合较好,文中所用的方法不仅减少了计算量,而且能扩展到其它类型的流动,为经济有效的数值求解边界层稳定性问题提供了一种途径.     

多体非连续变形系统的计算力学模型及其应用

针对由不同特性物体所组成的多体系统,探讨了能够涵盖各种变形状态和运动形式的广义有限单元模式及其插值函数形式。对于多体接触问题,发展了能够合理描述界面特性的接触力元模型,即采用某种应力插值函数将界面上的相互作用力由接触对上的接触应力来表达,并将接触对上的接触应力当作需满足界面上屈服准则与流动法则等状态控制条件的参变量,将其作为约束条件加入系统控制方程。根据非连续变形系统的分区参变量最小势能变分原理,联立变分驻值条件与参变量的状态控制条件建立了多体系统非连续变形计算力学分析的基本控制方程,将问题最终归结为一个含有自由变量和等式约束条件的线性互补问题,对此发展了数值解法,并进行了多个算例的数值分析。计算结果表明该模型不仅能够对多体系统进行静、动力耦合分析,而且还能够模拟多体系统的变形与应力及接触界面上的接触应力和相对运动等复杂的非线性过程。     

基于VOF/PLIC界面追踪方法的RTM工艺充模过程数值模拟

针对基于Darcy定律的树脂传递模塑(RTM)工艺的充模过程数值模拟的局限性,将纤维预制体内的充填流动作为两相流(树脂相和空气相)处理,在动量方程中考虑了惯性项和粘性项,采用有限体积方法(FVM)离散控制方程,并与VOF/PLIC界面追踪方法相结合,发展了求解树脂在纤维预制体内非稳态流动问题的数值模拟方法.在此基础上开发了RTM工艺的充模过程数值模拟程序,其算例的数值模拟结果与解析解或实验结果吻合良好,验证了此数值模拟方法的有效性和可靠性.     

饱和多孔介质不同动力耦合形式数值分析

Biot饱和多孔介质波动行为的数值模拟在众多工程领域中具有重要的意义和作用,由于固相与液相耦合方程难以解耦,使该问题的数值模拟难度较大。针对饱和多孔介质中部分耦合u-p及全耦合u-p-U方程形式的特征,推导了相应动力耦合控制方程的有限元弱形式,并引入不同耦合形式的饱和多孔介质时域黏性边界,综合利用Comsol Multiphysics提供的偏微分方程应用模式进行二次开发求解,通过一维饱和多孔介质动力响应的解析解和数值解验证了模型求解技术的合理性和可行性,基于u-p-U耦合形式探讨了冲击荷载作用下干砂饱和砂地基动力固结中应力波传播特性。计算结果表明慢纵波对动力固结的影响比较显著,合理的冲击荷载持续时间有利于固结效果的改善。     

运动界面追踪的CVOFLS方法

针对VOF方法中界面法向量计算不准确和level set方法中界面质量守恒性差的缺陷,本文提出了一种新的界面捕捉方法—CVOFLS(coupled volume of fluid and level set)方法.该方法吸收了VOF和level set两种方法的优点,在每个时间步内,首先求解level set输运方程,根据level set符号距离函数计算界面法向量;然后求解VOF输运方程,根据流体体积分数重构界面并校正流体质量.Zalesak圆盘旋转及圆面剪切的数值模拟表明,该方法能准确捕捉复杂的界面演化过程,具有良好的质量守恒性,并可以提高计算效率.     

Insertion of triangulated surfaces into a meccano tetrahedral discretization by means of mesh refinement and optimization procedures

In this paper, we present a new method for inserting several triangulated surfaces into an existing tetrahedral mesh generated by the meccano method. The result is a conformal mesh where each inserted surface is approximated by a set of faces of the final tetrahedral mesh. First, the tetrahedral mesh is refined around the inserted surfaces to capture their geometric features. Second, each immersed surface is approximated by a set of faces from the tetrahedral mesh. Third, following a novel approach, the nodes of the approximated surfaces are mapped to the corresponding immersed surface. Fourth, we untangle and smooth the mesh by optimizing a regularized shape distortion measure for tetrahedral elements in which we move all the nodes of the mesh, restricting the movement of the edge and surface nodes along the corresponding entity they belong to. The refining process allows approximating the immersed surface for any initial meccano tetrahedral mesh. Moreover, the proposed projection method avoids computational expensive geometric projections. Finally, the applied simultaneous untangling and smoothing process delivers a high‐quality mesh and ensures that the immersed surfaces are interpolated. Several examples are presented to assess the properties of the proposed method.     

Certification Institutions Expected to Update and Transfer Information on Present Standards

Standards are the basis for production enterprises to organize production, ex-factory inspection, trade (delivery) and technical exchanges, product certification, quality arbitration and supervision.……     

The Research of Anti-Wear Properties of Organomolybdenum Compounds in the Engine Oils

A four-ball tester was used to evaluate the anti-wear performance of three kinds of organomolybdemun compounds in the engine oils, i. e., molybdenum dialkyldithiophosphate (MoDDP), molybdenum dialkyldithiocarbamate ( MoDTC), and sulphur and phosphorus freeorganomolybdeum (Molybdate). The results indicate that a low concentration of MoDDP doesn' t improve the anti-wear properties of the commercial engine oils, but a high concentration of MoDDP can obviously improve the anti-wear properties and the load-carrying capacity of the engine oils. MoDTC doesn' t improve the antiwear properties of the engine oils, but worsens the anti-wear properties of the oils. Signifi can timprove ment of frictional and wear characteristics is obtained with Molybdate added in the commercial engine oils and the formulated oils.     

Flow calorimeter for enthalpy increment measurements on condensed gases

A flow calorimeter for enthalpy increment measurements on condensed gases is presented. A better knowledge of the properties of the liquefied natural gas is needed, and therefore a liquid loop has been designed for our flow calorimeter. The fluid loop in the calorimeter is designed in order to avoid the two-phase region, since two phases would give compositional disturbances in the measurements. The avoidance of the two-phase region is made possible by increasing the pressure of the test fluid after the measurement section, then heating the fluid at super-critical pressure past the critical point. Finally, the fluid is throttled to the low-pressure gas state at the inlet condition of the compressor that circulates the fluid. To perform the pressure increase, a new cryogenic pump has been designed. To evaluate the new equipment, measurements were taken on liquid ethane over the temperature range 146–256 K at pressure between 0.9 and 5.1 MPa.     

Refrigerators Enter a New Era of Health Preservation Frestech initiates the establishment of China's first standard for household sterilization refrigerators

On November 30, 2007, the China Association for Standardization (CAS) held a press conference at Beijing Diaoyutai State Guest House. Leaders from the China Household Electric Appliance Research Institute, the China Household Electric Appliance Association, and the China Consumers' Association attended and made speeches.……     

Optimum-Simulation von Ernteprozessen

Zusammenfassung Bei der Gestaltung von Ernteprozessen — beispielsweise des Prozesses Mähdrusch, Körnertransport und -abnahme — als transportverbundene Fließarbeitsverfahren mit mindestens zwei verschiedenen Arbeitsmitteln gibt es einen großen Optimierungsspielraum. Solche Prozesse sind mehrstufige Bedienungsprozesse mit gemischter Anordnung der Bedienungseinrichtungen. Außerdem besitzen sie stark ausgeprägte räumliche Aspekte. Es wird vor allem gezeigt, welche Merkmale solcher Prozesse in Simulationsmodellen abgebildet sein müssen, welche Zielgrößen in Abhängigkeit von welchen Einflußgrößen unter Berücksichtigung welcher Prozeßbedingungen ermittelt werden sollten und wie der Bewertungsprozeß zur Ermittlung der gewünschten Informationen gestaltet werden muß.     

A study of centrifugal cast high boron high speed steel

In the present study a high‐boron high speed steel (HSS) roll material was designed. Many expensive alloy elements have been substituted by cheap boron alloy, and high‐boron high speed steel roll has been manufactured by centrifugal casting method. The microstructures, mechanical properties and wear resistance of centrifugal casting high‐boron high speed steel roll have been investigated by optical microscopy (OM), scanning electron microscopy (SEM), and X‐ray diffraction (XRD) analysis, hardness test, impact test and wear test. The results indicated that the solidification microstructures of high‐boron high speed steel roll consisted of M₂(B,C), (W,Mo)₂(B,C), M₃(B,C), M₂₃(B,C)₆ type borocarbides and martensite, a small amount of retained austenite. Borocarbides were continuously distributed over the grain boundary. After quenching from 1050 °C, local broken network appeared in partial borocarbides, and fine secondary borocarbide precipitated from the matrix. After tempering from 525 °C, the amount of precipitated borocarbide increased significantly. After heat treatment, the hardness of high‐boron high speed steel roll excelled 60 HRC, and its impact toughness excelled 8.0 J/cm². The single groove steel rolling amount of high‐boron high speed steel rolls increases by 500% than that of bainite cast iron roll, when the rolls are used in K₁ mill housing of bar mill.     

Certification Becomes Important for the I Quality; and Safety Supervision of the Chinese Liquor Sector

The end of 2007, over 200 unit products of more than 80 Chinese firms have passed the quality grade certification for liquor products. These products involve distilled spirits, beer, wine, yellow wine, fruit syrup wine and others, and cover over 80% of the national top-branded liquor products.……     

Deposition and Coating Properties on CVD Tungsten

Surface characterization and microstructure studies are performed on chemical vapor deposited (CVD) tungsten coating. There is about 2 μm thickness diffusion layer of tungsten in the molybdenum substrate. The thermal shock test shows tungsten coating has good adhesion with molybdenum substrate, but the elements of oxygen and carbon in the tungsten coating have the bad affection to the adhesion. The result of high-temperature diffusion experiment is the diffusion rate from molybdenum substrate to tungsten coating is faster.