抛物方程的多重非均匀网格模型及其应用

提出了抛物方程的多重非均匀网格模型,以准确求解三维空间存在多辐射源的电波传播问题。通过对不同辐射源建立不同的坐标系,并对其仿真空间采用不同的非均匀网格划分,构建了抛物方程的多重非均匀网格模型。在此基础上,实现了三维多辐射源问题的并行计算。实例仿真了空间存在四个辐射源的电波传播特性。结果表明,抛物方程的多重非均匀网格模型能够准确求解多源的空间电磁场分布特性,且在该算例中,并行技术使得抛物方程的计算速度提升了2.41倍,极大地提高了抛物方程对三维多源问题的求解效率。     

森林环境电波传播抛物方程模型的改进研究

将森林看成空气和植物组成的混合物, 应用两相混合物折射模型求解了森林的等效介电常数, 通过与实验结果的对比, 验证了该模型的正确性. 将该森林介电常数求解方法引入到抛物方程的森林模型中, 改进了抛物方程的森林模型. 相对于传统的森林环境电波传播模型, 该模型能考虑森林各组成要素对电波传播的影响, 更适合于实际不同地区、不同种类植物分布的森林环境中电波传播特性的求解. 此外, 引入了非均匀网格技术, 有效提高了大区域森林环境中电波传播问题的求解效率.最后基于该模型仿真分析了森林的植物体积含量、重量含水量等要素对电波传播特性的影响. 关键词： 森林 抛物方程 非均匀网格 电波传播     

海洋声学中三维抛物方程非均匀网格模型

在海洋声学中,三维抛物方程模型可以有效考虑三维空间的声传播效应。然而,采用三维抛物方程模型分析三维空间内的声传播问题时,计算时间较长,并且需要消耗较大的计算机内存,因此给远距离声场的快速精确计算带来了很大困难。为此,将非均匀网格Galerkin离散化方法用于三维直角坐标系下的水声抛物方程模型中,深度算子和水平算子Galerkin离散方式由均匀网格变为非均匀网格。仿真结果表明,三维直角坐标系下非均匀网格离散的抛物方程模型,在保持计算精度、提高计算速度的同时,可以实现远距离声场的快速预报。另外,针对远距离局部海底地形与距离有关的三维声传播问题,给出了声场快速计算方法;在海底保持水平的区域,采用经典Kraken模型,重构抛物方程算法的初始场,随后依次递推求解地形与距离有关海底下的三维声场。采用改进模型,证明了远距离楔形波导声强增强效应。     

基于柱坐标系抛物方程和矩量法的电波传播混合算法

针对包含近源障碍物条件下的电波传播问题,提出了一种新颖的电波传播预测混合建模方法:矩量法(MOM)和圆柱坐标系抛物方程法(PEM)混合建模方法(MOM-PEM);MOM用于包含辐射源和近源障碍物的小圆柱区域内的电波传播建模,PEM用于MOM计算空间外的大区域范围内电波传播建模。MOM和PEM的计算过渡区域进行精细化网格剖分处理以避免场强数值传递的不兼容。仿真模拟了三类近源障碍物存在场景下的电波传播问题:有限开窗屏障碍物、立方体障碍物以及包含辐射源的半封闭空间障碍物,并将混合算法计算得到的结果和相同环境下采用全矩量法计算得到的结果进行了数值对比,结果表明混合算法和矩量法在精度上吻合较好。     

基于数字地图的二维电波传播问题仿真

为快速预测二维地理环境下的电磁环境特性,应用二维抛物方程模型对电磁环境进行仿真。为了得到真实有效的地形数据,研究了从GeoTiff中抽取网格点上地理信息的方法,并利用双线性插值法计算了任意位置处的高程值。同时研究了地球表面两点之间计算距离的方法,将该方法的计算结果与GIS软件结果进行对比,验证了该方法的可靠性。在标准大气环境下,利用二维抛物方程模型仿真分析了不同距离处不同高度的电波传播传播因子的变化情况,为预测真实地理环境中的电波传播特性提供了一种有效的方法。     

基于JASMIN框架的抛物方程有限差分解法并行计算及其应用

抛物方程有限差分解法的网格步长严格受波长限制,在求解城市小区电波传播问题时,计算速度明显变慢,为此,基于JASMIN框架研究了抛物方程有限差分解法的并行方法,通过将同一步进面划分成多个网格片,并分配到不同的处理器进行运算,实现了抛物方程有限差分解法的并行计算。与解析解的对比验证了并行程序的正确性,同时通过实例分析了并行程序的高效性,算例表明,抛物方程有限差分解法的求解效率得到了有效的提高。最后,模拟和分析了某一电信基站天线在包含9栋规则建筑物的城市小区环境中的电磁特性,结果表明,该方法能够得到基站在空间各处的信号覆盖强弱,可以为基站选址提供参考。     

基于数字地图的二维电波传播问题仿真

为快速预测二维地理环境下的电磁环境特性,应用二维抛物方程模型对电磁环境进行仿真。为了得到真实有效的地形数据,研究了从GeoTiff中抽取网格点上地理信息的方法,并利用双线性插值法计算了任意位置处的高程值。同时研究了地球表面两点之间计算距离的方法,将该方法的计算结果与GIS软件结果进行对比,验证了该方法的可靠性。在标准大气环境下,利用二维抛物方程模型仿真分析了不同距离处不同高度的电波传播传播因子的变化情况,为预测真实地理环境中的电波传播特性提供了一种有效的方法。     

基于JASMIN框架的抛物方程有限差分解法并行计算及其应用

