一种基于传输线方程的埋地电缆电磁脉冲耦合时域分析方法

 埋地电缆和其它金属长导体对核电磁脉冲的耦合分析对于确定合理的工程防护措施具有重要意义。采用Agrawal模型和时域有限差分算法对线缆耦合问题进行了分析，针对该模型中需要计算土壤阻抗与电流项卷积、土壤导纳与电压项卷积的问题，提出了一种土壤阻抗特性的线性拟合处理方法，可避免复杂的卷积计算。介绍了该方法的基本原理和试用结果。这一模型拟合法适用于确定无时域解析表达式的传输线模型参数，与传输线模型的时域有限差分方法相结合，可用于分析非均匀场作用下传输线的耦合。     

亚稳态氦原子束在激光驻波场中聚焦的衍射像差研究

用量子力学方法研究了一维激光驻波场对亚稳态氮原子束聚焦的衍射像差。建立了半经典模型，该模型在一定条件下和已有模型一致，提出了直接计算时间演化算符的数值方法，这种计算方法速度较快，收敛性较好，数值模拟结果表明在一定的激光强度和光束半径等参量下，亚稳态氮原子束被聚焦，焦点的半峰全宽为54nm，对比度16：1，焦长、焦深与入射原子束的速度有关，从而与动能有关，焦点的有限宽度主要是来自衍射效应的影响。     

螺旋型电磁超声辐射声场对缺陷检测影响的 计算机模拟*

电磁超声(EMAT)作为一种非接触超声检测技术，已得到了普遍的关注和研究。然而，EMAT目前只得到了一些有限的应用，尤其在缺陷探伤应用方面还不能替代压电超声，深入探究其原因并加以改进，对推动EMAT在缺陷检测方面的应用具有重要意义。利用有限元方法建立了螺旋型EMAT模型，分析了其辐射声场和声场特征；通过将有限元计算所得到的洛伦兹力分布简化为门函数，并结合格林函数法，得到了螺旋型EMAT辐射声场的简化解析结果。这些方法可用于计算螺旋型EMAT辐射声场指向性以及在缺陷检测中比较关心的一些特征参数，如扩散角、偏离角等参数。通过螺旋型EMAT辐射声场实验，验证了理论分析结果。通过与压电超声检测理论与实验的对比，研究认为：螺旋型EMAT辐射声场的具有横波和纵波等多种波模及多种模态转换的特点，以及相对复杂的异型声场几何分布等复杂特性，使其在缺陷检测应用方面受到局限。这些研究方法对改进EMAT设计，优化控制EMAT辐射声场，拓展EMAT在缺陷检测方面的应用具有参考价值。     

用有限电流元法识别HL-2A 装置的等离子体边界

在电流丝模型对HL-2A装置等离子体边界重建取得满意结果[3]、[5]的基础上，本文用有限电流元模型对HL-2A等离子体边界进行了重建研究。计算和实验结果表明，在通常情况下，有限电流元模型比固定电流丝模型的重建误差稍小，前者误差小于3mm，后者误差小于6mm。当部分有限电流元分布在等离子体边界之外时，用有限电流元模型仍然可以成功重建边界。有限电流元的位置在一定范围内变化时，重建的误差都很小。用普通奔腾4 2.4GHz PC机计算一组等离子体放电边界的时间不超过1ms。有限电流元法能准确而快速地识别等离子体的偏滤器位形，这对于HL-2A装置的实时位形控制是基本和重要的。     

脉冲强激光冲击力学效应的数值模拟

利用一维理想弹塑性流体动力学模型和有限差分方法，研究了激光辐照固态靶材时产生的形成和衰减过程以及动态断裂。     

基于高光谱图像信息融合的红提糖度无损检测

红提糖度是重要的内部品质衡量指标,传统的检测方法均为破坏性生化检测,本文基于高光谱成像技术,提出了一种基于高光谱信息融合的红提糖度含量无损检测方法。采集并提取260个红提样本的光谱信息和图像信息,对光谱信息分别利用SNV、S-G等光谱预处理方法建立PLSR模型,确定最好的光谱预处理方法,分别采用一次降维(GA、CARS、IRIV)算法和组合降维算法(CARS-SPA、IRIV-SPA、GA-SPA)共六种降维方法对光谱信息进行特征变量提取;获取灰度共生矩阵的图像纹理信息,结合图像的颜色信息(R、G、B、H、S、V、L、a、b),组成19个图像特征参数,采用PCA算法对图像信息进行降维,分别建立基于降维处理后的光谱信息、图像信息以及两者融合的红提糖度线性预测模型PLSR、非线性预测模型LSSVM,并对比分析模型的优劣。结果表明,若只利用光谱信息建模,IRIV-SPA可有效地提取红提糖度光谱信息的特征波长,提高模型的预测性能;若只利用图像信息进行建模,模型的预测性能不好,PCA降维有效地提高了模型的预测性能,但提高的性能有限;将IRIV-SPA特征波段提取后的光谱和经PCA降维后的图像信息进行融合,分别建立PLSR和LSSVM模型,红提糖度的最优PLSR模型的校正集和预测集相关系数分别为0.943,0.941;红提糖度的最优LSSVM模型的校正集和预测集相关系数分别为0.954,0.952。LSSVM所建模型的效果好于PLSR所建模型,但模型的运算时间较长。两种模型的精度均比单方面基于光谱或图像信息的模型都有较大的提高,表明融合高光谱图像的光谱与图像信息不仅可以提高模型的运算速度、简化模型,同时有效地提高了红提糖度预测性能,为红提糖度的检测找到了一种新的方法。     

