[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>