Finite difference time domain method forward simulation of complex geoelectricity ground penetrating radar model

The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.

Comparison of delay compensation methods for real-time hybrid simulation using frequency-domain evaluation index

The delay compensation method plays an essential role in maintaining the stability and achieving accurate real-time hybrid simulation results. The effectiveness of various compensation methods in different test scenarios, however, needs to be quantitatively evaluated. In this study, four compensation methods （i.e., the polynomial extrapolation, the linear acceleration extrapolation, the inverse compensation and the adaptive inverse compensation） are selected and compared experimentally using a frequency evaluation index （FEI） method. The effectiveness of the FEI method is first verified through comparison with the discrete transfer fimction approach for compensation methods assuming constant delay. Incomparable advantage is further demonstrated for the FEI method when applied to adaptive compensation methods, where the discrete transfer function approach is difficult to implement. Both numerical simulation and laboratory tests with predefined displacements are conducted using sinusoidal signals and random signals as inputs. Findings from numerical simulation and experimental results demonstrate that the FEI method is an efficient and effective approach to compare the performance of different compensation methods, especially for those requiring adaptation of compensation parameters.

Nonlinear Time-Domain Theory for the Simulation of Moored Floating Body Motion

Nonlinear wave loads can induce low-frequency and high-frequency resonance motions of a moored platform in deep water. For the analysis of the nonlinear response of an offshore platform under the action of irregular waves, the most widely used method in practice is the Cummins method, in which the second-order exciting forces in the time domain are computed by a two-term Volterra series model based on incident waves, first-order body motion response, and quadratic transfer functions(QTFs). QTFs are bichromatic waves acting on a body and are computed in the frequency domain in advance. For moving bodies, QTFs are related to the first-order body response, which is to be determined in the simulation process of body motion response but is unknown in the computation procedure of QTFs. In solving this problem, Teng and Cong(2017) proposed a method to divide the QTFs into different components,which are unrelated to the body response. With the application of the new QTF components, a modified Cummins method can be developed for the simulation of the nonlinear response of a moored floating platform. This paper presents a review of the theory.

NUMERICAL SIMULATIONS OF WAVE-INDUCED SHIP MOTIONS IN TIME DOMAIN BY A RANKINE PANEL METHOD

A time domain prediction of wave-induced ship motions by a Rankine panel method is investigated. Linear boundary conditions on free surface and mean wetted body surface are adopted, while the numerical damping method is used for the radiation conditions. The motions of two ships in regular head waves are computed by the present method. The related numerical results are compared with the experiment data and those from linear strip theory. The comparison shows satisfactory agreements for pitch and heave transfer functions.

A numerical simulation of surface wave excitation in a rectangular planar-type plasma source

The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.

Time domain simulations of radiation and diffraction by a Rankine panel method

The radiation and the diffraction of a ship with a forward speed are studied by using a time domain Rankine panel method. The free surface conditions are linearized onto an undisturbed free surface based on the double body flow. The linearized body boundary condition is applied on the mean wetted hull surface. The fluid domain boundary is discretized by a collection of quadric panels. The unknown quantities, including the free surface elevation, the normal flux over the free surface and the potential on the fluid domain boundary, are determined at each time step. The numerical results are compared with experimental data and other numerical solutions, showing satisfactory agreements.

Prediction of Non-Cavitation Propeller Noise in Time Domain

The blade frequency noise of non-cavitation propeller in a uniform flow is analyzed in time domain. The unsteady loading （dipole source） on the blade surface is calculated by a potential-based surface panel method. Then the time- dependent pressure data is used as the input for Ffowcs Williams-Hawkings formulation to predict the acoustics pressure. The integration of noise source is performed over the true blade surface rather than the nothickness blade surface, and the effect of hub can be considered. The noise characteristics of the non-cavitation propeller and the numerical discretization forms are discussed.

三维时域自由面格林函数法的时域绕射力计算研究

