GSvit documentation

open source FDTD solver with GPU support

User Tools

Site Tools


fdtd:boundary_conditions

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
fdtd:boundary_conditions [2018/01/24 08:14]
127.0.0.1 external edit
fdtd:boundary_conditions [2018/01/29 23:47]
pklapetek
Line 5: Line 5:
  
 There are two absorbing boundary conditions implemented in GSvit. Liao boundary conditions is highly effective second order boundary condition, using space and time derivatives of the values at the computational volume boundary. Convolutional perfectly matched layer (CPML) is one of the best available boundary conditions for FDTD at present, providing a layer of graded optical properties and stretched coordinates that lead to absorption of an electromagnetic wave with almost no backreflection. There are two absorbing boundary conditions implemented in GSvit. Liao boundary conditions is highly effective second order boundary condition, using space and time derivatives of the values at the computational volume boundary. Convolutional perfectly matched layer (CPML) is one of the best available boundary conditions for FDTD at present, providing a layer of graded optical properties and stretched coordinates that lead to absorption of an electromagnetic wave with almost no backreflection.
 +
 +----
 +//
 +{{ :​samples:​img_boundary_liao.png?​100|}}
 +Sample parameter file: {{samples:​boundary_liao.tar.gz|2nd order boundary condition}}.
 +\\
 +A 200x200x200 computational domain with single point source
 +//
 +----
 +//
 +{{ :​samples:​img_boundary_cpml.png?​100|}}
 +Sample parameter file: {{samples:​boundary_cpml.tar.gz|CPML boundary condition}}.
 +\\
 +A 200x200x200 computational domain with single point source.
 +//
 +----
  
 To compare the efficiency of Liao boundary condition and CPML we performed a simulation of a short pulse traveling in the computational volume, showing the backreflection from one side. As seen below, CPML provides much smaller reflection. The difference between both approaches is even more pronounced if we consider waves that are incident at high angles with respect to computational volume boundary normal. To compare the efficiency of Liao boundary condition and CPML we performed a simulation of a short pulse traveling in the computational volume, showing the backreflection from one side. As seen below, CPML provides much smaller reflection. The difference between both approaches is even more pronounced if we consider waves that are incident at high angles with respect to computational volume boundary normal.
  
 {{ :​fdtd:​cpmlvsliao.png?​400 |}} {{ :​fdtd:​cpmlvsliao.png?​400 |}}
 +
  
 A special case is periodic boundary condition. Here the fields from one side of computational domain are copied to the other side so the space becomes periodic. Note that this approach is not valid for some cases, e.g. for planar wave illuminating the computational volume from an angle in TSF source formulation. To overcome this a special algorithm is necessary, not yet implemented in GSvit. A special case is periodic boundary condition. Here the fields from one side of computational domain are copied to the other side so the space becomes periodic. Note that this approach is not valid for some cases, e.g. for planar wave illuminating the computational volume from an angle in TSF source formulation. To overcome this a special algorithm is necessary, not yet implemented in GSvit.
  
 The use of boundary conditions in GSvit can be independent for every face - so it is possible to use nonabsorbing boundary condition for some faces of the computational domain only or to use different boundary conditions for different faces. The use of boundary conditions in GSvit can be independent for every face - so it is possible to use nonabsorbing boundary condition for some faces of the computational domain only or to use different boundary conditions for different faces.
fdtd/boundary_conditions.txt ยท Last modified: 2018/01/29 23:47 by pklapetek