基于数据重建的三角网格模型简化优化方法

由CT切片数据重建得到的三角网格模型常存在数据量大、狭长三角形多等问题,针对这些问题,研究了一种保持特征的高质量三角网格模型的简化优化方法。该方法分为网格简化和网格优化两个阶段。首先,采用基于曲面变化的二次误差度量计算边折叠代价,并按代价值的大小进行迭代的折叠简化,可较好地保持模型表面的特征;其次,通过二阶加权伞算子对简化模型中局部存在的狭长三角形进行优化处理,改善三角网格模型的质量。实验结果表明,该方法能够较好地保持特征区域的细节信息,并可靠地生成高质量、低几何误差的简化模型。     

基于工业CT图像重构的网格模型的保特征简化方法

对基于工业CT图像重构的网格模型进行网格简化时,大多数现有网格模型简化算法会丢失特征,出现网格质量不好的问题。因此提出一种网格模型保特征简化方法,该方法用三角形折叠法对原始模型进行简化,当简化后模型的平均二面角角度误差达到允许误差后,再使用边折叠法对模型进行简化。在三角形折叠法中提出了利用被折叠三角形的法向量、各个顶点的高斯曲率及其在周边三角形上的投影确定该三角形的折叠点,利用局部体积误差与二面角角度误差的无因次化和确定折叠代价的方法;在边折叠法中提出了将二面角角度误差引入到二次误差测度(QEM)法的折叠代价中的改进QEM法。实验结果表明:与其他算法相比,该方法能够生成保特征、高质量、低几何误差的网格模型。     

保特征的自适应三角网格规范化算法

为了提高基于工业CT图像重构的三角网格质量，提出了顶点调整和特征因子相结合的网格规范化算法。引入模型特征因子和局部网格质量提高程度作为网格调整的控制条件，保留原始模型的局部细节特征；采用法矢量阈值按规律递增的方式自适应控制网格调整，实现不同曲率特征的自动识别。实验结果表明，与现有方法相比，该方法能更好地规范模型的三角网格，同时保留了原始模型的细节特征。     

基于约束Delaunay三角化的二维非结构网格生成方法

给出基于局部重构和边交换技术的两种约束Delaunay三角剖分方法并证明其收敛性.采用边界指示法恢复流场形状;在预设尺度的指导下融合流场边界曲率、中轴线、梯度限制等信息修正流场尺度;运用Spring方法布置边界点,通过符号面积函数和概率筛选法布置计算区域节点;运用Spring-Laplace方法优化节点位置,伴同边交换和边吞噬技术优化网格结构.该方法可自由进行局部自适应加密或稀疏,并应用于映射曲面网格生成和移动网格技术.     

基于工业CT图像重构的网格模型的全局优化方法

针对用工业CT切片图像直接重构得到的网格模型质量不高的问题,提出一种不受拓扑结构限制的隐式曲面重构全局优化方法。该方法将三维表面模型用隐式函数来表示,通过模型提供的点云信息计算出隐式函数,提取等值面,实现曲面重构。针对隐式曲面重构数据处理量大的问题,引入FFTW快速傅里叶变换来提高效率。实验结果表明,该方法能够同时实现三角网格模型的去噪、网格平滑、简化以及孔洞修补,与保特征的均匀化网格平滑算法相比,去噪效果更好,效率更高。     

基于相关法的哈特曼图像网格自动划分方法

在哈特曼波前测量过程中,需要对哈特曼图像进行网格划分用于计算波前斜率。传统的网格划分方法通过观察哈特曼图像手动调整网格位置和网格间距。该方法不但需要人工干预,还给波前实时测量带来困难。为实现哈特曼图像网格划分的自动化,研究了基于相关法的哈特曼图像网格自动划分方法。该方法首先对哈特曼图像进行自相关求取网格行列数以及网格间距,然后利用该网格建立相应的模拟光斑图像,最后对模拟光斑图像和哈特曼图像进行互相关求取网格最佳位置。该方法无需人工干预,耗时短,实现了网格实时自动划分。     

三维物理-热工耦合过程的网格优化