船舶工业CAE软件势流求解器的开发对势流计算方法的可靠性提出要求,本文针对三维时域面元法开展绕射力计算研究。首先引入三维时域自由面格林函数,其瞬时项的求解采用Hess&Smith方法,记忆项的求解采用Beck团队的方法;其次,推导出全浪向绕射力脉冲函数法编程所用的数学表达式;然后,采用分布源模型求解源强和速度势,并采用压力直接积分法求出时域绕射力;最后,针对Wigley I船型开展绕射力和FK力计算,并与公开的试验和计算结果进行对比,验证本文方法的可靠性。本文方法和代码可用于船舶工业CAE软件势流求解器和波浪中船舶非线性失稳运动预报。

时域间断伽辽金-粒子模拟方法研究进展

时域间断伽辽金-粒子模拟(discontinuous Galerkin time-domain particle-in-cell,DGTD-PIC)方法采用非结构网格能准确描述复杂结构,采用高阶基函数进一步提高了计算的精度,同时避免了求解大型逆矩阵,在等离子体和高功率电磁学相关数值模拟方面具有广泛的用途。近年来国内外报道了大量DGTD-PIC方面的研究工作,该方法在核心关键技术方面已取得了许多突破性进展。为促进DGTD-PIC方法的发展及应用,本文对该方法在修正势的DGTD电磁场计算方法、复频移完全匹配层边界、非结构网格中带电粒子定位与电荷和电流分配等关键技术方面所取得的研究成果进行了详细介绍。

基于非线性解耦方法的交直流配电系统暂态稳定性分析

随着柔性交直流配电系统的迅速发展,电力电子装置和恒功率负载等非线性元件的大规模接入使系统的暂态稳定性问题日趋严重。首先,建立了双环控制下AC/DC与DC/DC换流器并联的具有恒功率负载的交直流配电系统数学模型,并进行非线性泰勒展开。其次,根据非线性解耦理论将原始高阶系统解耦成一系列低阶子系统,确定主导解耦子系统并采用相平面法和逆轨迹法对其进行暂态稳定性分析,通过实物仿真验证理论分析的有效性。然后,基于非线性解耦原理引入非线性参与因子概念,对状态变量进行参与因子的计算,提出恒功率负载扰动下控制参数影响稳定性的分析方法,进一步给出控制参数的可行域。最后,利用PLECS软件搭建了恒功率负载下双环控制的交直流配电系统,仿真及实验验证了所提方法的有效性。

基于迁移学习的快速多边响应法

针对基于多边响应的快速时域仿真计算中,需要通过瞬态仿真得到所有边沿响应的时间花费问题,提出了基于迁移学习的多响应预测方法。首先,根据非线性链路的响应特点,将拖尾响应划为线性区和非线性区分别处理,以减少所需测量的响应长度;其次,以“0电平”响应作为源域数据构建预训练模型网络;然后,利用预训练模型的网络参数训练其他三种响应的迁移学习模型;最后,将所有非线性区响应绘制在一个单位间隔内,并将线性区的积累电压扩充至眼图轮廓。实验结果表明,该方法能以更少的响应数据,更全面且准确预测其余响应,从而提高快速多边响应法的计算速率。

核级截止阀抗冲击性能数值分析

为确保对核级截止阀的抗冲击性能进行全面而准确的评估,以核级截止阀为对象,分别采用动力学设计法(DDAM)和时域模拟法在频域和时域上对核级截止阀进行抗冲击计算分析。首先计算截止阀分别在纵向、垂向和横向三个方向上的频域和时域冲击载荷,然后对截止阀的抗冲击特性进行仿真计算,得到截止阀三个方向上的位移及应力分布。两种方法计算结果表明,截止阀的阀体两侧拐角部位是截止阀在冲击载荷下的薄弱部分。同时,截止阀对冲击载荷的响应与冲击加载方向有很大关系,截止阀响应的特性都以横向响应为主,要对横向抗冲击设计重点考虑。

