Although experimental studies on the impact of feed (FS) and draw solutions (DS) on the forward osmosis (FO) applications are reported in literature, systematic mathematical modeling considering the dynamic change in solution properties is lacking. In this study, asymmetric FO membrane simulation model was established using Aspen Plus-MATLAB subroutines algorithm to account for the effect of concentration polarization (CP), types of FS and DS and in their properties on FO performance. The developed model was validated by comparing the simulation with experimental results. The model successfully predict the performance of FO process under wide varieties of operational conditions, FS and DS flow rates and concentrations. The model showed that the variation of MCFDS concentration had a marked effect on water flux (WF) in contrast to flow rate. The WFs obtained from seawater (SW) increased from 5.28 L/m2.h to 42.08 L/m2.h as MCFDS changes from 150 g/L to 300 g/L which corresponding to 11.66% to 45.33% of water recovery. As for synthetic aquaculture wastewater (SAWW), 9.70 L/m2.h to 37.32 L/m2.h of WFs were exhibited with the increase of MCFDS concentration from 50 g/L to 200 g/L, respectively. The effect of concentrated external CP (CECP) was found to be significant in case of SW and negligible with SAWW. Whereas, increasing MCFDS concentration increases the severity effect of dilutive internal CP (DICP). The degree of DICP depends on the solute resistivity (KD) of porous layer, which were elevated (4.22-5.88 s/m) as MCFDS concentration increases (150-300 g/L). The study demonstrated the effectiveness and suitability of the developed Aspen Plus-MATLAB model simulating the FO process. Copyright © 2019 Elsevier B.V. All rights reserved.