Just add the appropriate boundary integral term in a similar way as you do a Dirichlet BC, e.g. something like
Fe(i) += JxW_face[qp] * phi_face[i][qp] * g; where g is the Neumann BC. Boundary IDs are useful if you only want to impose the BC on a part of the boundary... Note that you don't include a penalty term in this case. - Dave Ted Kord wrote: > Hi > > How do you apply a non-homogeneous Neumann boundary condition? All the > examples only deal with Dirichlet B.Cs. > > Thanks in advance. > > Ted Kord > ------------------------------------------------------------------------------ > Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day > trial. Simplify your report design, integration and deployment - and focus on > what you do best, core application coding. Discover what's new with > Crystal Reports now. http://p.sf.net/sfu/bobj-july > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users > ------------------------------------------------------------------------------ Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day trial. Simplify your report design, integration and deployment - and focus on what you do best, core application coding. Discover what's new with Crystal Reports now. http://p.sf.net/sfu/bobj-july _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
