[Getdp] Coupling Analysis (Thermo MagDyn)

[ABE Hiroshi]

Dear All,

I am working on a coupling analysis of magnetodynamics and thermal dynamics. Referring to  “indheat” sample, I build a model.
It uses A-V formulation regarding Magnetodynamics, and I would like to couple the electric current in the thermal formulation.

They are:

Formulation {

  { Name MagStaDyn_av_js0_3D ; Type FemEquation ;
    Quantity {
      { Name a  ; Type Local ; NameOfSpace HSpace ; }
      { Name v  ; Type Local ; NameOfSpace USpace ; }
    }

    Equation {
      Galerkin { [ nu[] * Dof{d a} , {d a} ] ;
        In Domain ; Jacobian Vol ; Integration II ; }
      Galerkin { DtDof[ sigma[] * Dof{a} , {a} ] ;
        In DomainC ; Jacobian Vol ; Integration II ; }
      Galerkin { [ sigma[] * Dof{d v} , {a} ] ;
        In DomainC ; Jacobian Vol ; Integration II ; }

      Galerkin { [ -js0[], {a} ] ;
        In  DomainS ; Jacobian Vol ; Integration II ; }


      Galerkin{ DtDof[ sigma[] * Dof{a}, {d v} ] ;
	In DomainC ; Jacobian Vol ; Integration II ; }
      Galerkin{ [ sigma[] * Dof{d v} , {d v} ] ;
	In DomainC ; Jacobian Vol ; Integration II ; }

    }
  }
  
  { Name Thermal; Type FemEquation;
    Quantity {
      { Name t; Type Local; NameOfSpace TSpace; }
      { Name a; Type Local; NameOfSpace HSpace; }
      { Name v; Type Local; NameOfSpace USpace; }
    }
    Equation {
      Galerkin { [ K[] * Dof{d t}, {d t} ];
	In DomainC; Integration II; Jacobian Vol; }
      Galerkin { DtDof [ rho[]*Cp[] * Dof{t}, {t} ];
	In DomainC; Integration II; Jacobian Vol; }
      Galerkin { [ -1./sigma[]*( 
	    Re[-(Dt[{a}]+{d v})]*
	    Re[-(Dt[{a}]+{d v})]+
	    Im[-(Dt[{a}]+{d v})]*
	    Im[-(Dt[{a}]+{d v})]), {t} ];
	In DomainC; Integration II; Jacobian Vol; }
      
      Galerkin { [ Ht[]*Dof{t}, {t} ];
	In Skin_ECore; Jacobian Sur ; Integration II ; }
      Galerkin { [-Ht[]*Tamb[], {t} ];
	In Skin_ECore; Jacobian Sur ; Integration II ; }
      Galerkin { [ sigma_sb*Ep[]*(Dof{t})^4, {t} ];
	In Skin_ECore; Jacobian Sur ; Integration II ; }
      Galerkin { [ -sigma_sb*Ep[]*(Tamb[])^4, {t} ];
	In Skin_ECore; Jacobian Sur ; Integration II ; }
      
    }
  }
}

The MagStaDyn_av_js0_3D formulation gives a resonable solution but Thermal gives weird results.

I know the “indheat” example uses T-O formulation for coupling analysis. Are there any reasons to take T-O formulation instead of A-V formulation?
Any ways for A-V formulation?

Thank you so much.

Best,

ABE Hiroshi
 from Tokorozawa, JAPAN