The Pressure Acoustics, Boundary Elements Interface
The Pressure Acoustics, Boundary Elements (pabe) interface (), found under the Pressure Acoustics branch () when adding a physics interface, is used to compute the pressure variation for the propagation of acoustic waves in gases and liquids. It is suited for all frequency-domain simulations with harmonic variations of the pressure field. The formulation is based on the boundary element method (BEM) and the interface is available in 2D and 3D. The physics interface solves Helmholtz equation for constant-valued material properties and uses the pressure as the dependent variable.
The interface is fully multiphysics enabled and can be coupled seamlessly with the physics interfaces that are based on the finite element method (FEM). This includes coupling to vibrating structures with the Acoustic–Structure Boundary multiphysics coupling, to FEM acoustic domains, with the Acoustic FEM–BEM Boundary multiphysics coupling, to thermoviscous acoustic domains with the Acoustic–Thermoviscous Acoustic Boundary multiphysics coupling, and to poroelastic waves domains with the Acoustic–Porous Boundary multiphysics coupling. This approach allows modeling in a FEM-BEM framework, using the strength of each formulation adequately. The BEM-based interface is especially well suited for radiation and scattering problems.
The advantage of the boundary element method is that only boundaries need to be meshed and the degrees of freedom (DOFs) solved for are restricted to the boundaries. This introduces some clear ease-of-use for handling complex geometries. However, the BEM technique results in fully populated or dense matrices that need dedicated numerical methods. The BEM method is so to speak more expensive per DOF than the FEM method but has fewer DOFs. Assembling and solving these can be very demanding. This means that when solving acoustic models of small and medium size, The Pressure Acoustics, Frequency Domain Interface will often be faster than solving the same problem with the BEM interface. The challenge for the FEM interface is to set up open boundaries, for example, using PMLs, in an efficient way. When the geometries are complex or two structures are far apart, large air domains need to be meshed. This costs a lot on the computational side as the frequency is increased.
For acoustically large models (problems that contain many wavelengths, at high frequency or for large domains) the stabilized formulation option (see Stabilization) ensures efficient convergence at the cost of some additional degrees of freedom. For low to medium frequencies (small to medium models), running without stabilization is more efficient. The stabilized formulation only gives a benefit in computing time for the acoustically large models.
Head and Torso HRTF Computation. Application Library path Acoustics_Module/Tutorials,_Pressure_Acoustics/head_torso_hrtf
Spherical Scatterer: BEM Benchmark. Application Library path Acoustics_Module/Verification_Examples/spherical_scatterer_bem_benchmark
Submarine Target Strength. Application Library path Acoustics_Module/Underwater_Acoustics/submarine_target_strength
The Bessel Panel. Application Library path Acoustics_Module/Tutorials,_Pressure_Acoustics/bessel_panel
The Loudspeaker Driver in a Vented Enclosure. Application Library path Acoustics_Module/Electroacoustic_Transducers/vented_loudspeaker_enclosure
The governing Helmholtz equation defined by the Pressure Acoustics, Boundary Element interface is given by:
(2-1)
where pt is the total acoustic pressure, keq is the wave number, ρc is the density, and cc is the speed of sound. The subscript “c” denotes that these can be complex-valued quantities in models with damping. The governing equations and boundary conditions are formulated using the total pressure pt with a so-called scattered field formulation. In the presence of a Background Pressure Field defining a background pressure wave pb (this could, for example, be a plane wave), the total acoustic pressure pt is the sum of the pressure solved for p (which is then equal to the scattered pressure ps) and the background pressure wave: pt = p+pb. The equations then contain the information about both the scattered field and the background pressure field.
When this physics interface is added, these default nodes are also added to the Model BuilderPressure Acoustics, Sound Hard Boundary (Wall), and Initial Values.
Then, from the Physics toolbar, add other nodes that implement boundary conditions. You can also right-click Pressure Acoustics, Boundary Elements to select physics features from the context menu. Infinite conditions like a symmetry plane or an infinite sound hard boundary are defined in the Symmetry/Infinite Boundary Condition section.
Physics Nodes — Equation Section in the COMSOL Multiphysics Reference Manual
Settings
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 physics interface in the model) is pabe.
Equation
Expand the Equation section to see the equations solved for with the Equation form specified. The default selection is Equation form is set to Study controlled. The available studies are selected under Show equations assuming. If Frequency domain is selected as the Equation form the Frequency of the study can be either taken From solver or User defined.
Physics Symbols
Select the Enable physics symbols check box to display the infinite boundary condition lines or planes in the geometry.
Sound Pressure Level Settings
The zero level on the dB scale varies with the type of fluid. That value is a reference pressure that corresponds to 0 dB. This variable occurs in calculations of the sound pressure level Lp based on the root mean square (rms) pressure prms, such that
where pref is the reference pressure and the star (*) represents the complex conjugate. This is an expression valid for the case of harmonically time-varying acoustic pressure p.
Select a Reference pressure for the sound pressure level based on the fluid type:
Use reference pressure for air to use a reference pressure of 2μPa (20·106 Pa).
Use reference pressure for water to use a reference pressure of 1 μPa (1·106 Pa).
User-defined reference pressure to enter a reference pressure pref, SPL (SI unit: Pa). The default value is the same as for air, 20 μPa.
Symmetry/Infinite Boundary Condition
In this section you can specify infinite boundary conditions for the boundary element problem. These are conditions that are not added on a boundary but apply to a line in 2D or a plane in 3D. The planes where the conditions are applied are visualized if the Enable physics symbols is turned on (the default) in the Physics Symbols section.
For pressure acoustics a symmetry condition is mathematically equivalent to a sound hard boundary, and an antisymmetry condition is equivalent to a sound soft boundary. In this way you can include a Symmetry/Infinite sound hard boundary or an Antisymmetry/Infinite sound soft boundary along “infinite lines” (in 2D) or “infinite planes” in (3D). Only one condition per Cartesian direction can be applied.
For the infinite conditions to be mathematically valid it is important that all sources and BEM boundaries are located on the same side of (or on) the infinite conditions. The condition can, for example, not be used as an infinite baffle with sound radiating through a hole. In this case, the propagation domain is on a different side than the source. The Exclude Boundary feature can be used to exclude boundaries from the BEM problem for advanced symmetry configurations.
Choose an option from the Condition for the x = x0 plane, Condition for the y = y0 plane, and Condition for the z = z0 plane lists (when applicable). In 2D, these are out-of-plane surfaces. Choose one of the following options:
Off, for no symmetry (the default)
Symmetric/Infinite sound hard boundary
Antisymmetric/Infinite sound soft boundary
Then enter the value for the plane location x0, y0, or z0 (the default is 0 m). This allows an offset of the infinite condition planes along the main coordinate axes.
There are important aspect on how results are visualized when symmetries are used:
First, when setting the Parameter Bounds in the grid dataset. Say the condition represents a symmetry plane then visualize the solution in the entire domain. On the other hand if the condition represents an infinite sound hard boundary, like the ground or an infinite baffle, then only visualize the solution above that plane by setting the Parameter Bounds in the Grid 3D or Grid 2D dataset in Results.
Far-Field Approximation
To display this section, click the Show More Options button () and select Advanced Physics Options.
These settings are used for matrix assembly and postprocessing. They allow characterization of interactions occurring in the boundary element method as near-field or far-field interactions. While the near-field interactions are represented explicitly, the far-field interactions can be represented in an approximate way. This approach results in considerable memory and performance improvements when used in combination with iterative solvers using matrix-free format or during postprocessing. The near-field part of the stiffness matrix is used as input by the Direct and Sparse Approximate Inverse preconditioners.
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.
Condition at Infinity
In this section, you specify the condition to apply at infinity for an unbounded problem. For the Helmholtz equation, solved here, choose to specify an Outgoing wave (the default) or an Incoming wave. In almost all cases the Outgoing wave option should be selected.
Quadrature
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
The quadrature settings control the details of the numerical integration performed when setting up the BEM problem.
Select Thin gap handling if narrow regions exist in the model. The Thin gap handling check box is used to increase the accuracy of the quadrature between almost-touching mesh element pairs, that is, pairs that are extremely close with respect to element size. The Threshold distance (relative to mesh size) controls the distance at which elements are considered almost-touching. A larger value should provide more accurate results, but at an increased computational cost. Under the Integration order settings, the field Integration order, very close elements (thin gap) becomes active (when Manual is selected), allowing a higher-order rule for such elements. In 3D, an additional variable transformation is applied to improve the accuracy even further.
The Integration orders 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. To see the used values or to set user-defined values, switch to Manual. Note that the values are reset to default when switching back and forth between the Automatic and Manual settings. The Integration order, very close elements (thin gap) is only editable when Thin gap handling is selected.
More details are found under Quadrature in The PDE, Boundary Elements Interface documentation in the COMSOL Multiphysics Reference Manual.
Stabilization
To display this section, click the Show More Options button () and select Stabilization in the Show More Options dialog box.
For acoustically large models (problems that contain many wavelengths, at high frequency or for large domains) enable the Stabilized formulation option to ensure efficient convergence at the cost of some additional degrees of freedom.
When Stabilized formulation is selected, a text field for the Stabilization parameter is enabled with the default value sqrt(abs(pabe.k[m])). This is a parameter that should scale inversely with the wavelength. The default gives good performance in most cases.
If the Stabilized formulation is enabled after a nonstabilized model has been solved (maybe showing slow convergence) it is important to reset the solver default. A different solver configuration is used for the stabilized BEM formulation.
For a model that uses the Stabilized formulation, see the Submarine Target Strength tutorial. Application Library path Acoustics_Module/Underwater_Acoustics/submarine_target_strength
Postprocessing Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
This setting controls the interpolation that is done close to boundaries (typically in the last mesh element) between variables on boundaries (these are the variables solved for in BEM based interfaces) and their value evaluated with the kernel in the domains. Numerically the kernel evaluation is prone to errors near the boundaries while the DOFs solved for are exact on the boundaries. Using an interpolation in the inner mesh element removes singularities in the results when postprocessing values in the domain near a boundary.
Enter an expression for the Interpolation distance from boundary. The entered value defines the distance beyond which the kernel evaluation is used. The default is 0.5*h which gives reasonable results. This default corresponds to half the mesh size h. If the value is set to 0 no interpolation is done.
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. 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.
Dependent Variables
This physics interface defines one dependent variable (field), the Pressure p. If required, edit the name, which changes both the field name and the dependent variable name. The new field name must not coincide with the name of a field of another type, or with a component name belonging to some other field.
To do a kernel evaluation of the BEM solution at a given coordinate use the at3_spatial() operator in 3D or the at2_spatial() operator in 2D. Both called with the minc argument, for mesh independent evaluations.
In a 3D model you can, for example, evaluate the sound pressure level in the point (x,y,z) = (1 m,0,0) by typing:
at3_spatial(1[m],0,0,pabe.Lp,’minc’).