Advective terms

Squared velocity

To begin with, we consider

1J¯¯Jhξ1u1ξ1δξ1u1ξ11J¯¯Jhξ2u2ξ1δξ2u1ξ21J¯¯Jhξ3u3ξ1δξ3u1ξ31Ju1¯δξ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ξ31Ju2¯δξ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ξ31Ju3¯δξ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

ujuixjuiujxj=ujuixj,

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:

N3k=1N2j=1N1+12i=12Ju1(1J¯¯Jhξ1u1ξ1δξ1u1ξ11J¯¯Jhξ2u2ξ1δξ2u1ξ21J¯¯Jhξ3u3ξ1δξ3u1ξ3),N3k=1N212j=12N1i=1Ju2(1J¯¯Jhξ1u1ξ2δξ1u2ξ11J¯¯Jhξ2u2ξ2δξ2u2ξ21J¯¯Jhξ3u3ξ2δξ3u2ξ3),N312k=12N2j=1N1i=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

N3k=1N2j=1N1+12i=12Ju1{1Jδξ1(¯Jhξ1u1ξ1¯u1ξ1)+1Jδξ2(¯Jhξ1u2ξ1¯u1ξ2)+1Jδξ3(¯Jhξ1u3ξ1¯u1ξ3)},N3k=1N212j=12N1i=1Ju2{1Jδξ1(¯Jhξ2u1ξ2¯u2ξ1)+1Jδξ2(¯Jhξ2u2ξ2¯u2ξ2)+1Jδξ3(¯Jhξ2u3ξ2¯u2ξ3)},N312k=12N2j=1N1i=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ξ31JT{δξ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:

N3k=1N2j=1N1i=1JT(1J¯Jhξ1u1δξ1Tξ11J¯Jhξ2u2δξ2Tξ21J¯Jhξ3u3δξ3Tξ3),

giving

N3k=1N2j=1N1i=1JT{1Jδξ1(Jhξ1u1¯Tξ1)+1Jδξ2(Jhξ2u2¯Tξ2)+1Jδξ3(Jhξ3u3¯Tξ3)}.

Thus the advective contribution vanishes.