The aim of this thesis is to implement multicomponent adsorption models in the custom OpenFOAM computational fluid dynamics solver adsorpFoam developed at Vienna University of Technology. This includes equilibrium and kinetics models. For this, two multicomponent equilibrium models, the Extended Langmuir Model ELM and the Ideal Adsorbed Solution Theory IAST, are used. They solely depend on single-component isotherm data. For interspecies-dependent kinetics, a diffusion-based approach is chosen. As the results of a zero-dimensional model show, the quality of prediction of equilibria is dependent on the chosen system of species. The model predictions are compared with experimental data of six multicomponent systems taken from literature. If experimental data are available, a simple extension to the ELM is possible. This is done by introducing empirical interaction coefficients to account for competitive adsorption, which improves the prediction of most systems. For this approach, data of multicomponent adsorption experiments have to be obtained. OpenFOAM, an open-source suite of CFD programs, is used in this thesis. The main reasons for this choice are its openness and extensibility. The implementation in OpenFOAM includes the adaptation of the governing equations, calculation of adsorption equilibrium loading and rate of adsorption. Additionally, the released heat of adsorption increases the temperature distribution of the adsorbing walls. If the calculated rate of adsorption leads to nonphysical results, e.g. more mass adsorbing in one cell than available, limiters are applied. At the end, a working multicomponent adsorption model was included in the solver adsorpFoam. It allows to define multiple adsorbing sites with different parameters per species and site. This implementation is a first step towards multicomponent mass transfer and serves as a basis for further work on the three-dimensional simulation of multicomponent adsorption.