General Variation of Integral Functional/Dependent on n Variables

Theorem
Let $ \mathbf x $ be an n-dimensional vector.

Let $ u = u \left ( { \mathbf x } \right ) $ be a real mapping.

Let $ J $ be a functional such that:


 * $ \displaystyle J \left [ { u } \right ] = \int_R F \left ( { \mathbf x, u, \frac{ \partial u }{ \partial \mathbf x } } \right ) \mathrm d x_1 \dots \mathrm d x_n $

Let $ \mathbf x^* $, $ u^* $ be such that:


 * $ \displaystyle \mathbf x^* = \boldsymbol \Phi \left ( { \mathbf x, u, \frac{ \partial u }{ \partial \mathbf x }; \epsilon } \right ) = \mathbf x + \epsilon \boldsymbol \phi \left ( { \mathbf x, u, \frac{ \partial u }{ \partial \mathbf x }  } \right ) + \mathcal O \left ( { \epsilon^{ 1 + 0_+ } } \right ) $


 * $ \displaystyle u^* = \Psi \left ( { \mathbf x, u, \frac{ \partial u }{ \partial \mathbf x }; \epsilon } \right ) = \mathbf x + \epsilon \psi \left ( { \mathbf x, u, \frac{ \partial u }{ \partial \mathbf x }  } \right ) + \mathcal O \left ( { \epsilon^{ 1 + 0_+ } } \right ) $

Then:


 * $ \displaystyle \delta J = \epsilon \int_R \left ( { F_u - \frac{ \partial F_{ u_{ \mathbf x } } }{ \partial \mathbf x } } \right ) \overline{ \psi } \mathrm d x_1 \dots \mathrm d x_n + \epsilon \int_R \frac{ \partial }{ \partial \mathbf x } \left ( { F_{ u_{ x } } \overline{ \boldsymbol \psi } +F \boldsymbol \phi } \right ) \mathrm d x_1 \dots \mathrm d x_n $

where


 * $ \displaystyle \overline{ \psi } = \psi - u_{ \mathbf x } \boldsymbol{ \phi } $

Proof
By definition:


 * $ \displaystyle \frac{ \partial x_i^* }{ \partial x_j } = \delta_i^j + \epsilon \frac{ \partial \phi_i }{ \partial x_j } + \mathcal O \left ( { \epsilon^2 } \right ) $


 * $ \displaystyle \det \left ( { J_{ \mathbf x^* } } \right ) = 1 + \epsilon \frac{ \partial \boldsymbol \phi }{ \partial \mathbf x } + \mathcal O  \left ( { \epsilon^2 } \right ) $


 * $ \displaystyle \Delta J = \int_R \left [ { F \left ( { \mathbf x^*, u^* \frac{ \partial u^* }{ \partial \mathbf x^*} } \right ) \left ( { 1 + \epsilon \frac{ \partial \boldsymbol \phi }{ \partial \mathbf x } } \right ) - F \left ( { \mathbf x, u, \frac{ \partial u }{ \partial \mathbf x }  } \right )  } \right ] \mathrm d x_1 \dots \mathrm d x_n + \mathcal O \left ( { \epsilon^2 } \right ) $


 * $ \displaystyle \delta J = \int_R \left [ { F_{ \mathbf x } \delta \mathbf x + F_u \delta u + F_{ u_{ \mathbf x } } \delta u_{ \mathbf x } + \epsilon F \frac{ \partial \boldsymbol \phi }{ \partial \mathbf x } } \right ] \mathrm d x_1 \dots \mathrm d x_n $


 * $ \displaystyle \delta J = \int_R \left [ { F_{ \mathbf x } \delta \mathbf x + F_u \overline{ \delta u } + F_u u_{ \mathbf x } \delta \mathbf x + F_{ u_{ \mathbf x } } \left ( { \overline{ \delta u } } \right )_{ \mathbf x } + \sum_{ i, j \mathop = 1 }^n F_{ u_{ x_i } } u_{ x_i x_j } \delta x_j + F \delta \left ( { \mathbf x } \right )_{ \mathbf x } } \right ] \mathrm d x_1 \dots \mathrm d x_n $


 * $ \displaystyle \frac{ \partial }{ \partial \mathbf x } \left ( { F \delta \mathbf x } \right ) = F_{ \mathbf x } \delta \mathbf x + F \left ( { \overline{ \delta \mathbf x } } \right )_{ \mathbf x } + F_u u_{ \mathbf x } \delta \mathbf x + \sum_{ i, j \mathop = 1 }^n F_{ u_{ x_i } } u_{ x_i x_j } \delta x_j $


 * $ \displaystyle F_{ u_{ \mathbf x } } \left ( { \overline{ \delta u } } \right )_{ \mathbf x } = \frac{ \partial }{ \partial \mathbf x } \left ( { F_{ u_{ \mathbf x } } \overline{ \delta u } } \right ) - \frac{ \partial F_{ u_{

\mathbf x } }}{ \partial \mathbf x } \delta u $


 * $ \displaystyle \delta J = \int_R \left ( { F_u - \frac{ \partial F_{ u_{ \mathbf x } }}{ \mathbf x } } \right ) \overline{ \delta u } \mathrm d x_1 \dots \mathrm d x_n + \int_R \frac{ \partial }{ \partial \mathbf x } \left ( { F_{ u_{ \mathbf x } } \overline{ \delta u } + F \delta \mathbf x } \right ) \mathrm d x_1 \dots \mathrm d x_n $


 * $ \overline{ \delta u } = \epsilon \overline \psi $


 * $ \delta \mathbf x = \epsilon \boldsymbol \phi $