A non-steady state mathematical model of amino acid (phenylalanine (Phe)) and mineral salt (NaCl) solution separation by neutralization dialysis (ND) carried out in a batch mode is proposed. The model takes into account the characteristics of membranes (thickness, ion-exchange capacity and conductivity) and solutions (concentration, composition), as well as the flow rate of the solution in the dialyzer compartments. As compared to previously developed models, the new one considers the local equilibrium of Phe protolysis reactions in solutions and membranes and, also, the transport of all the phenylalanine forms (zwitterionic, positively and negatively charged) through membranes. A series of experiments on ND demineralization of NaCl and Phe mixed solution was carried out. In order to minimize Phe losses, the solution pH in the desalination compartment was controlled by changing the concentrations of the solutions in the acid and alkali compartments of the ND cell. The validity of the model was verified by comparison of simulated and experimental time dependencies of solution electric conductivity and pH, as well as the concentration of Na+, Cl– ions and Phe species in the desalination compartment. Based on the simulation results, the role of Phe transport mechanisms in the losses of this amino acid during ND was discussed.