Anisotropic mesh adaptation: towards user‐independent,mesh‐independent and solver‐independent CFD. Part I: general principles

The present paper is the lead article in a three‐part series on anisotropic mesh adaptation and its applications to structured and unstructured meshes. A flexible approach is proposed and tested on two‐dimensional, inviscid and viscous, finite volume and finite element flow solvers, over a wide range of speeds. The directional properties of an interpolation‐based error estimate, extracted from the Hessian of the solution, are used to control the size and orientation of mesh edges. The approach is encapsulated into an edge‐based anisotropic mesh optimization methodology (MOM), which uses a judicious sequence of four local operations: refinement, coarsening, edge swapping and point movement, to equi‐distribute the error estimate along all edges, without any recourse to remeshing. The mesh adaptation convergence of the MOM loop is carefully studied for a wide variety of test cases. The mesh optimization generic coupling of MOM with finite volume and finite element flow solvers is shown to yield the same final mesh no matter what the starting point is. It is also shown that on such optimized meshes, the need for computational fluid dynamics (CFD) stabilization artifices, such as upwinding or artificial viscosity, are drastically reduced, if not altogether eliminated, in most well‐posed formulations. These two conclusions can be considered significant steps towards mesh‐independent and solver‐independent CFD. The structure of the three‐part series is thus, 1, general principles; 2, methodology and applications to structured and unstructured grids; 3, applications to three‐dimensional flows. Copyright © 2000 John Wiley & Sons, Ltd.     

计算含动边界非定常流动的无网格算法

在无网格算法中考虑了含动边界的流动问题,研究了可以计算处理包含一定位移及扭转动边界非定常流动的算法. 创建了无网格算法的动点法则,并引入抗扭方法对弹簧方法进行改进来处理离散点运动,提高了方法的可用度及精度. 发展了求解基于无网格的ALE方程组的算法,在点云离散的基础上采用曲面逼近计算空间导数及HLLC格式计算数值通量,运用四步龙格-库塔法进行时间推进. 在跨、超音速条件下,计算模拟了典型翼型简谐振动流场,计算结果与实验结果及文献对比吻合,验证了该算法的正确性.     

弹簧支撑自适应Z-pinch单层柱面丝阵负载研制

在靶室阴、阳电极呈水平放置的电磁内爆加速器中,抽真空时电极间距会收缩变短,电极间的丝阵负载无法维持预定的分布状态,为此,设计出了一套单层柱面状丝阵负载制备技术,通过弹簧预先压缩足够长度使负载能自动适应电极间距变化,并利用一套可以拆卸的辅助装置固定负载各部件。对制备出的多种单层柱面丝阵负载的测试及靶场实验结果表明,负载能自动适应电极间距变化,最大收缩量达3mm,负载丝仍保持预定的分布并且绷紧,制备技术也便于负载的装配、运输和安装。     

带弹簧支撑超静定结构位移的简化计算方法

针对带有弹性支撑的超静定结构的位移计算问题, 对利用虚功原理得到的位移计算一般公式进行了进一步推导, 提出了该类结构位移的简化计算方法. 文中计算公式表达形式简洁, 且与一般刚性支撑超静定结构的位移计算方法一致. 通过与有限元分析结果的对比,证明了该公式的正确性.     

Low band gap frequencies and multiplexing properties in 1D and 2D mass spring structures

This study reports on the propagation of elastic waves in 1D and 2D mass spring structures.An analytical and computation model is presented for the 1D and 2D mass spring systems with different examples.An enhancement in the band gap values was obtained by modeling the structures to obtain low frequency band gaps at small dimensions.Additionally,the evolution of the band gap as a function of mass value is discussed.Special attention is devoted to the local resonance property in frequency ranges within the gaps in the band structure for the corresponding infinite periodic lattice in the 1D and 2D mass spring system.A linear defect formed of a row of specific masses produces an elastic waveguide that transmits at the narrow pass band frequency.The frequency of the waveguides can be selected by adjusting the mass and stiffness coefficients of the materials constituting the waveguide.Moreover,we pay more attention to analyze the wave multiplexer and DE-multiplexer in the 2D mass spring system.We show that two of these tunable waveguides with alternating materials can be employed to filter and separate specific frequencies from a broad band input signal.The presented simulation data is validated through comparison with the published research,and can be extended in the development of resonators and MEMS verification.     

