MAX-DOAS仪器观测O4吸收对气溶胶参量的敏感性研究

研究了基于地基差分吸收光谱技术观测的O4吸收反演气溶胶光学参量的敏感性.利用大气辐射传输模型McArtim分析了不同波长、不同气溶胶光学参量(光学厚度、边界层高度、单次散射反照率、非对称因子)对O4吸收(大气质量因子)的影响.结果表明,大气质量因子对中心波长变化不敏感,气溶胶光学厚度和边界层高度对O4吸收具有重要影响,气溶胶光学厚度从0.1增加到1时,3°仰角测量的大气质量因子减少了28%,边界层高度从0.1 km增加到1 km,3°仰角测量的大气质量因子增加了9.2%.平纬圈观测模式下O4吸收对单次散射反照率和非对称因子具有较好的灵敏度,提供了一种基于地基MAX-DOAS观测O4吸收反演气溶胶光学参量的新方法.     

涡轮叶栅流动分离涡模拟验证研究

本文采用DDES模型分析透平叶栅内部亚音速和跨音速流动。DDES模型能有效且较精确地计算亚音速流动,但计算跨音速流动时,出现修正函数f_d"保护"边界层功能失效的问题,极大地降低了计算精度;人为设定RANS和LES计算域,保证边界层完全被RANS计算域覆盖,计算结果得到很大改善。     

基于MATLAB的黑体辐射量计算

介绍了黑体辐射理论,利用MATLAB求解了普朗克黑体辐射公式在有界域λ1～λ2内的积分,通过实例对梯形法和辛普森法两种数值积分方法进行了分析和比较。并在此基础上对距黑体一定距离处探测器表面的辐照度值E进行了研究。通过分析可知,利用梯形法和辛普森法能够方便求出黑体在有界域内的辐出度值。此外,利用MATLAB编制的相应的计算界面,能够实现辐出度和辐照度的快速计算。     

MAX-DOAS仪器观测O_4吸收对气溶胶参量的敏感性研究

研究了基于地基差分吸收光谱技术观测的O_4吸收反演气溶胶光学参量的敏感性.利用大气辐射传输模型McArtim分析了不同波长、不同气溶胶光学参量(光学厚度、边界层高度、单次散射反照率、非对称因子)对O_4吸收(大气质量因子)的影响.结果表明,大气质量因子对中心波长变化不敏感,气溶胶光学厚度和边界层高度对O_4吸收具有重要影响,气溶胶光学厚度从0.1增加到1时,3°仰角测量的大气质量因子减少了28%,边界层高度从0.1km增加到1km,3°仰角测量的大气质量因子增加了9.2%.平纬圈观测模式下O_4吸收对单次散射反照率和非对称因子具有较好的灵敏度,提供了一种基于地基MAX-DOAS观测O_4吸收反演气溶胶光学参量的新方法.     

轴流压气机环壁边界层的流量堵塞效应

本文讨论了轴流压气机环壁边界层的计算,它对流量的堵塞效应和在S_2流面计算中对此进行合理修正的方法.首先介绍了一种适合于叶轮机械设计和分析用的边界层的计算方法.在此基础上进行了主流与边界层的迭代计算。对三种迭代方案进行了比较.本文给出了一台单级和一台十级轴流压气机的算例以及计算与激光测量结果的比较.     

三维边界层内诱导横流失稳模态的感受性机理

边界层感受性问题是层流向湍流转捩的初始阶段,在转捩过程中起关键性作用,尤其是三维边界层流动.因此,研究三维边界层感受性问题对进一步理解层流向湍流转捩机理以及湍流成因具有重要的理论意义.采用数值方法研究自由来流湍流与三维壁面局部粗糙相互作用下三维边界层的感受性问题,确定是否能在三维边界层内寻找一种新的横流失稳模态;确定在何种条件下三维边界层内能诱导出定常、非定常的横流失稳模态;探索自由来流湍流的强度、展向波数和法向波数以及三维壁面局部粗糙的大小和结构类型等因素在自由来流湍流与三维壁面局部粗糙作用下三维边界层内被激发出的感受性过程中有何影响,并确定何种横流失稳模态在三维边界层感受性过程中占据何种地位.对自由来流湍流与三维壁面局部粗糙作用激发三维边界层内感受性问题的深入研究,将有助于完善流动稳定性与湍流理论,为层流向湍流转捩过程的预测与控制提供合理的理论依据.     

