粒子输运方程的线性间断有限元方法

将空间线性间断有限元方法应用于动态粒子输运方程的求解.数值算例表明,空间线性间断有限元方法在网格边界的数值精度方面明显高于指数格式和菱形格式,并且通量在时间上的微分曲线相对光滑,避免了指数格式、菱形格式数值解的非物理振荡现象.     

一种高精度差分格式

在微分－Ｔｈｏｍｐｓｏｎ变换结合时域有限差分技术计算复杂目标的电磁散射特性的方法中，差分格式的构造和选取与解的精度存在着密切关系。提出一种新的高精度差分格式，其数值实现进一步证实了该方法能精确模拟任意形状目标的电磁磁射过程。     

飞秒激光诱导自组织纳米光栅偏振散射特性研究

实验证实飞秒激光在透明材料内部诱导的周期性纳米条纹具有特殊的光学特性。采用时域有限差分方法(FDTD)分别对单元纳米光栅和纳米光栅阵列的散射特性进行研究,分析结果显示纳米光栅的散射对入射激光的偏振方向敏感;散射强度与入射激光波长相关,波长越长、散射越弱。对于纳米光栅阵列,平行于纳米光栅的散射强度比垂直于纳米光栅的散射强度小两万倍。最后从理论上验证了Ⅱ类波导和纳米光栅偏振导光是由纳米光栅的散射引起的而不是双折射效应。     

粒子模拟激光斜入射边界条件的处理

讨论了电磁场模拟中的非反射边界条件,详细推导了斜入射情况非反射边界条件的有限差分格式,并在粒子模拟程序LARED-P上将其实现,同时给出了激光斜入射时的数值模拟结果.     

多模激光场中电子-原子散射

利用微扰理论对多模激光场中电子-原子散射进行了初步研究. 在激光场的极化方向平行于入射电子方向这种特殊的散射模式下, 计算得到了电子-原子微分散射截面与激光场的相移、交换光子数之间的关系. 结论和其它相关理论是一致的.     

一维球几何菱形差分中子输运方程的扩散综合加速解法

本文导出了一维球几何定态中子输运方程菱形格式的扩散综合加速方程,并给出了差分公式。所给出的加速方法可以加速菱形格式的输运方程的迭代求解。并给出了部分模型的数值计算结果。     

基于数值方法的激光海面漫反射通信能量分布特性

 针对激光点对点通信方式的不足,提出了利用海面作为激光漫反射媒介进行组网通信的设想,采用前后向迭代的数值方法结合Green函数谱积分加速算法(FBM/SAA)对激光海面漫反射通信的能量分布特性进行了研究。通过对激光光束入射海面后产生的散射场的分析、计算及实验验证,得出了较为准确的2维激光海面双站散射系数,并对激光光束入射海面后的散射场进行了分析。研究结果表明：入射激光束经过粗糙海面散射后,能量大部分集中在前向散射区域上,而后向散射强度很弱,且在散射场的边缘处能量迅速衰减,说明了激光海面漫反射组网通信方法的可行性。通过与前后向迭代法和Kirchhoff近似方法的计算结果的比较,说明了FBM/SAA是一种高效、准确的计算方法。     

脉冲激光入射下介质薄膜双频互相关散射截面及功率展宽

 结合特殊的镀金聚酯薄膜表面的粗糙度和半球反射率,依据基尔霍夫近似及粗糙面脉冲散射互相关函数,数值计算了脉冲激光（1.06 μm）在不同入射角照射下该材料薄膜双频互相关散射截面随相干带宽频差和散射角的变化情况。并给出δ脉冲和高斯脉冲波入射下,其散射功率随时间和散射角的波形分布。计算结果表明：该薄膜材料激光双频互相关散射截面在镜反射方向有最大的散射峰值,在非镜反射方向上,其散射值随相干带宽频差的增大迅速减小,当窄脉冲垂直入射时,粗糙面散射功率展宽现象明显。     

二维三温能量方程的差分格式及平面靶数值模拟

通过统一二维直角坐标与柱坐标下三温能量方程的差分格式,实现了LARED-H(laser radiationelectron dynamic holehurm)程序既可以模拟正入射激光平面靶过程,也可模拟斜入射平面靶问题,从而扩展了LARED-H程序的功能.使用发展后的LARED-H程序进行了直角坐标系下线聚焦打靶耦合过程数值模拟,给出了不同入射角度下靶等离子体的膨胀过程.     

纺织纤维光散射特性的模拟与实验

