Дальше, наверное, надо вводить метрику. И вот с метрикой правило Лейбница развернётся во всей красе.
До введения метрики можно рассказать ещё про:
1) Картановскую внешнюю дифференциальную 2-форму кривизны, и вот на ней-то показать тождества Бьянки.
Правда сначала придётся грубо вводить в дифференциальные формы в том случае если ранее у читателя с ними опыта работы не было.2) Производную Ли от всяческих тензорных полей, правило Лейбница для неё. Показать что в производной Ли частные производные можно заменить на ковариантные (связности

всё равно сокращаются).
3) Инвариантное дифференцирование тензорных полей по параметрам. Там тоже используется производная Ли, правда, не многие об этом знают.
(Если в двух словах, то...)
Если в двух словах, то есть значит

координат

и

параметров

. В координатном пространстве есть связность

и кривизна

. И в пространстве параметров, вообще говоря, тоже есть какая-то своя связность

и своя кривизна

. Покуда координаты от параметров не зависят, эти два пространства живут сами по себе. Смотрите что же происходит когда координаты (да и все тензорные поля в

-пространстве) начинают зависеть от параметров

. Допустимые преобразования координат и параметров таковы:


То есть параметры

преобразуются только сами через себя, а вот координаты

преобразуются сами через себя и ещё зависят от параметров. Из-за таких преобразований возникают специфические поля

(они навроде связностей). По трансформационным свойствам поля

ведут себя подобно величинам

. Инвариантное дифференцирование тензорных полей по параметрам

:



В этих выражениях, что замечательно, частные производные по

могут быть заменены на ковариантные (это потому, что связности

всё равно сократятся). В этом нет ничего удивительного, так как там стоят как раз таки производные Ли вдоль полей

, а производные Ли такой симметрией обладают.
Инвариантные производные по параметрам некоммутативны

, соответственно дополнительно к

и

появляется ещё одна кривизна -- "кросс-кривизна" (у неё два индекса тензорные из пространства X, и два индекса тензорные из пространства параметров). В общем, богатенькая получается математическая структура.
Частный случай

, единственный параметр

имеет следующий физический смысл: инвариантная производная Эйлера по времени

при наличии поля скоростей

. Здесь же где-то рядом лежит и кориолисово ускорение.