A control system for controlling an operation of a heating ventilation and air conditioning, HVAC, system configured to condition an indoor environment by pushing air to the environment through a set of inlets arranged at a set of locations on one or multiple walls of the environment, wherein a thermal state of the air pushed to the environment at an inlet location includes one or combination of a temperature, a velocity and a humidity of the air, the control system comprising: an input interface configured to accept data indicative of a target distribution of thermal state in the environment, wherein a thermal state at a location in the environment includes one or combination of a temperature, a velocity, and humidity of the air; a memory configured to store an airflow dynamics model, ADM, defining a distribution of a thermal state in the environment subject to boundary conditions for thermal state of the air at the walls of the environment; and an HVAC model defining dynamics of the HVAC system; a processor configured to inverse the ADM to estimate values of the boundary conditions for the inlet locations defining target thermal state at the inlet locations that result in the target distribution of thermal state in the environment; determine, using the HVAC model, target control parameters of actuators of the HVAC system resulting in the target thermal state at the inlet locations; and submit control commands to the HVAC system to operate the actuators of the HVAC system according to the control parameters, wherein different combinations of values of target thermal state at the inlet locations result in the target distribution of thermal state in the environment, and wherein the processor selects a combination of thermal state based on a metric of performance of the HVAC system, wherein the metric of performance is defined by a multi-objective cost function of a combination of a cost of operation of the HVAC system and a difference between the target thermal distribution and a corresponding current thermal distribution, such that the processor is configured to determine the target thermal state by minimizing the multi-objective cost function, wherein the processor iteratively minimizes the multi-objective cost function until a termination condition is met, wherein, for performing an iteration, the processor is configured to: determine a sensitivity of the cost function to an update of the boundary conditions for the inlet locations; update the boundary conditions at the inlet locations in a direction of the sensitivity; determine the current distribution of thermal state according to the ADM (224, 1006) with the updated boundary conditions; and determine the cost of operation of the HVAC system resulting in the updated boundary conditions at the inlet locations. an input interface configured to accept data indicative of a target distribution of thermal state in the environment, wherein a thermal state at a location in the environment includes one or combination of a temperature, a velocity, and humidity of the air; a memory configured to store an airflow dynamics model, ADM, defining a distribution of a thermal state in the environment subject to boundary conditions for thermal state of the air at the walls of the environment; and an HVAC model defining dynamics of the HVAC system; a processor configured to inverse the ADM to estimate values of the boundary conditions for the inlet locations defining target thermal state at the inlet locations that result in the target distribution of thermal state in the environment; determine, using the HVAC model, target control parameters of actuators of the HVAC system resulting in the target thermal state at the inlet locations; and submit control commands to the HVAC system to operate the actuators of the HVAC system according to the control parameters, inverse the ADM to estimate values of the boundary conditions for the inlet locations defining target thermal state at the inlet locations that result in the target distribution of thermal state in the environment; determine, using the HVAC model, target control parameters of actuators of the HVAC system resulting in the target thermal state at the inlet locations; and submit control commands to the HVAC system to operate the actuators of the HVAC system according to the control parameters, determine a sensitivity of the cost function to an update of the boundary conditions for the inlet locations; update the boundary conditions at the inlet locations in a direction of the sensitivity; determine the current distribution of thermal state according to the ADM (224, 1006) with the updated boundary conditions; and determine the cost of operation of the HVAC system resulting in the updated boundary conditions at the inlet locations.