Menu

Discussion Forum

New Discussion
Filtrer

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.

Limiting variables in Monolithic solver

Posted 27 août 2024, 18:57 UTC−4 Studies & Solvers, Structural Mechanics, Material Models 4 Replies

Aida Bouaziz
Send Private Message Flag post as spam

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!

4 Replies Last Post 4 sept. 2024, 05:29 UTC−4
Adam Wahlsten COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 1 year ago 28 août 2024, 04:48 UTC−4

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
Aida Bouaziz
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 year ago 30 août 2024, 10:38 UTC−4

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 year ago 2 sept. 2024, 03:42 UTC−4

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

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
Aida Bouaziz
Send Private Message Flag post as spam

Please login with a confirmed email address before reporting spam

Posted: 1 year ago 4 sept. 2024, 05:29 UTC−4

Hi Magnus,

Thank you for the suggestion! I tried to add the Weak Contribution subnodes as needed, and it worked perfectly. Thanks again for your help!

Best regards, Aida

Hi Magnus, Thank you for the suggestion! I tried to add the Weak Contribution subnodes as needed, and it worked perfectly. Thanks again for your help! Best regards, Aida

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.

Suggested Content