Vegetable oils are among the fastest expanding commodities in the world, comprising a wide range of industrial applications, especially in the food and pharmaceutical sectors, thus playing a significant role in the human diet. In order to make them edible, vegetable oils are submitted to a refining process, in which undesirable compounds are removed. Those compounds are responsible for unpleasant odours, tastes, and inadequate colour. Palm oil is the most produced vegetable oil in the world, being mainly used for food purposes, but also for the production of biodiesel, and hence issues related to sustainability, have been raised in recent years. In the literature, there is a gap of works dealing with multicomponent mixtures in parastripping, as is the case of vegetable oils, which are composed of large molar mass molecules, such as tri-, diand monoacylglycerols, in addition to free fatty acids and nutraceutical compounds. In this context, this work aimed to model and simulate the multicomponent parastripping process applied to physical deacidification/deodorisation processes of palm oil. Parastripping is an adaptation to the conventional process (stripping) and consists of dividing the steam flow at the bottom of the column into two ascending and parallel streams that come into contact with a single descending liquid stream. This structural rearrangement can increase the number of theoretical ideal stages per unit of column height, reducing operational and design costs. The simulations were carried out by implementing the subroutines for calculating the physical properties for each of the oil and steam components, as well as an adaptation of the algorithm proposed by Naphtali & Sandholm, in the open-source software Python® using its libraries for solving systems of non-linear equations. The model palm oil was represented by a multicomponent mixture containing triacylglycerols, as major compounds, and diacylglycerols, monoacylglycerols, and tocopherol, as minor compounds. The influence of operational parameters (i.e., temperature, pressure, stripping steam flow rate) was investigated regarding the compositions of the outlet streams (i.e., final oil acidity, neutral oil loss, tocopherol recovery), in addition to a comparison with the results obtained by the conventional process. For the same column height and degree of deacidification, the parastripping column achieved similar results as the stripping column but with 14% less stripping steam. For the same amount of stripping steam and the same degree of deacidification, it was possible to reduce the height of the parastripping column in 18.2%. There were few variations in the neutral oil losses and tocopherol recoveries, comparing the two processes, under the same operating conditions. It can be thus concluded that it is possible to use a lower height column, consuming the same amount of stripping steam, or a same height column, with a lower consumption of stripping steam. That means the parastripping process confers operational and construction advantages, when co...