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.

Nonlinear term in the Coefficient Form PDE solver

Please login with a confirmed email address before reporting spam

I would like to solve the following equation

div ( u grad(u) ) = 0

The problem is I am trying to set the 'c' coefficient to be

c = u

But it doesn't work. The result will not converge. Can this be done in the PDE solver?

1 Reply Last Post 25 lug 2017, 15:33 GMT-4
Jeff Hiller COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 25 lug 2017, 15:33 GMT-4
Updated: 7 years ago 26 lug 2017, 08:56 GMT-4
There's no conceptual problem with setting up coefficients that depend on the solution itself. But obviously nonlinear problems require more caution than linear ones. With nonlinear equations, you have to be especially careful not to set up a problem that does not have any solution, or has several, for instance (The same concerns exist also with linear equations, of course, but such situations are much easier to detect and remedy). For nonlinear equations, it's also particularly important to set up suitable initial values for the variables that are being solved for. In your case, the default initial value (zero) is particularly poorly suited, since it results in a singular stiffness matrix.
If these pointers do not help, consider posting your model file.
Best,
Jeff
There's no conceptual problem with setting up coefficients that depend on the solution itself. But obviously nonlinear problems require more caution than linear ones. With nonlinear equations, you have to be especially careful not to set up a problem that does not have any solution, or has several, for instance (The same concerns exist also with linear equations, of course, but such situations are much easier to detect and remedy). For nonlinear equations, it's also particularly important to set up suitable initial values for the variables that are being solved for. In your case, the default initial value (zero) is particularly poorly suited, since it results in a singular stiffness matrix. If these pointers do not help, consider posting your model file. Best, Jeff

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.