From: Bao Kai <paeanball@gm...>  20120129 19:47:55

Hi, I am just starting to lean Libmesh. I have a question on ex13. I want to know what the Newton term means in ex13. The code is listed as follows, for (unsigned int i=0; i<n_u_dofs; i++) { Fu(i) += JxW[qp]*(u_old*phi[i][qp]  // massmatrix term (1.theta)*dt*(U_old*grad_u_old)*phi[i][qp] + // convection term (1.theta)*dt*p_old*dphi[i][qp](0)  // pressure term on rhs (1.theta)*dt*(grad_u_old*dphi[i][qp]) + // diffusion term on rhs theta*dt*(U*grad_u)*phi[i][qp]); // Newton term Fv(i) += JxW[qp]*(v_old*phi[i][qp]  // massmatrix term (1.theta)*dt*(U_old*grad_v_old)*phi[i][qp] + // convection term (1.theta)*dt*p_old*dphi[i][qp](1)  // pressure term on rhs (1.theta)*dt*(grad_v_old*dphi[i][qp]) + // diffusion term on rhs theta*dt*(U*grad_v)*phi[i][qp]); // Newton term Thank you very much. Best Regards, Kai 