The hydrodynamics of the Expanded Bed Adsorption process is studied through simulations combining Computational Fluid Dynamics and the Discrete Element Method. A representative base case is defined, based on process design parameters commonly encountered in literature. Then, 19 other cases are defined, each representing a singular adjustment to the column design, material properties, or operating conditions. The parameters that are varied are the expansion factor, liquid viscosity, bed aspect ratio, mean particle density, width of the particle density distribution, width of the particle size distribution, column taper angle, and column alignment angle. The impact of each adjustment on the bed behaviour is discussed, using the local particle size distribution and solids dispersion coefficient as main indicators of bed stability. Optimal performance was found for an expansion factor of two to three, and the combination of particle size distribution and particle density distribution was found to greatly improve bed stability. The mixing process of the liquid and solid phases is concluded to be of highly complex nature, and cannot simply be predicted from the liquid flow velocity.