Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Errors of singularities on meshing after integrating the derivative of the solution

Please login with a confirmed email address before reporting spam

Hi,

We are calculating the optical force within metal scatterers in a PEC bounded structure 2D model. The model is bounded by PEC with PMLs on the left and right sides. We define the background field Eb = E0exp(-jk*x) in the y direction and solve for the scattered fields (Please see the attached mph file).

The problem here is that the size of the mesh element can have significant impact on the integral of the force density. Please see the attached jpeg file. The x-axis indicates the maximum mesh element size in the model, which is the wavelength (633nm) over a number # (You may change the parameter, mesh_d, in the mph to have the same figure). The surface integral 1 calculates the gradient term (integral of the derivative of the solutions) and the surface integral 2 gives the Poynting term. The total is equal to the sum of the gradient and Poynting terms. We can see that the gradient term switches from negative value to positive value when varying the maximum meshing size.

https://www.comsol.com/blogs/how-identify-resolve-singularities-model-meshing/

This article explains that some issues may happen when singularities exist. In the conclusion, point 3 says that it is acceptable when calculating the integral of some quantity including the derivative of the solutions. However, we still have significant errors when applying different meshing. We also removed the sharp corners by filleting the scatterers, but the problem remains, especially for the gradient term.

Do you have any solutions toward this issue? Thanks.



0 Replies Last Post 21 févr. 2018, 17:05 UTC−5
COMSOL Moderator

Hello Li-Fan Yang

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.

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.