Dear Jean-François,
The simplest way of coupling a 1D to a 3D domain is by means of the
"Interpolate" function in the general weak form assembly language in
GetFEM. You can add an identity interpolation between the 1D mesh,
representing your lines, and the 3D mesh representing your volume, with
add_interpolate_transformation_from_expression("3Dto1D", *mesh3D*, *mesh1D*,
*"X"*)
then, in order to couple the two domains, you can use the interpolate name
"3Dto1D", or whatever you call it, in any assembly expression integrated on
your 1D domain, .
Best regards
Kostas
On Tue, Nov 10, 2020 at 4:29 PM Jean-François Barthélémy <
[email protected]> wrote:
> Dear Getfem users,
>
> I have a 3D mesh in which physical lines are identified (ie set of edges).
> These lines may be on the boundary or in the interior of the domain.
>
> I would like to build a model incorporating (linear and even nonlinear)
> terms on these lines (ie integration domain of dimension m.dim()-2). This
> could be used for example to model a concentrated line of current injection
> in an electric ohmic problem (rhs term) or a lineic domain of high
> conductivity (contributing then to the matrix term). Is there any way to do
> so in Getfem?
>
> Thank you in advance
>
> Best regards
>
> Jean-François Barthélémy
>
>
