[en] Numerical simulation method was examined for chemical reactions of actinide elements U, Pu, Np, and Tc etc. in an aqueous nitric acid solution. It is known that the numerical calculation for the Purex process with chemical reactions and liquid flow becomes stiff, because time constant for the chemical reactions is two to three order of magnitude smaller due to the very fast reactions than that of mass transfer or of reaching distribution equilibrium. Recently in order to increase a time step Δt the partial equilibrium (P.E.) model, in which some very fast reactions are treated by the equilibrium law whereas other reactions are by the rate law, has been proposed. In the present study concentration change of the solutes in an aqueous solution with 30 chemical reactions, of which 4 are expressed by equilibrium equations, has been calculated. Description of the P.E. model and the comparison of the results and cpu time between the kinetic and the P.E. models are given. (author)