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.

What's the formulation of "mechanical energy flux" in COMSOL?

Please login with a confirmed email address before reporting spam

There are "mechanical energy flux," "complex mechanical energy flux"
which I can derive or plot in the result.
They are cool but still I want to know the definitions of them, how it calculates, and I did not find them in the user manuals.
Could somebody tell me? Thanks.

8 Replies Last Post 23 août 2016, 10:25 UTC−4
COMSOL Moderator

Hello Ting-Wei Liu

Your Discussion has gone 30 days without a reply. If you still need help with COMSOL and have an on-subscription license, please visit our Support Center for help.

If you do not hold an on-subscription license, you may find an answer in another Discussion or in the Knowledge Base.


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 23 mai 2013, 10:12 UTC−4
Hi,

Nobody answer you? because i'm also interested to kwow more about this indicator and to what is the formulation difference between «mechanical energy flux» and «complex mechanical energy flux».
Hi, Nobody answer you? because i'm also interested to kwow more about this indicator and to what is the formulation difference between «mechanical energy flux» and «complex mechanical energy flux».

Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 24 mai 2013, 02:29 UTC−4
Hi,

The mechanical energy flux is a vector formed by the multiplication of the stress tensor and the velocity vector. Mathematically:



The reason to have the minus sign in the definition is that if you put a pressure on an external boundary, and it moves in the direction of the load, then a positive power input in the direction of the load is obtained.

In time domain, the expression above is sufficient, but in frequency domain, there are two flavors:

The "Complex mechanical energy flux" is now the complex vector formed by the operation above, but with the velocity replaced by its complex conjugate.



The "Mechanical energy flux" in frequency domain is a real quantity; the cycle average of the "Complex mechanical energy flux"



In the appended screenshots you can see the definitions as displayed in Equation View in time and frequency domain.

Thanks for pointing out that this is not properly documented. It will be improved in upcoming versions.

Regards,
Henrik
Hi, The mechanical energy flux is a vector formed by the multiplication of the stress tensor and the velocity vector. Mathematically: [math] I_i =- \sigma_{ij} v_j [/math] The reason to have the minus sign in the definition is that if you put a pressure on an external boundary, and it moves in the direction of the load, then a positive power input in the direction of the load is obtained. In time domain, the expression above is sufficient, but in frequency domain, there are two flavors: The "Complex mechanical energy flux" is now the complex vector formed by the operation above, but with the velocity replaced by its complex conjugate. [math] I_i =- \sigma_{ij} v_j^* [/math] The "Mechanical energy flux" in frequency domain is a real quantity; the cycle average of the "Complex mechanical energy flux" [math] I_i =\frac{1}{2}Re(- \sigma_{ij} v_j^*) [/math] In the appended screenshots you can see the definitions as displayed in Equation View in time and frequency domain. Thanks for pointing out that this is not properly documented. It will be improved in upcoming versions. Regards, Henrik


Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 24 mai 2013, 03:54 UTC−4
Thank you for the explanation, the definition (stress dot velocity) looks as people's exception for mechanical energy flux.
But in the COMSOL of my computer (COMSOL 4.3), the (complex) mechanical energy flux has the unit of W/m^3, where I saw it using "point evaluation," and it was W/m for surface integration!?
Shouldn't they be W/m^2 & W in the sense of energy flux (and also the definition)?
Where did I miss it?
Thank you!

Liu
Thank you for the explanation, the definition (stress dot velocity) looks as people's exception for mechanical energy flux. But in the COMSOL of my computer (COMSOL 4.3), the (complex) mechanical energy flux has the unit of W/m^3, where I saw it using "point evaluation," and it was W/m for surface integration!? Shouldn't they be W/m^2 & W in the sense of energy flux (and also the definition)? Where did I miss it? Thank you! Liu

Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 24 mai 2013, 07:45 UTC−4


But in the COMSOL of my computer (COMSOL 4.3), the (complex) mechanical energy flux has the unit of W/m^3, where I saw it using "point evaluation," and it was W/m for surface integration!?
Shouldn't they be W/m^2 & W in the sense of energy flux (and also the definition)?



You are right. The unit for this quantity was wrong in versions prior to 4.3b. This does not affect the values (only the unit display) as long as results are displayed in the unit system of the model, but unit conversions will be wrong. To get the automatic unit conversion to work you can, as an example, evaluate solid.IX*1[m] instead of just solid.IX (assuming that the SI system is used).

Regards,
Henrik
[QUOTE] But in the COMSOL of my computer (COMSOL 4.3), the (complex) mechanical energy flux has the unit of W/m^3, where I saw it using "point evaluation," and it was W/m for surface integration!? Shouldn't they be W/m^2 & W in the sense of energy flux (and also the definition)? [/QUOTE] You are right. The unit for this quantity was wrong in versions prior to 4.3b. This does not affect the values (only the unit display) as long as results are displayed in the unit system of the model, but unit conversions will be wrong. To get the automatic unit conversion to work you can, as an example, evaluate solid.IX*1[m] instead of just solid.IX (assuming that the SI system is used). Regards, Henrik

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 24 mai 2013, 08:30 UTC−4
Thank you for these explanations!
Thank you for these explanations!

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 25 mai 2013, 00:32 UTC−4
Thank you very much!!
Thank you very much!!

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 14 août 2016, 10:01 UTC−4
Hi Henrik,
In my case I also need to calculate the energy flux through a face in piezoelectric module, but why the energy flux through a free boundary condition face is not zero?
Thank you!
Hi Henrik, In my case I also need to calculate the energy flux through a face in piezoelectric module, but why the energy flux through a free boundary condition face is not zero? Thank you!

Henrik Sönnerlind COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 8 years ago 23 août 2016, 10:25 UTC−4
The stresses are not exactly zero on a free boundary since they are a result computed from the gradients of the displacements.

As a consequence, the energy flux may also be non-zero. It should be small when compared to the flux at other places in the model, though.

How close to zero the stresses are on a free boundary depends on the mesh size.

Regards,
Henrik
The stresses are not exactly zero on a free boundary since they are a result computed from the gradients of the displacements. As a consequence, the energy flux may also be non-zero. It should be small when compared to the flux at other places in the model, though. How close to zero the stresses are on a free boundary depends on the mesh size. Regards, Henrik

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.