在反应堆物理-热工耦合过程中,网格划分尺度会影响计算精度和计算时间。利用蒙特卡罗程序和FLUENT程序,对压水堆单棒模型进行不同尺度的网格划分,评估网格划分尺度对耦合结果的影响,得到单个网格中密度差值、温度差值对有效增殖因子和功率分布引入的误差。研究表明当燃料温度差值小于50 K,慢化剂密度差值3 kg/m3时,有效增殖因子相对误差小于10-4,功率相对误差小于1%。使用该规律,对典型的压水堆单棒模型和33通道模型进行网格划分并进行耦合计算。结果表明,单棒模型网格总数减少至1/100,计算时间减少至1/4,33通道模型网格总数减少至1/50,计算时间减少至1/10,但其结果仍然精确有效。     

流动数值模拟中一种并行自适应有限元算法

给出了一种流动数值模拟中的基于误差估算的并行网格自适应有限元算法.首先,以初网格上获得的当地事后误差估算值为权,应用递归谱对剖分方法划分初网格,使各子域上总体误差近似相等,以解决负载平衡问题.然后以误差值为判据对各子域内网格进行独立的自适应处理.最后应用基于粘接元的区域分裂法在非匹配的网格上求解N-S方程.区域分裂情形下N-S方程有限元解的误差估算则是广义Stokes问题误差估算方法的推广.为验证方法的可靠性,给出了不可压流经典算例的数值结果.     

多点择优推进阵面法生成曲面三角形网格

在现有曲面非结构网格生成法的基础上,提出了一种新的曲面网格生成法——多点择优推进阵面法.它可在曲面上直接进行三角形网格划分,克服了映射法的网格变形问题,并且可以在网格生成结束后,对曲面网格直接进行Laplace格点松弛光顺.该方法使用简单,不受曲面块类型的限制,且网格质量高,可以为三维非结构网格生成提供高质量的初始阵面,并给出了若干个算例.     

水平集方程在四边形网格上的数值离散方法

在多介质大变形数值计算中，正交的笛卡儿网格很难满足计算要求，经常要在一般四边形网格上进行数值计算，因此需要在一般四边形网格上实现波阵面的捕捉算法。水平集方法不需显式地构造界面，容易处理界面拓扑结构的变化，因此选用水平集方法来捕捉波阵面。无曲率的水平集方程、重新初始化方程是齐次Hamilton-Jacobi方程，可以直接运用正格式进行离散。文中只考虑曲率项是曲率的线性函数的情况，与曲率相关的水平集方程含有水平集函数的二阶导数，不能采用正格式离散，采用伽辽金有限元方法进行离散。     

偏振成像在光学曲率半径测量中的应用

为了保证整个光学系统的质量,对光学曲率半径的准确测量与检验至关重要.本文将机械球径仪法与光学投影法结合,采用光电图像法对矢高进行判读,经计算可以间接测得曲率半径,具体分析了被测元件边缘特征的准确性对矢高及曲率半径测量精度的影响.本文分别采用偏振成像与普通成像的方法进行了测量与对比,发现偏振图像具有更好的边缘细节特征,矢...     

Accuracy and analysis of long-radius measurement with long trace profiler

The long trace profiler(LTP) is proposed to measure radius of curvature(R) and surface figure of a longradius spherical surface in an optical shop.Equipped with a motorized rotary stage and a two-dimensional tilt stage,the LTP scans the full aperture and calculates the absolute radius of curvature of each scanning line based on the least square method.Nonlinear error and manufacture error difference between center and the edge are obtained by comparing R results.The R-limit is validated and expressed as D/R,where D is the aperture of the mirror under test.A full-aperture three-dimensional figure is also reconstructed based on triangle interpolation.     

Reduced surface point selection options for efficient mesh deformation using radial basis functions

Previous work by the authors has developed an efficient method for using radial basis functions (RBFs) to achieve high quality mesh deformation for large meshes. For volume mesh deformation driven by surface motion, the RBF system can become impractical for large meshes due to the large number of surface (control) points, and so a particularly effective data reduction scheme has been developed to vastly reduce the number of surface points used. The method uses a chosen error function on the surface mesh to select a reduced subset of the surface points; this subset contains a sufficiently small number of points so as to make the volume deformation fast, and a correction function is used to correct those surface points not included. Hence, the scheme is split such that both parts are working on appropriate problems. RBFs are an excellent way of finding smooth orthogonality preserving global deformations, but are less suitable for enforcing an exact geometry for a large number of points, while a simpler approach is ideal for diffusing small changes evenly but has quality (and possibly expense) drawbacks if used for the entire volume. However, alternatives exist for the error function used to select the reduced data set, so here a comparison is made between three different options: the surface error function, the unit function and the power function. Tests run on structured and unstructured meshes show that the surface error function gives the lowest errors, but this also requires a deformed surface shape to be known in advance of the simulation. The unit and power functions both avoid the need for a deformed surface, and the unit function is shown to be superior.     

