Maxwell's theory of multicomponent diffusion and subsequent extensions are based on systems of mass and momentum conservation equations. The partial stress tensor, which is involved in these equations, is expressed in terms of the gradients of velocity fields by statistical and continuum mechanical methods. We propose a method for the solution of Maxwell's equations of diffusion coupled with Muller's expression for the partial stress tensor. The proposed method consists in a singular perturbation process, followed by a weak (finite element) analysis of the resulting PDE systems. The singularity involved in the obtained equations was treated by a special technique, by which lower-order systems were supplemented by proper combinations of higher-order equations. The method proved particularly efficient for the solution of the Maxwell-Muller system, eventually reducing the number of unknown fields to that of the classical Navier-Stokes/Fick system. It was applied to the classical Stefan tube problem and the Hagen-Poiseuille flow in a hollow-fiber membrane tube. Numerical results for these problems are presented, and compared with the Navier-Stokes/Fick approximation. It is shown that the 0-th order term of the Maxwell-Muller equations differs from a properly formulated Navier-Stokes/Fick system, by a numerically insignificant amount. Numerical results for 1st-order terms indicate a good agreement of the classical approximation (with properly formulated Navier-Stokes and Fick's equations) with the Maxwell-Muller system, in the studied cases. We propose a method for the solution of Maxwell's equations of diffusion with Muller's expression for the stress tensor.The method is particularly efficient for the solution of the Maxwell-Muller system for low Mach numbers.It reduces the number of unknown fields to that of the classical Navier-Stokes/Fick system.It is applied to Stefan tube and Hagen-Poiseuille flow in a hollow-fiber membrane tube.Numerical results indicate good agreement with classical approximation.