The PDE, Boundary Elements interface (), found under the Mathematics>PDE Interfaces branch () when adding an interface, solves scalar stationary source-free PDEs in an unbounded domain.

When you add this interface, these default nodes are also added to the Model Builder: PDE, Boundary Elements; Zero Flux; and Initial Values. Then, from the Physics toolbar, add other nodes that implement, for example, boundary conditions. You can also right-click the main PDE, Boundary Elements node to select features from the context menu.

Use the PDE, Boundary Elements for a first-order scalar Laplace equation. Assuming a dependent variable u, these problems take the form

c is the only coefficient. It must be constant within each modeling region.

The plots can plot the boundary element field (default name: pdebe.u) and the normal boundary flux (default name: pdebe.bemflux) on boundaries and also the boundary element field in all domains and voids.

The Label is the default physics interface name.

The Name is used primarily as a scope prefix for variables defined by the physics interface. Refer to such physics interface variables in expressions using the pattern <name>.<variable_name>. In order to distinguish between variables belonging to different physics interfaces, the name string must be unique. Only letters, numbers, and underscores (_) are permitted in the Name field. The first character must be a letter.

The default Name (for the first PDE, Boundary Elements interface in the model) is pdebe.

Define the domains and voids in which the physics interface should be active. The list can include all solid domains and finite voids. By default, it also contains the infinite void. From the Selection list, you can choose All domains, All voids, or All domains and voids. All voids include all finite voids and the infinite void.

Select the Enable physics symbols check box to display the symmetry lines or planes in the geometry.

In 2D, choose a symmetry option from the Symmetry in line perpendicular to x and Symmetry in line perpendicular to y lists. In 3D; choose a symmetry option from the Symmetry in the yz plane, Symmetry in the xz plane, or Symmetry in the xy plane lists. Choose one of the following options:

Select a Dependent variable quantity. Dimensionless (1) is the default and there are many other options to choose from. Select a Source term quantity. None is the default and there are many other options to choose from. Enter a Unit.

For the condition, choose None for no condition or Laplace equation (the default), from the Condition type list.

For the Laplace equation, choose to specify either a Total flux through boundary (the default) or an Asymptotic value at infinity (3D only), or Value at reference distance (2D only) from the Condition at infinity list. Depending on the condition, enter a value in the Total flux through boundary, Asymptotic value, or Reference value field to define the condition at infinity.

To display this section, click the Show button () and select Advanced Physics Options.

The Use far-field approximation check box is selected by default in order to accelerate the solution process. If the check box is cleared, the solution will be slightly more accurate but the computational time and memory consumption may become prohibitively high.

The Approximation type can be either ACA+ or ACA. These alternatives correspond to two different versions of the adaptive-cross-approximation (ACA) method, which is a fast matrix multiplication method based on far-field approximations.

When the Use far-field approximation check box is selected, the Stationary Solver step in the study creates an octree structure, which is a tree data structure that divides the model into 2-by-2-by-2 blocks recursively until each of the smallest boxes contains at most the number of degrees of freedom specified in the Box size splitting limit field. The Far-field minimum relative distance decides if the interaction between two boxes occurs in the near field or far field. Boxes that interact in the near field at one level may also interact with a far-field approximation at a level with smaller boxes. Boxes at the smallest level that do not fulfill the far-field minimum relative distance criteria are considered to interact in the near field. For such boxes the system matrix that is defined by the integral equation and the elements in the two boxes are explicitly computed using no approximation. For two boxes that have been classified as interacting in the far field, an approximation of the resulting matrix is computed. The algorithm for computing the approximation rewrites the matrix defined by the interaction between the two boxes using a low-rank matrix approximation. The matrix rank for the approximation is chosen so that the relative error between the approximation and the actual matrix is estimated to be smaller than the value in the Relative tolerance field.

The ACA+ and ACA algorithms differ in the implementation of the fast matrix multiplication. The ACA+ algorithm is more robust but slightly more computationally expensive as compared to the ACA algorithm. If the Use SVD compression check box is selected, then after the ACA+ or ACA far-field approximation, a further approximation and data compression is made based on a singular value decomposition (SVD) algorithm. This additional data compression step reduces memory usage but increases computation time.

The damping parameter is related to the near field matrix that the preconditioner sees. By default, the Use damping check box is selected. The default value in the Damping parameter field is 1, and increasing this value increases the numerical damping. The parameter may impact the convergence of the solver but does not change the solution it converges to.

To display this section, click the Show button () and select Discretization.

The quadrature settings are by default set to Automatic. This means that the quadrature integration order values will follow the element order selection in the Discretization section. Higher element orders automatically generate higher values for the quadrature integration orders.

The quadrature integration order settings determine the level of accuracy of the computation of double integrals over the singular integral kernels that need to be evaluated for solving the BEM integral equations. For contributions from overlapping mesh element pairs (identical elements, elements sharing common edge, or elements with common vertex), specialized regularization transformations are applied before proceeding with numerical integration. For evaluation of contributions from nonoverlapping mesh element pairs, standard numerical integration is used. Due to the presence of singularities in the double integrals, higher integration order is needed for close mesh element pairs than for distant mesh element pairs. The following integration order settings are available when you select Manual:

Integration order, distant elements: Integration order for evaluating standard double integral BEM contributions from mesh pairs separated by a large distance relative to their size.

Integration order, close elements: Integration order for evaluating standard double integral BEM contributions from mesh pairs separated by a short distance relative to their size.

Integration order, elements with common vertex: Integration order for evaluating regularized double integral BEM contributions from mesh pairs with a common vertex.

Integration order, elements with common edge: Integration order for evaluating regularized double integral BEM contributions from mesh pairs with a common edge.

Integration order, same element: Integration order for evaluating regularized double integral BEM contributions from identical mesh pairs.

Integration order, weak contribution: Integration order for evaluating standard single integral BEM contributions. This integration does not contain any singularities.

To display this section, click the Show button () and select Discretization.

From the Dependent variable/Normal boundary flux list, choose from predefined options for the boundary element discretization order for the dependent variable and the normal boundary flux, respectively. The predefined options represent the suitable combinations of element orders such as Quadratic/Linear (the default).

The settings under Value types when using splitting of complex variables are important for sensitivity and optimization computations. See the description of the built-in operators fsens and fsensimag.

Define the name of the dependent variable. The default name is u.