Circular Dependency Error: Modeling T-dependent BH Saturation in Frequency-Transient Study

Please login with a confirmed email address before reporting spam

Hello everyone,

I am working on an Induction Heating simulation using the Frequency-Transient study step. My goal is to model a ferromagnetic workpiece heated above its Curie temperature (Tc), capturing two nonlinear effects simultaneously:

1. Magnetic Saturation: Dependent on ∣H∣ (B-H curve).
2. Phase Transition: The drop in permeability to unity as T>Tc.

I attempted to combine these effects by creating an analytic function for relative permeability:

mu_r_combined(T, H) = (mu_r_T(T) - 1) * mu_r_H(H) + 1

mu_r_T(T): Interpolation of relative permeability vs. Temperature.
mu_r_H(H): Interpolation of relative permeability derived from B-H curve data.
The logic is to scale the saturation curve by the temperature factor, ensuring μr→1 when T>Tc.

The Problem: When I define the material property using this function (passing the magnetic field norm, e.g., mf.normH or similar, as the argument), the solver fails with a Circular Dependency Error.

Is there a standard workaround for this error in the Frequency Domain? Or is there a better way to implement "Temperature-dependent Effective BH curve" to avoid manually defining μr?

Any advice on the correct setup to handle this nonlinearity would be appreciated.

Thank you!


Reply

Please read the discussion forum rules before posting.

Please log in to post a reply.

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.