Решение

При решении задачи я буду принебрегать зависимостью скорости звука $c$ от температуры и давления. Ось $x$ направим вдоль земли по направлению ветру, а ось $z$ в атмосферу. Поскольку скорость ветра $V(z)$ заметно уменьшается при приближении к поверхности земли, то градиент скорости ветра $\frac{d V}{dz}\gt 0$. Скорость ветра зависит только от координаты $z$.

Источник звука находится в точке (0,0). Атмосферу можно представить в виде бесконечного числа тонких горизонтальных слоев. Чтобы волновой фронт (гребень волны) при переходе из слоя в слой оставался единым целым и не разрывался, скорость движения точки пересечения фронта с горизонтальной плоскостью (горизонтальная фазовая скорость $\upsilon_{x}$) должна быть строго одинаковой на любой высоте:

$$\upsilon_{x}=const$$

Ещё раз: если бы в каком-то слое гребень волны бежал по горизонтали быстрее, а в соседнем — медленнее, то волновой фронт просто разорвался бы. Это наш инвариант.

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

Пусть угол между направлением скорости звука и вертикальной осью $z$ равен $\theta$. На рисунке $\textbf{а}$: $BC$ - это фронт волны, $AC$ это расстояниие, которое проходит фронт до горизонтального слоя, точка $B$ это точка пересечения фронта с горизонтальным слоем, $AB$ это расстояние, которое пройдёт точка $B$. Из прямоугольного $\Delta ACB$:

