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.

Volume Integration in COMSOL using different limit

Please login with a confirmed email address before reporting spam

Dear COMSOL Users, I am integrating a geometrical model in COMSOL using volume integration in derived value tab. It integrate whole geometry at once through total height (h) but I want to take ratio of integration of total height to integration of 1/3rd of the height. Now I don't know How I can select 1/3 limit in COMSOL. Can anyone give me an idea about integration this single geomtry in limit from 0 to h/3 and 0 to h ? It should not consist of two domains as doing so will affect my simulation result.

Thank you


3 Replies Last Post 9 gen 2018, 11:02 GMT-5
Magnus Ringh COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 9 gen 2018, 03:36 GMT-5
Updated: 7 years ago 10 gen 2018, 10:08 GMT-5

Hi,

You can define an Integration coupling operator under Component>Definitions>Component Couplings (say, intop1), valid in all domains, and then evaluate it using a Global Evaluation node.

If you want to integrate the entire volume, then use

intop1(1)

If you want to integrate to another limit (in the z direction), then use

intop1(1*(z < zlim))

where zlim is the limit in the z direction, and the geometry starts at z = 0.

With 1 as the integrand, the intop1 operator computes the volume. You can use any other quantity to integrate instead, if desired.

You can also insert an "artificial boundary" at z = zlim so that you get separate domains for the top and bottom parts of the geometry. It should not affect your simulation results (other than enforcing the mesh to respect the boundary between the domains).

Best regards, Magnus Ringh, COMSOL

Hi, You can define an Integration coupling operator under Component>Definitions>Component Couplings (say, intop1), valid in all domains, and then evaluate it using a Global Evaluation node. If you want to integrate the entire volume, then use intop1(1) If you want to integrate to another limit (in the z direction), then use intop1(1*(z < zlim)) where zlim is the limit in the z direction, and the geometry starts at z = 0. With 1 as the integrand, the intop1 operator computes the volume. You can use any other quantity to integrate instead, if desired. You can also insert an "artificial boundary" at z = zlim so that you get separate domains for the top and bottom parts of the geometry. It should not affect your simulation results (other than enforcing the mesh to respect the boundary between the domains). Best regards, Magnus Ringh, COMSOL

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 9 gen 2018, 07:09 GMT-5

Dear Magnus Ringh, Thank you for your response. I am integrating a composite column. I did create the intop1, but in global evolution it is giving error. If you could help me here. I am using Solid Mechanics, Eigenfrequency module and want to find elestic strian energy density (solid.Ws) through total height of composite column and different height value like(H/2,H/3,H/4 etc) so that I can take the ratio of these two integrated strain energy density values which I calculated. If you could comment on this. I did have tried dividing my column in two domains but first it is very lengthy process and secondly result obtained is not accurate, it contains error due to interface b/w two domains. COMSOL may take that interface as weak zone. Which affect my results.

Thank you

Dear Magnus Ringh, Thank you for your response. I am integrating a composite column. I did create the intop1, but in global evolution it is giving error. If you could help me here. I am using Solid Mechanics, Eigenfrequency module and want to find elestic strian energy density (solid.Ws) through total height of composite column and different height value like(H/2,H/3,H/4 etc) so that I can take the ratio of these two integrated strain energy density values which I calculated. If you could comment on this. I did have tried dividing my column in two domains but first it is very lengthy process and secondly result obtained is not accurate, it contains error due to interface b/w two domains. COMSOL may take that interface as weak zone. Which affect my results. Thank you

Magnus Ringh COMSOL Employee

Please login with a confirmed email address before reporting spam

Posted: 7 years ago 9 gen 2018, 11:02 GMT-5
Updated: 7 years ago 10 gen 2018, 10:09 GMT-5

Hi,

If you have defined the integration coupling operator so that it is active in the geometry, an expression such as intop1(solid.Ws) should work to integrate the elastic strain energy using a Global Evaluation node. For a ratio of integrated values, for example, you can evaluate an expression such as

intop1(solid.Ws*(z < H/2))/intop1(solid.Ws*(z < H/3))

if that is what you are interested in. To investigate why an added interior boundary would affect your simulation, we would need access to the model file. If this does not help you, I suggest that you contact COMSOL Multiphysics Support: COMSOL Multiphysics Support Center.

Best regards, Magnus Ringh, COMSOL

Hi, If you have defined the integration coupling operator so that it is active in the geometry, an expression such as intop1(solid.Ws) should work to integrate the elastic strain energy using a Global Evaluation node. For a ratio of integrated values, for example, you can evaluate an expression such as intop1(solid.Ws\*(z < H/2))/intop1(solid.Ws\*(z < H/3)) if that is what you are interested in. To investigate why an added interior boundary would affect your simulation, we would need access to the model file. If this does not help you, I suggest that you contact COMSOL Multiphysics Support: [COMSOL Multiphysics Support Center](http://www.comsol.com/support/). Best regards, Magnus Ringh, COMSOL

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.