In this paper, a numerical approach to operability analysis is developed. Based on the concept of the steady-state region of attraction, which indicates the initial operating conditions that can be driven to a given setpoint, this analysis helps design engineers assess whether a process can be effectively controlled by linear control systems. Using steady-state operating data that are generated with process simulators (Such as ASPEN PLUS(R)), the steady-state region of attraction are calculated by solving the proposed constrained nonlinear optimization problem. The objective function for the optimization problem is developed using the ellipsoid theory while the constraint functions are formulated from the complement set of the initially defined constraint set and the boundary of the available input space (AIS). Furthermore, as the solution of a nonlinear optimization problem generally depends on the initial starting point, a method of generating a set of ellipsoids in the hyperspace to start the calculation is also developed. The proposed constrained nonlinear optimization problem can be solved using numerical tools such as MATILAB(R). This analysis can be performed based on process flowsheets with process simulators and thus is useful in early stages of process design. The use of the proposed numerical framework is illustrated by a case study of a methyl acetate reactive distillation column. (C) 2008 Elsevier Ltd. All rights reserved.