TITLE 'Fraktalola' COORDINATES cartesian3 VARIABLES Temp SELECT { method controls } ERRLIM = 0.001 ngrid = 30 THREADS = 2 DEFINITIONS { parameter definitions } ! geometry z0 = 0 z1 = 0.01 z2 = 0.03 z3 = 0.02 z4 = 0.03 z5 = 0.04 x1 = 0.052/2 area1 = 4*x1*x1 x2 = 0.05/2 area2 = 4*x2*x2 d1 = 0.01/2 ! physical parameters k_alu = 150 k_pu = 0.0001 k_cu = 400 k_specimen = 0.9 k_sonde = 15 ! Temperatursensor Q_sonde = 2.5/(d1*d1*3.1415*z2) temp_cold = 10 temp_hot = 40 ! standard values k = k_cu Q = 0 ! test parameter f_z4 = z4 ! INITIAL VALUES EQUATIONS { PDE's, one for each variable } k*div(-grad(Temp)) = Q EXTRUSION ! small z at first, you get merging errors if you define it the other way round (yalla...) SURFACE "S0" Z = z0 LAYER "L0" SURFACE "S1" Z = z0+z1 LAYER "L1" SURFACE "S2" Z = z0+z1+z2 LAYER "L2" SURFACE "S3" Z = z0+z1+z2+z3 LAYER "L3" SURFACE "S4" Z = z0+z1+z2+z3+z4 LAYER "L4" SURFACE "S5" Z = z0+z1+z2+z3+z4+z5 ! CONSTRAINTS { Integral constraints } BOUNDARIES SURFACE "S5" VALUE(Temp) = temp_cold SURFACE "S0" VALUE(Temp) = temp_hot REGION 1 "R1" k = k_cu Q = 0 START(-x1,-x1) periodic(x,-y) LINE TO (x1,-x1) line TO (x1,x1) nobc(Temp) line TO (-x1,x1) periodic(-x,y) line TO CLOSE Limited REGION 2 "R2" Layer "L3" k = k_specimen Q = 0 START(-x1,-x1) LINE TO (x1,-x1) TO (x1,x1) TO (-x1,x1) TO CLOSE Limited REGION 3 "R3" k = k_pu Q = 0 Layer "L0" Layer "L1" Layer "L2" START(-x1,-x1) LINE TO (x1,-x1) TO (x1,x1) TO (-x1,x1) TO CLOSE Limited REGION 4 "R4" k = k_alu Q = 0 Layer "L1" Layer "L2" START(-x2,-x2) LINE TO (x2,-x2) TO (x2,x2) TO (-x2,x2) TO CLOSE Limited REGION 5 "R5" k = k_sonde Q = Q_sonde Layer "L1" Start(0,0+d1) Natural(Temp)=0 Arc(Center=0,0) angle=360 Feature ! TIME 0 TO 1 { if time dependent } PLOTS { save result displays } contour(Temp) as 'temp' on x = 0 painted contour(Temp) as 'temp' on y = 0 painted ! contour(Temp) as 'temp' on z = z0+z1+z2+z3 painted contour(k) as 'wlf' on x = 0 painted levels = 0, 0.9, 9.9, 99.9, 299, 499 Summary EXPORT FILE="Simulation Info.txt" !http://www.pdesolutions.com/help/integralsinthreedimensions.html?zoom_highlightsub=integral Report(val(temp,0,0,z0+z1+z2+z3)) as "Temp on AluSurface" Report(sintegral(temp/area2,"S3","R4","L3")) as "temp_S3R4L3" Report(val(temp,0,0,z0+z1+z2+z3+z4)) as "Temp on cold Surface" Report(sintegral(temp/area1,"S4", "R2","L3")) as "temp_S4R2L3" Report("") Report(sintegral(-normal(k*grad(Temp))/abs(sintegral(temp/area2,"S3","R4","L3")-sintegral(temp/area1,"S4", "R2","L3"))*z4/area2,"S3","R4","L3")) as "int_S3R4L3" Report(sintegral(-normal(k*grad(temp))/abs(sintegral(temp/area2,"S3","R4","L3")-sintegral(temp/area1,"S4", "R2","L3"))*z4/area1,"S4", "R2","L4")) as "int_S4R2L4" END