Infrasound can be modeled using either ray tracing for geometric acoustics, parabolic equation simulations for amplitude fields or normal mode / full waveform approaches for wave pulse propagation.
Furthermore characterizing the atmospheric background by a combination of climatologies (HWM, MSISE) and meteorological (forecast/analysis) data (e.g. ECMWF, G2S) is important for realistic infrasound propagation.
Free software packages including most of these approaches are "NCPA-PROP" or "GeoAc". Other, older codes are "HARPA", "WASP-3d" or "Tau-P" for ray-tracing and "InfraMAP" for PE. There are also a lot of individual / institutional software packages solving wave equation and propagation by finite differences / parabolic equation / etc methods or by using adaptations from seismology, audible acoustics, geometric optics, hydrodynamics or hydroacoustics.