On Wed, Dec 8, 2010 at 9:53 AM, robert <[email protected]> wrote: > Hi, > > I have tried to modify example code 9. I am using a starting > concentration of one everywhere but want to set the condition for one > (or more) boundary to zero. > > for (unsigned int qp=0; qp<qface.n_points(); qp++) > { > //const Number value = exact_solution > (qface_points[qp](0), > // > qface_points[qp](1), > // time); > > // RHS contribution > for (unsigned int i=0; i<psi.size(); i++) > Fe(i) += 0.;//penalty*JxW_face[qp]*value*psi[i][qp]; > > // Matrix contribution > for (unsigned int i=0; i<psi.size(); i++) > for (unsigned int j=0; j<psi.size(); j++) > Ke(i,j) += > penalty*JxW_face[qp]*psi[i][qp]*psi[j][qp]; > } > > > Well, at the boundaries I get the values which I want to have. However, > I also get strong oscillations. > I have also tried to just put the initial concentration to zero for a > small region within the domain. Again, I got oscillations (but less > intense). > What am I doing wrong here?
You got oscillations... after a solve? In the initial data? -- John ------------------------------------------------------------------------------ This SF Dev2Dev email is sponsored by: WikiLeaks The End of the Free Internet http://p.sf.net/sfu/therealnews-com _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
