Acoustic radiation force analysis using finite difference time domain method

Acoustic radiation force exerted by standing waves on particles is analyzed using a finite difference time domain Lagrangian method. This method allows the acoustic radiation force to be obtained directly from the solution of nonlinear fluid equations, without any assumptions on size or geometry of the particles, boundary conditions, or acoustic field amplitude. The model converges to analytical results in the limit of small particle radii and low field amplitudes, where assumptions within the analytical models apply. Good agreement with analytical and numerical models based on solutions of linear scattering problems is observed for compressible particles, whereas some disagreement is detected when the compressibility of the particles decreases.     

Antireflection subwavelength structures analyzed by using the finite difference time domain method

Antireflection subwavelength structures (ASSs) are analyzed by using the finite difference time domain (FDTD) method in the visible light spectrum. Low reflectance can be obtained by both the conical and pyramidal shapes over a broadband range. Comparing the reflectance of different structure shapes and aspect ratios by the FDTD method, it shows that the antireflection efficiencies of the pyramidal structures are better than that of the conical structures when the aspect ratio is up to 0.8. It is found that, for the conical structure surface, the average transmittance increases gradually with the aspect ratio and the average transmittance is about 99.6% with the aspect ratio of 2.0. However, for the pyramidal structure with the aspect ratio ranging from 1.0 to 2.0, the average transmittance is up to 99.7%.     

嵌入式薄片模型在时域有限差分算法中的应用

将一等效薄片模型嵌入到时域有限差分算法（FDTD）中,以快速而有效地解决复合材料薄片在电磁计算中的多尺度问题。在该嵌入式薄片模型中,薄片不需要被剖分网格,而是被嵌入到相邻的网格间,从而可以使用相对较大的网格剖分周围物体,进而可节省大量的计算资源。在这一模型中,薄片被等效为一段传输线,并用其频域的导纳矩阵代替。使用数字滤波器理论以及逆Z变换可将频域的导纳矩阵转换到时域,并将其嵌入到时域有限差分算法中。该嵌入式薄片模型被用来计算一单层碳纤维复合材料薄片的反射以及透射性能,并与其解析解进行对比,从而验证该模型的准确性、收敛性以及高效性。该模型被用来计算三种具有不同电参数的单层碳纤维复合材料薄片的屏蔽性能,以研究各电参数对其屏蔽性能的影响。     

时域有限差分法的混凝土超声波损伤分析

混凝土质量的好坏直接影响到建筑物的安全性,选用方便实用的方法对其内部裂纹进行检测具有显著意义。该文首先借助时域有限差分法,建立了混凝土声场数学模型。之后对不同内部缺陷尺寸的混凝土模型进行了超声波数值模拟,分别从时域和频域对接收到的超声信号进行了分析。结果发现,当内部存在缺陷时,接收信号时域幅值明显变小,接收信号的幅值与缺陷尺寸大小呈指数关系,且接收信号传播时间与缺陷宽度呈线性关系。最后以初次到达信号的频域能量为特征指标,得出其与混凝土内部开裂缺陷尺寸呈指数关系。该方法为混凝土超声无损检测建立合适的数学模型提供了参考。     

Computation of effective propagation parameters in the optical domain with a finite difference time domain method

The study of optical scattering by heterogeneous media is a complex topic where homogenization is very helpful and rigorous methods are useful. Finite difference time domain (FDTD) method coupled with Monte Carlo process is used to compute the effective parameters of heterogeneous media. Effective parameters based on the coherent field propagation of a beam in the medium are determined in bidimensional geometry and for both polarizations. It is applied to media composed of small particles embedded in an host medium, for relevant ranges of particle sizes and optical constants for both binder and particles. The results are compared to the Maxwell–Garnett and Bruggeman mixing laws and the Foldy–Twersky and Keller perturbative approximations, leading to the assessment of their validity domain.     

Application of finite difference time domain method for the design of birdcage RF head coils using multi-port excitations

A three-dimensional finite difference time domain model was developed where the high pass birdcage coil and the imaged object are analysed as a single unit. A study was performed comparing linear, conventional quadrature, and four-port excitation at 64 MHz and 200 MHz for different coil loadings, namely muscle phantoms and an anatomically detailed human head model. A phase array concept was utilized to excite the birdcage coil in four ports. Two phase conditions were analyzed, the simple fixed phase and the variable phase. At 200 MHz, compared to the conventional quadrature drive, the four-port drive reduces the effects of the tissue-coil interactions leading to more uniform currents on the coil legs and consequently to a better B(1) field homogeneity. Also at 200 MHz, driving the coil in four ports provides an SAR distribution with peak values that are significantly less than those with linear or quadrature excitations.     