道路路基脱空病害雷达波场数值模拟与瞬时属性分析

针对道路路基脱空病害的快速识别,采用时间域有限差分法进行病害的雷达波长数值模拟,模拟了脱空病害不同大小、长度、位置及填充的差异,并通过图像增益和复信号分析技术提取雷达波场的振幅、相位和频率等瞬时属性特征,分析了不同脱空地质的雷达波场响应特征。研究结果表明:根据路基脱空病害的地质雷达瞬时振幅图像和单道波归一化幅值可推算出脱空病害上界面信息和脱空区域介质类型;根据充气脱空病害和充水脱空病害瞬时相位图像特征,可判断出路基结构层间信息和脱空病害区域和边界位置;根据充水脱空病害和充气脱空病害瞬时频率图像,可判断出高频和低频电磁波存在位置。此外,根据充水脱空病害单道波瞬时频率图,可得到正向电磁波存在位置,并作为脱空病害下界面位置判断依据。研究成果可为道路病害的智能识别和快速判识提供理论依据。

某舰用雷达抗冲击仿真研究

抗冲击设计是舰用雷达结构设计的重要组成部分。文中针对某舰用雷达抗冲击仿真分析问题,介绍了常用的舰用设备抗冲击数值仿真分析方法,剖析了主流方法在计算效率和适用场合上的不同之处,为仿真分析方法的选用提供依据,并建立了某舰用雷达有限元模型;针对不同仿真方法产生的抗冲击分析结果差异问题,分别使用动力学设计法(DDAM)和时域分析法进行抗冲击仿真分析,分析结果表明,两种仿真计算方法得到的结构应力分布趋势基本一致,得出的薄弱结构部位相同,DDAM计算效率更高,但由于DDAM算法本身的保守性,且没有考虑低频段密集模态问题,计算得到应力值结果偏大。因此舰用雷达抗冲击仿真推荐采用时域分析法。

基于FDTD算法的表面阻抗吸收边界的研究

吸收边界条件是时域有限差分(FDTD)算法解决电磁计算问题的一种重要条件。对4种吸收边界方法:卷积完美匹配层(CPML)、表面阻抗吸收边界条件(SIBC)和2种媒介吸收体的原理进行推导,通过Matlab软件对FDTD算法进行编程计算。将4种方法用于微带低通滤波器和分支耦合器2个微带结构中,计算其反射系数。通过改变空气层厚度和迭代步长,4种方法均能得到正确的结果。对比分析4种吸收边界方法在2个模型中的计算时间、内存占用、误差百分比,可知SIBC占用内存最少,计算效率更高,误差百分比为1.63%。

巷道电性源瞬变电磁响应三维数值模拟

近年来地面电性源得到成功应用,巷道电性源瞬变电磁探测将成为热点。本文采用三维时域有限差分法,对巷道底板以及巷道与巷道之间异常体的电性源瞬变电磁响应进行了模拟。模拟结果表明:对于巷道底板异常体,受异常体表面电荷影响,电场水平分量和垂直分量在异常体表面剧烈变化,水平电场对低阻体的分辨能力较强,对高阻体的分辨能力较弱,低阻体与高阻体的垂直电场响应的符号相反;低阻体和高阻体的∂B_(y)/∂t响应在早、晚期的符号相反,低阻体响应在早期的符号为负,晚期为正,而高阻体响应在早期的符号为正,晚期为负。对于巷道与巷道之间的异常体,水平电场对低阻体的分辨能力较强,对高阻体的分辨能力较弱;高阻体模型的∂B_(z)/∂t响应曲线高于背景曲线,低阻体模型的响应曲线低于背景曲线,说明磁场∂B_(z)/∂t响应对高阻体和低阻体均有分辨能力。

基于时域高频法的宽带近场回波仿真算法