纺织纤维的光散射特性在纺织材料微观结构、光学性质以及无损检测中至关重要。利用纺织纤维各向异性的结构特点改进现有的角谱展开法,得到纺织纤维对倾斜入射平面波的散射以及表征其光散射特性的Mueller矩阵。采用蒙特卡罗算法模拟了平行棉纤维束对倾斜入射偏振光的多次散射,并对其计算结果进行了实验验证。结果表明理论计算和实验测量的散射光斑在形状特征、光强分布上均一致,从而验证了蒙特卡洛方法模拟纺织纤维光散射的正确性。结果同时也表明蒙特卡罗方法在纺织材料光传播特性的理论研究中具有重要的指导意义。     

Evaluation of numerical scattering in finite volume method for solving radiative transfer equation by a central laser incidence model

The numerical scattering caused by spatial discretization in finite volume method is discussed. Based on an analysis of the generation process of numerical scattering, a physical model of central laser incidence to a two-dimensional rectangle containing semitransparent medium is established to validate the numerical scattering, with Monte Carlo method as benchmark, in which numerical scattering does not exist. Numerical scattering will be affected by spatial grid number, spatial differential schemes and spectral absorption coefficient. With the spatial grid number increasing, numerical scattering will be decreased. The accuracy of the diamond scheme is the highest, and the exponential scheme is a bit lower, the lowest accuracy of the three schemes is the step scheme. The tendency of numerical scattering is reverse, i.e., the step scheme produces minimum numerical scattering, and exponential scheme produces more, while the diamond scheme produces maximum among three methods. When the bias of absorption efficient is high, the numerical scattering cannot be eliminated only by increasing the grid number. If we set the direction of laser incidence as central axis, it can be seen that numerical scattering distributed symmetry along the axis, which can be called as symmetrical cross-scattering. All of the three schemes show symmetrical cross-scattering.     

A conservative numerical scheme for solving an autonomous Hamiltonian system

A new numerical scheme is proposed for solving Hamilton’s equations that possesses the properties of symplecticity. Just as in all symplectic schemes known to date, in this scheme the conservation laws of momentum and angular momentum are satisfied exactly. A property that distinguishes this scheme from known schemes is proved: in the new scheme, the energy conservation law is satisfied for a system of linear oscillators. The new numerical scheme is implicit and has the third order of accuracy with respect to the integration step. An algorithm is presented by which the accuracy of the scheme can be increased up to the fifth and higher orders. Exact and numerical solutions to the two-body problem, calculated by known schemes and by the scheme proposed here, are compared.     

Three-Step Predictor-Corrector of Exponential Fitting Method for Nonlinear Schroedinger Equations

We develop the three-step explicit and implicit schemes of exponential fitting methods. We use the three- step explicit exponential fitting scheme to predict an approximation, then use the three-step implicit exponential fitting scheme to correct this prediction. This combination is called the three-step predictor-corrector of exponential fitting method. The three-step predictor-corrector of exponential fitting method is applied to numerically compute the coupled nonlinear Schroedinger equation and the nonlinear Schroedinger equation with varying coefficients. The numerical results show that the scheme is highly accurate.     

Finite difference computation of acoustic scattering by small surface inhomogeneities and discontinuities

The use of finite difference schemes to compute the scattering of acoustic waves by surfaces made up of different materials with sharp surface discontinuities at the joints would, invariably, result in the generations of spurious reflected waves of numerical origin. Spurious scattered waves are produced even if a high-order scheme capable of resolving and supporting the propagation of the incident wave is used. This problem is of practical importance in jet engine duct acoustic computation. In this work, the basic reason for the generation of spurious numerical waves is first examined. It is known that when the governing partial differential equations of acoustics are discretized, one should only use the long waves of the computational scheme to represent or simulate the physical waves. The short waves of the computational scheme have entirely different propagation characteristics. They are the spurious numerical waves. A method by which high wave number components (short waves) in the wave scattering process is intentionally removed so as to minimize the scattering of spurious numerical waves is proposed. This method is implemented in several examples from computational aeroacoustics to illustrate its effectiveness, accuracy and efficiency. This method is also employed to compute the scattering of acoustic waves by scatterers, such as rigid wall acoustic liner splices, with width smaller than the computational mesh size. Good results are obtained when comparing with computed results using much smaller mesh size. The method is further extended for applications to computations of acoustic wave reflection and scattering by very small surface inhomogeneities with simple geometries.     

An analytic and numerical solution with spectral Green's function method for transport equation in spherical geometry

Since in many cases curvilinear geometry is more appropriate than cartesian geometry for precise modeling of the complex systems for reactor calculation, we have developed the spectral Green's function (SGF) method which is employed to obtain angular and scalar flux distributions in heterogeneous sphere geometry with isotropic scattering. In this study, we showed that the neutron transport problems of homogeneous spheres could be reduced to the solution of plane geometry equation.Finally, some results are discussed and compared with those already obtained by diamond difference scheme to test the accuracy of the results. The agreement is satisfactory. SGF method is very suitable for the numerical solution of the neutron transport equation with isotropic scattering.     

