An efficient numerical method to solve nonlinear sound propagation is presented. The frequency-domain Burgers equation, which includes nonlinear steepening and atmospheric absorption, is formulated in the form of the real and imaginary parts of the pressure. The new formulation effectively eliminates possible numerical issues associated with zero amplitude at higher frequencies occurring in a previous frequency-domain algorithm. In addition, to circumvent a high-frequency error that can occur in the truncated higher frequencies, a split algorithm is developed, in which the Burgers equation is solved below a cutoff frequency and a recursive analytic expression is used beyond the cutoff frequency. Finally, the Lanczos smoothing filter is incorporated to remove the Gibbs phenomenon. The new method is found to successfully eliminate high-frequency numerical oscillations and to provide excellent agreement with the exact solution for an initially sinusoidal signal with only a few harmonics. The new method is applied to a broad range of applications with a comparison to other methods to assess the robustness and numerical efficiency of the method. These include sonic boom, broadband supersonic jet engine noise, and helicopter high-speed impulsive noise. It is shown that the new method provides the fastest and most accurate predictions compared to the other methods for all the application problems. Nomenclature A = amplitude of the Fourier transform pressure, Pa s B n = nth harmonic amplitude c = blade chord length, m c 0 = speed of sound, m s 1 D j = jet exit diameter, m f = frequency of the acoustic signal, Hz f c = cutoff frequency, Hz f 0 = fundamental frequency of the sinusoidal signal, Hz k = wave number, m 1 M = Mach number M a = u a =c 0 , u a is the peak particle velocity at the source, M a p a = 0 c 2 0 m = parameter for plane, cylindrical, and spherical waves n = index of harmonic p = acoustic pressure, Pã p = Fourier transform of the pressure, Pa s p a = initial pressure amplitude of the sinusoidal wave, Pa Q = amplitude of the Fourier transform of the pressure square, Pa 2 s q = acoustic pressure squared, Pa 2 q = Fourier transform of the pressure square, Pa 2 s R = rotor blade radius, m r = propagation distance, m r f = final distance for the sound propagation, m r i = initial or starting distance for the sound propagation, m U = real part of the Fourier transform of the pressure square, Pa 2 s V = imaginary part of the Fourier transform of the pressure square, Pa 2 s X = real part of the Fourier transform of the pressure, Pa s x = shock formation distance, m Y = imaginary part of the Fourier transform of the pressure, Pa s = atmospheric absorption, neper m 1 0 = atmospheric absorption and dispersion, i d , m 1 = coefficient of nonlinearity d = atmospheric dispersion, m 1 = ratio of the specific heats = diffusivity of sound, m 2 s 1 "= nonlinear coefficients, = 0 c 3 0 , kg 1 s 0 = density, kg m 3 = dimensionless propagation distance, x= x = retarded time, s = phase of the Fourier transform of the pressure,...