基于交替隐式有限差分法的快速早期乳腺癌时域微波断层成像

采用时域微波断层成像技术进行早期乳腺癌检测能够准确地获得乳房的电参数分布,具有明确的物理解释和医学诊断价值.临床应用讲究即时性,为了提高检测的速度,本文将交替隐式有限差分法应用到乳腺癌检测中,基于正反演时间步进成像算法进行成像分析,结果显示在保证精度的前提下,采用交替隐式有限差分法的成像时间可缩短为传统时域有限差分法的23%,提高了微波断层成像技术的临床可应用性.     

侧入式导光板网点全自动优化设计研究

为解决目前侧入式导光板网点优化设计中存在经验式手动优化繁琐、且很难达到一个高均匀效果的问题,提出了一种基于模糊优化理论的网点自动优化设计方法。网点形状采用锥形结构,为使目标面上各处的光能分布主要由其正下方网点的散射光贡献,分析了其半顶角与网点位置的关系,并使得网点半顶角仅由其位置决定。以网点半径为主要优化参数,通过隶属度函数将网点结构模糊化,并自定义评价函数来解模糊化实现均匀照明。进而,采用动态数据交换(DDE)技术将Matlab与TracePro进行联立,通过Matlab语言与Scheme语言混编控制TracePro自动进行数据交换、光线追迹与模糊优化。导光板网点优化设计实例仿真结果表明,优化后均匀度达到92.17%,光能利用率达到70.37%,全程实现自动优化,无需任何手动调节。     

Shape reconstruction of metallic objects from intensity scattered field data only

We address an inverse scattering problem using only amplitude information of the scattered field. In particular, we are concerned with the reconstruction of the shape of a metallic planar obstacle, located in a known plane (aperture plane), starting from the knowledge of the intensity of the scattered field over two measurement planes in near zone and parallel to the aperture plane, when a single-frequency incident plane wave is exploited. The formulation of the inverse scattering problem is given under the physical optics approximation, and thus the resulting phase retrieval problem is quadratic. This allows us to apply some phase retrieval techniques already developed for antenna diagnostics. Reconstruction results with synthetic data indicate the feasibility of the proposed approach.     

Truss optimization on shape and sizing with frequency constraints based on orthogonal multi-gravitational search algorithm

Structural optimization on shape and sizing with frequency constraints is well-known as a highly nonlinear dynamic optimization problem with several local optimum solutions. Hence, efficient optimization algorithms should be utilized to solve this problem. In this study, orthogonal multi-gravitational search algorithm (OMGSA) as a meta-heuristic algorithm is introduced to solve truss optimization on shape and sizing with frequency constraints. The OMGSA is a hybrid approach based on a combination of multi-gravitational search algorithm (multi-GSA) and an orthogonal crossover (OC). In multi-GSA, the population is split into several sub-populations. Then, each sub-population is independently evaluated by an improved gravitational search algorithm (IGSA). Furthermore, the OC is used in the proposed OMGSA in order to find and exploit the global solution in the search space. The capability of OMGSA is demonstrated through six benchmark examples. Numerical results show that the proposed OMGSA outperform the other optimization techniques.     

一种强噪声背景下微弱超声信号提取方法研究

为解决在强噪声背景下获取超声信号的难题,基于粒子群优化算法和稀疏分解理论提出一种强噪声背景下微弱超声信号提取方法.该方法将降噪问题转换为在无穷大参数集上对函数进行优化的问题,首先以稀疏分解理论和超声信号的结构特点为依据构建了粒子群优化算法运行所需要的目标函数及去噪后信号的重构函数,从而将粒子群优化算法和超声信号降噪联系在一起;然后根据粒子群优化算法可以在连续参数空间寻优的特点建立了用于匹配超声信号的连续超完备字典,并采用改进的自适应粒子群优化算法在该字典中对目标函数进行优化;最后根据对目标函数在字典上的优化结果确定最优原子,并利用最优原子按照重构函数重构出降噪后的超声信号.通过对仿真超声信号和实测超声信号的处理,结果表明本文提出的方法可以有效提取信噪比低至-4 dB的强噪声背景下的微弱超声信号,且和基于自适应阈值的小波方法相比本文方法表现出更好的降噪性能.     

Electric double layer of anisotropic dielectric colloids under electric fields

Anisotropic colloidal particles constitute an important class of building blocks for self-assembly directed by electrical fields. The aggregation of these building blocks is driven by induced dipole moments, which arise from an interplay between dielectric effects and the electric double layer. For particles that are anisotropic in shape, charge distribution, and dielectric properties, calculation of the electric double layer requires coupling of the ionic dynamics to a Poisson solver. We apply recently proposed methods to solve this problem for experimentally employed colloids in static and time-dependent electric fields. This allows us to predict the effects of field strength and frequency on the colloidal properties.     

