Divergence on piezoelectric domain

Please login with a confirmed email address before reporting spam

Hello to all,

I'm trying to create a 2D model of an ultrasonic NDT using Elastic Waves, Time Explicit with Electrostatics. I used the tutorial "Angle Beam Nondestructive Test" (https://www.comsol.com/model/angle-beam-nondestructive-testing-78811) as starting point for my own model, but I'm having some trouble with the Piezoelectric domain.

I assigned an edge as Electric Potential boundary condition and to create a pulse I inserted an Analytic function as voltage input, just like the tutorial. The simulation starts fine, but after a while, the electric potential in the piezoelectric domain grows and I can't figure it out why. I attached a plot of the average electric potential on the edge, in this image is visible the input pulse and how exponentially the voltage grows after some point.

I thought that it could be a mesh problem, so I tried to refine it, but that didn't work. Anyone faced a similar problem?

Thank you in advance,

Rodrigo Pires



3 Replies Last Post May 21, 2024, 9:20 a.m. EDT
Mark Cops COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 weeks ago May 17, 2024, 3:22 p.m. EDT

The graph is indicating that the solution is diverging. Since there can be various reasons for this, you would need to submit your .mph model file in order to get some specific suggestions for your model.

-Mark

The graph is indicating that the solution is diverging. Since there can be various reasons for this, you would need to submit your .mph model file in order to get some specific suggestions for your model. -Mark

Please login with a confirmed email address before reporting spam

Posted: 4 weeks ago May 21, 2024, 8:58 a.m. EDT

Hello Mark, thanks for your reply!

I’m attaching the .mph file to this topic. I’m sorry, but I had to exclude the mesh and result to get the right file size for uploading.

Thanks in advance, Rodrigo

Hello Mark, thanks for your reply! I’m attaching the .mph file to this topic. I’m sorry, but I had to exclude the mesh and result to get the right file size for uploading. Thanks in advance, Rodrigo


Mark Cops COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 4 weeks ago May 21, 2024, 9:20 a.m. EDT
Updated: 4 weeks ago May 21, 2024, 9:29 a.m. EDT

It looks like there is no damping in any of the solids, even the "Damping material". It could help to add some Rayleigh damping - similar to what was done in the angle_beam_ndt model. Secondly, both piezo domains currently have the same coordinate system used for the poling direction. Is this correct? It would seem that due to the different angle, each piezo would need a separate base vector system to define its poling direction relative to the mount angle.

It looks like there is no damping in any of the solids, even the "Damping material". It could help to add some Rayleigh damping - similar to what was done in the angle_beam_ndt model. Secondly, both piezo domains currently have the same coordinate system used for the poling direction. Is this correct? It would seem that due to the different angle, each piezo would need a separate base vector system to define its poling direction relative to the mount angle.

Reply

Please read the discussion forum rules before posting.

Please log in to post a reply.

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.