Well,I think you can use the molecular dynamics to solve this problem.Set the Hamiltonian such as fpu model,phi-4model...,and then calculate the diffusion coefficient by fourier's law.BTW,if you do like this,you can find that abnormal diffusion is emerged in 2D and 1D.
needing more details you...