Efficient discretization of pressure‐correction equations on non‐orthogonal grids

An alternative discretization of pressure‐correction equations within pressure‐correction schemes for the solution of the incompressible Navier–Stokes equations is introduced, which improves the convergence and robustness properties of such schemes for non‐orthogonal grids. As against standard approaches, where the non‐orthogonal terms usually are just neglected, the approach allows for a simplification of the pressure‐correction equation to correspond to 5‐point or 7‐point computational molecules in two or three dimensions, respectively, but still incorporates the effects of non‐orthogonality. As a result a wide range (including rather high values) of underrelaxation factors can be used, resulting in an increased overall performance of the underlying pressure‐correction schemes. Within this context, a second issue of the paper is the investigation of the accuracy to which the pressure‐correction equation should be solved in each pressure‐correction iteration. The scheme is investigated for standard test cases and, in order to show its applicability to practical flow problems, for a more complex configuration of a micro heat exchanger. Copyright © 2003 John Wiley & Sons, Ltd.     

Thick grids with circular apertures: A comparison of theoretical and experimental performance

In this paper a comparison between a rigorous electromagnetic model for transmission through a hexagonal array of circular waveguides in a series of thick, metallic screens and experimental measurements in the far infrared is made. It is found that there is excellent agreement between theory and experiment when the frequency is below that where any diffracted orders propagate. The agreement is still very good above this frequency. Below a frequency approximately equal to the cut-off frequency of the circular waveguides little power is transmitted. As the thickness of the screen is increased, this decrease in transmission becomes more abrupt. Also, for thick screens, resonances appear in the transmission spectrum which are analogous to those which appear in the spectra of two grids separated by a distance comparable to the wavelength of the radiation being used.     

非结构子网格压力方法

纯拉氏流体力学计算失败的根源是砂漏运动和虚假涡旋,运用子网格拉氏质量与其相关的子网格密度得到子网格压力,子网格压力产生的扰动隅角力能消除这些虚假运动,以便消除困扰拉氏流体力学算法的网格畸变带来的数值计算困难.将包括子网格压力的相容流体算法推广到非结构网格,研制包括非结构子网格压力的相容流体程序,对Saltzman活塞问题进行数值模拟,从数值结果分析可知,非结构子网格压力方法能抑制虚假涡旋,消除砂漏畸变.     

CGNS API和FVM在非结构混合网格计算中的应用

用CGNS API（CFD General Notation System Application Programming Interface）作为非结构混合网格求解器的前处理和后处理,用FVM（Finite Volume Method）作为偏微分方程求解方法．在前处理过程中,用hash表法对内部网格面和边界网格面进行编号,并计算出相应的几何信息,以满足FVM求解器的需要．从FVM求解器计算出来的各种场信息可以写入原来的CGNS文件,该文件可以被许多专业商业后处理软件（如Tecplot,Fluent,CFX等）读取和进行可视化;对于求解器,用基于网格中心的FVM及SIMPLEC（Semi Implicit Methodfor Pressure Linked Equation Consistent）方法求解压力速度耦合．最后给出两个说明算例．     

Magnetic Coupling Constants of Self‐Assembled CuII [3×3] Grids: Alternative Spin Model from Theoretical Calculations

This paper reports a theoretical analysis of the electronic structure and magnetic properties of a ferromagnetic Cu^II [3×3] grid. A two‐step strategy, combining calculations on the whole grid and on binuclear fragments, has been employed to evaluate all the magnetic interactions in the grid. The calculations confirm an S=7/2 ground state, which is in accordance with the magnetisation versus field curve and the thermal dependence of the magnetic moment data. Only the first‐neighbour coupling terms present non‐negligible amplitudes, all of them in agreement with the structure and arrangement of the Cu 3d magnetic orbitals. The results indicate that the dominant interaction in the system is the antiferromagnetic coupling between the ring and the central Cu sites (J₃=J₄≈?31 cm^?1). In the ring two different interactions can be distinguished, J₁=4.6 cm^?1 and J₂=?0.1 cm^?1, in contrast to the single J model employed in the magnetic data fit. The calculated J values have been used to determine the energy level distribution of the Heisenberg magnetic states. The effective magnetic moment versus temperature plot resulting from this ab initio energy profile is in good agreement with the experimental curve and the fitting obtained with the simplified spin model, despite the differences between these two spin models. This study underlines the role that the theoretical evaluations of the coupling constants can play on the rationalisation of the magnetic properties of these complex polynuclear systems.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

TWO-SCALE FINITE ELEMENT DISCRETIZATIONS FOR PARTIAL DIFFERENTIAL EQUATIONS

Some two-scale finite element discretizations are introduced for a class of linear partialdifferential equations. Both boundary value and eigenvalue problems are studied. Basedon the two-scale error resolution techniques, several two-scale finite element algorithmsare proposed and analyzed. It is shown that this type of two-scale algorithms not onlysignificantly reduces the number of degrees of freedom but also produces very accurateapproximations.     

Unstructured Adaptive Grid Flow Simulations of Inert and Reactive Gas Mixtures

Unstructured adaptive grid flow simulation is applied to the calculation of high-speed compressible flows of inert and reactive gas mixtures. In the present case, the flowfield is simulated using the 2-D Euler equations, which are discretized in a cell-centered finite volume procedure on unstructured triangular meshes. Interface fluxes are calculated by a Liou flux vector splitting scheme which has been adapted to an unstructured grid context by the authors. Physicochemical properties are functions of the local mixture composition, temperature, and pressure, which are computed using the CHEMKIN-II subroutines. Computational results are presented for the case of premixed hydrogen–air supersonic flow over a 2-D wedge. In such a configuration, combustion may be triggered behind the oblique shock wave and transition to an oblique detonation wave is eventually obtained. It is shown that the solution adaptive procedure implemented is able to correctly define the important wave fronts. A parametric analysis of the influence of the adaptation parameters on the computed solution is performed.     

二维矢量边界推进非结构网格生成方法

提出二维矢量边界推进生成非结构三角形网格方法并证明其可行性.根据流场边界尺度布置边界节点并运用矢量边界推进方法生成背景网格,运用符号面积函数和概率筛选方法布置初始点阵,提出Spring-Laplace方法优化节点位置,同时利用边交换技术优化网格结构.该方法可包含任意点源、线源和内嵌边界,可自由进行局部自适应加密或稀疏,实现任意平面域内与尺度要求一致的高效光滑三角网格剖分.