[Top][All Lists]

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [Getfem-users] Jump conditions on a surface.

From: Yves Renard
Subject: Re: [Getfem-users] Jump conditions on a surface.
Date: Tue, 15 May 2007 16:20:44 +0200
User-agent: KMail/1.9.5

Dear Ronan,

You can use the level set strategy to describe S and define a finite element 
method which is discontinuous across this interface. But, there is in fact a 
little difficulty to compute the integral you need because on the surface 
itself, the value on the interface of the base function of the xfem (the 
mesh_fem_level_set object) is not guaranty. This is indeed a problem that we 
want to solve in a near future because we need to prescribe non elementary 
interface limit condition (such as contact condition or jump condition as you 

If the interface is static and the mesh is conformal to this interface, there 
is a tool to partition a finite element method into two with regard to this 
interface. You can partition a fem into two or more parts with the method

mf_u.set_dof_partition(cv, i)

where cv is the element number and i is the zone number. Doing so, the dof 
between two zones will not be identified and you can represent a 
discontinuity between these two zones. In order to compute the integral you 
need, you have to define two mesh regions which will describe one side of the 
surface S and the other one respectively. Then you can integrate the terms 
you need with the usual strategy.


Le mardi 15 mai 2007 09:16, vous avez écrit :
> Dear getfem experts,
> I would like to solve a (relatively simple) problem with jump conditions
> on an surface S.
> More precisely, if O denotes the domain, the variational form is the
> following :
> \int_O \sigma \grad u \grad v dx + \int_S \alpha [u][v] ds = <source terms>
> where [u] is the jump of u over S.
> Is there some elegant solution to implement this?
> If I have to use the level set strategy (NB: the surface is static and I
> can easily include it in the mesh), how can I translate the \int_S
> \alpha [u][v] ds integration?
> I will be glad to supply you further details.
> Thank you in advance for your answer,
> Best regards,
> Ronan Perrussel

  Yves Renard (address@hidden)       tel : (33)
  Pole de Mathematiques, INSA de Lyon          fax : (33)
  Institut Camille Jordan - CNRS UMR 5208
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE

reply via email to

[Prev in Thread] Current Thread [Next in Thread]