On the adjacencies of triangular meshes based on skeleton-regular partitions

For any 2D triangulation τ, the 1-skeleton mesh of τ is the wireframe mesh defined by the edges of τ, while that for any 3D triangulation τ, the 1-skeleton and the 2-skeleton meshes, respectively, correspond to the wireframe mesh formed by the edges of τ and the “surface” mesh defined by the triangular faces of τ. A skeleton-regular partition of a triangle or a tetrahedra, is a partition that globally applied over each element of a conforming mesh (where the intersection of adjacent elements is a vertex or a common face, or a common edge) produce both a refined conforming mesh and refined and conforming skeleton meshes. Such a partition divides all the edges (and all the faces) of an individual element in the same number of edges (faces). We prove that sequences of meshes constructed by applying a skeleton-regular partition over each element of the preceding mesh have an associated set of difference equations which relate the number of elements, faces, edges and vertices of the nth and (n−1)th meshes. By using these constitutive difference equations we prove that asymptotically the average number of adjacencies over these meshes (number of triangles by node and number of tetrahedra by vertex) is constant when n goes to infinity. We relate these results with the non-degeneracy properties of longest-edge based partitions in 2D and include empirical results which support the conjecture that analogous results hold in 3D.     

Domain decomposition preconditioners of Neumann-Neumann type for hp-approximations on boundary layer meshes in three dimensions

We develop and analyse Neumann–Neumann methods for hpfinite-element approximations of scalar elliptic problems ongeometrically refined boundary layer meshes in three dimensions.These are meshes that are highly anisotropic where the aspectratio typically grows exponentially with the polynomial degree.The condition number of our preconditioners is shown to be independentof the aspect ratio of the mesh and of potentially large jumpsof the coefficients. In addition, it only grows polylogarithmicallywith the polynomial degree, as in the case of p approximationson shape-regular meshes. This work generalizes our previousone on two-dimensional problems in Toselli & Vasseur (2003a,submitted to Numerische Mathematik, 2003c to appear in Comput.Methods Appl. Mech. Engng.) and the estimates derived here canbe employed to prove condition number bounds for certain typesof FETI methods.     

复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

The influence of patterned nanofiber meshes on human mesenchymal stem cell osteogenesis

A specially designed electroconductive collector enables the electrospinning of P-NFM composed of areas of parallel/uniaxially aligned fibers and areas of random/orthogonal nanofiber distribution. The biological relevance of P-NFM is demonstrated using hBMSCs as an autologous cell source. The structures induce cell orientation along the uniaxially aligned fibers, mainly during earlier culturing periods under basal and osteogenic differentiation conditions. The microtopography of the P-NFM also controls the deposition of mineralized extracellular matrix along the pre-defined fiber direction. Genotypic characterization confirms the successful differentiation into the osteogenic lineage.     

谱体积方法单元分割与问题单元标记方法的改进研究

谱体积方法是一种本质上解决网格依赖性的高精度CFD计算方法,本文研究了二维Euler方程的谱体积方法,提出一种基于切比雪夫多项式的单元分割方法,建立了基于WENO的变量限制器方法,并发展了结合谱体积和控制体的问题单元标记方法.采用15°超声速压缩拐角和NACA0012跨声速流动两个典型算例进行验证,结果表明,该分区方法具有更好的计算精度,标记方法可有效识别不连续区域,在较少的网格下即可获得与密网格传统有限体积法相当的计算精度.     

各向异性网格下抛物方程一个新的非协调混合元收敛性分析

本文将 Crouzeix-Raviart 型非协调线性三角形元应用到抛物方程,建立了一个新的混合元格式.在抛弃传统有限元分析的必要工具 Ritz 投影算子的前提下,直接利用单元的插值性质和导数转移技巧, 分别得到了各向异性剖分下关于原始变量u 的H^-1-模和积分意义下L²-模以及通量p=-▽u 在L²-模下的最优阶误差估计.数值结果与我们的理论分析是相吻合的.     

二维非结构网格有限体积格式的壁热现象和人工热流方法

基于Noh人工热流方法的思想,推导出任意多边形网格人工热流的格式,将Noh的结构网格人工热流方法推广到任意多边形非结构网格的拉氏有限体积方法中,抑制或消除了壁热现象.几个数值实验,包括Noh激波反射问题和两块钢金属碰撞问题,验证了方法的可行性.     

各向异性网格下特征值问题的有限元分析

主要目的是在各向异性网格下研究二阶椭圆特征值问题的两类非协调有限元—类Wilson矩形元和Carey三角形元—的收敛性分析.通过新的技巧和方法,得到了与传统有限元网格剖分下相同的特征对的最优误差估计.推广了已有的结果.     

Comparison of adjoint‐based and feature‐based grid adaptation for functional outputs

Adjoint‐based and feature‐based grid adaptive strategies are compared for their robustness and effectiveness in improving the accuracy of functional outputs such as lift and drag coefficients. The output‐based adjoint approach strives to improve the adjoint error estimates that relate the local residual errors to the global error in an output function via adjoint variables as weight functions. A conservative adaptive indicator that takes into account the residual errors in both the primal (flow) and dual (adjoint) solutions is implemented for the adjoint approach. The physics‐based feature approach strives to identify and resolve significant features of the flow to improve functional accuracy. Adaptive indicators that represent expansions and compressions in the flow direction and gradients normal to the flow direction are implemented for the feature approach. The adaptive approaches are compared for functional outputs of three‐dimensional arbitrary Mach number flow simulations on mixed‐element unstructured meshes. Grid adaptation is performed with h‐refinement and results are presented for inviscid, laminar and turbulent flows. Copyright © 2006 John Wiley & Sons, Ltd.