Squared velocity
To begin with, we consider
1J¯¯Jhξ1u1ξ1δξ1u1ξ11J¯¯Jhξ2u2ξ1δξ2u1ξ21J¯¯Jhξ3u3ξ1δξ3u1ξ3−1Ju1¯δξ1(Jhξ1u1)+δξ2(Jhξ2u2)+δξ3(Jhξ3u3)ξ1,1J¯¯Jhξ1u1ξ2δξ1u2ξ11J¯¯Jhξ2u2ξ2δξ2u2ξ21J¯¯Jhξ3u3ξ2δξ3u2ξ3−1Ju2¯δξ1(Jhξ1u1)+δξ2(Jhξ2u2)+δξ3(Jhξ3u3)ξ2,1J¯¯Jhξ1u1ξ3δξ1u3ξ11J¯¯Jhξ2u2ξ3δξ2u3ξ21J¯¯Jhξ3u3ξ3δξ3u3ξ3−1Ju3¯δξ1(Jhξ1u1)+δξ2(Jhξ2u2)+δξ3(Jhξ3u3)ξ3,
where the first three terms for each direction are the advective terms in the momentum balance, while the last terms are the product of the corresponding velocity and the averaged incompressibility constraint and thus zero.
With some algebra, we find that they are equal to
−1Jδξ1(¯Jhξ1u1ξ1¯u1ξ1)−1Jδξ2(¯Jhξ1u2ξ1¯u1ξ2)−1Jδξ3(¯Jhξ1u3ξ1¯u1ξ3),−1Jδξ1(¯Jhξ2u1ξ2¯u2ξ1)−1Jδξ2(¯Jhξ2u2ξ2¯u2ξ2)−1Jδξ3(¯Jhξ2u3ξ2¯u2ξ3),−1Jδξ1(¯Jhξ3u1ξ3¯u3ξ1)−1Jδξ2(¯Jhξ3u2ξ3¯u3ξ2)−1Jδξ3(¯Jhξ3u3ξ3¯u3ξ3),
respectively.
In short, the above relation is the discrete description of
−uj∂ui∂xj−ui∂uj∂xj=−∂ujui∂xj,
i.e., the advective and the divergence forms of the advective terms are equivalent.
Now the global energy balance attributed to the advective terms are considered:
N3∑k=1N2∑j=1N1+12∑i=12Ju1(1J¯¯Jhξ1u1ξ1δξ1u1ξ11J¯¯Jhξ2u2ξ1δξ2u1ξ21J¯¯Jhξ3u3ξ1δξ3u1ξ3),N3∑k=1N2−12∑j=12N1∑i=1Ju2(1J¯¯Jhξ1u1ξ2δξ1u2ξ11J¯¯Jhξ2u2ξ2δξ2u2ξ21J¯¯Jhξ3u3ξ2δξ3u2ξ3),N3−12∑k=12N2∑j=1N1∑i=1Ju3(1J¯¯Jhξ1u1ξ3δξ1u3ξ11J¯¯Jhξ2u2ξ3δξ2u3ξ21J¯¯Jhξ3u3ξ3δξ3u3ξ3).
Our objective is to prove they are all zero, indicating that the advective terms do not alter the net amount of the quadratic quantities.
By using the relations derived in the prerequisite, each term leads to
N3∑k=1N2∑j=1N1+12∑i=12Ju1{1Jδξ1(¯Jhξ1u1ξ1¯u1ξ1)+1Jδξ2(¯Jhξ1u2ξ1¯u1ξ2)+1Jδξ3(¯Jhξ1u3ξ1¯u1ξ3)},N3∑k=1N2−12∑j=12N1∑i=1Ju2{1Jδξ1(¯Jhξ2u1ξ2¯u2ξ1)+1Jδξ2(¯Jhξ2u2ξ2¯u2ξ2)+1Jδξ3(¯Jhξ2u3ξ2¯u2ξ3)},N3−12∑k=12N2∑j=1N1∑i=1Ju3{1Jδξ1(¯Jhξ3u1ξ3¯u3ξ1)+1Jδξ2(¯Jhξ3u2ξ3¯u3ξ2)+1Jδξ3(¯Jhξ3u3ξ3¯u3ξ3)}.
The resulting terms are product of the velocity and the advective terms of the divergence form whose sign is flipped, and as a result the aforementioned conclusion is obtained.
Squared Temperature
We consider
1J¯Jhξ1u1δξ1Tξ11J¯Jhξ2u2δξ2Tξ21J¯Jhξ3u3δξ3Tξ3−1JT{δξ1(Jhξ1u1)+δξ2(Jhξ2u2)+δξ3(Jhξ3u3)},
giving
−1Jδξ1(Jhξ1u1¯Tξ1)−1Jδξ2(Jhξ2u2¯Tξ2)−1Jδξ3(Jhξ3u3¯Tξ3).
Now we focus on the global balance of the quadratic quantity:
N3∑k=1N2∑j=1N1∑i=1JT(1J¯Jhξ1u1δξ1Tξ11J¯Jhξ2u2δξ2Tξ21J¯Jhξ3u3δξ3Tξ3),
giving
N3∑k=1N2∑j=1N1∑i=1JT{1Jδξ1(Jhξ1u1¯Tξ1)+1Jδξ2(Jhξ2u2¯Tξ2)+1Jδξ3(Jhξ3u3¯Tξ3)}.
Thus the advective contribution vanishes.