************************************************************* Tem3P ACE3P Codes Source Date: Mon Mar 6 23:26:49 2023 -0800 ACE3P Codes Source Branch: master ACE3P Codes Source Tag: 503652f41066a31de4a9b7088dd1fd7572ada43a Support Lib Source Date: Fri Dec 2 09:43:53 2022 -0800 Support Lib Source Branch: master Support Lib Source Tag: 2f7bd8bf8ec6eb3646dc05e32622a4475531a105 Compilation Date: Mon 06 Mar 2023 11:38:55 PM PST ************************************************************* Copyright 2023, Stanford University Authors make no representations or warranties, expressed or implied. By way of example, but not limitation, authors make no representatinos or warranties of merchantibility or fitness for any particular purpose or that the use of the software componentns of documentation will not infringe any patents, copyrights, trademarks or other rights. The Authors shall not be held liable for any liability nor for any direct, indirect or consequential damages with respect to any claim by users or any third party on account of or arising from the use of this software. ************************************************************* Advanced Computations Department SLAC National Accelerator Laboratory https://slacportal.slac.stanford.edu/sites/ard_public/acd/Pages/Default.aspx Contact: ace3p@slac.stanford.edu Thank you for citing ACE3P when publishing related results. ************************************************************* Starting master process on nid004422 Number of MPI processes: 128 Number of compute nodes: 1 Number of processes per node: 128 Data precision: 64 bits Compiler: 11.2.0 20210728 (Cray Inc.) Beginning run at Sat Mar 11 13:51:39 2023 --------------------------------- Entering ThermoElasticProblem::solve()... --------------------------------- In ThermoElasticProblem::solve(): begining tem3p iter 0 with 'runThermostaticProblem' = 1 with 'runThermoTransientProblem' = 0 with 'runElasticProblem' = 0 1 In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: instantiating ScalarDiscreteDomain... 2 In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting up ScalarDiscreteDomain... Read Mesh: 2D-TTF3-FPC-Body.ncdf Check Mesh: 2D-TTF3-FPC-Body.ncdf Time for read model: 0.04120367299992722 Time for setBasisOrder: 0.001406965000114724 Time for making right-orientation and calculating Adjacency: 0.09187002299995584 Partitioning Method: PARMETIS Time for partitioning mesh: 0.08737734100031957 Time for replicating remote meshes: 0.02352289999998902 Time for numberingHierarchical: 0.0483747560001575 *********************************************************** * Total Number of Elements read: 315091 * Total Number of Elements used: 315091 * Total Number of DOFs: 101857 *********************************************************** Time for numbering remote nodes: 0.005346340999949462 Time for numbering remote edges: 0.01691337000011117 Time for numbering remote faces: 0.02491872500013415 Time for numbering remote volumes: 0.0003098279998994258 Time for numbering DOF entities: 0.0959137190002366 Time for set model: 0.2996180279997134 Time for read midpoints: 0.06859559200029253 Time for creating higher-order elements: 0.001468224999825907 Shell number 1 Shell boundary id 6 Shell basis order 1 Shell layer 1 Shell thickness 1e-05 Shell material 221 Shell number 2 Shell boundary id 7 Shell basis order 1 Shell layer 1 Shell thickness 0.0001 Shell material 221 loading RF solution for TEM3P... Omega3P directory name: s3p_results Boundary conditions: 0 = INTERIOR 1 = IMPEDANCE 2 = MAGNETIC 3 = MAGNETIC 4 = WAVEGUIDE 5 = WAVEGUIDE Material 1 has the following properties: Eps_r: (9, -0.0009) Mu_r: (1, 0) Sigma: 0 Boundary 1 has surface conductivity of 57000000.00000001 S/m Surface impedance model Metal properties: Epsilon = 1 MU = 1 Sigma = 57000000 Frequency = 1000000000 Read Mesh: /pscratch/sd/l/liling/cw23/tem3p/SrfCoupler/2D-TTF3-FPC-Vacuum.ncdf Time for reading the model: 0.3262526749999779 Using curved quadratic tetrahedrons Setting global vector finite element basis order to p=2 Partitioning Method: parmetis *********************************************************** * Total Number of Elements read: 139754 * Total Number of Elements used: 139754 * Total Number of DOFs: 968476 *********************************************************** Time for setting up finite element framework: 3.378400882999813 /********************************/ /* input parameters, KVC syntax */ /********************************/ Mesh : { ReplicatedElementDistribution : { stddev : 0.0 total : 0 min : 0 max : 0 average : 0.0 } File : /pscratch/sd/l/liling/cw23/tem3p/SrfCoupler/2D-TTF3-FPC-Vacuum.ncdf ElementDistribution : { stddev : 0.0 total : 139754 min : 139754 max : 139754 average : 139754.0 } MeshCoords : 32081 } /********************************/ In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting map of boundaries... EM frequency : (1300000000,0) Powerinput : 437.5 : : 437.5 The scaling factor is 437.5 Volumetric Heating Boundary 5 dielectri constants de(T) = EM frequency : (1300000000,0) Powerinput : 437.5 : : 437.5 The scaling factor is 437.5 EM Heating Boundary 6 surface resis R(T) = EM frequency : (1300000000,0) Powerinput : 437.5 : : 437.5 The scaling factor is 437.5 EM Heating Boundary 7 surface resis R(T) = EM frequency : (1300000000,0) Powerinput : 437.5 : : 437.5 The scaling factor is 437.5 EM Heating Boundary 8 surface resis R(T) = EM frequency : (1300000000,0) Powerinput : 437.5 : : 437.5 The scaling factor is 437.5 EM Heating Boundary 9 surface resis R(T) = In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: map of boundaries set succesfully In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting linear solver... Leaving StaticProblem::readKvcParameters1() of ThermostaticProblem Volume id: 221 Volume id: 1 Volume id: 2 Volume id: 3 Volume id: 4 Volume id: 5 Entering StaticProblem::instantiateMatrix() of ThermostaticProblem Entering assembleMatrix(). Time for assembling thermal matrix: 0.07940910700017412 Entering ThermostaticProblem::solveWithNewton() for global iteration 0 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1049824970000373 ThermostaticProblem::computeRHS Initial norm : 2941.062976517171 DOFS norm : 4138.193204769444 Picard limit : 100 Max Nonlinear Iter : 500 Absolute Tol : 1e-09 Relative Tol : 1e-09 Nonlinear Newton iteration 0 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1093818939998528 ThermostaticProblem::computeRHS Newton residual norm : 12.80728725514524 Initial norm : 2941.062976517171 DOFS norm : 102091.7014813783 Nonlinear Newton iteration 1 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1025023980000697 ThermostaticProblem::computeRHS Newton residual norm : 0.5354560505086164 Initial norm : 2941.062976517171 DOFS norm : 87059.88012526993 Nonlinear Newton iteration 2 Entering assembleMatrix(). Time for assembling thermal matrix: 0.102281266999853 ThermostaticProblem::computeRHS Newton residual norm : 0.01844334888494861 Initial norm : 2941.062976517171 DOFS norm : 86461.01621543053 Nonlinear Newton iteration 3 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1021798719998515 ThermostaticProblem::computeRHS Newton residual norm : 0.003293372469854504 Initial norm : 2941.062976517171 DOFS norm : 86261.50081796799 Nonlinear Newton iteration 4 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1020625679998375 ThermostaticProblem::computeRHS Newton residual norm : 0.0006122174728530146 Initial norm : 2941.062976517171 DOFS norm : 86224.59992749234 Nonlinear Newton iteration 5 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1018799940002282 ThermostaticProblem::computeRHS Newton residual norm : 0.0001051732992730229 Initial norm : 2941.062976517171 DOFS norm : 86219.06944617006 Nonlinear Newton iteration 6 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1022953290002988 ThermostaticProblem::computeRHS Newton residual norm : 7.905809985234996e-06 Initial norm : 2941.062976517171 DOFS norm : 86219.14645450425 Nonlinear Newton iteration 7 Entering assembleMatrix(). Time for assembling thermal matrix: 0.1018686960001105 ThermostaticProblem::computeRHS Newton residual norm : 3.299021151355816e-06 Initial norm : 2941.062976517171 DOFS norm : 86219.18121534205 Nonlinear Newton iteration 8 Entering assembleMatrix(). Time for assembling thermal matrix: 0.103043555000113 ThermostaticProblem::computeRHS Newton residual norm : 5.887218503899363e-07 Initial norm : 2941.062976517171 DOFS norm : 86219.15635929929 Leaving ThermostaticProblem::solveWithNewton() for global iteration 0 Calculating the total power flow through every sideset... ===================================== Max T = 324.3520007076649 Min T = 2 P (SideSet 1) = 0.08725785613942286 P (SideSet 2) = 0.0005649793014200234 P (SideSet 3) = 0.03957574428881457 P (SideSet 4) = 0.7813392638376635 P (SideSet 5) = 0.1005412253626891 P (SideSet 6) = 0 P (SideSet 7) = 0 P (SideSet 8) = -0.03257071919747352 P (SideSet 9) = -0.002889115633092689 P (SideSet 10) = -0.1809940299856579 P (SideSet 11) = -0.2004387158296005 ===================================== Temperature output file : tem3p_results/Temperature.mod --------------------------------- Leaving ThermoElasticProblem::solve() --------------------------------- Leaving StaticProblem::destructor() of ThermostaticProblem Ending run at Sat Mar 11 14:05:29 2023 Total run time = 830 seconds