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.
time-dependent extension of a stationary problem
Posted 3 mar 2011, 20:06 GMT-5 Low-Frequency Electromagnetics, Studies & Solvers Version 5.1 9 Replies
Please login with a confirmed email address before reporting spam
As you know transient analysis is not available for 3D magnetic field analysis in comsol,
I have done my simulation for stationary case, however, I have to perform time-dependant calculations using sweeping parameter analysis, to achieve the steady-state response.
I would be grateful if you could inform me how to perform it, and introduce me an example inside the model gallery.
Best regards
Behrooz
Please login with a confirmed email address before reporting spam
indeed MF/MEF are not set up for time dependent solvers.
"My way" around, is to consider the physics and identify the time frames for the different items,
i.e. take a magnet falling down in a Cu tube:
The eddy currents induced in the tube will "hold" back and make the magnet fall very slowly.
How to model that in time dependent ?
If we consider the magnetic field you can solve it in stationary mode, extract the Eddy current forces related to a velocity displacement (note: in v4.1 the Lorentz velocity term is not yet set up in 2D, you must complete the BC with the equations, good exercise to write out the curl ;) and then add a global equation GE solving the velocity of the falling magnet expressed as derived from canceling out / minimising gravity forces with the magnetic Eddy current forces. Then you solve this independent ODE in the time domain and at each step you run a stationary for the MEF. You have hence a coupled problem and you solve it in the time domain, with MF/MEF
--
Good luck
Ivar
Please login with a confirmed email address before reporting spam
thanks for your helps,
i understand your explanation, it is similar to the magnetic brake example in the model gallery:
www.comsol.com/showroom/documentation/model/2014/
but it is done by version 3.5
behrooz
Please login with a confirmed email address before reporting spam
In version 4.2, should the Magnetic Fields module be capable of time dependent solutions in 3D? The specifications table (www.comsol.com/products/specifications/acdc/) seems to say that it is, but I haven't been able to get it to work.
I tried starting with a very simple problem and checked that the stationary solver could handle it, then included the variable t in a boundary condition. I got the error "Failed to find consistent initial values. Matrix has zero on diagonal. Last time step is not converged," but I don't see a reason for the error.
Thanks,
Phil
Please login with a confirmed email address before reporting spam
i think time-dependent solver should be selected,
similar to the following example in the model gallery:
www.comsol.com/showroom/documentation/model/2014/
Behrooz
Please login with a confirmed email address before reporting spam
The problem turned out to be that the conductivity of air was set to 0 (i.e., I didn't change it from the built-in value).
Ivar explained this in a different thread: www.comsol.com/community/forums/acdc-module/thread/21628/
The time-dependent solver does work in 3D.
-Phil
Please login with a confirmed email address before reporting spam
Please login with a confirmed email address before reporting spam
Hi
indeed MF/MEF are not set up for time dependent solvers.
"My way" around, is to consider the physics and identify the time frames for the different items,
i.e. take a magnet falling down in a Cu tube:
The eddy currents induced in the tube will "hold" back and make the magnet fall very slowly.
How to model that in time dependent ?
If we consider the magnetic field you can solve it in stationary mode, extract the Eddy current forces related to a velocity displacement (note: in v4.1 the Lorentz velocity term is not yet set up in 2D, you must complete the BC with the equations, good exercise to write out the curl ;) and then add a global equation GE solving the velocity of the falling magnet expressed as derived from canceling out / minimising gravity forces with the magnetic Eddy current forces. Then you solve this independent ODE in the time domain and at each step you run a stationary for the MEF. You have hence a coupled problem and you solve it in the time domain, with MF/MEF
--
Good luck
Ivar
Dear Ivar,
I am using a coupled 2D MF + Structural model (version 5.1) in order to simulate a pitch-catch setup of an electromagnetic acoustic (EMAT) transducer. I managed to get my model to partially work, i.e. I can see the Rayleigh wave generation and propagation in real time but I failed to correctly model the reception process. I am using a single-turn coil domain (meander-line coil) above the conductor and beneath a dc magnet domain and have added a velocity Lorentz term to my conductor domain where the components come from the structural mechanics module (solid.u_tX & solid.u_tY). Unfortunately, I cannot see the corresponding voltage across my RX coil which should have a profile similar to that of my Rayleigh wave ( see attached picture).
Having read your comment above, I am wondering if this is in fact due to the incapability of the MF module for the time dependent analysis, as you pointed out. If so, how can I work around this issue ?
Your advice is much appreciated in advance,
Regards,
Hamed
Attachments:
Please login with a confirmed email address before reporting spam
Regards,
Hamed
Please login with a confirmed email address before reporting spam
Hi
indeed MF/MEF are not set up for time dependent solvers.
"My way" around, is to consider the physics and identify the time frames for the different items,
i.e. take a magnet falling down in a Cu tube:
The eddy currents induced in the tube will "hold" back and make the magnet fall very slowly.
How to model that in time dependent ?
If we consider the magnetic field you can solve it in stationary mode, extract the Eddy current forces related to a velocity displacement (note: in v4.1 the Lorentz velocity term is not yet set up in 2D, you must complete the BC with the equations, good exercise to write out the curl ;) and then add a global equation GE solving the velocity of the falling magnet expressed as derived from canceling out / minimising gravity forces with the magnetic Eddy current forces. Then you solve this independent ODE in the time domain and at each step you run a stationary for the MEF. You have hence a coupled problem and you solve it in the time domain, with MF/MEF
--
Good luck
Ivar
Hello Ivar,
I have done all the geometries, replicating the 2d axisymmetric version of the falling_magnet.mph in 3d but I do not have a velocity, lorentz term that I can use for my copper tube domain.
But I do not know how to "complete the BC with the equations" .I know that the equation I have to add is : Ji=sigma*E + sigma*v x B where Ji is the induced current density.
How does one give such an equation and where? How are the vectors "v" and "B" defined? I have defined sigma as a parameter and hence, I have no problem with it.
Please tell me as to how I have to proceed, in a detailed manner.
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.
Suggested Content
- FORUM Generalized Plain Strain and out-of-plane mode extension
- KNOWLEDGE BASE Improving Convergence of Nonlinear Stationary Models
- KNOWLEDGE BASE Solving time dependent models with inconsistent initial values
- KNOWLEDGE BASE Solving Complicated Models
- BLOG How to Implement the Weak Form for Time-Dependent Equations