Re: [Getfem-users] Dirichlet condition on curved boundary

[Yves Renard]

Dear Roman,

Yes, for the moment, you can only use the generalized Dirichlet condition and 
precompute the unit normal field or use the 
asm_normal_component_dirichlet_constraints and use some explicit matrix bricks 
to add manually the constraint to your problem.

I will check if it is easy to build the brick corresponding to a Dirichlet 
condition on the normal component for the new brick system. 

Yves.




On jeudi 22 septembre 2011, Roman Putanowicz wrote:
> Dear All,
> 
> So far I have had no problems with straight boundaries and Dirichlet
> conditions. However now I have a curved boundary and I would like
> to prescribe displacement  in the direction normal to the boundary.
> 
> In the old brick system there is a calass:
>   getfem::mdbrick_normal_component_Dirichlet
> which seems to be designed exactly for the above task.
> 
> I do not see however a similar class in the new brick system.
> 
> There is a function add_generalized_Dirichlet_condition to set
> a condition of a form H*u=r but using it for my task would
> require to precompute the matrix field of boundary normals.
> 
> There is a function
> asm_normal_component_dirichlet_constraints
> but using it gets me out of the getfem::model framework -- I mean
> I cannot use bricks and have to deal explicitly with he tangent matrix,
> rhs vector and solving the system of equation.
> 
> Is the above picture correct or am I missing something?
> What would be the best approach in such situation ?
> 
> Regards,
> 
> Roman


-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------