确定大气边界层顶高度的新方法及数值实验

提出了一种确定大气边界层顶高度的数值微分新方法,该方法使用了正则化技术,把对弯角廓线求导数的数值微分问题转化为求目标泛函极小值的问题,采用双参数模型函数方法来选择正则化参数,最后利用最大梯度法确定边界层顶高度.首先通过两个数值实验验证了新方法的有效性,实验结果显示,随着掩星资料噪音的增多,由差分法和结合L曲线方案的数值微分方法得到的边界层顶高度波动增大,而通过双参数模型函数方法得到的高度很稳定,这说明新方法能够很好地过滤噪音,从而保留廓线的主要信息.随后基于2007-2011年1,4,7,10月的COSMIC弯角数据,利用新方法分析了全球海洋大气边界层顶高度的季节特征,并与用掩星资料自带的大气边界层顶高度数据zbalmax得到的季节分布进行对比.结果表明,两者的季节分布特征十分一致:海温相对周围海域高的区域,边界层顶高度较高,反之,边界层顶高度较低;在暖流经过的海域,边界层顶高度较高,在寒流经过的海域,边界层顶的高度相对较低.     

一类Marangoni对流边界层方程的近似解析解

利用Adomain解析拆分和Padé逼近方法对由Marangoni对流诱发的层流边界层问题进行了研究, 提供了一种求解边界层方程的解析分析方法. 得到了问题的近似解析解并对相应的流动及传热特性进行了探讨. 本文所提出的思想方法可以用于解决其他科学和工程技术问题. 关键词： Marangoni对流 非线性 Adomain拆分法 近似解析解     

前缘曲率变化对平板边界层感受性问题的影响

边界层感受性问题是层流向湍流转捩的初始阶段,是实现边界层转捩预测和控制的关键环节.目前已有的研究成果显示,在声波扰动或涡波扰动作用下前缘曲率变化对边界层感受性机制有着显著的影响.本文采用直接数值模拟方法,研究了在自由来流湍流作用下具有不同椭圆形前缘平板边界层感受性问题,揭示椭圆形前缘曲率变化对平板边界层内被激发出Tollmien-Schlichting (T-S)波波包的感受性机制以及波包向前传播群速度的影响;通过快速傅里叶分析方法从波包中提取获得了不同频率的T-S波,详细分析了前缘曲率变化对不同频率的T-S波的幅值、色散关系、增长率、相速度以及形状函数的作用;确定了前缘曲率在平板边界层内激发T-S波的感受性过程中所占据的地位.通过上述研究能够进一步认识和理解边界层感受性机制,从而丰富和完善了流动稳定性理论.     

谐和激励与有界噪声作用下具有同宿和异宿轨道的Duffing振子的混沌运动

研究了具有同宿轨道、异宿轨道的双势阱Duffing振子在谐和激励与有界噪声摄动下的混沌运动.基于同宿分叉和异宿分叉，由Melnikov理论推导了系统出现混沌运动的必要条件及出现分形边界的充分条件.结果表明：当Wiener过程的强度参数大于某一临界值时，噪声增大了诱发混沌运动的有界噪声的临界幅值，相应地缩小了参数空间的混沌域，且产生混沌运动的临界幅值随着噪声强度的增大而增大.同时数值计算了最大Lyapunov指数，由最大Lyapunov指数为零从另一角度得到了系统出现混沌运动的有界噪声的临界幅值，发现在Wi 关键词： 混沌 同宿和异宿分叉 随机Melnikov方法 最大Lyapunov指数     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

Near-field performance analysis of locally-conformal perfectly matched absorbers via Monte Carlo simulations