A planar rod model with flexible cross-section for the folding and the dynamic deployment of tape springs: Improvements and comparisons with experiments

A planar rod model with flexible cross-section has been recently proposed in literature (Guinot et al., 2012). This model is especially suitable for the modeling of tape springs, which develop localized folds due to the flattening of the cross-section. Starting from a complete non-linear elastic shell model, original kinematics assumptions (inspired from the elastica model) have been made to describe the important in-plane changes of the cross-section shape. In the present work, the choice of the position of the rod reference line is discussed. This choice plays an important role in the overall behavior because of the large changes of the cross-section shape. We show that the model published in Guinot et al. (2012) can be improved by considering the centerline as the rod reference line. This enhanced model is then validated through quantitative comparisons with experimental results of dynamic deployments taken from literature.     

屈服后表面裂纹与穿透裂纹的等效裂纹长度换算

在较小范围屈服条件下，将D-M模型和弹塑性线弹簧模型相结合，建立了屈服后一般扁壳表面裂纹的计算模型，导出了其控制方程和相应的COD计算公式。在此基础上，提出了屈服后一般扁壳表面裂纹与平板穿透裂纹的等COD换算方法。通过各种影响因素的分析并去除次要因素，最终给出了该等效换算关系的简单拟合公式。计算结果表明，该分析方法简便且有足够的工程精度，适应于在役压力容器等结构的现场快速断裂评定。     

应力梯度不均匀时平板表面裂纹的应力强度因子

本文采用考虑裂纹面上具有任意分布载荷的线弹簧模型,在Kirchhoff板弯曲理论的假设下,将含半椭圆型表面裂纹的平板问题化为一组耦合的积分方程组进行求解,对均匀拉伸和纯弯曲两种载荷作用下的应力强度因子数值解,同经典线弹簧模型和有限元解进行了比较,并给出了经典线弹簧模型不能得到的、裂纹面上承受幂次不均匀应力分布时应力强度因子的数值解.     

Thermo-elastic properties of heterogeneous materials with imperfect interfaces: Generalized Levin's formula and Hill's connections

We study the thermo-elastic properties of heterogeneous materials containing spherical particles or cylindrical fibres. The interface between the matrix and second-phase inhomogeneity is imperfect with either the displacement or the stress experiencing a jump across it. We relate the effective coefficient of thermal expansion (CTE) to the effective elastic moduli and thereby generalize Levin's formula, and reveal two connections among the effective elastic moduli, thereby generalizing Hill's connections. In contrast to the classical results, the effective CTE in the presence of an imperfect interface is strongly dependent on the size of the inhomogeneity, besides the interface elastic and thermo-elastic properties. This size dependence has been accurately captured by simple scaling laws.     

Use of amidoximated hydrogel for removal and recovery of U(VI) ion from water samples

Poly(acrylamidoxime-co-2-acrylamido-2-methylpropane sulfonic acid) (PAMSA) hydrogel was prepared by copolymerization of acrylonitrile and 2-acrylamido-2-methylpropane sulfonic acid as monomer, N,N′-methylenebis(acrylamide) as crosslinking agent and potassium peroxodisulfate as initiator. Amidoximated copolymer network was prepared by the reaction of copolymer network with hydroxylamine hydrochloride. A batch procedure was used for the determination of the characteristics of the U(VI) solid phase extraction from the amidoximated hydrogel. The determination of U(VI) was performed by spectrophotometric method using arsenazo-III as complexing agent. Optimal pH value for the quantitative preconcentration was 3, and full desorption was achieved with 3 mol L⁻¹ HClO₄. The adsorption process can be well described by the pseudo-second-order kinetic model, and the equilibrium adsorption isotherm was closely fitted with the Langmuir model. A preconcentration factor of 20 and the three sigma detection limit of 2.8 μg L⁻¹ (n = 20) were achieved for uranium(VI) ions. The PAMSA hydrogel was used for separating and preconcentrating the uranyl ion existing in sea water samples, thermal spring water samples and the certified reference materials (TMDA 64; fortified lake water sample).