A general finite difference time domain method for hybrid dispersive media model

The shift operator finite difference time domain (SO-FDTD) method is developed to deal with a hybrid dispersive media model. First, we prove that the complex permittivity of the hybrid dispersive media model can be described by a rational polynomial fraction with respect to jω. Then, the relation between electric displacement D and electric field strength E is derived in the time domain by introducing z_t as a shift operator. The constitutive relation in the discretized time domain and the recursive formulation of D and E available for FDTD computation are obtained. Finally, the reflection of the hybrid dispersive slab is computed. The computed results are in good agreement with that obtained by analytic method. This illustrates the generalization and the feasibility of the presented scheme.     

Subwavelength structured surfaces with a broadband antireflection function analyzed by using a finite difference time domain method

Subwavelength structured surfaces with a broadband antireflection function are analyzed by using a finite difference time domain method. The 2D periodic pure conical and pyramidal structures and several composite structured surfaces composing of conical and pyramidal structures are presented and discussed. It is found that the antireflection function will decrease for the same size of structures when the structure arrangement is discontinuous and the results of the structures with an aspect ratio of 2.0 are obviously better than the 1.0. It is also found that some composite structured surfaces have powerful antireflection function in the spectral ranges of 400-1400 nm within the 800-1200 nm spectral ranges.     

用传输线和时域有限差分法计算电缆X射线响应

 介绍了时域传输线和时域有限差分两种计算电缆X射线瞬态响应的方法,并利用电缆芯线总电荷以及介质层的电荷分布,结合实验结果,分别给出了时域传输线模型中的电流源和时域有限差分模型中电流密度的计算方法。作为例子,在60 keV单能X射线垂直于电缆轴向辐照的条件下,采用两种方法计算了电缆两端接匹配负载时负载上的电流响应;采用传输线方法计算了不同能量X射线辐照电缆时的芯线负载电流响应。计算结果表明,两种方法计算的电缆X射线瞬态响应是一致的;电缆介质层中存在静电场;X射线能量不同时,电缆芯线负载电流存在最大值。计算中,电缆与X射线的作用采用蒙特卡罗方法计算。     

Design of the sample cell in near-field surface-enhanced Raman scattering by finite difference time domain method

The rational design of the sample cell may improve the sensitivity of surface-enhanced Raman scattering(SERS) detection in a high degree. Finite difference time domain (FDTD) simulations of the configurationof Ag film-Ag particles illuminated by plane wave and evanescent wave are performed to provide physicalinsight for design of the sample cell. Numerical solutions indicate that the sample cell can provide more"hot spots" and the massive field intensity enhancement occurs in these "hot spots". More information onthe nanometer character of the sample can be got because of gradient-field Raman (GFR) of evanescentwave.     

Solving the Schrödinger equation for a charged particle in a magnetic field using the finite difference time domain method

We extend our finite difference time domain method for numerical solution of the Schrödinger equation to cases where eigenfunctions are complex-valued. Illustrative numerical results for an electron in two dimensions, subject to a confining potential V(x,y), in a constant perpendicular magnetic field demonstrate the accuracy of the method.     

一种处理色散介质问题的通用时域有限差分方法

色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系.这就给用时域有限差分方法计算色散介质中波的散射和传播带来了困难.现有算法往往要针对不同色散介质模型推导相应的递推公式,算法的通用性较差.本文完善和发展了移位算子-时域有限差分方法,使之成为一种处理色散介质电磁问题的通用方法.首先,证明了常见的三种色散介质模型(德拜模型、洛伦兹模型和德鲁模型)的介电系数均可以写成适于移位算子法计算的有理分式函数形式.然后,用/t代替jω,过渡到时域,再引入时域移位算子z_t代替时间微分算子来处理有理分式函数形式的介电系数,给出离散时域本构关系的表示式,进而导出时域有限差分方法当中电位移矢量和电场强度之间的关系.最后,计算了几种色散介质的电磁散射,数值结果表明了本文方法和程序的通用性和正确有效性. 关键词： 时域有限差分方法 色散介质 移位算子     