光子图像统计处理

讨论了一种对有限光子图像累加的统计处理方法。在超微弱发光的研究中（例如细胞的超微弱光），由于发光强度较弱，需要用像增强器对超微弱发光图像进行增强得到可视光子图像，超微弱发光的光子图像不可避免地受到像增强系统暗噪声及背景噪声的影响，使光子图像湮没在噪声中。为从原始图像中检验出信号，常将采集到的有限光子图像进行累加以增加样本数，根据信号光子和噪声光子的不同统计分布，运用似然比检验方法得到有效判断累加图像中某像素是否含有信号光子，由此剔除图像中的噪声光子，得到信噪比改善的光子图像。     

F－L模型的状态方程与退禁闭相变图象

用温度场论方法计算了Ｆ－Ｌ模型在有限温度和密度下的状态方程，分析了压强对净重子数密度的等温线．结果表明，在平均场近似下，Ｆ－Ｌ模型所给出的退禁闭相变为一级相变．     

超宽带微波检测早期乳腺肿瘤三维仿真

超宽带微波成像技术可作为有效的方法用于乳房早期乳腺肿瘤的检测. 该方法基于乳房组织和肿瘤之间较大的电学特性差异的特点进行成像, 能提供足够的分辨率以及足够的穿透深度. 本文采用时域有限差分方法建立超宽带微波信号在三维乳房组织中传播的模型, 并采用单极德拜模型完成了生物组织色散特性的模拟. 利用8发9收的天线阵列实现微波肿瘤探测, 利用共焦成像算法对乳房组织进行图像重构并进行肿瘤定位. 仿真结果显示共焦成像算法能够实现最小直径3 mm的肿瘤的检测, 同时证实了超宽带微波成像技术应用于早期乳腺肿瘤检测的有效性. 关键词： 微波成像 肿瘤检测 时域有限差分 共焦成像算法     

惯性约束聚变实验中空间干涉检测与预警系统设计

针对惯性约束聚变实验中多诊断系统在有限空间同时作业引发碰撞风险的问题，设计了基于碰撞目标检测和空间距离测算方法的空间干涉检测与预警系统，分析了真空靶室特定环境下不同实现方法的适用条件和可行性。实验室条件下模拟验证结果表明，基于结构光深度图像获取方法可获得较好的碰撞距离测量精度，在靶室半径探测距离内空间测距误差不大于3cm，为靶室空间干涉检测与预警提供了一种可行的解决方案。     

Free Mesh Method: fundamental conception, algorithms and accuracy study

The finite element method (FEM) has been commonly employed in a variety of fields as a computer simulation method to solve such problems as solid, fluid, electro-magnetic phenomena and so on. However, creation of a quality mesh for the problem domain is a prerequisite when using FEM, which becomes a major part of the cost of a simulation. It is natural that the concept of meshless method has evolved. The free mesh method (FMM) is among the typical meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, especially on parallel processors. FMM is an efficient node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm for the finite element calculations. In this paper, FMM and its variation are reviewed focusing on their fundamental conception, algorithms and accuracy.     

介质光波导的数值计算及近似方法

综述介质光波导中模式分析方法的最新进展。着重分析严格的数值分析方法及近似方法最近的研究进展。依据方法的计算精确度、计算效率和计算复杂度等因素．对有限元法和有限差分法这两种数值分析方法以及有效指数法和频谱指数法等近似方法进行讨论。论述磁光波导的分析方法。     

A Fast Finite Element Analysis of Millimeter Wave Waveguide Filters

The multifrontal method is applied for solving a large system of linear equations resulting from the use of edge-based finite element method (FEM). The finite element method combined with the perfectly matched layers (PML) is given for simulation of microwave filters and the algorithm of multifrontal method is described in detail. The reflection and insertion losses of millimeter wave filter are analyzed as the examples and the obtained results are compared with those obtained from literature. In order to demonstrate the efficiency of the multifrontal method, the computational time is compared with that of both symmetric successive overrelaxation (SSOR) preconditioned conjugate gradient (PCG) and conjugate gradient methods (CG) for the thick-iris waveguide bandpass filer.     