抛物方程有限差分解法的网格步长严格受波长限制,在求解城市小区电波传播问题时,计算速度明显变慢,为此,基于JASMIN框架研究了抛物方程有限差分解法的并行方法,通过将同一步进面划分成多个网格片,并分配到不同的处理器进行运算,实现了抛物方程有限差分解法的并行计算。与解析解的对比验证了并行程序的正确性,同时通过实例分析了并行程序的高效性,算例表明,抛物方程有限差分解法的求解效率得到了有效的提高。最后,模拟和分析了某一电信基站天线在包含9栋规则建筑物的城市小区环境中的电磁特性,结果表明,该方法能够得到基站在空间各处的信号覆盖强弱,可以为基站选址提供参考。     

复杂地理环境的图像分割及电波传播特性

复杂地理环境是电波传播不可避免的传播环境,不仅不规则地形会对电波传播产生影响,不同的地表媒质对电波传播也会产生不同的影响。为了使得电波传播特性的预测结果更加地准确有效,通过图像分割算法实现地表环境的简单分类,同时对不同的媒质赋予不同的电磁参数,并结合数字高程模型（DEM)数据,实现了既具有地形起伏信息又具有地表电磁环境参数的复杂地理环境建模。在此基础上,对地表电磁环境信息做了网格剖分处理,利用抛物方程（PE）模型对复杂地理环境下的电波传播特性进行了预测。     

耦合辐射的三维热化学非平衡流数值模拟

发展耦合辐射的三维热化学非平衡流场计算方法,可用于非结构网格.采用Jameson有限体积法求解耦合辐射源项的三维N-S方程.辐射源项通过求解辐射输运方程(Radiative Transport Equation RTE)获得.在空间方向上离散后,采用有限体积法求解辐射输运方程.化学模型包含11个组元,20个化学反应.采用该数值方法计算MUSES-C模型在速度为11.6 km·s^-1时的绕流流场及前驻点处的辐射热流密度.并通过对比,分析热辐射对流场的影响.     

New approaches to nonlinear diffractive field propagation

In many domains of acoustic field propagation, such as medical ultrasound imaging, lithotripsy shock treatment, and underwater sonar, a realistic calculation of beam patterns requires treatment of the effects of diffraction from finite sources. Also, the mechanisms of loss and nonlinear effects within the medium are typically nonnegligible. The combination of diffraction, attenuation, and nonlinear effects has been treated by a number of formulations and numerical techniques. A novel model that incrementally propagates the field of baffled planar sources with substeps that account for the physics of diffraction, attenuation, and nonlinearity is presented. The model accounts for the effect of refraction and reflection (but not multiple reflections) in the case of propagation through multiple, parallel layers of fluid medium. An implementation of the model for axis symmetric sources has been developed. In one substep of the implementation, a new discrete Hankel transform is used with spatial transform techniques to propagate the field over a short distance with diffraction and attenuation. In the other substep, the temporal frequency domain solution to Burgers' equation is implemented to account for the nonlinear accretion and depletion of harmonics. This approach yields a computationally efficient procedure for calculating beam patterns from a baffled planar, axially symmetric source under conditions ranging from quasilinear through shock. The model is not restricted by the usual parabolic wave approximation and the field's directionality is explicitly accounted for at each point. Useage of a harmonic-limiting scheme allows the model to propagate some previously intractable high-intensity nonlinear fields. Results of the model are shown to be in excellent agreement with measurements performed on the nonlinear field of an unfocused 2.25-MHz piston source, even in the near field where the established parabolic wave approximation model fails. Next, the model is used to compare the water path and in situ fields of a medical ultrasound device. Finally, the model is used to calculate the spatial heating rate associated with a nonlinear field and to simulate the phenomenon of saturation-induced beam broadening.     

A large-scale wave-current coupled module with wave diffraction effect on unstructured meshes

Based on the extended mild-slope equation,a large-scale wave module is developed.By combining the eikonal equation and the modified wave action equation,the wave model can account for diffraction in most situations such as in the lee of islands and breakwaters,and using unstructured meshes provides great flexibility for modelling the wave in the complex geomorphology of barriers and islands,also allowing for refinement of the grid resolution within computationally important domains.The numerical implementation of the module is based on the explicit second-order upwind finite-volume schemes in geographic space,the Flux-Corrected Transport(FCT)algorithm in frequency space and the implicit Crank-Nicolson method in directional space.The three-dimensional hydrodynamic module is then modified to couple with the wave model,where the wave readily provides the depth-dependent radiation stress and the wave-induced turbulence coefficient for the current fields,and the wave propagation takes into account the current-induced advection,refraction and diffraction of wave energy and the effect of water level.The applicability of the proposed model to calculate Snell’s Law,wave transformation over the breakwaters and the elliptic shoal,wave propagation over the rip current field and the undertow on a sloping beach is evaluated.Numerical results show that the present model makes better predictions of the near-shore wave propagation and complex three-dimensional(3D)near-shore circulation driven by the waves,considering analytical solutions and experimental values.     

Two-frequency mutual coherence function for electromagnetic pulse propagation over rough surfaces

The propagation of a transient electromagnetic pulse over irregular terrain is considered. We model the wave propagation using the parabolic wave equation, which is valid for near-horizontal propagation. We model the effect of scattering from the rough terrain by introducing a surface-flattening coordinate transform. This coordinate transform simplifies the boundary condition of our problem, and introduces an effective refractive index into our wave equation. As a result, the problem of propagation over an irregular surface becomes equivalent to the problem of propagation through random media. The parabolic equation is solved analytically using the path integral method. Both vertically polarized and horizontally polarized signals are treated. Cumulant expansion is introduced to obtain an approximate expression for the two-frequency mutual coherence function. From the mutual coherence function, spatial and temporal dependence of the propagating signal can be determined. It can be shown that scattering from the irregular surface can cause broadening of the transient signal. This can have a significant impact on the performance of radio communication systems.