A method for controlling an operation of an air-conditioning system generating airflow in a conditioned environment, comprising: determining values of the airflow measured in the conditioned environment during the operation of the air-conditioning system; updating a model of airflow dynamics connecting values of flow and temperature of air conditioned during the operation of the air-conditioning system, wherein the updating iteratively reduces an error between values of the airflow determined according to the model and the measured values of the airflow, characterized in that the model includes a reduced order model having a number of parameters less than a physical model of the airflow according to a Boussinesq equation, wherein the Boussinesq equation is a partial differential equation, PDE, and wherein the reduced order model is an ordinary differential equation, ODE, and the model further includes a stability parameter representing a difference between the reduced order model and the physical model, and wherein the updating comprises updating only the stability parameter of the model in response to detecting the error by minimizing a cost function of the error subject to stability constraints; wherein the model is x ̇ r = b + Ax r + x r T Bx r + F K x. wherein b, A, B are constants, x is a vector of states of the airflow, x r is a vector of the states with a reduced dimension r, F is the stability parameter represented as a function of its arguments, and K is a vector of stability coefficients ensuring a stability of the model of the airflow dynamics the method further comprising: modeling the airflow using the updated model; and controlling the operation of the air-conditioning system using the modeled airflow, wherein steps of the method are performed using at least one processor of a controller. determining values of the airflow measured in the conditioned environment during the operation of the air-conditioning system; updating a model of airflow dynamics connecting values of flow and temperature of air conditioned during the operation of the air-conditioning system, wherein the updating iteratively reduces an error between values of the airflow determined according to the model and the measured values of the airflow, characterized in that the model includes a reduced order model having a number of parameters less than a physical model of the airflow according to a Boussinesq equation, wherein the Boussinesq equation is a partial differential equation, PDE, and wherein the reduced order model is an ordinary differential equation, ODE, and the model further includes a stability parameter representing a difference between the reduced order model and the physical model, and wherein the updating comprises updating only the stability parameter of the model in response to detecting the error by minimizing a cost function of the error subject to stability constraints; wherein the model is x ̇ r = b + Ax r + x r T Bx r + F K x. wherein b, A, B are constants, x is a vector of states of the airflow, x r is a vector of the states with a reduced dimension r, F is the stability parameter represented as a function of its arguments, and K is a vector of stability coefficients ensuring a stability of the model of the airflow dynamics the method further comprising: modeling the airflow using the updated model; and controlling the operation of the air-conditioning system using the modeled airflow, wherein steps of the method are performed using at least one processor of a controller.