The implementation of multicomponent diffusion in porous catalyst particles into resolved-particle fixed bed computational fluid dynamics (CFD) presents challenges. Usually an isotropic effective diffusivity is employed for each species, which includes molecular diffusion, Knudsen diffusion and viscous flow. To obtain this quantity, constant and isotropic molar flux ratios are assumed. In fixed bed models, pellets composed of isotropic materials are situated in non-isotropic fluid-phase gradients in the bed, which call the assumptions into question. In this paper we present an extended model using a non-isotropic effective diffusivity tensor, which allows evaluation of the flux assumptions. The method is illustrated for resolved particle steady-state CFD simulations of ethylene oxidation in small test beds of spheres and four-hole cylinders. Comparisons to the original constant molar flux ratio model are made, which show that the flux changes affect reaction, species and temperature profiles in the particles, and are stronger in particles with longer diffusion paths. Larger differences between results from the two diffusion models are also seen when the effective diffusivity is decreased by reducing porosity and increasing tortuosity, typical of a reaction system where deactivation increases by carbon deposition. (C) 2018 Elsevier Ltd. All rights reserved.