基于偏折光谱的非线性迭代层析数值模拟.pdfVIP

第29卷，第2期 光谱学与光谱分析 200 and 9年2月 SpectralAnalysis February，2009 Spectroscopy 基于偏折光谱的非线性迭代层析数值模拟 宋一中1’2，赵志敏1 1．南京航空航天大学理学院。江苏南京210016 2．德州学院物理系，山东德州253000 摘要模拟流场，校验自研发的偏折层析系统。根据光线偏折的光学原理，以及层析的数学、物理意义， self-correlative 研究开发了偏折层析程序。结合自研发的简单自相关代数迭代重建技术(simple algebraic reconstruction technique：SSART)，开发了一种新的层析系统，文章称“基于偏折光谱的非线性迭代层析系 统”。用Gauss函数模拟复杂流场的一个切片，根据偏折层析算法计算该切片的偏折光谱，从偏折光谱中提 取投影。采用SSART迭代重建算法重建截面。比较模拟重建结果与原场的差别，采用均方差(mean-square error：MSE)和峰差(peakerror：PE)指数标定重建效果，以检验算法的有效性。结果发现，应用开发的偏折 09～o．000 11范围，PE稳定 层析程序，采用SSART迭代算法，能够精确重建模拟场，MSE稳定在0．000 在0．007～0．013之间。因此，基于偏折光谱的非线性迭代层析系统能够根据偏折信息精确重建复杂流场。 关键词层析；复杂流场；迭代方法；数值模拟；误差 中图分类号：029 文献标识码：A 1000-0593{2009)02-0561-04 DOI：10．3964／j．issrL 拟复杂三维流场的一个切片 引 言 应用光学计算机层析(OCT)技术[1’2]，将流场层析并用 (y一0．4)2)) (1) 图像显示出来，实现流场可视化[3“]。这是一种重要的流场 函数图像如图1所示。从图l可见，模拟流场截面含有 可视化解决方案，人们已做了大量的研究[5一]，但尚未开发 两个正的Gauss峰和一个负的Gauss峰。 出适合等离子体等高温流场检测的有效层析算法[7’8]。应用 该技术显示复杂的高温流场时，需要一种收敛速度快、稳定 性好、精度高，而且能重建含遮挡物、有限采样角等实际工 况下复杂流场的重建算法阳J…。为了开发这样的算法，许多 研究人员已经做了大量探索[6．11]。本课题组也从投影、数据 提取到重建层析的各个环节，都做了大量研究[12’”]。在前人 工作的基础上，本文开发了非线性偏折层析程序专门用于高 温复杂流场的层析。该程序的实用性首先需要通过数值模拟 的方法进行检验[14,15]，而复杂流场是较好的一种检验模型。 本文用自开发的偏折层析程序，采用ssART迭代方法 重建复杂流场，检验偏折层析程序的有效性。 1三维流场模型 层析是重建三维场的切片口，151。根据数值模拟层析算 法，本文用一个非对称、含有正负相的三Gauss峰函数来模 收稿日期：2007—12—25。修订日期：2008-03—26 作者简介：宋一中，1971年生，南京航空航天大学物理学院博士后e-mail：yizhongsong@126．啪 万方数据 582 光谱学与光谱分析