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 to get the current from current density

Please login with a confirmed email address before reporting spam

Hi Everyone,
I am simulating a piezo electric transducer in COMSOL Multiphysics 4.0a. In order to get the impedance of the transducer at different frequencies, I need the current flowing at the boundary of piezo ceramic. My model is 2D axisymmetric. Following are my problems:

1. For getting current(I) from the variable current density(J), I have to integrate the expression:
int(J*2*pi*r)dr over the boundary
However, I came to know from my friend that multiplication by (2 *pi *r) is applicable for COMSOL 3.5 but COMSOL 4.0a does this automatically through internal calculations for axisymmetric models. So there is no need to evaluate int(J*2*pi*r)dr. Rather int(J)dr is appropriate. Please confirm if it is correct.

2. How can I get the integration of current density. Using probe option, model coupling etc. Expedite if anyone have an idea

3. If I get the current by above means, then how to plot this impedance (Z=V/I), versus frequency. I mean how to access current value and corresponding frequencies to plot line graph.

Hope to hear anyone soon, Thanking you in anticipation
--
User COMSOL

5 Replies Last Post 8 août 2013, 06:20 UTC−4

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 27 janv. 2011, 04:28 UTC−5
Does anybody have any idea :(
--
User COMSOL
Does anybody have any idea :( -- User COMSOL

Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 28 janv. 2011, 11:09 UTC−5
Hi

I'm no longer sure about the 2*pi*r, I though it was only the 2*pi, but is it still there in 4.1 ? Ibelieve not that its back to 3.5a

The best as there are so many "small" things like this to remeber, is to check by estimating the volume from the 2D-axy symmetry:

integrate "1" then "r" and "2*pi*r" over a circle and check analytically.
Compare also the integration of "1" with the the option "Integration settings - Compute volume integral (axial symmetry)

Impedances are typically obtained either through energy calculation or via the lumped ports, or by the ratio imaginary versus real parts, all depends on your physics and what you are doing. Check the doc and the examples of the model library

--
Good luck
Ivar
Hi I'm no longer sure about the 2*pi*r, I though it was only the 2*pi, but is it still there in 4.1 ? Ibelieve not that its back to 3.5a The best as there are so many "small" things like this to remeber, is to check by estimating the volume from the 2D-axy symmetry: integrate "1" then "r" and "2*pi*r" over a circle and check analytically. Compare also the integration of "1" with the the option "Integration settings - Compute volume integral (axial symmetry) Impedances are typically obtained either through energy calculation or via the lumped ports, or by the ratio imaginary versus real parts, all depends on your physics and what you are doing. Check the doc and the examples of the model library -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 31 janv. 2011, 04:47 UTC−5
Thanks, this was really a good idea to check and I got it...
--
User COMSOL
Thanks, this was really a good idea to check and I got it... -- User COMSOL

Ivar KJELBERG COMSOL Multiphysics(r) fan, retired, former "Senior Expert" at CSEM SA (CH)

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 31 janv. 2011, 08:25 UTC−5
Hi

indeed, these checks are so simple and quickly done, that if you have any doubts, its much faster to check it. Just be careful with the default values i.e. 2D thickness = 1[m] mostly then the division by 1 induces you often in error. Or that you have the origin at the CoG bay design, most of the cases, then you have once a part far off and you forget to take the difference of the true CoG i.e. for an inertia evaluation

--
Good luck
Ivar
Hi indeed, these checks are so simple and quickly done, that if you have any doubts, its much faster to check it. Just be careful with the default values i.e. 2D thickness = 1[m] mostly then the division by 1 induces you often in error. Or that you have the origin at the CoG bay design, most of the cases, then you have once a part far off and you forget to take the difference of the true CoG i.e. for an inertia evaluation -- Good luck Ivar

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 8 août 2013, 06:20 UTC−4
What is correct expression for your module ? Could you please let me know ?
What is correct expression for your module ? Could you please let me know ?

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.