针对宽带引信弹目交会多普勒回波的快速仿真需求,提出了一种基于时域高频法的宽带近场回波快速仿真算法和仿真流程。该算法将时域物理光学(TDPO)近场积分方程和时域面积分方程相结合,利用远场近似,推导得到导体目标的时域物理光学近场散射闭合表达式,避免了计算复杂的面积分,提高了仿真速度。同时,提出了基于时域高频法的宽带近场动态多普勒回波仿真流程,经一次计算便能够得到宽带的时域回波。仿真结果表明,采用该仿真算法和仿真流程可实现宽带引信动态交会多普勒回波的快速仿真,且通过与多普勒频率理论值的对比,验证了所提仿真算法和仿真流程的正确性。

船用典型动力设备抗冲击性能评估研究

介绍了舰船设备抗冲击计算的动态设计计算法(DDAM)和时域模拟法,根据舰船动力机组各组成设备的结构特点,分别采用DDAM方法和时域加载法在频域和时域上对其进行抗冲击计算分析。计算分析结果表明,不同设备对冲击载荷的响应与设备本身的结构特点有关,且与冲击加载方向有很大关系。对于主汽轮机组和主减速齿轮装置,由于各个构件之间存在连接关系,螺栓、销钉及其连接结构的动态响应更为敏感。所得分析结论可以对开展舰船设备抗冲击计算分析以及新型舰船设备的抗冲击设计提供一定的参考。

电力系统暂态稳定计算的在线应用

暂态稳定计算中直接法计算速度快,能对暂态稳定进行定量分析,但是在模型的详细程度、准确度和可靠性方面都无法和时域仿真法相比,而时域仿真法由于考虑较复杂的模型,计算时间较长。文章首先应用快速暂态计算求取系统暂态稳定裕度,对系统故障进行排序,对筛选出的故障做考虑复杂模型的时域仿真法校验计算,从而实现了在线稳定分析。文中所论述的在线稳定计算方法在甘肃电网中得到了应用,取得了较好的效果,同时对实际电网中出现的问题,在方法上进行了改进,增强了方法的实用性。

探地雷达GPR正演模拟的时域有限差分实现(英文)

近年来,随着数字处理技术和电子技术的飞速发展,探地雷达(GPR)的实际应用范围迅速扩大,现已覆盖考古、矿产资源勘探、水文地质调查、岩土勘查、无损检测、工程建筑物结构调查、军事等众多领域,解决了很多工程实际问题,成为浅层勘探的有力工具.而探地雷达的理论研究与实际的应用相比,具有明显的滞后性.但是解释人员要达到精确地对探地雷达实际资料的进行解析,必须事先了解地质体的雷达反射剖面的特征,所以作为反演与解释基础的复杂地电模型的探地雷达正演模拟技术,就成了探地雷达理论研究的主要内容之一.本文以麦克斯韦两个旋度方程为基本出发点,运用K.S.Yee的空间网格模型理论和时域有限差分法的基本原理,推导出二维空间的探地雷达正演方程组,并详细地分析了差分格式中半空间步长与半时间步长的实现方法,及其雷达波电场与磁场分量在计算机上相互关系的C程序实现.然后讨论了数值频散关系及其产生原因,通过同时考虑时域有限差分法及Yee氏网格的特点,推导出了符合探地雷达实际传播规律的理想频散关系,作者自制了探地雷达正演程序,并分别计算了Mur超吸收边界条件及无边界条件下的雷达地电模型,通过对比可知,超吸收边界条件可利用,大大地减少截断边界处的干扰波,达到用有限区域达到在无限空间传播的效果.最后作者利用自制程序,对“V”字形和同一斜面上的五个圆的两个典型的探地雷达地电模型进行了正演模拟,得到了正演剖面图,消除了边界反射后的雷达剖面能很好地指导工作人员对雷达实测剖面的地质解释,同时使正演研究更符合实际的地质情况.