Hello,
In Steps 13 and 14, for primal problem, the boundary condition is calculate
based on an analytical solution, therefore it is easy to set different
boundary indicators and apply the BC on different boundary indicators.
In my problem I have a cylinder (boundary indicator of R_in is 0 and
boundary indicator of R_out is 1).
How can I implement two different BC on these BCs?
for example on R_in the displacement U_x = x + y and on R_out displacement
U_x = x^2 + y.
I don't want to destroy the extensibility of the code, therefore I would
like to modify template <int dim> double
CurvedRidges<dim>::BoundaryValues::value (const Point<dim> &p, const
unsigned int /*component*/) const {...} part NOT the
Solver<dim>::assemble_linear_system part.
I was thinking of passing boundary indicator to the boundary function or
something similar; like:
VectorTools::interpolate_boundary_values (dof_handler,
0,
**boundary_values(0),*
boundary_value_map);
and then relate the BC in CurvedRidges<dim>::BoundaryValues::value(...) to
boundary indicator, but I am not sure if it is possible and which parts
should I modify.
Thanks
Elena
--
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.