An element-free Galerkin(EFG) method for numerical solution of the coupled Schrdinger-KdV equations

The present paper deals with the numerical solution of the coupled Schrdinger-KdV equations using the elementfree Galerkin(EFG) method which is based on the moving least-square approximation.Instead of traditional mesh oriented methods such as the finite difference method(FDM) and the finite element method(FEM),this method needs only scattered nodes in the domain.For this scheme,a variational method is used to obtain discrete equations and the essential boundary conditions are enforced by the penalty method.In numerical experiments,the results are presented and compared with the findings of the finite element method,the radial basis functions method,and an analytical solution to confirm the good accuracy of the presented scheme.     

Numerical simulation of continuum models for fluid-fluid interface dynamics

This paper is concerned with numerical methods for two-phase incompressible flows assuming a sharp interface model for interfacial stresses. Standard continuum models for the fluid dynamics in the bulk phases, for mass transport of a solute between the phases and for surfactant transport on the interface are given. We review some recently developed finite element methods for the appropriate discretization of such models, e.?g., a pressure extended finite element (XFE) space which is suitable to represent the pressure jump, a space-time extended finite element discretization for the mass transport equation of a solute and a surface finite element method (SurFEM) for surfactant transport. Numerical experiments based on level set interface capturing and adaptive multilevel finite element discretization are presented for rising droplets with a clean interface model and a spherical droplet in a Poisseuille flow with a Boussinesq-Scriven interface model.     

求解辐射传递的非结构混合有限体积／有限元法

本文给了一种适用于任意非结构网格的有限体积/有限元法的混合算法用于求解多维半透明吸收、发射、散射性灰矩形介质内的辐射传递.该方法使用有限元法进行角度离散,有限体积法进行空间离散.与基于辐射传递离散坐标方程的方法不同的是,该方法在迭代求解的过程中,针对每一个空间体元,所有角度方向的辐射强度同时耦合求出.通过两个算例验证了该解法的正确性.     

相变温度场中无网格伽辽金法的应用

简述了无单元法的基础理论,推导出相变温度场的无单元法计算公式,采用罚函数法引入了第一类边界条件,编制了相应的计算程序.通过经典相变的应用例子,和有限元计算结果及解析解的比较,说明了无单元法应用于相变温度场具有连续性好,精度高,前后处理简单等优点.     

The finite element method for micro-scale modeling of ultrasound propagation in cancellous bone

Quantitative ultrasound for bone assessment is based on the correlations between ultrasonic parameters and the properties (mechanical and physical) of cancellous bone. To elucidate the correlations, understanding the physics of ultrasound in cancellous bone is demanded. Micro-scale modeling of ultrasound propagation in cancellous bone using the finite-difference time-domain (FDTD) method has been so far utilized as one of the approaches in this regard. However, the FDTD method accompanies two disadvantages: staircase sampling of cancellous bone by finite difference grids leads to generation of wave artifacts at the solid–fluid interface inside the bone; additionally, this method cannot explicitly satisfy the needed perfect-slip conditions at the interface. To overcome these disadvantages, the finite element method (FEM) is proposed in this study. Three-dimensional finite element models of six water-saturated cancellous bone samples with different bone volume were created. The values of speed of sound (SOS) and broadband ultrasound attenuation (BUA) were calculated through the finite element simulations of ultrasound propagation in each sample. Comparing the results with other experimental and simulation studies demonstrated the capabilities of the FEM for micro-scale modeling of ultrasound in water-saturated cancellous bone.     

Simulation of Millimeter Wave Circulators by FEM with Unifrontal/Multifrontal Technique

The combination of unifrontal and multifrontal techniques is applied for solving a large, sparse, and unsymmetric system of linear equations resulting from the use of edge-based finite element method (FEM) for millimeter wave ferrite devices. The formulation of finite element method combined with the perfectly matched layers (PML) is given and the combined algorithm of unifrontal and multifrontal methods is described. The lower data movement of unifrontal method and the lower fill-in of normal multifrontal methods are combined for an algorithm with their features. The performance of typical waveguide junction circulators is analyzed and the calculated results are compared with those obtained from the literature.     

Computational performance of Free Mesh Method applied to continuum mechanics problems

The free mesh method (FMM) is a kind of the meshless methods intended for particle-like finite element analysis of problems that are difficult to handle using global mesh generation, or a node-based finite element method that employs a local mesh generation technique and a node-by-node algorithm. The aim of the present paper is to review some unique numerical solutions of fluid and solid mechanics by employing FMM as well as the Enriched Free Mesh Method (EFMM), which is a new version of FMM, including compressible flow and sounding mechanism in air-reed instruments as applications to fluid mechanics, and automatic remeshing for slow crack growth, dynamic behavior of solid as well as large-scale Eigen-frequency of engine block as applications to solid mechanics.