Condensed mode operation of gas-phase ethylene polymerization adds more complexity to the hydrodynamics of fluidized bed reactors (FBRs). This study aims at the development of an industrial ethylene polymerization FBR model operating in condensed mode, and in this regard, the classical two-phase theory was modified. Two modeling strategies were investigated: a simple twophase model (model I) and a more complex one taking the presence of a wet zone (model II) into consideration. To access the predictive capabilities of the two models, model predictions were compared with plant data. The absolute relative deviation percentages for model I in prediction of the production rate and catalyst flow rate were 10.4 and 14.8, respectively. However, these errors for model II were 0.80 and 10.36, respectively. It was observed that the temperature profile predicted by model I was more consistent with that of plant data when the bubble diameter was considered to be small because of enhanced heat transfer rates between the bubble and emulsion phases. However, this assumption was not necessary for model II, and assuming a maximum stable bubble diameter for bubbles in the reactor, which is a more realistic assumption, the predicted temperature profile was reasonable. According to the results of model II, two zones, wet (including liquid droplets) and dry (including evaporated liquid), are formed in the condensed-mode operation of FBRs. It was found that for low-condensing agent flow rates (70 ton/h in this study), the difference between process conditions (concentration and temperature) of wet and dry zones increases (ΔC = 22.16 mol/m 3 and ΔT = 40.19 K), which can produce bimodal polyethylene (PE). Furthermore, a search strategy was developed based on the imperialist competitive algorithm and model II to find polymerization conditions needed to produce a PE grade with a specified microstructure. The results showed that the proposed algorithm is efficient for such a purpose. The sensitivity analysis of model I proved that in low gas velocities, the emulsion phase temperature is very sensitive to the catalyst flow rate. Moreover, the results of model I demonstrate that the dependency of temperature and monomer concentration of the emulsion phase on particle elutriation is negligible, especially under normal operating conditions.