一种基于压缩感知的三维导体目标电磁散射问题的快速求解方法

提出了一种将压缩感知和特征基函数结合的方法来计算三维导体目标的雷达散射截面.利用压缩感知理论,将随机选择的矩量法阻抗矩阵作为测量矩阵,将激励电压视为测量值,然后再用恢复算法可实现二维或二维半目标感应电流的求解.对于三维导体目标,使用Rao-Wilton-Glisson基函数表示的感应电流在常用的离散余弦变换基、小波基等稀疏基上不稀疏.为此,本文将计算出的目标特征基函数作为稀疏基,用广义正交匹配追踪算法作为恢复算法来加速恢复过程,并应用到三维导体目标的雷达散射截面计算中.数值结果证明了本文方法的准确性与高效性.     

Numerical Analysis of Electromagnetic Fields in Multiscale Model

Modeling technique for electromagnetic fields excited by antennas is an important topic in computational electromagnetics, which is concerned with the numerical solution of Maxwell's equations. In this paper, a novel hybrid technique that combines method of moments(MoM) with finite-difference time-domain(FDTD) method is presented to handle the problem. This approach employed Huygen's principle to realize the hybridization of the two classical numerical algorithms. For wideband electromagnetic data, the interpolation scheme is used in the MoM based on the dyadic Green's function. On the other hand, with the help of equivalence principle, the scattered electric and magnetic fields on the Huygen's surface calculated by MoM are taken as the sources for FDTD. Therefore, the electromagnetic fields in the environment can be obtained by employing finite-difference time-domain method. Finally, numerical results show the validity of the proposed technique by analyzing two canonical samples.     

An overlooked electrostatic force that acts on a non-charged asymmetric conductor in a symmetric (parallel) electric field

Usually, when a material that has charge Q is placed in an electric field E, an electrostatic force F = QE acts on the material. This force does not act on a non-charged material. Nevertheless, when a non-charged material is placed in a convergent field, another electrostatic force acts. This force is called the gradient force. If the material is small and the shape is a sphere, the gradient force can be calculated by an approximate formula, but it cannot be calculated for other shapes. In this paper the gradient force that acts on a symmetric rod conductor in a convergent (asymmetric) field was simulated by an axis symmetry finite difference method.Under same simulation conditions without the next two points, the shape of the conductor and the form of the field were reversed. The shape of the conductor was changed into an asymmetric shape (e.g. bat shape), and the form of the field was changed into a symmetric (parallel) one. The electrostatic force that acts on the asymmetric conductor in the symmetric (parallel) field was simulated. It was found that approximately the same intensity force as in the first simulation also acts on this conductor. This force is thought to be an overlooked electrostatic force. I provisionally call it the asymmetric force in this paper.The asymmetric force with differently shaped conductors was simulated and it was found that the asymmetric force was maximized for a cup shaped conductor.Finally, the asymmetric force with the cup shaped conductor in normal and reversed parallel (symmetric) fields was simulated, and it was confirmed that the asymmetric force remains the same in both fields.     

Two-dimensional fracture analysis of piezoelectric material based on the scaled boundary node method

A scaled boundary node method(SBNM) is developed for two-dimensional fracture analysis of piezoelectric material,which allows the stress and electric displacement intensity factors to be calculated directly and accurately. As a boundarytype meshless method, the SBNM employs the moving Kriging(MK) interpolation technique to an approximate unknown field in the circumferential direction and therefore only a set of scattered nodes are required to discretize the boundary. As the shape functions satisfy Kronecker delta property, no special techniques are required to impose the essential boundary conditions. In the radial direction, the SBNM seeks analytical solutions by making use of analytical techniques available to solve ordinary differential equations. Numerical examples are investigated and satisfactory solutions are obtained, which validates the accuracy and simplicity of the proposed approach.     

Semi-empirical bone model for determination of trabecular structure properties from backscattered ultrasound

A novel semi-empirical scattering model of trabecular bone facilitating its characterization and allowing optimization of the interrogating pulse-echo transducer performance was developed. The model accounts for spatial density distribution of the trabeculae and includes measurement conditions such as pressure–time waveform of the probing ultrasound wave, the emitted field structure, and the transfer function and limited bandwidth of the acoustic source operating in pulse-echo mode. These measurement conditions are of importance as they modify the scattered echoes, which in turn are linked to the micro-architecture of the bone. The bone was modeled by a random distribution of long and thin cylindrical scatterers having randomly varying diameters and mechanical properties, and oriented perpendicularly to the ultrasound beam axis. To mimic clinically encountered conditions the relevant empirical data obtained at 1 MHz were input to the model. The data included pulse-echo source pressure field distribution in the focal zone and the above mentioned transfer function. With these data the model allowed frequency dependent backscattering coefficient of the simulated bone structure and its statistical properties to be determined. The results obtained indicated that the computer simulation is of particular relevance in studying scattering properties of the cancellous bone and holds promise as a tool to determine the relationship between the physical dimensions and shape of the scatterers and for monitoring of osteoporosis. The results of simulations also indicated that the new bone model proposed is well suited to mimic clinically relevant conditions. In contrast to the existing bone models, which usually assume scatterers to be randomly distributed as infinitely long identical cylinders with a cross-section much smaller than the probing ultrasound wave, the new model includes two populations of scatterers having different physical dimensions and also allows the mechanical properties of the scatterers to be varied.     

