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.

density change due to deformation

Please login with a confirmed email address before reporting spam

I have multiple custom PDE's coupled that are modeling wood pyrolysis. They are all solving for the densities of the constituents of the pyrolyzing material. The equations them selves are basically just reaction equations. I now have also coupled a structural mechanics node with a LinElMat with thermal expansion. I would like the densities to change with the deformation but i cant find a way to do so during the simulation. In post thats not a problem since you can just multiply the results with the deformation gradient (everything is calculated in the spatial frame). But since the reaction rates depend on the densities that seems unsatisfying to me. Any help is greatly appreciated, thanks.


2 Replies Last Post 25 ago 2021, 03:27 GMT-4
Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 3 years ago 25 ago 2021, 02:53 GMT-4

You can use the multiplication with determinant of the deformation gradient (solid.J) also in your own equations.

This said, having changing densities is always tricky. Is mass conservation to be valid in your calculation domain, or are you losing mass to the surrondings? Note also that the thermal coefficent of expansion in general is not imdependent of density, since it describes the change in volume of a piece of material containg a certain number of molecules.

-------------------
Henrik Sönnerlind
COMSOL
You can use the multiplication with determinant of the deformation gradient (solid.J) also in your own equations. This said, having changing densities is always tricky. Is mass conservation to be valid in your calculation domain, or are you losing mass to the surrondings? Note also that the thermal coefficent of expansion in general is not imdependent of density, since it describes the change in volume of a piece of material containg a certain number of molecules.

Please login with a confirmed email address before reporting spam

Posted: 3 years ago 25 ago 2021, 03:27 GMT-4
Updated: 3 years ago 25 ago 2021, 03:28 GMT-4

Thanks for the fast reply!

I assume you mean the derivative of solid.J? Because thats what i came up with as well (as an absorption coefficient) and it seems to be working.

As I have set it up right now I shouldnt loose any mass but when I integrate all densities over the domain the total mass does change and I dont know why.

The density dependence of the thermal expansion shouldnt be too much of a concern in the model since the pointwise transitions happen relatively fast.

Thanks for the fast reply! I assume you mean the derivative of solid.J? Because thats what i came up with as well (as an absorption coefficient) and it seems to be working. As I have set it up right now I shouldnt loose any mass but when I integrate all densities over the domain the total mass does change and I dont know why. The density dependence of the thermal expansion shouldnt be too much of a concern in the model since the pointwise transitions happen relatively fast.

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.