Dear Mates,

I was running a simulation, in which I applied acoustic pressure from the top and set the peripheral boundaries on Perfectly Matched Boundary for considering the domain as free air between the sound source and the membrane which is supposed to be excited by acoustic waves. I faced with an error as below:

- Feature: Stationary Solver 1 (sol1/s1)

Division by zero.

- Function: /

Failed to evaluate variable.

- Variable: comp1.acpr.pmb1.thiFactor

- Defined as: 1/(1.7320508075688771932*(((comp1.acpr.nx*comp1.acpr.pmb1.ndirx)+(comp1.acpr.ny*comp1.acpr.pmb1.ndiry))+(comp1.acpr.nz*comp1.acpr.pmb1.ndirz)))

Failed to evaluate expression.

- Expression: comp1.acpr.pmb1.thiFactor

Failed to evaluate variable.

- Variable: comp1.acpr.pmb1.thiFactor

- Defined as: 1/(1.7320508075688771932*(((comp1.acpr.nx*comp1.acpr.pmb1.ndirx)+(comp1.acpr.ny*comp1.acpr.pmb1.ndiry))+(comp1.acpr.nz*comp1.acpr.pmb1.ndirz)))

Failed to evaluate expression.

- Expression: comp1.acpr.pmb1.thiFactor

Failed to evaluate expression.

- Expression: comp1.atonly(comp1.acpr.pmb1.thiFactor)

Failed to evaluate variable.

- Variable: comp1.acpr.pmb1.xn

- Defined as: ((comp1.acpr.pmb1.thiFactor*comp1.acpr.pmb1.scaleFactor_p)*((xd1/unit_m_cf)^comp1.acpr.pmb1.gamma_p))*(1-i)

Failed to evaluate variable.

- Variable: comp1.acpr.pmb1.T23

- Defined as: d(comp1.acpr.pmb1.coordvary,xi2)+((d(comp1.acpr.pmb1.ndiry,xi2)*comp1.acpr.pmb1.typicalWavelength)*comp1.acpr.pmb1.xn)

Failed to evaluate variable.

- Variable: comp1.acpr.pmb1.invT21

- Defined as: ((comp1.acpr.pmb1.T23*comp1.acpr.pmb1.T31)-(comp1.acpr.pmb1.T33*comp1.acpr.pmb1.T21))/((((((comp1.acpr.pmb1.T33*(comp1.acpr.pmb1.T22*comp1.acpr.pmb1.T11))+((comp1.acpr.pmb1.T23*comp1.acpr.pmb1.T12)*comp1.acpr.pmb1.T31))+(comp1.acpr.pmb1.T32*(comp1.acpr.pmb1.T13*comp1.acpr.pmb1.T21)))-(comp1.acpr.pmb1.T32*(comp1.acpr.pmb1.T23*comp1.acpr.pmb1.T11)))-(comp1.acpr.pmb1.T33*(comp1.acpr.pmb1.T12*comp1.acpr.pmb1.T21)))-(comp1.acpr.pmb1.T31*(comp1.acpr.pmb1.T22*comp1.acpr.pmb1.T13)))

Failed to evaluate expression.

- Expression: comp1.acpr.pmb1.invT21

Failed to evaluate Jacobian of expression.

- Expression: ((comp1.acpr.pmb1.invT11*comp1.acpr.pmb1.pT1)+(comp1.acpr.pmb1.invT21*comp1.acpr.pmb1.pT2))+(comp1.acpr.pmb1.invT31*comp1.acpr.pmb1.pT3)

Failed to evaluate Jacobian of expression.

- Expression: ((-test(comp1.acpr.pmb1.invT11*comp1.acpr.pmb1.pT1+comp1.acpr.pmb1.invT21*comp1.acpr.pmb1.pT2+comp1.acpr.pmb1.invT31*comp1.acpr.pmb1.pT3)*(comp1.acpr.pmb1.invT11*comp1.acpr.pmb1.pT1+comp1.acpr.pmb1.invT21*comp1.acpr.pmb1.pT2+comp1.acpr.pmb1.invT31*comp1.acpr.pmb1.pT3)*comp1.atonly(comp1.acpr.delta)*comp1.acpr.pmb1.detT*comp1.atonly(comp1.acpr.pmb1.detSign)/(comp1.atonly(comp1.acpr.rho_c)*dvol))*(acpr_pmb1_xdim.atonly(dvol)))*(dvol_spatial)

When the geometry was set to two cylinders over each other, the model worked correctly, but this error turned out when I changed the geometry to a trapezoidal shape. Do you have any idea how to fix it?