<div dir="ltr">Hi everyone,<div>I am trying to solve a 3D heat conduction problem . when I open my *.pro file I get an error :"Redefinition of JacobianMethod" .</div><div>Here is my heat.pro file which is included in another *.pro file :</div><div><br></div><div><br></div><div><div>DefineConstant[ Flag_AnalysisType = 0 ];</div><div><br></div><div>Jacobian {</div><div>  { Name JVol ;</div><div>    Case {</div><div>      { Region All ; Jacobian Vol ; }</div><div>    }</div><div>  }</div><div>  { Name JSur ;</div><div>    Case {</div><div>      { Region All ; Jacobian Sur ; }</div><div>    }</div><div>  }</div><div>}</div><div><br></div><div>Integration {</div><div>  { Name I1 ;</div><div>    Case {</div><div>      { Type Gauss ;</div><div>        Case {</div><div><span class="Apple-tab-span" style="white-space:pre">   </span>  { GeoElement Point       ; NumberOfPoints  1 ; }</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>  { GeoElement Line        ; NumberOfPoints  3 ; }</div><div><span class="Apple-tab-span" style="white-space:pre">     </span>  { GeoElement Triangle    ; NumberOfPoints  4 ; }</div><div><span class="Apple-tab-span" style="white-space:pre">       </span>  { GeoElement Quadrangle  ; NumberOfPoints  4 ; }</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>  { GeoElement Tetrahedron ; NumberOfPoints  4 ; }</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>  { GeoElement Hexahedron  ; NumberOfPoints  6 ; }</div><div><span class="Apple-tab-span" style="white-space:pre">        </span>  { GeoElement Prism       ; NumberOfPoints  6 ; }</div><div><span class="Apple-tab-span" style="white-space:pre">      </span>}</div><div>      }</div><div>    }</div><div>  }</div><div>}</div><div><br></div><div>FunctionSpace {</div><div>  { Name Hgrad_T; Type Form0;</div><div>    BasisFunction {</div><div>      { Name sn; NameOfCoef Tn; Function BF_Node; Support Tot_The;</div><div>        Entity NodesOf[All]; }</div><div>    }</div><div>    Constraint {</div><div>      { NameOfCoef Tn; EntityType NodesOf ; NameOfConstraint Temperature; }</div><div>    }</div><div>  }</div><div><br></div><div>}</div><div><br></div><div>Function{</div><div>  DefineFunction[Flux,qVol,h,hr,TConv];</div><div>}</div><div>Group{</div><div>  DefineGroup[SurConv_The,SurRad_The];</div><div>}</div><div><br></div><div>Formulation {</div><div><br></div><div>  { Name The_T ; Type FemEquation;</div><div>    Quantity {</div><div>      { Name T;  Type Local; NameOfSpace Hgrad_T; }</div><div>    }</div><div>    Equation {</div><div>      Galerkin { [ k[] * Dof{d T} , {d T} ];</div><div>                 In Vol_The; Integration I1; Jacobian JVol;  }</div><div><br></div><div>      Galerkin { DtDof [ rhoc[] * Dof{T} , {T} ];</div><div>                 In Vol_The; Integration I1; Jacobian JVol;  }</div><div><br></div><div>      Galerkin { [ -qVol[] , {T} ];</div><div>                 In Vol_The; Integration I1; Jacobian JVol;  }</div><div><br></div><div>      Galerkin { [ -Flux[] , {T} ]; // - sign for incoming flux</div><div>                 In Sur_The; Integration I1; Jacobian JSur;  }</div><div><br></div><div>      Galerkin { [ h[] * Dof{T} , {T} ] ;</div><div>                 In SurConv_The ; Integration I1; Jacobian JSur;  }</div><div><br></div><div>      Galerkin { [ -h[] * TConv[] , {T} ] ;</div><div>                 In SurConv_The ; Integration I1; Jacobian JSur;  }</div><div><br></div><div>      Galerkin { [ hr[{T}] * (({T}+273.)^4-(TConv[]+273.)^4) , {T} ] ;</div><div>                 In SurRad_The ; Integration I1; Jacobian JSur;  }</div><div>    }</div><div>  }</div><div><br></div><div>}</div><div><br></div><div>Resolution {</div><div>  { Name analysis;</div><div>    System {</div><div>      { Name T; NameOfFormulation The_T; }</div><div>    }</div><div>    Operation {</div><div>      If(Flag_AnalysisType == 0) // steady state</div><div>        Generate[T] ; Solve T ; SaveSolution T;</div><div>      EndIf</div><div>      If(Flag_AnalysisType == 1) // transient general</div><div>        InitSolution[T] ; SaveSolution[T] ;</div><div>        TimeLoopTheta [t0, t1, dt, 1.0] {</div><div>          Generate[T] ; Solve[T];</div><div>          Test[SaveFct[]] {</div><div>            SaveSolution[T];</div><div>          }</div><div>        }</div><div>      EndIf</div><div>      If(Flag_AnalysisType == 2) // transient linear fast</div><div>        InitSolution[T] ;  SaveSolution[T] ;</div><div>        GenerateSeparate[T] ;</div><div>        TimeLoopTheta [t0, t1, dt, 1.0] {</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>  Update[T, TimeFct[]] ;</div><div>          SolveAgain[T] ;</div><div><span class="Apple-tab-span" style="white-space:pre">       </span>  Test[SaveFct[]] {</div><div>            SaveSolution[T];</div><div>          }</div><div><span class="Apple-tab-span" style="white-space:pre"> </span>}</div><div>      EndIf</div><div>    }</div><div>  }</div><div>}</div><div><br></div><div><br></div><div>PostProcessing {</div><div>  { Name The; NameOfFormulation The_T;</div><div>    Quantity {</div><div>      { Name T; Value{ Local{ [ {T} ] ; In Vol_The; Jacobian JVol; } } }</div><div>    }</div><div>  }</div><div><br></div><div>}</div></div><div><br></div></div>