A Newton-like finite element scheme for compressible gas flows

Semi-implicit and Newton-like finite element methods are developed for the stationary compressible Euler equations. The Galerkin discretization of the inviscid fluxes is potentially oscillatory and unstable. To suppress numerical oscillations, the spatial discretization is performed by a high-resolution finite element scheme based on algebraic flux correction. A multidimensional limiter of TVD type is employed. An important goal is the efficient computation of stationary solutions in a wide range of Mach numbers, which is a challenging task due to oscillatory correction factors associated with TVD-type flux limiters. A semi-implicit scheme is derived by a time-lagged linearization of the nonlinear residual, and a Newton-like method is obtained in the limit of infinite CFL numbers. Special emphasis is laid on the numerical treatment of weakly imposed characteristic boundary conditions. Numerical evidence for unconditional stability is presented. It is shown that the proposed approach offers higher accuracy and better convergence behavior than algorithms in which the boundary conditions are implemented in a strong sense.     

二维欧拉程序中一种新的分界面跟踪技术

１．引言 二维流体力学计算方法按其采用Ｌａｇｒａｎｇｅ坐标系还是Ｅｕｌｅｒ坐标系,分为Ｌａｓｒａｎｇｅ方法和 Ｅｕｌｅｒ方法两大类． Ｌａｇｒａｎｇｅ坐标是随材料质团运动的,由于每个有限网格点总是与初始时同一部分流体相联系,所以拉氏程序能够清晰显示解域内部多种物质的分界向和自由界向,而且描述材料分界面的运动是自动完成的．然而在二维计算中,当材料发生大变形时,拉氏网格发生扭曲,甚至使相邻网格彼此相交和重叠,使计算过程无法进行下去．因此,Ｌａｇｒａｎｇｅ方法在二维计算中有其局限性,只能用了变形较小的模…     

共轴圆线圈和圆板的互感系数计算

本文给出了同轴圆线圈和圆板间互感系数的近似函数表达式。并用Matlab绘制了该互感系数随线板间距离和圆板大小的变化曲线。此外,根据磁感应强度的解析表达式,计算了载流线圈的磁力线分布图,分析了其磁感应强度的分布,并给出了它的径向分布曲线。     

Flux limiters fourth and sixth order high resolution method for hyperbolic conservation laws

This paper presents new fourth and six order accurate explicit finite difference scheme for hyperbolic conservation laws. The technique of obtaining high resolution, total variation non increasing oscillation free and explicit scalar difference scheme is derived by adding multi-limiters of antidiffusive flux to a first order scheme.     

Energy flux characteristics of seismic waves at the interface between soil layers with different saturations

Based on the multiphase poroelasticity theory describing the propagation of waves in the unsaturated fluid-saturated porous medium, the reflection and transmission coefficients of the seismic waves at the interface between soil layers with different saturations are obtained. Our unsaturated model consists of a deformable skeleton in which two compressible and viscous fluids (i.e., water and gas) flow in the interstices. Three compressional waves (i.e., P1, P2, and P3 waves) and one shear (i.e., S wave) wave exist in the unsaturated soils. The expressions for the energy ratios of the various reflected and transmitted waves at the interface during the transmission and reflection processes are presented in explicit forms accordingly. At last, numerical computations are performed and the results obtained are respectively depicted graphically. The variation of the energy ratios with the incident angle, wave frequency and saturation degrees of the upper and lower soil layers is illustrated in detail. The calculation results show that the allocation of incident seismic waves at the interface is influenced not only by the angle and frequency of the incident seismic waves, but also by the saturations of the upper and lower soil layers. It is also verified that, at the interface, the sum of energy ratios of the reflected and transmitted waves is approximately equal to unity as was expected. This study is of importance to several fields such as geotechnical engineering, seismology, and geophysics.     

锅炉过热器/再热器管道热流密度的评估方法

通过建立数值模型进行编程计算,分析了锅炉过热器/再热器管道在不同的管道几何尺寸、烟气温度、蒸汽流量和温度下,各个界面热流密度的变化情况。通过回归分析得到了热流密度与管道内壁侧氧化层厚度和运行时间的关系表达式,即热流密度的评估关系式。     

一种矩形微米流道的加工与流量测定

介绍了一种比较成熟的矩形微米流道(微槽道)加工工艺,利用厚胶光刻工艺在硅片上制备胶模板,进而在该模板上用高分子聚合物PDMS(Polydimethylsiloxane,聚二甲基硅氧烷)制作了微米量级的微槽道(高、宽范围20~100μm),对制备的微槽道的几何精确度进行了20×100倍放大显微观测,对其流量特性进行了实验测定.实验结果与泊肃叶理论公式(NS方程的管流表达式)相对比,可以看到在微米量级的矩形微流道中,泊肃叶公式仍然适用.     

低温日用光泽透明釉的研制

研制了一种以硅灰石、锂辉石为主要助熔剂，烧结温度在1150℃的低温日用光泽透明釉，性能达到了国家标准。文章通过对其烧成过程和显微结构观测，探讨了其烧成机理，以及助熔剂在成釉过程中的作用。     

新型横向磁场电机的标量磁位三维磁场分析与参数计算

提出了一种新型横向磁场电机（TFM），定子采用E型铁心，可以大大降低加工难度和制造成本。但是，该电机磁路复杂，需采用三维场场进行分析。在三维有限元磁场中，如果对电流区域进行适当处理，采用标量磁位进行分析，与采用矢量磁位相比，可大大提高计算速度。对电流区域在三维磁场计算中的处理方法进行了推导，根据所推导的原理编制三维有限元软件，利用该软件对样机磁场进行计算。试验结果表明了该方法的正确性。     

基于两相漂移流理论的蒸汽发生器计算程序(CCM)的应用

介绍了一个建立在两相漂移流和相边界可连续移动模型基础上的、可较为全面描述单相及两相流的传热和流动的程序CCM.该程序通过"修正的有限元方法"对控制体节点之间建立联系,保证模拟效果的同时,节省了计算时间.以CCM为基础的UTSG-3程序可较好地模拟U型管蒸汽发生器的稳态和瞬态过程,可得到更为广泛的应用.通过不同运行参数的扰动研究了大亚湾核电站蒸汽发生器动态响应过程,结果显示UTSG-3程序具有模拟U型管蒸汽发生器动态行为的能力.