基于全变分最小化和交替方向法的康普顿散射成像重建算法

康普顿散射成像技术利用射线与物质作用后的散射光子信息对物质的电子密度进行成像.与传统的透射成像方式相比,康普顿散射成像具有系统结构灵活、成像对比度高、辐射剂量低等优势,在无损检测、医疗诊断、安全检查等领域有着广阔的应用前景.但其重建问题是一个非线性的逆问题,通常是不适定的,其解对噪声和测量误差非常敏感.为解决此问题,本文结合全变分最小化正则化方法和交替方向法提出了一种新的康普顿散射成像重建算法.该算法首先将问题对应的TV模型转化为与之等价的带约束的优化问题,然后利用增广拉格朗日乘子法将优化问题分解为两个具有解析解的子问题,并通过交替求解子问题使增广拉格朗日函数达到最小,进而得到重建的图像.在仿真实验中,通过与主流的ASD-POCS方法进行对比,证明了该算法在重建精度和重建效率方面的优势.     

Scattering from a coated nihility circular cylinder placed in chiral metamaterial

The scattering of electromagnetic plane wave from a coated nihility circular cylinder placed in chiral metamaterial is investigated. The nihility cylinder is coated with chiral or chiral-nihility metamaterial and the medium hosting the coated cylinder is also considered as chiral or chiral-nihility metamaterial. The vector wave equation is used to find out the solution for the fields scattered from the concentric chiral ciruclar cylinders immersed in chiral metamaterial. Inner chiral cylinder is reduced to nihility cylinder taking permittivity and permeability approaching to zero with chirality parameter equal to zero. Coated perfect electric conductor (PEC)/perfect magnetic conductor (PMC) are obtained by taking very large value of permittivity/permeability of the inner cylinder. Numerical results for fields reflected from coated nihility cylinder are compared with those obtained for coated PEC/PMC cylinder. Our results for the special case of dielectric coated nihility/PEC cylinders in free space are shown to be in agreement with already published results thus validating the scheme of analysis.     

基于自适应粒子群算法的多峰谱线分离方法研究

电感耦合等离子体原子发射光谱分析法（ICP-AES）已成为一种常规的元素分析方法,但在ICP-AES分析过程中,大多元素的分析谱线会受到背景或其他谱线的重叠干扰,形成的光谱干扰严重影响了谱线分析的准确性,所以在元素的分析过程中,需要通过适当的光谱干扰校正方法才能得到合适的元素分析线。根据光谱强度具有叠加性的特征,利用谱图将谱线形状表示为Voigt线型函数加和的多峰谱线叠加模型,以多峰谱线叠加模型与目标谱线的均方根误差构建多元函数作为评价函数的数学模型,设计自适应粒子群优化（APSO）算法寻找分离谱线特征参数的最优解,APSO算法在标准PSO算法的基础上,引入压缩因子同时使得种群参数惯性权重根据粒子个体适应度值自适应变化以及学习因子线性变化,在算法迭代过程中协调粒子种群内全局搜索能力和局部开发能力,保证算法有效且迅速收敛,实现多峰谱线分离,减少干扰谱线的影响从而得到更精准的元素分析线。以ICP-AES检测器返回的含Pr元素溶液特征波长为390.844 nm和汞灯特征波长为313.183 nm两条谱线的光强AD采样值作为两组实测数据,以两个Voigt线型近似函数构成的三种不同重叠程度的叠加合成曲线作为三组模拟数据,在数据曲线上分别选取50个能够包含曲线全部特征参数信息的点作为数据点,通过对上述五组目标数据点进行APSO算法处理,结果表明APSO算法得到的多峰谱线叠加模型相关参数能够较准确地拟合出相应的目标数据曲线,目标数据点与拟合曲线函数值相对误差较低,算法表明能够有效扣除谱线重叠干扰,同组目标数据经过多次算法处理,选择最小的最优适应度值相应的特征参数向量作为Voigt线型函数相关参数,以此拟合出的多峰谱线叠加模型曲线精准度越高、相对误差越小。这种算法具有良好的收敛性和适应性,可应用于ICP-AES在元素定性、定量方面的分析研究。     

