The equimaterial approach for the numerical solution of the one dimensional transient diffusion equation with zero flux boundary conditions

Abstract

A new approach is proposed for the grid motion for the numerical solution of a general transient diffusion equation in one spatial dimension with zero flux boundary conditions. The new criterion for grid motion is that the solute amount contained in each discretization section should be a pre-described fraction of the total solute amount at each time step. This requirement is not explicitly enforced to the solution technique but it is implicitly included in the equation through the appropriate variable transformation. The results showed that although the technique leads to the required grid motion the numerical results are of pure quality due to the appearance of singularities during the variable transformation procedure. Nevertheless, it is shown that by appropriate numerical handling of the solution at the singularity region the technique can lead to accurate results and potentially can replace the existing moving grid algorithms at least for the particular problem at hand.