摘要

A hybrid ADI-FDTD sub-gridding method is proposed to improve strong dispersion medium GPR forward simulation, which uses fine mesh ADI-FDTD in the partial region of dramatic physical parameter changes, and coarse mesh FDTD in other regions, because the ADI-FDTD is not limited by CFL stability condition and we can choose larger time step according to coarse mesh to improve computational efficiency. This paper mainly discussed the field value exchange of the combined coarse Mesh and fine mesh and gave out the calculation procedure of the algorithm, through studying the radar wave based on Debye equation, then deduced the coarse mesh FDTD and fine mesh ADI-FDTD difference scheme of frequency dispersion media. Choosing a thin layer model as an example, we applied coarse mesh, fine grid and hybrid ADI-FDTD subgridding separately in forward simulation; the computing cost and the simulation accuracy clearly show the advantages of the hybrid ADI-FDTD subgridding method. Finally, we established the dispersion media and non-dispersion media composition model, applied three methods in these two models, and compared the advantages and disadvantages of three methods to analyze the waveform characteristics of the radar profiles in the non-dispersion and dispersion media and effectively guided the precise interpretation of radar data.