In the numerical solution of some boundary value problems by the finite element method (FEM), the unbounded domain must be truncated by an artificial absorbing boundary or layer to have a bounded computational domain. The perfectly matched layer (PML) approach is based on the truncation of the computational domain by a reflectionless artificial layer which absorbs outgoing waves regardless of their frequency and angle of incidence. In this paper, we present the near-field numerical performance analysis of our new PML approach, which we call as locally-conformal PML, using Monte Carlo simulations. The locally-conformal PML method is an easily implementable conformal PML implementation, to the problem of mesh truncation in the FEM. The most distinguished feature of the method is its simplicity and flexibility to design conformal PMLs over challenging geometries, especially those with curvature discontinuities, in a straightforward way without using artificial absorbers. The method is based on a special complex coordinate transformation which is ‘locally-defined’ for each point inside the PML region. The method can be implemented in an existing FEM software by just replacing the nodal coordinates inside the PML region by their complex counterparts obtained via complex coordinate transformation. We first introduce the analytical derivation of the locally-conformal PML method for the FEM solution of the two-dimensional scalar Helmholtz equation arising in the mathematical modeling of various steady-state (or, time-harmonic) wave phenomena. Then, we carry out its numerical performance analysis by means of some Monte Carlo simulations which consider both the problem of constructing the two-dimensional Green’s function, and some specific cases of electromagnetic scattering.     

Eigenmodes of index-modulated layers with lateral PMLs

Maxwell equations are solved in a layer comprising a finite number of homogeneous isotropic dielectric regions ended by anisotropic perfectly matched layers (PMLs). The boundary-value problem is solved and the dispersion relation inside the PML is derived. The general expression of the eigenvalues equation for an arbitrary number of regions in each layer is obtained, and both polarization modes are considered. The modal functions of a single layer ended by PMLs are found, and their orthogonality relation is derived. The present method is useful to simulate scattering problems from dielectric objects as well as propagation in planar slab waveguides. Its potential to deal with more complex problems such as the scattering from an object with arbitrary cross section in open space using the multilayer modal method is briefly discussed.     

层状介质时域有限差分方法斜入射平面波引入新方式

应用二维时域有限差分方法分析层状介质中的目标散射时,在总场-散射场边界斜入射平面波源用常规方法难以引入,因为在总场-散射场边界处设置的入射波实际上包含了入射脉冲以及各分层界面的反射和多次反射.为解决这个问题,提出了斜入射平面波的混合引入方式,即对总场-散射场的四个边界面采取不同的处理方式.对于总场-散射场的纵向侧边界,用含有斜入射角度的修正一维时域有限差分方法,只要在自由空间位置加入入射脉冲就会自行产生由各分层界面形成的反射波,包括多次反射.同时,把纵向总场-散射场侧边界向下延伸,使得总场-散射场下边界位于完全匹配层内,这样透射波和散射波均为外向行波而被吸收.对于总场-散射场的上边界,由于完全位于自由空间中,边界上各点的入射波将是总场-散射场纵向边界角点处入射波的带有时间延迟的复制.数值模拟结果表明了本文所提出方法的正确性和有效性. 关键词： 时域有限差分 层状介质 斜入射平面波 修正一维麦克斯韦方程     

一种新时域交替隐式差分算法在散射问题中的应用

提出一种新时域交替隐式有限差分(ADI-FDTD)算法格式. 传统ADI-FDTD算法的 场量步进方程涉及周围若干网格的较多场量，导致两个区域的步进方程处理较困难：一个是 邻近完全匹配层(PML)和散射场交界区，另一个是邻近连接边界区. 特别是后者，考虑入射 波影响需对场量所在区域判断，根据不同情况对原有方程进行修正，一维和二维散射问题 相对简单，可三维问题修正式有数十种之多而几乎无法完成. 本方法基于分裂场形式的ADI- FDTD技术，使得散射场区和PML吸收层区的表达形式完全一致，从而忽略两者差别.另 关键词： 时域交替隐式有限差分算法 电磁散射     

An Iterative Two-Grid Method of a Finite Element PML Approximation for the Two Dimensional Maxwell Problem

