Limiting variables in Monolithic solver

Please login with a confirmed email address before reporting spam

Hello COMSOL Community,

In the typical segregated approach, it’s possible to limit a variable between an upper and lower bound. In the case of, the phase field variable needs to remain between 0 and 1. However, when using the monolithic or fully coupled solver, I haven’t found a way to enforce these limits.

My question is: How can I ensure that the phase field variable stays within the range of 0 to 1 when calculating energy, which depends on this variable, while using the monolithic solver?

I would appreciate any insights or suggestions on how to approach this issue.

Thank you!


3 Replies Last Post 2024年9月2日 GMT+2 09:42
Adam Wahlsten COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 6 days ago 2024年8月28日 GMT+2 10:48

Hi Aida,

If possible, it is usually better to enforce the limit on the phase field from the governing equations instead of setting a limit in the solver, since the latter will disregard the equations where your bound is violated. This means that you can get an arbitrarily large residual.

In version 6.2, with the Phase Field in Solids interface, you can use the Bounds subnode to a Phase Field Model node to specify these kind of inequality constraints.

If you are working with a phase field damage model like AT1, which is unbounded, an alternative approach to adding inequality constraints is to analytically derive the minimum energy that is necessary for a zero solution and add this as an additional source term in the phase field equation.

Best regards,
Adam Wahlsten
COMSOL

Hi Aida, If possible, it is usually better to enforce the limit on the phase field from the governing equations instead of setting a limit in the solver, since the latter will disregard the equations where your bound is violated. This means that you can get an arbitrarily large residual. In version 6.2, with the Phase Field in Solids interface, you can use the Bounds subnode to a Phase Field Model node to specify these kind of inequality constraints. If you are working with a phase field damage model like AT1, which is unbounded, an alternative approach to adding inequality constraints is to analytically derive the minimum energy that is necessary for a zero solution and add this as an additional source term in the phase field equation. Best regards, Adam Wahlsten COMSOL

Please login with a confirmed email address before reporting spam

Posted: 4 days ago 2024年8月30日 GMT+2 16:38

Thank you for your response, Mr Adam
However, I've noticed that the "Weak Contribution" sub-node doesn't exist under the " Phase Field in Solids" node. How can I resolve this issue?

Thank you for your response, Mr Adam However, I've noticed that the "Weak Contribution" sub-node doesn't exist under the " Phase Field in Solids" node. How can I resolve this issue?

Magnus Ringh COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 1 day ago 2024年9月2日 GMT+2 09:42

Hi Aida,

Could it be that you have not yet activated theEquation-Based Contributionscheckbox in theShow More Optionswindow (which you open from the top of theModel Builderwindow). Then you should be able to addWeak Contributionsubnodes.

Best regards,

Magnus

Hi Aida, Could it be that you have not yet activated the **Equation-Based Contributions** checkbox in the **Show More Options** window (which you open from the top of the **Model Builder** window). Then you should be able to add **Weak Contribution** subnodes. Best regards, Magnus

Reply

Please read thediscussion forum rulesbefore posting.

Pleaselog into 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 theSupport Centerfor a more comprehensive response from the Technical Support team.

Baidu
map