From: Bao K. <pae...@gm...> - 2012-01-29 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] - // mass-matrix 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] - // mass-matrix 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 |
From: John P. <jwp...@gm...> - 2012-01-29 21:27:57
|
On Sun, Jan 29, 2012 at 12:47 PM, Bao Kai <pae...@gm...> wrote: > 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. > theta*dt*(U*grad_u)*phi[i][qp]); // Newton term The Newton method of Example 13 is atypical: rather than solving for the "update" du_{k+1} = u_{k+1} - u_{k}, it solves for u_{k+1} directly. So the "standard" Newton iteration: J_{k} du_{k+1} = -F_{k} is instead written as: J_{k} u_{k+1} = -F_{k} + J_{k} u_{k} This modified RHS can be computed by hand; it involves the term you asked about above. -- John |
From: Bao K. <pae...@gm...> - 2012-01-30 17:59:18
|
Hi, John, Thank you for your helpful response. Since the lid driven cavity problem is actually a steady problem, it seems that always the steady result will be obtained after the non-linear iteration, which also means we get the final result in the first time step. Is it right? Regards, Kai On Mon, Jan 30, 2012 at 12:27 AM, John Peterson <jwp...@gm...> wrote: > On Sun, Jan 29, 2012 at 12:47 PM, Bao Kai <pae...@gm...> wrote: >> 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. > >> theta*dt*(U*grad_u)*phi[i][qp]); // Newton term > > The Newton method of Example 13 is atypical: rather than solving for > the "update" du_{k+1} = u_{k+1} - u_{k}, it solves for u_{k+1} > directly. > > So the "standard" Newton iteration: > > J_{k} du_{k+1} = -F_{k} > > is instead written as: > > J_{k} u_{k+1} = -F_{k} + J_{k} u_{k} > > This modified RHS can be computed by hand; it involves the term you > asked about above. > > -- > John |
From: John P. <jwp...@gm...> - 2012-01-30 18:04:29
|
On Mon, Jan 30, 2012 at 10:59 AM, Bao Kai <pae...@gm...> wrote: > Hi, John, > > Thank you for your helpful response. > > Since the lid driven cavity problem is actually a steady problem, it > seems that always the steady result will be obtained after the > non-linear iteration, which also means we get the final result in the > first time step. Is it right? Not in general, no. The lid driven cavity has a steady state solution (for some Re) but we solve it in ex13 using a transient algorithm and time derivative terms are present in the PDE. The parameters (delta_t and Re) used in ex13 cause it to converge to steady state quite quickly, but certainly not in 1 timestep. -- John |