In this paper, we propose an iterative two-grid method for the edge finite element discretizations (a saddle-point system) of Perfectly Matched Layer (PML) equations to the Maxwell scattering problem in two dimensions. Firstly, we use a fine space to solve a discrete saddle-point system of $H(grad)$ variational problems, denoted by auxiliary system 1. Secondly, we use a coarse space to solve the original saddle-point system. Then, we use a fine space again to solve a discrete$\boldsymbol{H}(curl)$-elliptic variational problems, denoted by auxiliary system 2. Furthermore, we develop a regularization diagonal block preconditioner for auxiliary system 1 and use $H$-$X$ preconditioner for auxiliary system 2. Hence we essentially transform the original problem in a fine space to a corresponding (but much smaller) problem on a coarse space, due to the fact that the above two preconditioners are efficient and stable. Compared with some existing iterative methods for solving saddle-point systems, such as PMinres, numerical experiments show the competitive performance of our iterative two-grid method.     

偏振-米散射激光雷达的研制

研制的偏振-米散射激光雷达（PML）,可用于探测卷云和沙尘气溶胶的后向散射光退偏振比以及研究流层大气气溶胶的消光特性。采用窄带滤光片和光阑,将接收到的激光大气回波信号谱线（米散射和瑞利散射光谱）从天空太阳背景噪声中分离出来,以提高系统的白天探测能力。介绍了偏振-米散射激光雷达的结构、技术参数、测量方法和数据处理方法。对偏振-米散射激光雷达的性能参数进行了测定,并对测定结果进行了分析与讨论,给出了偏振-米散射激光雷达对合肥市地区(117.16°E, 31.90°N)上空大气气溶胶的消光特性和卷云的结构、退偏振比垂直廓线以及光学厚度的典型探测结果,对这些结果进行了分析和讨论。结果表明：研制的偏振-米散射激光雷达性能可靠,能对大气气溶胶和卷云的物理和光学特性进行有效的探测。     

A computationally efficient finite element model with perfectly matched layers applied to scattering from axially symmetric objects

A frequency-domain finite-element (FE) technique for computing the radiation and scattering from axially symmetric fluid-loaded structures subject to a nonsymmetric forcing field is presented. The Berenger perfectly matched layer (PML), applied directly at the fluid-structure interface, makes it possible to emulate the Sommerfeld radiation condition using FE meshes of minimal size. For those cases where the acoustic field is computed over a band of frequencies, the meshing process is simplified by the use of a wavelength-dependent rescaling of the PML coordinates. Quantitative geometry discretization guidelines are obtained from a priori estimates of small-scale structural wavelengths, which dominate the acoustic field at low to mid frequencies. One particularly useful feature of the PML is that it can be applied across the interface between different fluids. This makes it possible to use the present tool to solve problems where the radiating or scattering objects are located inside a layered fluid medium. The proposed technique is verified by comparison with analytical solutions and with validated numerical models. The solutions presented show close agreement for a set of test problems ranging from scattering to underwater propagation.     

基于逆散射理论的地震波速度正则化反演

本文基于逆散射理论利用正则化有限差分对比源反演算法对地震波传播速度进行反演, 该方法是基于波动方程的频率域波形反演算法, 利用非线性共轭梯度法, 通过最小化目标优化函数不断迭代更新速度模型. 由于地球物理反演问题的病态性和不稳定性, 通过基于反演参数总变差的正则化处理, 使反演问题变为良性问题且算法具有较强的抗噪声干扰能力. 反演过程中使用了频率-空间域9点差分正演算子以及PML吸收边界条件. 与其他反演算法相比, 由于背景模型在反演迭代过程中保持不变, 可以避免在每次迭代过程中重新构造正演算子及矩阵分解等相关计算过程, 使得该算法非常适合于大规模三维反演计算. 此外, 本文采用基于MPI的并行计算, 进一步提高了反演计算的效率. 二维CSEG模型反演结果表明该方法可以反演得到高分辨率的地震波速度重建结果, 为地震勘探数据处理及解释提供准确的速度信息.     

Optical modelling of photonic crystals and VCSELs using eigenmode expansion and perfectly matched layers

We present a modelling approach for photonic crystal structures and vertical-cavity surface-emitting lasers (VCSELs). This method is based on vectorial eigenmode expansion combined with perfectly matched layer (PML) boundary conditions. Compared to other methods, a relatively small computational effort is required, while at the same time accurate results are obtained, even in the presence of strong scattering and diffraction losses.