/* InductionTestBox.pro */ Group { Air = Region[101]; Coil = Region[102]; // Susceptor = Region[3]; Ground1 = Region[8]; Ground2 = Region[10]; OuterBoundary = Region[{1,2,3,4,5,6}]; /* conducting regions = DomainC_Mag */ DomainC_Mag = Region[Coil]; /* non-conducting regions = DomainCC_Mag */ DomainCC_Mag = Region[Air]; /* entire domain */ Domain_Mag = Region[{DomainC_Mag, DomainCC_Mag}]; /* domain with source current */ DomainS_Mag = Region[Coil]; } Function { epsr[Air] = 1.; epsr[Coil] = 1.; //epsr[Susceptor] = 1.; } Function { sigma[Coil]=1.e6; sigma[Air]=0.; //sigma[Susceptor]=0.; } Function { mu0 = 4.e-7 * Pi; nu [ Domain_Mag] = 1. / mu0; } Function { js[Coil] = Vector[0 ,1000,0]; } Constraint { { Name ElectricScalarPotential; Type Assign; Case { { Region Ground1; Value 0.; } { Region Ground2; Value 1.; } } } { Name MagneticVectorPotential_3D; Type Assign; Case { { Region OuterBoundary ; Value 0.; } } } } Include "Jacobian_Lib.pro" Include "Integration_Lib.pro" Include "MagDyn_av_3D.pro" PostOperation { { Name MapE; NameOfPostProcessing EleSta_v; Operation { Print [ v, OnElementsOf DomainC_Mag, File "Coil_v.pos" ]; Print [ e, OnElementsOf DomainC_Mag, File "Coil_e.pos" ]; Print [ j, OnElementsOf DomainC_Mag, File "Coil_j.pos" ]; } } } PostOperation { { Name MapA; NameOfPostProcessing MagSta; Operation { Print [ a, OnElementsOf Domain_Mag, File "MagSta_a.pos" ]; Print [ h, OnElementsOf Domain_Mag, File "MagSta_h.pos" ]; } } } PostOperation { { Name MapDyn; NameOfPostProcessing MagDyn_avf; Operation { Print [ a, OnElementsOf Domain_Mag, File "MagDyn_a.pos" ]; Print [ h, OnElementsOf Domain_Mag, File "MagDyn_h.pos" ]; Print [ j, OnElementsOf DomainC_Mag, File "MagDyn_j.pos" ]; } } }