$$\frac{c\cdot dt}{\sin\theta}=dS\Longrightarrow \frac{dS}{dt}=\upsilon^{'}=\frac{c}{\sin\theta}$$

Добавляя скорость движения ветра, получим:

$$\upsilon_{x}=\frac{c}{\sin\theta(z)}+V(z)=C_{0}=const$$

Зададим начальные условия на поверхности земли $z=0$, где скорость ветра $V(0)=0$, а угол выхода луча из источника звука равен $\theta_{0}$:

$$C_{0}=\frac{c}{\sin\theta_{0}}$$

Итоговый инвариант связывает угол луча на любой высоте с начальным углом вылета:

$$\frac{c}{\sin\theta(z)}+V(z)=\frac{c}{\sin\theta_{0}} \tag{1}$$

Точка на гребне волны движется со скоростью, равной векторной сумме скорости звука и скорости ветра. Смотри рисунок $\textbf{б}$. Кинематические уравнения движения этой точки (звукового луча) в проекциях на оси $x$ и $z$ имеют вид:

$$\frac{dx}{dt}=c\sin\theta(z)+V(z)$$
$$\frac{dz}{dt}=c\cos\theta(z)$$

Разделим первое уравнение на второе и получим дифференциальное уравнение траектории луча:

$$\frac{dx}{dz}=\tan\theta(z)+\frac{V(z)}{c\cos\theta(z)} \tag{2}$$

Выразим тригонометрические функции угла $\theta(z)$ из нашего инварианта (1):

$$\sin\theta(z)=\frac{\sin\theta_{0}}{1-\frac{V(z)}{c}\sin\theta_{0}}$$

Используя основное тригонометрическое тождество $\cos\theta=\sqrt{1-\sin^{2}\theta}$:

$$\cos\theta(z)=\frac{\sqrt{\left( 1-\frac{V(z)}{c}\sin\theta_{0} \right)^{2}-\sin^{2}\theta_{0}}}{1-\frac{V(z)}{c}\sin\theta_{0}}$$

Подставим полученные функции в дифференциальное уравнение траектории луча (2):

$$\frac{dx}{dz}=\frac{\sin\theta_{0}+\frac{V(z)}{c}\left( 1-\frac{V(z)}{c}\sin\theta_{0} \right)}{\sqrt{\left( 1-\frac{V(z)}{c}\sin\theta_{0} \right)^{2}-\sin^{2}\theta_{0}}}$$

Интегрируя это выражение по высоте от 0 до $z$, получаем общее уравнение траектории для произвольного профиля ветра $V(z)$:

$$x(z)=\int_{0}^{z}\frac{\sin\theta_{0}+\frac{V(\zeta)}{c}\left( 1-\frac{V(\zeta)}{c}\sin\theta_{0} \right)}{\sqrt{\left( 1-\frac{V(\zeta)}{c}\sin\theta_{0} \right)^{2}-\sin^{2}\theta_{0}}}d\zeta \tag{3}$$

Рассмотрим приземное приближение для слабого линейного ветра:

$$V(z)\ll c, \quad V(z)=\beta z$$

В приземном слое скорость ветра много меньше скорости звука, поэтому отношение $\epsilon=\frac{V(\zeta)}{c}=\frac{\beta\zeta}{c}\ll 1$ является малым параметром. Разложим подынтегральное выражение в ряд Тейлора по $\epsilon$, удерживая только линейные члены. Так как выражения весьма громоздки, я сделаю это отдельно для числителя и отдельно для знаменателя.

Разложение знаменателя (под корнем):

$$(1-\epsilon\sin\theta_{0})^{2}-\sin^{2}\theta_{0}=1-2\epsilon\sin\theta_{0}+\epsilon^{2}\sin^{2}\theta_{0}-\sin^{2}\theta_{0}\approx \cos^{2}\theta_{0}-2\epsilon\sin\theta_{0}$$

Используя приближение $(1-x)^{-1/2}\approx 1+\frac{x}{2}$, запишем:

$$\frac{1}{\sqrt{\cos^{2}\theta_{0}-2\epsilon\sin\theta_{0}}}=\frac{1}{\cos\theta_{0}}\left( 1-2\epsilon\frac{\sin\theta_{0}}{\cos^{2}\theta_{0}} \right)^{-1/2}\approx \frac{1}{\cos\theta_{0}}\left( 1+\epsilon\frac{\sin\theta_{0}}{\cos^{2}\theta_{0}} \right)$$

Разложение числителя:

$$\sin\theta_{0}+\epsilon(1-\epsilon\sin\theta_{0})\approx \sin\theta_{0}+\epsilon$$

Тогда подинтегральное выражение:

$$\frac{\sin\theta_{0}+\epsilon}{\cos\theta_{0}}\left( 1+\epsilon\frac{\sin\theta_{0}}{\cos^{2}\theta_{0}} \right)\approx \frac{\sin\theta_{0}}{\cos\theta_{0}}+\frac{\epsilon}{\cos\theta_{0}}+\epsilon\frac{\sin^{2}\theta_{0}}{\cos^{3}\theta_{0}}=\tan\theta_{0}+\frac{\epsilon}{\cos^{3}\theta_{0}}$$

Вернем замену $\epsilon=\frac{\beta\zeta}{c}$ и подставим упрощенное выражение под интеграл (3):

$$x(z)\approx \int_{0}^{z}\left( \tan\theta_{0}+\frac{\beta}{c\cdot\cos^{3}\theta_{0}}\zeta \right)d\zeta$$

Интегрируя, получаем уравнение параболических траекторий звуковых лучей:

$$x(z)\approx z\cdot\tan\theta_{0}+\frac{\beta z^{2}}{2c\cdot\cos^{3}\theta_{0}} \tag{4}$$

Проанализируем полученное уравнение. Поскольку источник находится на земле, физически лучи могут излучаться только в верхнее полупространство, что ограничивает начальный угол пределами $\theta_{0}\in (-\frac{\pi}{2};\frac{\pi}{2})$.

Диапазон по ветру: $\theta_{0}\in (0;\frac{\pi}{2})$. Оба коэффициента уравнения параболы положительны. Луч выходит вправо-вверх, и с ростом высоты его наклон к горизонту уменьшается, так как $\frac{dx}{dz}$ растет. Волновой фронт опрокидывается вперед, прижимая звуковую энергию к земле. Это зона отличной слышимости.

Диапазон против ветра: $\theta_{0}\in (-\frac{\pi}{2};0)$. Линейный член отрицателен, квадратичный — положителен:

$$x(z)=-z\cdot\left| \tan\theta_{0} \right|+\frac{\beta z^{2}}{2c\cdot\cos^{3}\theta_{0}}$$

Луч сначала идет влево-вверх (против ветра), но квадратичный член неуклонно тянет траекторию вправо. Парабола имеет точку разворота по горизонтали. Достигнув некоторого максимального удаления против ветра, луч разворачивается и уносится ветром в обратную сторону, уходя высоко в атмосферу.

Траектория лучей представлена на рисунке:

На этом задача полностью решена, но очень хочется найти геометрическую границу области акустической тени - места, куда звук не будет доходить от источника звука.

ПРОДОЛЖЕНИЕ СЛЕДУЕТ ...