Modeling of non-equilibrium solidification in multi-component alloys is of singular importance in microstructure control, which however owing to the complex systems with complex additional constraints is still an open problem. In this work, the thermodynamic extremal principle was applied to solve the complex additional constraints self-consistently in thermodynamics. Consequently, short-range solute redistribution and long-range solute diffusion that share the same mobility are integrated naturally into the solute diffusion equations, thus avoiding the introduction of additional kinetic coefficients (e.g. interface permeability) to describe solute redistribution. Application to the non-equilibrium solidification of Al-Si-Cu alloys shows that anomalous solute trapping and anomalous solute profiles within the diffuse interface could occur, thus highlighting the important effect of the interaction among the component elements on the interface kinetics. The current phase-field model might be preferred for simulations not only because of its simplest form of evolution equations but also its feasibility to increase the simulation efficiency by the "thin interface limit" analysis. (C) 2017 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.