************************************************************* 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 nid006472 Number of MPI processes: 16 Number of compute nodes: 1 Number of processes per node: 16 Data precision: 64 bits Compiler: 11.2.0 20210728 (Cray Inc.) Beginning run at Sat Mar 11 14:53:25 2023 --------------------------------- Entering ThermoElasticProblem::solve()... --------------------------------- In ThermoElasticProblem::solve(): begining tem3p iter 0 with 'runThermostaticProblem' = 1 with 'runThermoTransientProblem' = 0 with 'runElasticProblem' = 1 1 In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: instantiating ScalarDiscreteDomain... 2 In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting up ScalarDiscreteDomain... Read Mesh: RfGunBody.ncdf Check Mesh: RfGunBody.ncdf Time for read model: 0.04009315199982666 Time for setBasisOrder: 0.007313972000247304 Time for making right-orientation and calculating Adjacency: 0.08057344100006958 Partitioning Method: PARMETIS Time for partitioning mesh: 0.07195291299967721 Time for replicating remote meshes: 0.01989548599976843 Time for numberingHierarchical: 0.1655464119999124 *********************************************************** * Total Number of Elements read: 229035 * Total Number of Elements used: 229035 * Total Number of DOFs: 337804 *********************************************************** Time for numbering remote nodes: 0.003500390000226616 Time for numbering remote edges: 0.008915444999729516 Time for numbering remote faces: 0.01336482100032299 Time for numbering remote volumes: 0.0009816209999371495 Time for numbering DOF entities: 0.1924188809998668 Time for set model: 0.3690928759997405 Time for read midpoints: 0.00795454500030246 Time for creating higher-order elements: 0.01025465300017459 loading RF solution for TEM3P... Omega3P directory name: omega3p_results Boundary conditions: 0 = INTERIOR 1 = ELECTRIC 2 = MAGNETIC 3 = MAGNETIC 4 = MAGNETIC 5 = ELECTRIC 6 = EXTERIOR Read Mesh: /pscratch/sd/l/liling/cw23/tem3p/RfGun-Coupler/RfGunVacuum.ncdf Time for reading the model: 0.09417208300010316 Using curved quadratic tetrahedrons Setting global vector finite element basis order to p=2 Partitioning Method: parmetis *********************************************************** * Total Number of Elements read: 83124 * Total Number of Elements used: 83124 * Total Number of DOFs: 512188 *********************************************************** Time for setting up finite element framework: 1.223409560000164 /********************************/ /* input parameters, KVC syntax */ /********************************/ Mesh : { File : /pscratch/sd/l/liling/cw23/tem3p/RfGun-Coupler/RfGunVacuum.ncdf ReplicatedElementDistribution : { total : 0 min : 0 average : 0.0 stddev : 0.0 max : 0 } MeshCoords : 16207 ElementDistribution : { total : 83124 min : 83124 average : 83124.0 stddev : 0.0 max : 83124 } } /********************************/ In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting map of boundaries... EM frequency : (2853388398.852401,0) Parameter used in computing gradient: /********************************/ /* input parameters, KVC syntax */ /********************************/ SymmetryFactor : 1.0 Trajectory : { EndPoint : 0.0001 , 0.0001 , 0.017018 TrajectoryMethod : 1 NumberOfPointsPerSegment : 2000 StartPoint : 0.0001 , 0.0001 , -0.017018 } /********************************/ Number of sampled points: 2001 Vr, Vi, and Voltage on the line1: -3.202134019277217 5.130237752476776 6.047561630471444 Scaling Info : Scaling Factor : Target Gradient : EM gradient : Effective Length : Duty Factor : Symmety Scaling Info : 41050781.14681887 : 60000000 : 177.6813265504597 : 0.034036 : 0.00036 : 0 The scaling factor is 41050781.14681887 EM Heating Boundary 6 sigmaSqrt = 7615.773105863908 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 Entering StaticProblem::instantiateMatrix() of ThermostaticProblem Entering assembleMatrix(). Time for assembling thermal matrix: 0.2783120640001471 ThermostaticProblem::computeRHS In ThermostaticProblem::solveLinearSystem(), tem3p iter 0: solving Ax = b with ||b||_2 = 26.15985519308913 In ThermostaticProblem::solveLinearSystem(), tem3p iter 0: solved Ax = b with ||x||_2 = 14323.44863842977, after 1234 iterations 1 In StaticProblem::readKvcParameters1() of ElasticProblem, rank 0: instantiating ScalarDiscreteDomain... 2 In StaticProblem::readKvcParameters1() of ElasticProblem, rank 0: setting up ScalarDiscreteDomain... Read Mesh: RfGunBody.ncdf Time for read model: 0.03054097100039144 Time for setBasisOrder: 0.006184748000123363 Time for making right-orientation and calculating Adjacency: 0.09763487899999745 Partitioning Method: PARMETIS Time for partitioning mesh: 0.08246546299960755 Time for replicating remote meshes: 0.02270723899982841 Time for numberingHierarchical: 0.159927413999867 *********************************************************** * Total Number of Elements read: 229035 * Total Number of Elements used: 229035 * Total Number of DOFs: 1013412 *********************************************************** Time for numbering remote nodes: 0.002147114000308648 Time for numbering remote edges: 0.007858971000132442 Time for numbering remote faces: 0.01142940500039913 Time for numbering remote volumes: 0.0007405559999824618 Time for numbering DOF entities: 0.1821941349999179 Time for set model: 0.3892272819998652 Time for read midpoints: 0.007164934999764228 Time for creating higher-order elements: 0.01150815699975283 loading RF solution for TEM3P... In StaticProblem::readKvcParameters1() of ElasticProblem, rank 0: setting map of boundaries... box 1 : 0.0001 box 0.0001 : 0.0001 box 0.0001 : 0.0001 box -0.017018 : 0.017018 Parameter used in computing gradient: /********************************/ /* input parameters, KVC syntax */ /********************************/ Trajectory : { EndPoint : 0.0001 , 0.0001 , 0.017018 TrajectoryMethod : 1 NumberOfPointsPerSegment : 2000 StartPoint : 0.0001 , 0.0001 , -0.017018 } GradientDirection : 0.0 , 0.0 , 1.0 SymmetryFactor : 1.0 /********************************/ Number of sampled points: 2001 Vr, Vi, and Voltage on the line1: -3.202134019277217 5.130237752476776 6.047561630471444 Scaling Info : Gradient Scaling Factor : Target Gradient : EM gradient : Effective Length : Duty Factor Scaling Info : 337683.2060231193 : 60000000 : 177.6813265504597 : 0.034036 : 1 The LF scaling factor is 114029947630.0524 mixed boundary count is 3 In StaticProblem::readKvcParameters1() of ElasticProblem, rank 0: map of boundaries set succesfully In StaticProblem::readKvcParameters1() of ElasticProblem, rank 0: setting linear solver... Leaving StaticProblem::readKvcParameters1() of ElasticProblem EM Mesh input directory: omega3p_results Entering StaticProblem::instantiateMatrix() of ElasticProblem Entering StaticProblem::setSolver() of ElasticProblem Memory usage: used mem per MPI process: min: 373.6 MB, max: 407.8 MB, avg: 390.5 MB, stddev: 9.148 MB, total: 6248 MB used mem per node in GB : min: 32.44 GB, max: 32.44 GB, avg: 32.44 GB, stddev: 0 GB, total: 32.44 GB used mem per node in % : min: 6.447 %, max: 6.447 %, avg: 6.447 %, stddev: 2.752e-15 % Leaving StaticProblem::setSolver() of ElasticProblem Deformed Vacuum Mesh Calculation ... Boundary conditions: 0 = INTERIOR 1 = ELECTRIC 2 = MAGNETIC 3 = MAGNETIC 4 = MAGNETIC 5 = ELECTRIC 6 = EXTERIOR Read Mesh: /pscratch/sd/l/liling/cw23/tem3p/RfGun-Coupler/RfGunVacuum.ncdf Time for reading the model: 0.02153 Using curved quadratic tetrahedrons Setting global vector finite element basis order to p=2 Partitioning Method: parmetis *********************************************************** * Total Number of Elements read: 83124 * Total Number of Elements used: 83124 * Total Number of DOFs: 512188 *********************************************************** Time for setting up finite element framework: 0.1493 /********************************/ /* input parameters, KVC syntax */ /********************************/ Mesh : { File : /pscratch/sd/l/liling/cw23/tem3p/RfGun-Coupler/RfGunVacuum.ncdf ReplicatedElementDistribution : { total : 31438 min : 762 average : 1964.875 stddev : 624.18970139961 max : 3114 } MeshCoords : 16207 ElementDistribution : { total : 83124 min : 4633 average : 5195.25 stddev : 208.09597144907 max : 5329 } } /********************************/ Vacuum Mesh Exterior ID : 6 Size of Multiphysics surface coord list : 4065 Size of Multiphysics surface mid coord list : 6337 Calculating the total power flow through every sideset... ===================================== Max T = 28.6 Min T = 23.12 P (SideSet 1) = 2.408e-07 P (SideSet 2) = 0.001142 P (SideSet 3) = -0.01299 P (SideSet 4) = -0.05743 P (SideSet 5) = 213.1 P (SideSet 6) = -215.6 ===================================== Temperature output file : tem3p_results/Temperature.mod Strain/Stress Max : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : -1e+09 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : Strain/Stress Min : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 1e+09 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : 0 : writing strain file tem3p_results/Strain.mod writing stress file tem3p_results/Stress.mod ===================================== Min u_x = -1.182e-06 Max u_x = 4.556e-07 Min u_y = -7.828e-07 Max u_y = 2.647e-07 Min u_z = 0 Max u_z = 3.311e-06 Max u_mag = 3.312e-06 Max u_mag block id = 1 ===================================== Displacement output file : tem3p_results/Displacement.mod --------------------------------- Leaving ThermoElasticProblem::solve() --------------------------------- Leaving StaticProblem::destructor() of ElasticProblem Leaving StaticProblem::destructor() of ThermostaticProblem Ending run at Sat Mar 11 14:54:21 2023 Total run time = 56 seconds