Я предлагаю свести Вашу задачу к задаче Коши на всей оси. Тогда будет применима формула Даламбера. Для этого нам надо так подобрать начальные данные при

и

, чтобы решение задачи Коши совпало с решением смешанной задачи в цилиндре

. Для этого всего-то и надо - обеспечить условия

.
Давайте начнем с более простого примера. Пусть данные уже заданы на полуоси

. Вопрос. Как подобрать нач. данные при

так, чтобы решение получившейся задачи Коши равнялось 0 в точке

? Идею, как это сделать я уже озвучивал. Таким образом мы "избавляемся" от стенки

. Аналогично можно разбираться и с

. Отмечу, что нам даже не надо продолжать нач. данные на всю ось. Достаточно продолжить на отрезок

. Но это уже прям предельный случай. На самом деле легко и просто получается продолжение на отрезок

.
Можно подойти к этой задаче и с другой стороны.
В чем суть метода характеристик? В том, что вдоль характеристик сохраняется либо

либо

(в зависимости от характеристики). А значит, выпуская характеристики из начальных данных, Вы можете получить значения

и

на стенках. А там

, значит и

, значит на стенках можно вычислить

. Теперь можно выпускать характеристики и со стенок. По сути это анализ отражения волн от стенок. Для

достаточно всего одного отражения ...
Есть еще один способ решить Вашу задачу. Сосчитать ряд Фурье ... Это сделать можно. Заметим, что

После этого остается сосчитать ряд

С помощью тригонометрии он сводится к ряду

А это просто дельта-функция ( данный ряд - разложение по косинусам)
Три метода - более чем достаточно. Выбор за Вами.