Note: For Task 1, the thermal boundary condition at the wall (excluding inlet, outlet, and bottom plate) should be just "thermally insulated", i.e., heat flux = 0. This is actually the default already. Further background related to Boussinesq setting: In case it might be useful, the following is the link to a paper that provides the empirically determined relation of ρ as a function of T for water. (See Eq. (1) or Eq. (3) in the paper.) NIST 1992 formula for ρ(T) By taking the derivative of ρ(T) one can readily obtain the thermal expansion coefficient. Note that the thermal expansion coefficient (as used in Fluent) is defined by β ≡ α−1(∂α/∂T) where α ≡ 1/ρ is specific volume. Alternatively, giving the operating temperature, the operating density and thermal expansion coefficient can also be obtained from look-up tables in standard engineering handbooks (some are available online).Discussion of solution
profile data file, flyingsaucer2DH.txt (For reference only. Matlab code that generates flyingsaucer2DH.txt)Discussion of solution