Study on diffraction characteristics of a planar diamond waveguide

A numerical simulation is performed to study the far-field diffraction properties of planar diamond waveguides.The far-field intensity distributions of a planar air waveguide and a diamond waveguide with different distances are given by numerical calculations.In the experiment,the diffraction patterns on the screen with different distances are recorded using a He-Ne laser as the light source,wherein the laser beam is coupled with and propagates in the diamond waveguide.The simulation results are found to be consistent with the experimental ones.     

A rapid boundary perturbation algorithm for scattering by families of rough surfaces

In this paper we describe a novel algorithm for the computation of scattering returns by families of rough surfaces. This algorithm makes explicit use of the fact that some scattering profiles of engineering interest (e.g., traveling ocean waves) come in branches parameterized analytically by a bifurcation quantity. Our approach delivers recursions which not only can be implemented to yield a rapid, robust and high-order numerical scheme, but also give a new proof of analyticity of scattering quantities with respect to the bifurcation parameter of the surface family. The real advantage of this new approach is that it computes, in one step, the scattered field for all possible members of the family of surfaces. By contrast, other state-of-the-art schemes must restart when the returns from a new surface are desired, so that the cost of our new approach is greatly advantaged when the number of samples of the family reaches even modest values. Numerical results which verify the accuracy of our approach and demonstrate their utility in computing grating efficiencies scattered by traveling surface ocean waves are presented.     

新型上风格式(AUSM+)的性能分析及其应用

研究了新型上风格式AUSM+的分辨率、效率等性能,并用它与Roe、vanLeer上风格式数值模拟了前向台阶激波反射流动,通过对激波、膨胀波、接触间断及其间相互干扰的复杂波系的模拟对比,分析探讨了AUSM+格式的低数值耗散、间断高分辨率等特性。     

Construction of low dissipative high-order well-balanced filter schemes for non-equilibrium flows

The goal of this paper is to generalize the well-balanced approach for non-equilibrium flow studied by Wang et al. (2009) [29] to a class of low dissipative high-order shock-capturing filter schemes and to explore more advantages of well-balanced schemes in reacting flows. More general 1D and 2D reacting flow models and new examples of shock turbulence interactions are provided to demonstrate the advantage of well-balanced schemes. The class of filter schemes developed by Yee et al. (1999) [33], Sjögreen and Yee (2004) [27] and Yee and Sjögreen (2007) [38] consist of two steps, a full time step of spatially high-order non-dissipative base scheme and an adaptive non-linear filter containing shock-capturing dissipation. A good property of the filter scheme is that the base scheme and the filter are stand-alone modules in designing. Therefore, the idea of designing a well-balanced filter scheme is straightforward, i.e. choosing a well-balanced base scheme with a well-balanced filter (both with high-order accuracy). A typical class of these schemes shown in this paper is the high-order central difference schemes/predictor–corrector (PC) schemes with a high-order well-balanced WENO filter. The new filter scheme with the well-balanced property will gather the features of both filter methods and well-balanced properties: it can preserve certain steady-state solutions exactly; it is able to capture small perturbations, e.g. turbulence fluctuations; and it adaptively controls numerical dissipation. Thus it shows high accuracy, efficiency and stability in shock/turbulence interactions. Numerical examples containing 1D and 2D smooth problems, 1D stationary contact discontinuity problem and 1D turbulence/shock interactions are included to verify the improved accuracy, in addition to the well-balanced behavior.     

On large time step TVD scheme for hyperbolic conservation laws and its efficiency evaluation

A large time step (LTS) TVD scheme originally proposed by Harten is modified and further developed in the present paper and applied to Euler equations in multidimensional problems. By firstly revealing the drawbacks of Harten’s original LTS TVD scheme, and reasoning the occurrence of the spurious oscillations, a modified formulation of its characteristic transformation is proposed and a high resolution, strongly robust LTS TVD scheme is formulated. The modified scheme is proven to be capable of taking larger number of time steps than the original one. Following the modified strategy, the LTS TVD schemes for Yee’s upwind TVD scheme and Yee–Roe–Davis’s symmetric TVD scheme are constructed. The family of the LTS schemes is then extended to multidimensional by time splitting procedure, and the associated boundary condition treatment suitable for the LTS scheme is also imposed. The numerical experiments on Sod’s shock tube problem, inviscid flows over NACA0012 airfoil and ONERA M6 wing are performed to validate the developed schemes. Computational efficiencies for the respective schemes under different CFL numbers are also evaluated and compared. The results reveal that the improvement is sizable as compared to the respective single time step schemes, especially for the CFL number ranging from 1.0 to 4.0.