================================================================ === FlexPDE Version 7.04/W64 16:06:44 Apr 11 2017 © 1996-2017 PDE Solutions Inc. === C:\Users\USER\Desktop\Neuer Ordner\Standard_Descriptor.pde 20:38:23 4/20/17 ================================================================ 1| TITLE 'Random dispersed spheres distribution test' { the problem identification } 2| 3| COORDINATES cartesian2 { coordinate system, 1D,2D,3D, etc } 4| 5| VARIABLES { system variables } 6| Temp 7| 8| SELECT 9| 10| NGRID = 101 !Crashes with 100, but not with 101 or 99 11| 12| DEFINITIONS 13| k 14| kMatrix = 0.01 15| kInlay = 1 16| Q = 0 17| 18| xLen = 1.0 19| yLen = 0.1 20| tHi = 50 21| tLo = 10 22| delT = abs(tHi - tLo) 23| 24| 25| {Modifikation Beginn} 26| stagenum = 7 27| number = 11 28| xOff = array(-0.46082523239,-0.220955057384,0.849612894687,0.351706571907,-0.84846261354,0.537128201088,0.131703067401,0.663488049349,-0.618867468533,-0.0565239326678,-0.336797341775) 29| yOff = array(0.0070427289739,0.028320953656,-0.0178591230337,-0.0068958166222,-0.0276488104398,0.0065735254396,-0.0182951924,-0.0174390360647,0.0192620623572,0.0209511601287,-0.0266391810126) 30| volFrac = 0.3 31| run = 36 32| {Modifikation Ende} 33| 34| 35| 36| meshArea = 4 * xLen * yLen 37| totInlayArea = meshArea * volFrac / number 38| inRad = sqrt(totInlayArea / 3.14152926) 39| 40| 41| 42| EQUATIONS 43| k*div(-grad(Temp)) = Q 44| 45| BOUNDARIES 46| 47| REGION 1 'Matrixebene' 48| k = kMatrix 49| START(-xLen , -yLen) 50| VALUE(Temp) = tHi LINE TO (xLen , -yLen) 51| NATURAL(Temp) = 0 LINE TO (xLen , yLen) 52| VALUE(Temp) = tLo LINE TO (-xLen , yLen) 53| NATURAL(Temp) = 0 LINE TO CLOSE 54| 55| REGION 2 'Kugelebene' 56| k = kInlay 57| repeat i = 1 by 1 to number 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 58| START(xOff[i] - inRad, yOff[i]) ARC(CENTER = xOff[i], yOff[i]) ANGLE = 360 TO CLOSE 59| endrepeat 60| 61| FEATURE 62| START 'SurfKalt' ( -xLen , -yLen) LINE TO (xLen , -yLen) 63| START 'SurfWarm' (-xLen , yLen) LINE TO (xLen , yLen) 64| 65| MONITORS 66| 67| SUMMARY EXPORT FILE = "Sim_Info_" + $stagenum + "_avgrun=" + $run + ".txt" 68| REPORT(" ") 69| REPORT(" ") 70| REPORT("Report data") 71| REPORT('-----') 72| REPORT (LINE_INTEGRAL((normal(-k*grad(Temp))), 'SurfWarm') + LINE_INTEGRAL(normal(k*grad(Temp)), 'SurfKalt'))/2/delT*yLen/xLen as 'Mittlere Waermeleitfaehigkeit des Compound-Materials' 73| REPORT volFrac as 'Volume Fraction' 74| REPORT inRad as 'Inlay Radius' 75| REPORT number as '# Inlays/filler particles' 76| 77| END Threads=4 Order=3 Speed=1 Stability=0 E6=1 E7=0 Ej=1 FlexPDE Professional Trial 7.04/W64 3D (22 days) © 1996-2017 PDE Solutions Inc. FlexPDE Professional Trial 7.04/W64 3D (22 days) © 1996-2017 PDE Solutions Inc. SUMMARIZE 18 PARAMETERS (NO COMMON SUBEXPRESSION ELIMINATION) Param Compartment Complexity Memory -------------------------------------------------------------- K 1 1 80 K 2 1 80 KMATRIX 0 1 80 KINLAY 0 1 80 Q 0 1 80 XLEN 0 1 80 YLEN 0 1 80 THI 0 1 80 TLO 0 1 80 DELT 0 1 80 STAGENUM 0 1 80 NUMBER 0 1 80 XOFF 0 1 80 YOFF 0 1 80 VOLFRAC 0 1 80 RUN 0 1 80 MESHAREA 0 1 80 TOTINLAYAREA 0 1 80 INRAD 0 1 80 -------------------------- SUMMARIZE EQUATIONS -------------------------- DEFINES DEPENDS BASE CSE EQUATION VAR COMPARTMENT VAR FUNCTION COMPLEXITY COMPLEXITY MEMORY ------------------------------------------------------------------------------------------ 1 TEMP 1 Source 1 2 160 1 TEMP 1 Galerkin 9 9 720 1 TEMP 1 TEMP Jacobian 9 9 720 1 TEMP 1 Surf Galerkin 10 10 800 1 TEMP 1 TEMP Surf Jacobian 10 10 800 1 TEMP 2 Source 1 2 160 1 TEMP 2 Galerkin 5 5 400 1 TEMP 2 TEMP Jacobian 5 5 400 1 TEMP 2 Surf Galerkin 6 6 480 1 TEMP 2 TEMP Surf Jacobian 6 6 480 EQUATION 1 PARAMETER DEPENDENCIES: K KMATRIX KINLAY Q BC 1 REGION 1 BC TIME DERIVATIVE | TEMP.T ---------------------------------------- BC 2 REGION 1 BC TIME DERIVATIVE | ZERO ---------------------------------------- BC 3 REGION 1 BC TIME DERIVATIVE | TEMP.T ---------------------------------------- BOUNDARY CONDITION TAGS with associated BC INDEX BC TAG 1: val(1) BC TAG 2: nat(2) BC TAG 3: val(3) BOUNDARY CONDITION FORMULAS [BCINDEX] BC INDEX 1 VAL(TEMP) | THI ---------------------------------------- REGIONS 1 | 1 BC INDEX 2 NAT(TEMP) | 0.00000 ---------------------------------------- REGIONS 1 | 1 BC INDEX 3 VAL(TEMP) | TLO ---------------------------------------- REGIONS 1 | 1 -------------------------- SUMMARIZE BOUNDARY CONDITIONS -------------------------- BC DEFINES COMPARTMENT DEPENDS FUNCTION COMPLEXITY MEMORY ---------------------------------------------------------------------- 1 TEMP 1 Surf Galerkin 3 240 1 TEMP 1 TEMP Surf Jacobian 1 80 2 TEMP 1 Surf Galerkin 1 80 2 TEMP 1 TEMP Surf Jacobian 1 80 3 TEMP 1 Surf Galerkin 3 240 3 TEMP 1 TEMP Surf Jacobian 1 80 -------------------------------- SUMMARIZE FITS AND INTEGRALS --------------------------------- REGION/ BASE CSE FIT KIND COMPARTMENT COMPLEXITY COMPLEXITY MEMORY DEPENDENCY ----------------------------------------------------------------------------------------- 1 Boundary Integal 1 12 12 960 TEMP 2 Boundary Integal 1 11 11 880 TEMP ------------------------------- SUMMARIZE PLOT EQUATIONS ------------------------------- BASE CSE PLOT ARGUMENT COMPARTMENT COMPLEXITY COMPLEXITY MEMORY ------------------------------------------------------------------ ==UNUSED== PARAMETER DEFINITIONS: XOFF YOFF Lanczos/DInv FINAL ITERATION 11 of 1000 RNORM=0 XNORM=1.24181e-015 Flag=0 MS=1 Orthomin/ILU FINAL ITERATION 52 of 1000 RNORM=1.16284e-007 XNORM=4.66729 Flag=0 MS=64 Orthomin/ILU FINAL ITERATION 45 of 1000 RNORM=1.28276e-012 XNORM=1.22904e-005 Flag=0 MS=57 Orthomin/ILU FINAL ITERATION 52 of 1000 RNORM=2.29173e-018 XNORM=8.40759e-011 Flag=0 MS=64 RMS Spatial Errors 5.268e-4 DO_NEXTSTAGE pde 000000000490B950 MAXSTAGE=0