MathOverflow Asked by Amir Sagiv on December 17, 2020

The Leibniz integral rule, in its multivariate form, deals with differentiation of the following sort:

$$ frac{partial}{partial t} int_{D(t)} F({bf x}, t) , d{bf x} , , qquad D(t)in mathbb{R}^d , .$$

I am looking for a fully rigorous formulation of this theorem, as well as a proper proof. So far, I could only find:

- The one-dimensional case (see e.g., Courant calculus book).
- Physics-flavored proofs, where the normal speed $v_n ({bf x})$ is not well defined (see e.g., here). While insightful, this is not what I need.
- Differential-geometry and form-based proofs (see the above paper). While these are valid, they are far more general than what I need.

**My question:** Are you familiar with a multivariate-calculus version and proof of this theorem?

I searched most of the textbooks mentioned here and here, but still no luck. Any ideas? I believe it is not usually covered in graduate curses of advanced multivariate analysis, and so it is pretty difficult to find it in the textbooks.

(This is cross posted from math.se, since a week has passed and it recieved little-to-no-answer. Also, this is not to be confused with the Leibniz rule $(fg)’=f’g +fg’$, as in this post.)

This is probably just a more streamlined version of the proof in the article you linked to, mentioning somewhat less differential forms. I don't know a more elementary proof, but hope this is of use anyway.

I'll write $D_t subset mathbb{R}^d$ for the domain of integration and I'll assume that the map $tmapsto D_t$ is sufficiently regular in the sense that for sufficiently short times, the motion of the domain can be described by the flow of an ambient vector field $X$ on $mathbb{R}^d$. I.e. for each $tin mathbb{R}$ there should exist a vector field $X$ on $mathbb{R}^d$ with flow $Phi_h:mathbb{R}^d to mathbb{R}^d$ such that $$ D_{t+h}=Phi_{h}(D_{t}) qquad forall h in (-epsilon,epsilon). $$

(The field $X$ is not unique, one may add any vector field that's tangent to the boundary of $D_t$.)

With such a vector field $X$ given we compute begin{align} frac{d}{dt} int_{D_t}F(mathbf{x},t)dmathbf{x} &= lim_{hto 0}frac{1}{h}left( int_{D_{t+h}}F(mathbf{x},t+h)dmathbf{x} - int_{D_{t}}F(mathbf{x},t)dmathbf{x}right)\ &=lim_{hto 0}frac{1}{h}left( int_{Phi_h(D_{t})}F(mathbf{x},t+h)dmathbf{x} - int_{D_{t}}F(mathbf{x},t)dmathbf{x}right)\ &=lim_{hto 0}frac{1}{h}left( int_{D_{t}}Phi^*_h(F(mathbf{x},t+h)dmathbf{x}) - int_{D_{t}}F(mathbf{x},t)dmathbf{x}right)\ end{align}

(The previous equality is just the substitution rule for integrals in several variables. So $Phi^*_hleft(F(mathbf{x},t+h)dmathbf{x}right)$ denotes the pullback of the differential form $F(mathbf{x},t+h)dmathbf{x}$ along $Phi_h$. It's what you get by replacing every $mathbf{x}$ with $Phi_h(mathbf{x})$ in the expression $F(mathbf{x},t+h)dmathbf{x}$. We continue by swapping limit and integral)

begin{align} &=int_{D_{t}}lim_{hto 0} frac{Phi^*_h(F(mathbf{x},t+h)dmathbf{x}) - F(mathbf{x},t)dmathbf{x}}{h} end{align}

(To go on use Taylor $F(mathbf{x},t+h)=F(mathbf{x},t)+frac{partial F(mathbf{x},t)}{partial t}cdot h+O(h^2)$ as well as linearity of $Phi^*$ and $Phi^*_0=text{id}_{mathbb{R}^n}$.) $$ =int_{D_{t}} left(frac{partial F(mathbf{x},t)}{partial t}dmathbf{x} +lim_{hto 0}frac{Phi^*_h(F(mathbf{x},t)dmathbf{x}) - F(mathbf{x},t)dmathbf{x}}{h}right). $$ Now that limit in the integral is per definition $L_X(F(mathbf{x},t)dmathbf{x})$, the Lie derivative of the form $omega=F(mathbf{x},t)dmathbf{x}$ with respect to $X$. Using Cartan's magic formula $L_x omega = d(i_X omega)+i_X(domega)$, which in our case simplifies to $L_x omega = d(i_X omega)$ since $F(mathbf{x},t)dmathbf{x}$ is of top degree, we get $$ =int_{D_{t}}frac{partial F(mathbf{x},t)}{partial t}dmathbf{x} +int_{D_{t}} d left(i_X(F(mathbf{x},t)dmathbf{x})right) $$ Applying Stokes theorem to the second integral we finally arrive at $$ frac{d}{dt} int_{D_t}F(mathbf{x},t)dmathbf{x} = int_{D_{t}}frac{partial F(mathbf{x},t)}{partial t}dmathbf{x} +int_{partial D_{t}} i_X(F(mathbf{x},t)dmathbf{x}). $$ For those who want to get rid of that insertion in the second integral, a direct computation shows that $int_Gamma i_X( f(mathbf{x})dmathbf{x})=int_Gamma f(mathbf{x}) Xcdot dvec{A}$ where $dvec{A}$ denotes the usual "surface element" with respect to the standard metric. So you can write the Leibniz integral rule as $$ frac{d}{dt} int_{D_t}F(mathbf{x},t)dmathbf{x} = int_{D_{t}}frac{partial F(mathbf{x},t)}{partial t}dmathbf{x} +int_{partial D_{t}} F(mathbf{x},t), Xcdot dvec{A}. $$

Correct answer by Michael Bächtold on December 17, 2020

Get help from others!

Recent Questions

- How can I transform graph image into a tikzpicture LaTeX code?
- How Do I Get The Ifruit App Off Of Gta 5 / Grand Theft Auto 5
- Iv’e designed a space elevator using a series of lasers. do you know anybody i could submit the designs too that could manufacture the concept and put it to use
- Need help finding a book. Female OP protagonist, magic
- Why is the WWF pending games (“Your turn”) area replaced w/ a column of “Bonus & Reward”gift boxes?

Recent Answers

- Peter Machado on Why fry rice before boiling?
- Lex on Does Google Analytics track 404 page responses as valid page views?
- haakon.io on Why fry rice before boiling?
- Joshua Engel on Why fry rice before boiling?
- Jon Church on Why fry rice before boiling?

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP