Port
The Port boundary condition is used to excite and absorb acoustic waves that enter or leave waveguide structures, like small ducts or channels, in an acoustic model. The thermoviscous port formulation ensures that the non-trivial mode shapes of the acoustic velocity and thermal fields are captured correctly.
A given port condition supports one specific propagating mode. To provide the full acoustic description, combine several port conditions on the same boundary. Typically, only the plane wave mode is propagating in small structures where the thermoviscous representation is necessary. The port condition provides a superior nonreflecting or radiation condition for waveguides compared to a simple impedance condition or a perfectly matched layer (PML) configuration. The same port boundary condition feature should not be applied to several waveguide inlets/outlets. The port condition supports S-parameter (scattering parameter) calculation but it can also be used as a source to just excite a system.
The Port boundary condition exists for 3D, 2D, and 2D axisymmetric models.
On a given boundary, a combination of ports will define the total acoustic fields (sum of incident and outgoing pressure, temperature, and velocity waves) as
where the summation “i” is over all ports on the given boundary “bnd”, Sij is the scattering parameter, Ain is the amplitude of the incident field (at port “j”), and pi, ui, and Ti are the mode shape of the i-th port. The mode shape is normalized to have a unit maximum amplitude for the pressure pi. This definition means that the scattering parameter Sij defines the amplitude of mode i when a system is exited at port j (with mode j). This corresponds to a multi-mode expansion of the solution on the given port. The scattering parameters are automatically calculated when an acoustic model is set up with just one port exciting the system. To get the full scattering matrix The Port Sweep Functionality can be used.
Port Properties
Enter a unique Port name. Only nonnegative integer numbers can be used as Port name as it is used to define the elements of the S-parameter matrix. The numeric port names are also required for the port sweep functionality. The port name is automatically incremented by one every time a port condition is added.
Select a Type of port: User defined (the default), Numeric (plane wave mode), Circular (plane wave mode), or Slit (plane wave mode). Depending on the selection different options appear in the Port Mode Settings section (see below). Use the Circular (plane wave mode) for a port with a circular cross section in 3D or 2D axisymmetry and the Slit (plane wave mode) option in 2D. If the port has a different cross section then either of these, use the User defined option or the Numeric (plane wave mode) port.
Port Mode Settings
Depending on the option selected in the Type of port (see above):
For User defined enter user defined expressions for the Mode shape pn, un, Tn, and the Mode wave number kn (SI unit: rad/m). The modes shape will automatically be scaled before it is used in the port condition. Use the user defined option to enter a known analytical expression or to use the solution from The Thermoviscous Acoustics, Boundary Mode Interface. The solutions from the boundary mode analysis can be referenced using the withsol() operator.
The Numeric (plane wave mode) port options is used for waveguides of arbitrary cross sections. In this case the shape of the propagating plane wave mode (0,0) is solved on the port face assuming no-slip and isothermal wall conditions.
For this option it is necessary to make a few changes to the solver to get full performance. The port mode shape should be solved before the domain solution. Generate the default solver, then expand the Stationary Solver node, right click and add a Segregated solver. On the segregated node change the Termination technique to Iterations. Expand the segregated node, right click, and add a Segregated Step. The first segregated step node should only contain the port variables (psi, Psi_th, Psi_v, and vip) and the second should contain the remaining variables (p, u, T, and Sparam1). For the first step set the Linear solver to the default direct solver and for the second step select the direct solver or one of the iterative suggestions.
The Circular (plane wave mode) port option is used for waveguides of circular cross section. The analytical mode is a plane wave mode (0,0) given by a constant cross section pressure, no-slip condition for the velocity, and isothermal condition for the temperature. An example of the propagating mode shape in a cylindrical waveguide is seen below.
Figure 6-1: Plane wave mode for a circular duct of 1 mm in diameter at 250 Hz.
Select how the Circle radius of the cross section is defined, either Automatic (the default) or User defined. The latter option can for some geometry configurations increase the numerical precision of the computed mode.
The Slit (plane wave mode) port option only exists in 2D on a boundary. In 2D, the geometry is assumed infinite in the out-of-plane direction and represents a slit. The analytical mode is a plane wave mode (0,0) given by a constant cross section pressure, no-slip condition for the velocity, and isothermal condition for the temperature.
Incident Mode Settings
Activate if the given port is excited by an incident wave of the given mode shape. For the first Port condition added in a model, the Incident wave excitation at this port is set to On. For subsequent conditions added the excitation is set to Off per default. If more than one port in a model is excited, the S-parameter calculation is not performed.
When the Incident wave excitation at this port is set to On, then select how to define the incident wave. Select the Define incident wave: Amplitude (the only option).
For Amplitude enter the amplitude Ain (SI unit: Pa) of the incident wave. This is in general defined as the maximum amplitude for a given mode shape.
Enter the phase φin (SI unit: rad) of the incident wave. This phase contribution is multiplied with the amplitude defined by the above options. Note, that the Amplitude input can be a complex number.
Note, that when the Activate port sweep option is selected at the physics level, the options in the Incident Mode Settings section are deactivated. This is because this option automatically sends in a mode of unit amplitude, sweeping through one port at the time.
For the Circular and Slit options make sure to only select modes that are actually symmetric according to the symmetry planes.
When postprocessing, remember that absolute values like, for example, the outgoing power at port 1, ta.port1.P_out, needs to be multiplied with an appropriate factor. Multiplication with two if one symmetry plane is used etc.
Constraint Settings
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box.
Excluded Edges/Points
To display this section, click the Show More Options button () and select Advanced Physics Options in the Show More Options dialog box. See Suppressing Constraints on Lower Dimensions for details
The Port Sweep Functionality
The port sweep functionality is used to reconstruct the full scattering matrix Sij by automatically sweeping the port excitation through all the ports included in the model. When the port sweep is activated, the options in the Incident Mode Settings in the port conditions are deactivated and COMSOL controls which port is excited with an incident mode.
The port sweep functionality is activated at the main physics interface level by selecting Activate port sweep in the Port Sweep Settings section. Enter the Sweep parameter name, the default is PortName. Create a parameter with the same name under Global Definitions>Parameters 1. This is the name of the parameter to be used in a parametric sweep; it should represent the Port name integer values (defined when adding the port conditions). Add a parametric sweep study step and run the sweep over the PortName parameter with an integer number of values representing all the ports in the model. Once the model is solved, the full scattering matrix can be evaluated using the defined global variables ta.S11, ta.S21, ta.S12 etc. The transmission loss (TL) between two given ports is also computed, for example, the variable for the TL loss from port 1 to 2 is given by ta.TL_12.
Use the Global Matrix Evaluation under Derived Values to evaluate the full scattering matrix ta.S.
If only two ports are added to the thermoviscous model, COMSOL also automatically computes the transfer matrix of the system (variables ta.T11, ta.T12, ta.T21, ta.T22) and the impedance matrix of the system (ta.Z11, ta.Z12, ta.Z21, ta.Z22). These expressions are only true if plane wave modes are used. This is nearly the case in all configurations when working with microacoustic systems. For ports in thermoviscous acoustics, the Circular, Slit, and Numeric options are for plane waves only. Higher order modes can only be introduced with the User defined option. The transfer matrix representation is often used in electroacoustic modeling.
Wax Guard Acoustics: Transfer Matrix Computation: Application Library path Acoustics_Module/Tutorials,_Thermoviscous_Acoustics/wax_guard_acoustics