Interieur = 11; Exterieur = 12; Adiabatique = 13; Zone1 = 1; Zone2 = 2; Zone3 = 3; Zone4 = 4; Group { Interieur = Region[{Interieur}]; Exterieur = Region[{Exterieur}]; Adiabatique = Region[{Adiabatique}]; Ambiance = Region[{Exterieur,Interieur}]; AllGamma = Region[{Ambiance,Adiabatique}]; Zone1 = Region[{Zone1}]; Zone2 = Region[{Zone2}]; Zone3 = Region[{Zone3}]; Zone4 = Region[{Zone4}]; Omega = Region[{Zone1, Zone2, Zone3, Zone4}]; AllDomain = Region[{Omega, AllGamma}]; } Function { lambda[Zone1] = 1.15; lambda[Zone2] = 0.12; lambda[Zone3] = 0.029; lambda[Zone4] = 230; Rs[Interieur]= 0.06; Rs[Exterieur]= 0.11; Ts[Interieur]= 293.15; Ts[Exterieur]= 273.15; } FunctionSpace{ { Name Vh; Type Form0; BasisFunction{ {Name wn; NameOfCoef vn; Function BF_Node; Support AllDomain; Entity NodesOf[All];} } } } Jacobian { { Name JVol ; Case { { Region All ; Jacobian Vol ; } } } { Name JSur ; Case { { Region All ; Jacobian Sur ; } } } { Name JLin ; Case { { Region All ; Jacobian Lin ; } } } } Integration { { Name I1 ; Case { { Type Gauss ; Case { { GeoElement Point ; NumberOfPoints 1 ; } { GeoElement Line ; NumberOfPoints 4 ; } { GeoElement Triangle ; NumberOfPoints 6 ; } { GeoElement Quadrangle ; NumberOfPoints 7 ; } { GeoElement Tetrahedron ; NumberOfPoints 15 ; } { GeoElement Hexahedron ; NumberOfPoints 34 ; } } } } } } Formulation{ {Name Fourier; Type FemEquation; Quantity{ {Name u; Type Local; NameOfSpace Vh;} } Equation{ Galerkin{ [lambda[] * Dof{Grad u}, {Grad u}]; In Omega; Jacobian JVol; Integration I1;} Galerkin{ [ (1/Rs[])* Dof{u}, {u}]; In Ambiance; Jacobian JSur; Integration I1;} Galerkin{ [ (-1/(Rs[])*Ts[]), {u}]; In Ambiance; Jacobian JSur; Integration I1;} } } } Resolution{ {Name Fourier; System{ {Name Syst; NameOfFormulation Fourier;} } Operation{ Generate[Syst]; Solve[Syst]; SaveSolution[Syst]; } } } PostProcessing{ {Name Fourier; NameOfFormulation Fourier; Quantity{ { Name temperature; Value {Local{ [{u}-273.15]; In AllDomain; Jacobian JVol;}}} { Name flux_vector; Value {Local{ [-lambda[]*{Grad u}]; In AllDomain; Jacobian JVol;}}} { Name flux_norm; Value {Local{ [Norm[-lambda[]*{Grad u}]]; In AllDomain; Jacobian JVol;}}} } } } PostOperation{ {Name Fourier; NameOfPostProcessing Fourier; Operation{ Print[temperature, OnElementsOf Omega, File "temperature.pos"]; Print[flux_vector, OnElementsOf Omega, File "flux_vector.pos"]; Print[flux_norm, OnElementsOf Omega, File "flux_norm.pos"]; } } }