A Four-Zone Model and Nonlinear Dynamic Analysis of Solution Multiplicity of Buoyancy Ventilation in Underground Building. (31st August 2020)