[Getdp] coupling electrokinetic and Magnetostatic problems

[Marco Antolovic]

 Hi François,
thanks for clarifying!
Marco 
    Il domenica 15 marzo 2020, 11:12:10 CET, François Henrotte <francois.henrotte at uclouvain.be> ha scritto:  
 
 Hello Marco,

Yes.The first three terms in the formulation belowdescribe equations associated with the test functionsof the edge unknowns (or node unknowns in 2D) of the magnetic vector potential field ‘a’in the whole domain of computation ‘Domain’.They express the weak form fo the Euler-Lagrange equation:
curl (nu curl a) + sigma( dadt + grad u) = 0 (curl h = j)

The next 2 terms are associated with the test functionsof the node unknowns of the scalar electric potential ‘v’,only for nodes in the conducting domain ‘DomainC’. They express the weak form fo the Euler-Lagrange equation:
div ( sigma( dadt + grad u) ) = 0 (div j=0)
By analyzing the ‘Quantity’ section of the ‘Formulation’ field,which I reproduce below
    Quantity {
      { Name a  ; Type Local ; NameOfSpace Hcurl_a_3D ; }

      { Name v  ; Type Local  ; NameOfSpace Hregion_u_3D ; }
      { Name U  ; Type Global ; NameOfSpace Hregion_u_3D [U] ; }
      { Name I  ; Type Global ; NameOfSpace Hregion_u_3D [I] ; }
    }

together with the informations contained in the « FunctionSpace » section and the declared « Constraints »  acting on it,GetDP is able to decide which Dofs are fixedand which are unknowns (and then associated with a test function).This is all automatic.That is the reason why it is enough in the Formulation descriptionto give the weak form of the equations you want to solve. 
Regards, 
Fr. 


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

      Galerkin { DtDof[ sigma[] * Dof{a} , {d v} ] ;         // <====== A1       In DomainC ; Jacobian Vol ; Integration I1 ; }
      Galerkin { [ sigma[] * Dof{d v} , {d v} ] ;               // <====== A2        In DomainC ; Jacobian Vol ; Integration I1 ; }
      GlobalTerm { [ Dof{I} , {U} ] ; In SurfaceElecWithI ; }
      Galerkin { [ - js0[] , {a} ] ; In DomainS ;        Jacobian Vol ; Integration I1 ; }    } I don't understand how A1 and A2 can be part of the same integral. I mean, based on the micro strip tutorial, every Galerkin{....} term belonging to Equation{...} is added up and this sum is then made equal to zero. And this equation should hold for a single test function.But here A1 and A2 have a test function different from the other terms. Does this mean that GetDP generates automatically two equations? One for test functions a' and one for test functions v’?








I'm lost, help please :D
Marco_______________________________________________
getdp mailing list
getdp at onelab.info
https://eur03.safelinks.protection.outlook.com/?url=http%3A%2F%2Fonelab.info%2Fmailman%2Flistinfo%2Fgetdp&data=02%7C01%7Cfrancois.henrotte%40uclouvain.be%7Cefa4dbe094fe44e3665d08d7a83200b0%7C7ab090d4fa2e4ecfbc7c4127b4d582ec%7C0%7C0%7C637162805466833982&sdata=xGSf9biUZB6zyCJRp3JikPklFdzb%2BYMznDFJydncgG0%3D&reserved=0


--François Henrotte Dr Ir - francois.henrotte at uclouvain.be - francois.henrotte at uliege.be UCLouvain - Bât. Euler a.217 - Av. G. Lemaître 4-6 , B-1348 Louvain-la-Neuve - +32(0)10 47 23 64 ULiège - Institut Montefiore I154 - Allée de la Découverte 10, B-4000 Liège - +32(0)4 366 37 36 





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