May be it would be better to edit my question and remove the "cylinder" 
part, but I do not have access to modify the question.
Please consider my question for a general case of a geometry with multiple 
boundary indicators and different BCs on them.

Thanks.
Elena

On Tuesday, May 16, 2017 at 1:10:53 PM UTC+2, Elena Greco wrote:
>
> Dear Uwe
>
> Thanks for you prompt reply.
>
> But if I import a mesh from gmesh and for a not regular geometry I would 
> face a problem.
> Do you have any suggestion for this case?
>
> Thanks.
> Elena
>
> On Tuesday, May 16, 2017 at 12:50:36 PM UTC+2, Elena Greco wrote:
>>
>> 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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

Reply via email to