Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

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.

How is Dirichlet boundary conditions implemented in Comsol?

Please login with a confirmed email address before reporting spam

So I'm building a model with the Acoustics Module trying to model a piezoelectric disc, and when documenting the model the question of how the Dirichlet boundary conditions like voltage on either electrode or prescribed displacements is implemented came up.

The manual seems to be quite vague on this point, like is it using the method of Lagrange multipliers, or is it modifying the matrices directly to enforce the conditions, or something completely different?
So maybe someone here would know anything about it? =)

Thanks in advance,
Andreas Hagen

2 Replies Last Post 7 avr. 2017, 09:22 UTC−4
Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 7 avr. 2017, 07:24 UTC−4
Hi,

Most Dirichlet conditions (there are hundreds of them) are as a default implemented in an 'exact' (removing the equations) manner. This is what is called 'pointwise constraint' in our terminology.

In most cases you can choose a Lagrange multiplier approach instead. This is the 'Use weak constraints' check box in the settings for the individual constraint. This option is shown if you have 'Advanced Physics Options' on.

There are also a few cases (particularly in Multibody Dynamics) when a penalty implementation is available.

Finally, you can always inspect what is actually used in Equation View.

Regards,
Henrik
Hi, Most Dirichlet conditions (there are hundreds of them) are as a default implemented in an 'exact' (removing the equations) manner. This is what is called 'pointwise constraint' in our terminology. In most cases you can choose a Lagrange multiplier approach instead. This is the 'Use weak constraints' check box in the settings for the individual constraint. This option is shown if you have 'Advanced Physics Options' on. There are also a few cases (particularly in Multibody Dynamics) when a penalty implementation is available. Finally, you can always inspect what is actually used in Equation View. Regards, Henrik

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 7 avr. 2017, 09:22 UTC−4
Thanks for the reply =)

So how I read the documentation for the pointwise constraint, it looks to use Lagrange multipliers implicitly, but instead of leaving them in the final matrix, it just eliminates them from the matrix straight away before solving.

Is that correct?

Cheers
Andreas Hagen
Thanks for the reply =) So how I read the documentation for the pointwise constraint, it looks to use Lagrange multipliers implicitly, but instead of leaving them in the final matrix, it just eliminates them from the matrix straight away before solving. Is that correct? Cheers Andreas Hagen

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.