A second-order accurate, operator splitting scheme for reaction-diffusion systems in an energetic variational formulation.