Validation of adaptive unstructured hexahedral mesh computations of flow around a wind turbine airfoil

In this paper we investigate local adaptive refinement of unstructured hexahedral meshes for computations of the flow around the DU91 wind turbine airfoil. This is a 25% thick airfoil, found at the mid‐span section of a wind turbine blade. Wind turbine applications typically involve unsteady flows due to changes in the angle of attack and to unsteady flow separation at high angles of attack. In order to obtain reasonably accurate results for all these conditions one should use a mesh which is refined in many regions, which is not computationally efficient. Our solution is to apply an automated mesh adaptation technique. In this paper we test an adaptive refinement strategy developed for unstructured hexahedral meshes for steady flow conditions. The automated mesh adaptation is based on local flow sensors for pressure, velocity, density or a combination of these flow variables. This way the mesh is refined only in those regions necessary for high accuracy, retaining computational efficiency. A validation study is performed for two cases: attached flow at an angle of 6° and separated flow at 12°. The results obtained using our adaptive mesh strategy are compared with experimental data and with results obtained with an equally sized non‐adapted mesh. From these computations it can be concluded that for a given computing time, adapted meshes result in solutions closer to the experimental data compared to non‐adapted meshes for attached flow. Finally, we show results for unsteady computations. Copyright © 2005 John Wiley & Sons, Ltd.     

黏性不可压流体的自适应网格技术和基本特性方程分离算法的联合分析

组合基本特性方程分离算法和自适应网格技术,分析二维黏性不可压流体.该方法使用3节点三角单元,对速度分量和压力等变量分析,使用等阶次的插值函数.组合解法的主要优点在于,在自适应网格技术中,对解梯度变化大的区域,通过耦合误差估计生成小的单元,利于提高解的精度,在其它区域生成大单元,可以节省时间.最后,通过对一个黏性流体圆柱体绕流问题的瞬态和稳态特性分析,给出了组合解法性能的评价.     

高精度TVD有限差分格式的自由梯度条件

1自由流条件守恒性气动方程组的高精度TVD有限差分格式一般通过使用通量限制器加入非线性的人工粘性来构成。在曲线坐标系中，通量限制器比较的两个因子与各自所在单元的界面雅可比转换行列式相关。为了消除畸变网格单元界面雅可比转换行列式的变化对通量限制器比较结果的影响，本文提出了高精度TVD格式的自由梯度性质，并分别就Osher和HartenTVD格式讨论了同时满足自由流和自由梯度的度量张量计算方法，使高精度TVD格式在畸变网格上得以实现。任意曲线坐标系中守恒型二维欧拉方程表示为其中，E、F为笛卡尔坐标系中的通量；J为坐标变…     

Surface triangulation for polygonal models based on CAD data

This paper presents an approach to the generation of unstructured surface meshes for Computer‐Aided Design (CAD) surface models represented as lists of polygons with minimum user interventions. Stereolithography (STL) data are adopted as an interface between a CAD system and the surface grid generator. STL files often include problems such as overlapping surfaces, gaps, and intersections. They have to be revised quickly and automatically before the surface models are used for the background grid of the surface grid generation. In this paper, we describe an automatic revision method for use with STL‐defined surface models. The advancing front method using geometric features is adopted directly on the modified STL surfaces. The capability of the method is demonstrated for several 3D surface models. Copyright © 2002 John Wiley & Sons, Ltd.     

Anisotropic mesh adaptation for numerical solution of boundary value problems

We present an efficient mesh adaptation algorithm that can be successfully applied to numerical solutions of a wide range of 2D problems of physics and engineering described by partial differential equations. We are interested in the numerical solution of a general boundary value problem discretized on triangular grids. We formulate a necessary condition for properties of the triangulation on which the discretization error is below the prescribed tolerance and control this necessary condition by the interpolation error. For a sufficiently smooth function, we recall the strategy how to construct the mesh on which the interpolation error is below the prescribed tolerance. Solving the boundary value problem we apply this strategy to the smoothed approximate solution. The novelty of the method lies in the smoothing procedure that, followed by the anisotropic mesh adaptation (AMA) algorithm, leads to the significant improvement of numerical results. We apply AMA to the numerical solution of an elliptic equation where the exact solution is known and demonstrate practical aspects of the adaptation procedure: how to control the ratio between the longest and the shortest edge of the triangulation and how to control the transition of the coarsest part of the mesh to the finest one if the two length scales of all the triangles are clearly different. An example of the use of AMA for the physically relevant numerical simulation of a geometrically challenging industrial problem (inviscid transonic flow around NACA0012 profile) is presented. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004.     

一类无结构三角网上抛物方程的有限差分区域分解算法

本文讨论了一类在无结构三角网上数值求解二维热传导方程的有限差分区域分解算法．在这个算法中,将通过引进两类不同类型的内界点,将求解区域分裂成若干子区域．一旦内界点处的值被计算出来,其余子区域上的计算可完全并行．本文得到了稳定性条件和最大模误差估计,它表明我们的格式有令人满意的稳定性和较高的收敛阶．     

A novel finite point method for flow simulation

A novel finite point method is developed to simulate flow problems. The mashes in the traditional numerical methods are supplanted by the distribution of points in the calculation domain. A local interpolation based on the properties of Taylor series expansion is used to construct an approximation for unknown functions and their derivatives. An upwind‐dominated scheme is proposed to efficiently handle the non‐linear convection. Comparison with the finite difference solutions for the two‐dimensional driven cavity flow and the experimental results for flow around a cylinder shows that the present method is capable of satisfactorily predicting the flow separation characteristic. The present algorithm is simple and flexible for complex geometric boundary. The influence of the point distribution on computation time and accuracy of results is included. Copyright © 2002 John Wiley & Sons, Ltd.     

AN ACCURATE NUMERICAL SOLUTION OF A TWO DIMENSIONAL HEAT TRANSFER PROBLEM WITH A PARABOLIC BOUNDARY LAYER

1.IntroductionSingularlyperturbeddifferentialequationsarecharacterizedbythepresenceofasmallparameterEmultiplyingthehighestorderderivatives.Suchproblemsariseinmanyareasofappliedmathematics.Thesolutionsofsingularlyperturbeddifferentialequationstypicallyhave…     

QUADRILATERAL MESH

IntroductionQuadrilateral mesh Is widely used In the finite eleme     

A higher-order hybrid spline difference method on adaptive mesh for solving singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions

We propose a hybrid numerical scheme to discretize a class of singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions on an equidistributed grid. The hybrid difference scheme is developed by using a modified backward difference scheme in time, a combination of the cubic spline and exponential spline difference scheme in space. The proposed scheme uses a cubic spline difference scheme for the discretization of robin-boundary conditions. For the time discretization of the problem, we use the standard uniform mesh while a layer adapted equidistributed grid is generated for the spatial discretization. By equidistributing a curvature-based monitor function, the spatial adaptive grid is able to capture the presence of parabolic boundary layers without using any prior information about the solution. Parameter uniform error estimates are derived to illustrate an optimal convergence of first-order in time and second-order in space for the proposed discretization. The accuracy of the proposed scheme is confirmed by the numerical experiments that underpin the theoretical analysis.