Аннотация:Рассмотрена задача ускорения алгоритма расчета нелинейных эффектов при моделировании высокоинтенсивных ультразвуковых пучков на основе однонаправленного уравнения Вестервельта.При построении численного решения для сильно искаженных волн с ударными фронтами необходимо учитывать большое число гармоник (до 1000) на пространственных сетках с размером матрицпорядка 10000 на 10000, что требует обработки больших объемов данных и длительного временирасчетов. В данной работе реализация оператора нелинейности проводится во временном представлении с использованием удароулавливающей схемы типа Годунова, которая позволяет моделировать нелинейные волны с ударными фронтами с небольшим (3) количеством узлов сетки на ударном фронте. В работе проводится сравнение эффективности использования данного метода при егореализации на центральном процессоре (CPU) и графических ускорителях (GPU) по сравнению соспектральным методом, реализованным ранее для квазилинейного распространения волны. Проводится анализ скорости выполнения алгоритмов на CPU и GPU в зависимости от размеров массивов входных данных.