In this thesis, modeling, simulation and control of a reactive distillation tower for production of Ethyl acetate were studied. Simulation of the tower was performed based on mass and energy balance equations using the equilibrium staged model. In order to evaluate effects of simultaneous chemical reaction and phase equilibrium, by using of k_value approach and introducing of correlation factor z, the phase - reaction equilibrium was reduced to quasi-phase equilibrium and then effect of chemical reaction was examined. Wilson activity coefficient was used to consider non-ideal behavior of the liquid phase. Applications of the correction factor results in a faster convergence for solution of equations in simulation of the tower’s behavior toward the next steady state. The behavior of the developed dynamic model was validated through comparison of the data with experimental data. The average of errors in liquid phase mole fraction was 7.7%. Then, the dynamic model was used to study effects of operational parameters such as reflux ratio, reboiler duty, reaction rate, liquid hold-up, temperature and feed composition on operation of the tower. The gained results show that increasing of reflux ratio enhances top product concentration, but the conversion is reduced and consequently concentrations of reactants on top of tower increase. This is unlike to conventional distillations that increasing of reflux ratio improves distillate concentration. Furthermore, by study of the reboiler duty effect on concentration of ethyl acetate in top product, the optimal value of reboiler duty for was obtained. In order to control top product concentration and according to non-linear behavior of the system, an adaptive control system based gain scheduling method was developed to adjust controller parameters continuously to accommodate changes in process dynamics. Also for finding the best manipulated variable, slope criterion, sensitivity criterion and SVD analysis was used. In the control system, PID controller was implemented and its parameters were tuned initially by Cohen-Coon method. In the next step, these parameters were optimized by integral methods through numerical methods and PSO algorithm. By comparison of the system responses, time response and overshoot value were improved by 20 and 1.6%, respectively. Keywords : Reactive distillation, Ethyl acetate, Modeling, Simulation, Control, Gain scheduling, Optimization, PSO algorithm.