When I replaced
eqnFEM = Table[D[Xcomp[t, x][[i]], t] + Inactive[Div][Sum[((-D[[i,j]]*IdentityMatrix[1]).Inactive[Grad][Xcomp[[j]],{x}],{j,1,3}],{x}]==0, {i, 1, 3}];
bound1FEM =
Table[DirichletCondition[Xcomp[t, x][[i]] == list[[i]],
x == xbegin], {i, 1, 3}];
bound2FEM =
Table[DirichletCondition[Xcomp[t, x][[i]] == list[[i + 3]],
x == xend], {i, 1, 3}];
by
eqnFEM = Table[D[Xcomp[t, x][[i]], t] + Sum[D[-D[[i,j]],x]*D[Xcomp[[j]],x],{j,1,3}], {i, 1, 3}]+Sum[-D[[i,j]]*D[Xcomp[[j]],{x,2}],{j,1,3}]=NeumannValue[0,x==xend], {i, 1, 3}];
bound1FEM =
Table[DirichletCondition[Xcomp[t, x][[i]] == list[[i]],
x == xbegin], {i, 1, 3}];
I fixed the problem. I happy that it works, but I do not understand yet why it fixed the problem. Guide from Wolfram recommends first form, rather than second.