室内声学时域有限差分模拟中的边界条件

给出了时域有限差分法用于室内声学问题模拟中的边界条件,结合声波方程的基本差分格式,模拟并分析了高斯脉冲在4m×4m房间中的波动过程和脉冲响应;模拟了一9m×6m×4m房间的简正频率,并与经典理论计算值进行了对比;模拟了一12m×5m×4m水平地面房间中的坐席吸声低谷效应,并与Joe LoVtri的模拟结果进行了对比;模拟并实际测量了一10.6m×5.8m×3.4m房间在几个受声点的脉冲响应和早期衰变时间EDT,将模拟结果与实际测量结果进行了对比分析,计算程序是用Metlab语言编写的。模拟与经典理论、相关研究、实际测量几方面的对比分析,验证了本边界条件的可靠性。     

Image recovery by backward iteration of the non-standard finite difference time domain algorithm

The time-reversed finite difference time domain (FDTD) algorithm is a useful tool for image recovery. By using a negative time step, a scattered wave field can propagated backwards to extract information about the size and shape of a scatterer, but errors, in both the FDTD algorithm and noise in scattered field data, reduce the quality of the recovered image. Using a new high accuracy FDTD algorithm derived from a nonstandard (NS) finite difference model, we can greatly reduce FDTD error. Using NS-FDTD we have investigated image recovery as a function of noise in numerical experiments.     

高精度有限差分法模拟Kelvin-Helmholtz不稳定性

WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。     

高精度有限差分法模拟Kelvin-Helmholtz不稳定性

 WENO有限差分格式有较高的分辨精度,适合复杂流场的计算,在国际上被广泛采用。本文利用WENO有限差分格式求解2维守恒型欧拉方程,实现了对无粘流体中Kelvin-Helmholtz不稳定性的数值模拟。速度剪切方向采用周期边界条件;扰动增长方向采用嵌边出流边界条件,一个不稳定波长分布64个网格。数值模拟给出的扰动幅值线性增长率与线性稳定性分析给出的结果很好符合,显示了该格式的有效性和精度。数值模拟给出了清晰的密度等值线,表明该方法还具有较好的界面变形捕捉能力。     

雪层覆盖的粗糙地面与上方矩形截面柱复合电磁散射的时域有限差分法研究

采用土壤介电常数的四成分模型和雪的介电常数模型分别表示实际的地面和雪层的介电特性, 应用指数型分布粗糙面模型和Monte Carlo方法模拟实际的粗糙地面, 运用时域有限差分方法研究了雪层覆盖的粗糙地面与上方矩形截面柱的复合电磁散射问题. 得出了复合电磁散射系数的角分布曲线, 计算了双站复合散射系数随土壤与雪层粗糙度参数、介电参数、 矩形截面柱几何参数、介电参数等的变化情况, 并做了详细分析与讨论. 得到了雪层覆盖的粗糙地面与上方矩形截面柱复合电磁散射特性.     

Non-destructive imaging using the time domain topological energy method

The time domain topological gradient has recently been presented as a promising technique for imaging of complex structures [N. Dominguez, et al., Wave Motion 42 (1) (2005) 31-52], showing results obtained on simulated data. The method relies on the adequate combination of two computed ultrasonic fields, one forward and one adjoint. The adjoint field, that carries the information about the defects, is analogous to a time reversal operation. The forward field plays the role of a “photographic developer”. The great interest of the method is that the time reversal operation is not done experimentally but numerically, allowing implementation of low cost systems. The expected low sensitivity of the method to the presence of measurement or structure noise had been discussed in a former work. This paper shows results obtained on experimental data in the context of non-destructive testing, showing the efficiency of the method even in complex configuration such that composite material testing.     

Free vibration analysis of stiffened plates using finite difference method

Free vibration characteristics of rectangular stiffened plates having a single stiffener have been examined by using the finite difference method. A variational technique has been used to minimize the total energy of the stiffened plate and the derivatives appearing in the energy functional are replaced by finite difference equations. The energy functional is minimized with respect to discretized displacement components and natural frequencies and mode shapes of the stiffened plate have been determined as the solutions of a linear algebraic eigenvalue problem. The analysis takes into consideration inplane deformation of the plate and the stiffener and the effect of inplane inertia on the natural frequencies and mode shapes. The effect of the ratio of stiffener depth to plate thickness on the natural frequencies of the stiffened plate has also been examined.