[Getdp] Question magestadyn_av_js0_3d.pro

[Ruth Vazquez Sabariego]

Dear Frederic,
When the core is non-conducting, there is no difference at all between what you wrote and what follows:
// Magnetic vector potential a (b = curl a)
 { Name Hcurl_a_3D ; Type Form1 ;
   BasisFunction {// a = a_e * s_e
     { Name se ; NameOfCoef ae ; Function BF_Edge ;
       Support Domain ; Entity EdgesOf[ All ] ; }
   }

Only when the core is conducting, it makes sense to split the basis function in two parts and that for allowing different constraints.
Indeed, the tree-cotree gauging will be kept outside the conducting region, so the gauge constraint will be applied only on the first part, excluding the boundary of the conductors.
See hereafter.

 { Name Hcurl_a_3D ; Type Form1 ;
    BasisFunction {// a = a_e * s_e
      { Name se ; NameOfCoef ae ; Function BF_Edge ;
        Support Domain ; Entity EdgesOf[ All, Not SkinDomainC ] ; }
      { Name se2 ; NameOfCoef ae2 ; Function BF_Edge ;
        Support Domain ; Entity EdgesOf[ SkinDomainC ] ; }
    }
    Constraint {
      { NameOfCoef ae  ; EntityType EdgesOf ; NameOfConstraint MVP_3D ; }
      { NameOfCoef ae2 ; EntityType EdgesOf ; NameOfConstraint MVP_3D ; }

      { NameOfCoef ae  ; EntityType EdgesOfTreeIn ; EntitySubType StartingOn ;
          NameOfConstraint GaugeCondition_a ; }
    }
  }

Best regards,
Ruth Sabariego


—
Prof. Ruth V. Sabariego
KU Leuven, Dept. Electrical Engineering ESAT/Electa, EnergyVille
http://www.esat.kuleuven.be/electa
http://www.energyville.be








On 13 Sep 2019, at 02:41, Frederic Trillaud <ftrillaudp at gmail.com<mailto:ftrillaudp at gmail.com>> wrote:

Dear all,

I have quick question, what does the following piece of code mean?, checking the code, when the core is not conductive it seems that SkinDomainC is empty, when it is conducting, it is the surface of the Core domain. I am not clear why the skin of the core should be defined apart.

 // Magnetic vector potential a (b = curl a)
 { Name Hcurl_a_3D ; Type Form1 ;
   BasisFunction {// a = a_e * s_e
     { Name se ; NameOfCoef ae ; Function BF_Edge ;
       Support Domain ; Entity EdgesOf[ All, Not SkinDomainC ] ; }
     { Name se2 ; NameOfCoef ae2 ; Function BF_Edge ;
       Support Domain ; Entity EdgesOf[ SkinDomainC ] ; }
   }

Best,

Frederic Trillaud



_______________________________________________
getdp mailing list
getdp at onelab.info<mailto:getdp at onelab.info>
http://onelab.info/mailman/listinfo/getdp

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://onelab.info/pipermail/getdp/attachments/20190913/e14b3b57/attachment.html>