[Getdp] "Jacobian method error"

Ali Masoudi masoudi.mech at gmail.com
Mon Sep 14 19:02:25 CEST 2015


Hi everyone,
I am trying to solve a 3D heat conduction problem . when I open my *.pro
file I get an error :"Redefinition of JacobianMethod" .
Here is my heat.pro file which is included in another 3D.pro file :


DefineConstant[ Flag_AnalysisType = 0 ];

Jacobian {
  { Name JVol ;
    Case {
      { Region All ; Jacobian Vol ; }
    }
  }
  { Name JSur ;
    Case {
      { Region All ; Jacobian Sur ; }
    }
  }
}

Integration {
  { Name I1 ;
    Case {
      { Type Gauss ;
        Case {
  { GeoElement Point       ; NumberOfPoints  1 ; }
  { GeoElement Line        ; NumberOfPoints  3 ; }
  { GeoElement Triangle    ; NumberOfPoints  4 ; }
  { GeoElement Quadrangle  ; NumberOfPoints  4 ; }
  { GeoElement Tetrahedron ; NumberOfPoints  4 ; }
  { GeoElement Hexahedron  ; NumberOfPoints  6 ; }
  { GeoElement Prism       ; NumberOfPoints  6 ; }
}
      }
    }
  }
}

FunctionSpace {
  { Name Hgrad_T; Type Form0;
    BasisFunction {
      { Name sn; NameOfCoef Tn; Function BF_Node; Support Tot_The;
        Entity NodesOf[All]; }
    }
    Constraint {
      { NameOfCoef Tn; EntityType NodesOf ; NameOfConstraint Temperature; }
    }
  }

}

Function{
  DefineFunction[Flux,qVol,h,hr,TConv];
}
Group{
  DefineGroup[SurConv_The,SurRad_The];
}

Formulation {

  { Name The_T ; Type FemEquation;
    Quantity {
      { Name T;  Type Local; NameOfSpace Hgrad_T; }
    }
    Equation {
      Galerkin { [ k[] * Dof{d T} , {d T} ];
                 In Vol_The; Integration I1; Jacobian JVol;  }

      Galerkin { DtDof [ rhoc[] * Dof{T} , {T} ];
                 In Vol_The; Integration I1; Jacobian JVol;  }

      Galerkin { [ -qVol[] , {T} ];
                 In Vol_The; Integration I1; Jacobian JVol;  }

      Galerkin { [ -Flux[] , {T} ]; // - sign for incoming flux
                 In Sur_The; Integration I1; Jacobian JSur;  }

      Galerkin { [ h[] * Dof{T} , {T} ] ;
                 In SurConv_The ; Integration I1; Jacobian JSur;  }

      Galerkin { [ -h[] * TConv[] , {T} ] ;
                 In SurConv_The ; Integration I1; Jacobian JSur;  }

      Galerkin { [ hr[{T}] * (({T}+273.)^4-(TConv[]+273.)^4) , {T} ] ;
                 In SurRad_The ; Integration I1; Jacobian JSur;  }
    }
  }

}

Resolution {
  { Name analysis;
    System {
      { Name T; NameOfFormulation The_T; }
    }
    Operation {
      If(Flag_AnalysisType == 0) // steady state
        Generate[T] ; Solve T ; SaveSolution T;
      EndIf
      If(Flag_AnalysisType == 1) // transient general
        InitSolution[T] ; SaveSolution[T] ;
        TimeLoopTheta [t0, t1, dt, 1.0] {
          Generate[T] ; Solve[T];
          Test[SaveFct[]] {
            SaveSolution[T];
          }
        }
      EndIf
      If(Flag_AnalysisType == 2) // transient linear fast
        InitSolution[T] ;  SaveSolution[T] ;
        GenerateSeparate[T] ;
        TimeLoopTheta [t0, t1, dt, 1.0] {
  Update[T, TimeFct[]] ;
          SolveAgain[T] ;
  Test[SaveFct[]] {
            SaveSolution[T];
          }
}
      EndIf
    }
  }
}


PostProcessing {
  { Name The; NameOfFormulation The_T;
    Quantity {
      { Name T; Value{ Local{ [ {T} ] ; In Vol_The; Jacobian JVol; } } }
    }
  }

}


I have also attached all of my files.

-- 
Ali Masoudi, PhD Candidate
University of Houston
Department of Mechanical Engineering
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150914/b7856015/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 3D.geo
Type: application/octet-stream
Size: 4702 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150914/b7856015/attachment-0001.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 3D.pro
Type: application/octet-stream
Size: 11913 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150914/b7856015/attachment-0003.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: heat.pro
Type: application/octet-stream
Size: 2692 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150914/b7856015/attachment-0004.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Thermal.pro
Type: application/octet-stream
Size: 3098 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150914/b7856015/attachment-0005.pro>