Smoothing of Doppler spectra line shapes in electrophoretic analysis of biological cell suspensions

The study of cell and particle electrophoresis by laser Doppler techniques requires the application of an electric field to a conducting solution. Steady d.c. electric fields are usually not possible because of electrode surface reactions. Hence, it has been common practice to apply temporally varying electric fields, either pulses of alternating polarity or square-wave fields of constant amplitude. Oscillating particle motion leads to harmonic structure (multiples of the square-wave electric field frequency) in the scattered light spectra. This harmonic structure may be considerably smoothed by a suitable variation of the period of the applied square-wave field. The application of a pseudo-random period oscillating electric field yields line shapes that approach the line shapes to be expected from the application of a continuous d.c. electric field.     

均匀外电场中无限长介质圆柱壳的场分布的计算和讨论

用圆柱坐标系中的分离变量法计算了位于均匀外电场中的无限长介质圆柱壳各区域的电势和电场,由计算结果分析了无限长介质圆柱壳对外电场的屏蔽效果,并指出均匀外电场中的无限长导体圆柱壳、无限长介质圆柱体、无限长导体圆柱及无限大均匀电介质中开有一无限长的圆柱形空腔的电势和电场都可以由均匀电场中的无限长介质圆柱壳电势及电场给出.     

Optimized basis function for spectral reflectance recovery from tristimulus values

Spectral reflectance is a significant physical property of materials. It plays an important role in color constancy, illumination modeling, and color reproduction. Spectral reflectance basis functions are the most important impact factors for spectral reflectance recovery. Previous methods mainly calculated basis functions for the reflectance spectra data sets by employing the principle component analysis (PCA) and its improved methods. In this paper, we present a new method to solve this problem. Specifically, we propose a new cost function and some constraint conditions to convert the problem into an optimization problem by minimizing the cost function. Unlike the PCA method which yields the orthogonal basis functions for approximating the reflectance spectra, our method yields the nonorthogonal basis functions for better approximating the reflectance spectra. Experimental results show that our derived basis functions are better than those obtained by the PCA method for reflectance spectra recovery.     

Shape optimization of 3D curved slots and its application to the squirrel-cage elastic support design

The squirrel-cage elastic support is one of the most important components of an aero-engine rotor system.A proper structural design will favor the static and dynamic performances of the system.In view of the deficiency of the current shape optimization techniques,a new mapping approach is proposed to define shape design variables based on the parametric equations of 3D curves and surfaces.It is then applied for the slot shape optimization of a squirrel-cage elastic support.To this end,an automatic design procedure that integrates the Genetic Algorithm (GA) is developed to solve the problem.Two typical examples with different shape constraints are considered.Numerical results provide reasonable optimum designs for the improvement of stiffness and strength of the squirrel-cage elastic support.     

逆向傅里叶衍射定理快速预报水下目标三维声散射指向分布

提出逆向运用傅里叶衍射定理预报水下弱散射目标三维声散射指向分布的快速计算方法。依据目标形状、周围介质的密度和声速构建三维声场图像模型,建立散射远场积分结果与图像频域幅值的关系式,提取频域中半径为水中波数k;的球型表面上的幅值,获得精细化的宽带、全方位散射声压指向特性。数值计算表明:将傅里叶衍射定理逆向运用于解决声学正问题,适用于分层的、不均匀的、非规则及多体弱散射目标散射声场的求解。通过插值提取频域样本获得远场声压的方法,避免了有限元法(3D-FEM)所必须的大规模的网格划分和迭代运算,可以有效地减少计算成本并拓展散射频率响应的带宽。在水池中完成两种具有不同声学参数和形状目标指向性测试实验,得到散射声压指向性幅度函数与理论预报相一致。      

Continuous and Discrete Adjoint Approach Based on Lattice Boltzmann Method in Aerodynamic Optimization Part I: Mathematical Derivation of Adjoint Lattice Boltzmann Equations

The significance of flow optimization utilizing the lattice Boltzmann (LB) method becomes obvious regarding its advantages as a novel flow field solution method compared to the other conventional computational fluid dynamics techniques. These unique characteristics of the LB method form the main idea of its application to optimization problems. In this research, for the first time, both continuous and discrete adjoint equations were extracted based on the LB method using a general procedure with low implementation cost. The proposed approach could be performed similarly for any optimization problem with the corresponding cost function and design variables vector. Moreover, this approach was not limited to flow fields and could be employed for steady as well as unsteady flows. Initially, the continuous and discrete adjoint LB equations and the cost function gradient vector were derived mathematically in detail using the continuous and discrete LB equations in space and time, respectively. Meanwhile, new adjoint concepts in lattice space were introduced. Finally, the analytical evaluation of the adjoint distribution functions and the cost function gradients was carried out.