<div dir="ltr"><div><div><div><div><div>Dear List Members,<br><br></div>I am at the very beginning of learning how to use getdp. In fact I bumped into a problem at my first model, which is about an axisymmetric 2D coil in time domain. I use vector potential formalism with excitation current perpendicular to the plane. Therefore this is a nodal FEM problem because the vector potential has only phi component as the current density. <br><br></div>I did my best to tackle this problem but something is wrong with my *.pro file. However the example files run well. <br><br></div>I use getdp version 2.6.1 64 bit on Linux and get this message:<br><br><span style="font-family:monospace,monospace">getdp <a href="http://coil.pro">coil.pro</a> -solve '#1' -pos '#1' -msh coil.msh<br>Info    : Running 'getdp <a href="http://coil.pro">coil.pro</a> -solve #1 -pos #1 -msh coil.msh' [GetDP 2.6.1, 1 node]<br>Info    : Started (Fri Aug  7 00:13:51 2015, CPU = 0.076s, Mem = 255.719Mb)<br>Info    : Initializing Gmsh<br>Info    : Loading problem definition '<a href="http://coil.pro">coil.pro</a>'<br>Info    : Loading problem definition 'param.geo'<br>Info    : Selected Resolution 'Solution'<br>Info    : Loading Geometric data 'coil.msh'<br>Info    : System 'Syst' : Real<br>P r e - P r o c e s s i n g . . .<br>Info    : Treatment Formulation 'VectorPotentialFormulation'<br>make: *** [solve] Segmentation fault (core dumped)</span><br><br><br></div>Could you take a look at my project file what's wrong with it?  This is really my very first project in getDP. So any clue is welcome.<br><br></div><div>Thank you for your help in advance.<br><br></div><div>Peter<br><br>-----------<br></div><div><a href="http://coil.pro">coil.pro</a>:<br>------------<br><span style="font-family:monospace,monospace">// <a href="http://coil.pro">coil.pro</a>: weak formulation<br><br>Include "param.geo";<br><br>Group{ // OK<br>    Omega_air = Region[{Omega_air}];<br>    Omega_coil = Region[{Omega_coil}];<br>    Omega = Region[ All ];<br>    Gama_D = Region[{Gama_D}]; // Dirichlet for A=0<br>}<br><br>Function{ // OK<br>    InitialState[Omega] = 0.;<br>    mu0 = 4e-7*Pi;<br>    Freq = 50;<br>    t0 = 0;<br>    tmax = 2.;<br>    dt = 0.1;<br>    theta = 1.; // backward Euler<br>    source_current_t = {0, 1, 2};<br>    source_current_j = {0, 1, 1};<br>    Js[] = InterpolationLinear[$1]{ListAlt[source_current_t, source_current_j]};<br>//    Jext[] = Js[$Time];<br>    Jext[] = 50*(1-Exp[-0.5*$Time]);<br>    sigma[Omega_air] = 1;<br>    sigma[Omega_coil] = 6e7;<br>}<br><br>Jacobian { // OK<br> { Name JVol; // for 2D regions<br>   Case {<br>     {Region All; Jacobian VolAxi;}<br>   }<br> }<br> { Name JSur; // for 1D regions<br>   Case {<br>     {Region All; Jacobian SurAxi;}<br>   }<br> }<br>}<br><br>Integration {<br>  { Name I1 ;<br>    Case {<br>      { Type Gauss ;<br>        Case {<br>          { GeoElement Point ; NumberOfPoints  1 ; }<br>          { GeoElement Line ; NumberOfPoints  4 ; }<br>          { GeoElement Triangle ; NumberOfPoints  6 ; }<br>          { GeoElement Quadrangle ; NumberOfPoints 7 ; }<br>          { GeoElement Tetrahedron ; NumberOfPoints 15 ; }<br>          { GeoElement Hexahedron ; NumberOfPoints 34 ; }<br>        }<br>      }<br>    }<br>  }<br>}<br><br>Constraint{ // OK<br>    // Initial State<br>    {Name InitialData; Type Init;<br>        Case {<br>            {Region Omega; Value InitialState[];}<br>        }<br>    }<br>    // Dirichlet in Inf.<br>    {Name DirichletBC; Type Assign;<br>        Case {<br>            {Region Gama_D; Value 0;}<br>        }<br>    }<br>    // No constraint for the axis<br>}<br><br>FunctionSpace{ // OK<br>  { Name Afs; Type Form0;<br>    BasisFunction{<br>      {Name Ni; NameOfCoef Ai; Function BF_Node;<br>       Support Omega; Entity NodesOf[All];}<br>    }<br>    Constraint {<br>        {NameOfCoef Ai; EntityType NodesOf;<br>         NameOfConstraint DirichletBC;}<br>        {NameOfCoef Ai; EntityType NodesOf;<br>         NameOfConstraint InitialData;}<br>    }<br>  }<br>}<br><br>Formulation{ // OK<br>  {Name VectorPotentialFormulation; Type FemEquation;<br>   Quantity{<br>     {Name a; Type Local; NameOfSpace Afs;}<br>   }<br>   Equation{<br>    Galerkin{ [1/mu0*Dof{Curl a},{Curl a}];<br>      In Omega; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ DtDof[sigma[]*Dof{a},{a}];<br>      In Omega; Jacobian JVol; Integration I1;}<br><br>    Galerkin{ [-Jext[],{a}];<br>      In Omega_coil; Jacobian JVol; Integration I1;}<br>   }<br>  }<br>}<br><br>Resolution{<br>  {Name Solution;<br>    System{<br>      {Name Syst; NameOfFormulation VectorPotentialFormulation;}<br>    }<br>    Operation{<br>        InitSolution[Syst];SaveSolution[Syst];<br>        TimeLoopTheta[t0,tmax,dt,theta] {<br>            Generate[Syst]; Solve[Syst]; SaveSolution[Syst];<br>        }<br>    }<br>  }<br>}<br><br>PostProcessing{<br>  {Name VectorPotentialPostProc; NameOfFormulation VectorPotentialFormulation;<br>    Quantity {<br>      {Name vec_pot; Value {Local{[{a}];In Omega;Jacobian JVol;}}}<br>    }<br>  }<br>}<br><br>PostOperation{<br>  {Name Map_Vectorpotential; NameOfPostProcessing VectorPotentialPostProc;<br>    Operation{<br>      Print[vec_pot, OnElementsOf Omega, File "A.pos"];<br>    }<br>  }<br>}<br></span><br>---------------------<br></div><div>param.geo<br></div><div>---------------------<br><span style="font-family:monospace,monospace">// param.geo: the auxiliary file<br><br>// PhysicalIDs<br>Omega_coil = 700;<br>Omega_air = 1000;<br>Gama_D = 500;<br><br>// Geometry<br>r_enc = 60e-3;<br>z_enc = 40e-3;<br>inner_radius = 59e-3/2;<br>thick = 100e-6;<br>width = 4e-3;<br><br>// Mesh<br>lc_air = 0.008;<br>lc_coil = 0.5e-3;<br></span><br><br>---------------------<br></div><div>coil.geo<br></div><div>---------------------<br><br><span style="font-family:monospace,monospace">// coil.geo: creation of the geometry with GMSH<br>// <br><br>Include "param.geo";<br><br>// Enclosure points<br>p1 = newp; Point(p1) = {0,-z_enc,0,lc_air};<br>p2 = newp; Point(p2) = {r_enc,-z_enc,0,lc_air};<br>p3 = newp; Point(p3) = {r_enc,z_enc,0,lc_air};<br>p4 = newp; Point(p4) = {0,z_enc,0,lc_air};<br><br>// First turn points<br>p5 = newp; Point(p5) = {inner_radius,-width/2,0,lc_coil};<br>p6 = newp; Point(p6) = {inner_radius,width/2,0,lc_coil};<br>p7 = newp; Point(p7) = {inner_radius+thick,width/2,0,lc_coil};<br>p8 = newp; Point(p8) = {inner_radius+thick,-width/2,0,lc_coil};<br><br>// Enclosure lines<br>L1 = newreg; Line(L1) = {p1,p2};<br>L2 = newreg; Line(L2) = {p2,p3};<br>L3 = newreg; Line(L3) = {p3,p4};<br>L4 = newreg; Line(L4) = {p4,p1};<br><br>// First turn lines<br>L5 = newreg; Line(L5) = {p5,</span>p6};<br><span style="font-family:monospace,monospace">L6 = newreg; Line(L6) = {p6,p7};<br>L7 = newreg; Line(L7) = {p7,p8};<br>L8 = newreg; Line(L8) = {p8,p5};<br><br>// Line loops<br>Bound = newreg; Line Loop(Bound) = {L1,L2,L3,L4};<br>Turn = newreg; Line Loop(Turn) = {L5,L6,L7,L8};<br><br>SurfaceOmega_air = newreg; Plane Surface(SurfaceOmega_air) = {Bound,Turn};<br>SurfaceOmega_coil = newreg; Plane Surface(SurfaceOmega_coil) = {Turn};<br><br>Physical Surface(Omega_air) = {SurfaceOmega_air};<br>Physical Surface(Omega_coil) = {SurfaceOmega_coil};<br><br>// For Dirichlet bnd cnd<br>Physical Line(Gama_D) = {L1,L2,L3};<br><br>Mesh.Algorithm = 1;<br>//Mesh.SurfaceFaces = 1;<br>//Mesh.SurfaceNumbers = 1;<br>//Mesh.LabelType = 2; // Show Phys ID<br>//Mesh 2;</span><br></div></div>