Reducing the line curvature error of mechanically ruled gratings by interferometric control

Line curvature error greatly influences the quality of the diffraction wave fronts of machine-ruling gratings. To reduce the line curvature error, we propose a correction method that uses interferometric control. This method uses diffraction wave fronts of symmetrical orders to compute the mean line curvature error of the ruled grating, taking the mean line curvature error as the system line curvature error. To minimize the line curvature error of the grating, a dual-frequency laser interferometer is used as a real-time position feedback for the grating ruling stage, along with using a piezoelectric actuator to adjust the stage positioning to compensate the line curvature error. Our experiments show that the proposed method effectively reduced the peak-to-valley value of the line curvature error, improving the quality of the grating diffraction wave front.     

Metric construction by length distribution tensor and edge based error for anisotropic adaptive meshing

Metric tensors play a key role to control the generation of unstructured anisotropic meshes. In practice, the most well established error analysis enables to calculate a metric tensor on an element basis. In this paper, we propose to build a metric field directly at the nodes of the mesh for a direct use in the meshing tools. First, the unit mesh metric is defined and well justified on a node basis, by using the statistical concept of length distribution tensors. Then, the interpolation error analysis is performed on the projected approximate scalar field along the edges. The error estimate is established on each edge whatever the dimension is. It enables to calculate a stretching factor providing a new edge length distribution, its associated tensor and the corresponding metric. The optimal stretching factor field is obtained by solving an optimization problem under the constraint of a fixed number of edges in the mesh. Several examples of interpolation error are proposed as well as preliminary results of anisotropic adaptation for interface and free surface problem using a level set method.     

Temperature,roughness and depth resolution in ion sputter profiles

The depth resolutions of evaporated silver layers on polished polycrystalline copper substrates are studied at various temperatures during argon ion sputtemg with AES. A detailed analysis of the profile shape at the interface reveals the nature of contributions to the terms governing interface resolution. The profiles are all accurately described by an error function leading edge followed by an exponential trailing edge. The characteristic of the trailing edge is governed by the development of roughness which depends on the depth, z, sputtered as z^0.875. The roughness is reduced at elevated temperatures by the action of surface diffusion such that the maximum reduction occurs at the highest temperature and the lowest sputtering rate. The characteristic of the leading edge is composed of three terms added in quadrature, (i) defect-enhanced diffusion of copper into the silver film, (ii) roughness as above, and (iii) a constant term due to film nucleation. In the present samples the increased interdiffusion at elevated temperatures and low sputtering rates largely offsets the improvements in topography so that, overall, the depth resolution appears to be a very weak function of temperature. However, in other systems where interdiffusion is small, the resolution could be greatly enhanced.     

Efficient mesh motion using radial basis functions with data reduction algorithms

Mesh motion using radial basis functions has been demonstrated previously by the authors to produce high quality meshes suitable for use within unsteady and aeroelastic CFD codes. In the aeroelastic case the structural mesh may be used as the set of control points governing the deformation, which is efficient since the structural mesh is usually small. However, as a stand alone mesh motion tool, where the surface mesh points control the motion, radial basis functions may be restricted by the size of the surface mesh, as an update of a single volume point depends on all surface points. In this paper a method is presented that allows an arbitrary deformation to be represented to within a desired tolerance by using a significantly reduced set of surface points intelligently identified in a fashion that minimises the error in the interpolated surface. This method may be used on much larger cases and is successfully demonstrated here for a 10⁶ cell mesh, where the initial solve phase cost reduces by a factor of eight with the new scheme and the mesh update by a factor of 55. It has also been shown that the number of surface points required to represent the surface is only geometry dependent (i.e. grid size independent), and so this reduction factor actually increases for larger meshes.