Dear group members,
I am trying to apply the Neumann (force)
BC on the edge using the right_hand_side() function which I modified as
following:
void right_hand_side (const std::vector<Point<dim> > &points, // the
active cell is also imported when right hand side is called to induce some
Neumann BC
std::vector<Tensor<1, dim> > &values,
typename DoFHandler<dim>::active_cell_iterator &
cell_in_rhs)
for (unsigned int f=0; f < GeometryInfo<dim> :: faces_per_cell; ++f)
{
if (cell_in_rhs->face(f)->at_boundary() &&
cell_in_rhs->face(f)->boundary_id() == 3)
{
std::cout << "points.size(): "<<
points.size()<<" :"<<std::endl;
values[0][1] = 2.0;
values[2][1] = 2.0;
}
}
So currently I know that the quadrature points on the edge number 3 upon
which I shall assign values are "0" and "2" in my case. But I want to
generalize my code i.e. it detects the nearest or valid point near the edge
number 3.
Any help in this regard would be appreciated. I also want to use this
method later on to apply heat flux condition upon the same indexed edge in
future. Thank you!
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
