f ( r , u , t = t i ) = ρ i ( r ) δ D [ u - u i ( r ) ] ,    ( 1 )