------------------------------------------------------------------------------------- (rtest): Running Test 9 (2D, viscous, thermal nonequilibrium flow over a cylinder with adiabatic noncatalytic wall) -------------------------------------------------------------------------------------  Initializing $Id: physics.C 23608 2011-09-09 21:14:34Z roystgnr $ Initializing Properties Databases ----------------------------- | Flowfield Mesh Information: | ----------------------------- Mesh Information: mesh_dimension()=2 spatial_dimension()=3 n_nodes()=1891 n_local_nodes()=258 n_elem()=1800 n_local_elem()=225 n_active_elem()=1800 n_subdomains()=1 n_processors()=8 processor_id()=0 -------------------------- | Boundary ID Information: | -------------------------- Side Boundary conditions: ------------------------- (ID, number of sides) (0, 60) (1, 60) (4, 60) Initializing $Id: comp_ns.C 22870 2011-08-15 12:10:42Z roystgnr $ --------------------------- | Boundary Mesh Information: | --------------------------- Mesh Information: mesh_dimension()=1 spatial_dimension()=3 n_nodes()=0 n_local_nodes()=0 n_elem()=0 n_local_elem()=0 n_active_elem()=0 n_subdomains()=0 n_processors()=8 processor_id()=0 ModelDimensions::Thermochemistry::n_species() = 5 ModelDimensions::Thermochemistry::n_energies() = 2 Constructing LaminarFlowTurbulenceModel... Using libMesh NewtonSolver >>> Restarting from restart.xdr assembling boundary projection system Writing output 0 Running $Id: physics.C 23608 2011-09-09 21:14:34Z roystgnr $ ----- time step 0, dt=100, time=100 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 91228.3 Continuing... NL conv: |u|_oo = 4.434835e+04 |du|_2 = 2.268772e+01 |du|_oo = 3.146486e+00 |du/dt|_oo = 3.146486e-02 ~~~ Cumulative Walltime: 6.180430e-01 ~~~ ----- time step 1, dt=1.000000e+02, time=2.000000e+02 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 1.759659e+04 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434845e+04 |du|_2 = 5.013129e+00 |du|_oo = 5.045877e-01 |du/dt|_oo = 5.045877e-03 ~~~ Cumulative Walltime: 1.229435e+00 ~~~ ----- time step 2, dt=2.000000e+02, time=4.000000e+02 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 1.153243e+04 Continuing... NL conv: |u|_oo = 4.434844e+04 |du|_2 = 2.100006e+00 |du|_oo = 3.648564e-01 |du/dt|_oo = 1.824282e-03 ~~~ Cumulative Walltime: 1.841689e+00 ~~~ ----- time step 3, dt=2.000000e+02, time=6.000000e+02 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 1.478200e+04 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434845e+04 |du|_2 = 1.506919e+00 |du|_oo = 2.619569e-01 |du/dt|_oo = 1.309784e-03 ~~~ Cumulative Walltime: 2.453496e+00 ~~~ ----- time step 4, dt=4.000000e+02, time=1.000000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 9.551323e+03 Continuing... NL conv: |u|_oo = 4.434845e+04 |du|_2 = 8.615058e-01 |du|_oo = 1.397905e-01 |du/dt|_oo = 3.494762e-04 ~~~ Cumulative Walltime: 3.069896e+00 ~~~ assembling boundary projection system Writing output 1 ----- time step 5, dt=4.000000e+02, time=1.400000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 1.072224e+04 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434846e+04 |du|_2 = 5.850528e-01 |du|_oo = 9.455416e-02 |du/dt|_oo = 2.363854e-04 ~~~ Cumulative Walltime: 3.946831e+00 ~~~ ----- time step 6, dt=8.000000e+02, time=2.200000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 8.190157e+03 Continuing... NL conv: |u|_oo = 4.434846e+04 |du|_2 = 3.497105e-01 |du|_oo = 5.302859e-02 |du/dt|_oo = 6.628574e-05 ~~~ Cumulative Walltime: 4.563271e+00 ~~~ ----- time step 7, dt=8.000000e+02, time=3.000000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 7.944678e+03 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434846e+04 |du|_2 = 2.319322e-01 |du|_oo = 3.547225e-02 |du/dt|_oo = 4.434032e-05 ~~~ Cumulative Walltime: 5.178224e+00 ~~~ ----- time step 8, dt=1.600000e+03, time=4.600000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 6.967778e+03 Continuing... NL conv: |u|_oo = 4.434846e+04 |du|_2 = 1.394548e-01 |du|_oo = 2.058217e-02 |du/dt|_oo = 1.286386e-05 ~~~ Cumulative Walltime: 5.795655e+00 ~~~ ----- time step 9, dt=1.600000e+03, time=6.200000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 5.892667e+03 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434847e+04 |du|_2 = 9.053161e-02 |du|_oo = 1.353506e-02 |du/dt|_oo = 8.459410e-06 ~~~ Cumulative Walltime: 6.410695e+00 ~~~ assembling boundary projection system Writing output 2 ----- time step 10, dt=3.200000e+03, time=9.400000e+03 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 5.767654e+03 Continuing... NL conv: |u|_oo = 4.434847e+04 |du|_2 = 5.498902e-02 |du|_oo = 8.045197e-03 |du/dt|_oo = 2.514124e-06 ~~~ Cumulative Walltime: 7.293969e+00 ~~~ ----- time step 11, dt=3.200000e+03, time=1.260000e+04 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 4.422258e+03 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434847e+04 |du|_2 = 3.524969e-02 |du|_oo = 5.216150e-03 |du/dt|_oo = 1.630047e-06 ~~~ Cumulative Walltime: 7.907696e+00 ~~~ ----- time step 12, dt=6.400000e+03, time=1.900000e+04 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 4.680372e+03 Continuing... NL conv: |u|_oo = 4.434847e+04 |du|_2 = 2.158940e-02 |du|_oo = 3.149231e-03 |du/dt|_oo = 4.920674e-07 ~~~ Cumulative Walltime: 8.522687e+00 ~~~ ----- time step 13, dt=6.400000e+03, time=2.540000e+04 ----- Nonlinear solver reached maximum step 1, latest evaluated residual 3.356089e+03 Continuing... ~~~ Scaling new time step by 2.000000e+00 ~~~ NL conv: |u|_oo = 4.434847e+04 |du|_2 = 1.370593e-02 |du|_oo = 2.019286e-03 |du/dt|_oo = 3.155134e-07 ~~~ Cumulative Walltime: 9.136192e+00 ~~~ --> User-specified time convergence criteria satisfied... ------------------------------------------------------------------- | Processor id: 0 | | Num Processors: 8 | | Time: Tue Oct 4 09:27:52 2011 | | OS: Linux | | HostName: navier.pecos | | OS Release: 2.6.18-164.11.1.el5_lustre.1.8.3 | | OS Version: #1 SMP Fri Apr 9 18:00:39 MDT 2010 | | Machine: x86_64 | | Username: roystgnr | | Configuration: ./configure run on Mon Oct 3 11:04:17 CDT 2011 | ------------------------------------------------------------------- ----------------------------------------------------------------------------------------------------------- | Physics Performance: Alive time=9.13944, Active time=9.13619 | ----------------------------------------------------------------------------------------------------------- | Event nCalls Total Time Avg Time Total Time Avg Time % of Active Time | | w/o Sub w/o Sub With Sub With Sub w/o S With S | |-----------------------------------------------------------------------------------------------------------| | | | run() 1 9.1362 9.136194 9.1362 9.136194 100.00 100.00 | ----------------------------------------------------------------------------------------------------------- | Totals: 1 9.1362 100.00 | ----------------------------------------------------------------------------------------------------------- Finalizing $Id: physics.C 23608 2011-09-09 21:14:34Z roystgnr $ ---------------------------------------------------------------------------------------------------------------- | libMesh Performance: Alive time=10.2802, Active time=9.8875 | ---------------------------------------------------------------------------------------------------------------- | Event nCalls Total Time Avg Time Total Time Avg Time % of Active Time | | w/o Sub w/o Sub With Sub With Sub w/o S With S | |----------------------------------------------------------------------------------------------------------------| | | | Mesh Setup 1 0.0019 0.001890 0.0234 0.023388 0.02 0.24 | | output() 3 0.0562 0.018730 0.6577 0.219226 0.57 6.65 | | solve() 14 0.0009 0.000064 8.5995 0.614252 0.01 86.97 | | | | BoundaryInfo | | sync() 1 0.0003 0.000251 0.0006 0.000610 0.00 0.01 | | | | DofMap | | add_neighbors_to_send_list() 4 0.0008 0.000189 0.0010 0.000244 0.01 0.01 | | compute_sparsity() 2 0.0143 0.007132 0.0164 0.008225 0.14 0.17 | | create_dof_constraints() 4 0.0023 0.000563 0.0023 0.000563 0.02 0.02 | | distribute_dofs() 4 0.0033 0.000816 0.0182 0.004552 0.03 0.18 | | dof_indices() 46563 0.0428 0.000001 0.0428 0.000001 0.43 0.43 | | prepare_send_list() 4 0.0001 0.000026 0.0001 0.000026 0.00 0.00 | | reinit() 4 0.0127 0.003187 0.0127 0.003187 0.13 0.13 | | | | EquationSystems | | build_solution_vector() 3 0.0154 0.005134 0.0286 0.009528 0.16 0.29 | | read() 1 0.4437 0.443731 0.4739 0.473930 4.49 4.79 | | update() 1 0.0002 0.000170 0.0002 0.000170 0.00 0.00 | | write() 3 0.4945 0.164825 0.5094 0.169788 5.00 5.15 | | | | ExodusII_IO | | write_nodal_data() 3 0.0125 0.004177 0.0125 0.004177 0.13 0.13 | | | | FE | | compute_face_map() 29638 0.1014 0.000003 0.2083 0.000007 1.03 2.11 | | compute_map() 36613 0.0337 0.000001 0.0337 0.000001 0.34 0.34 | | compute_shape_functions() 36613 0.0174 0.000000 0.0174 0.000000 0.18 0.18 | | init_face_shape_functions() 14 0.0000 0.000003 0.0000 0.000003 0.00 0.00 | | init_shape_functions() 29669 0.0987 0.000003 0.0987 0.000003 1.00 1.00 | | inverse_map() 59276 0.1002 0.000002 0.1002 0.000002 1.01 1.01 | | | | Matrix Assembly | | closing 14 0.0338 0.002417 0.0338 0.002417 0.34 0.34 | | | | Mesh | | find_neighbors() 2 0.0024 0.001224 0.0026 0.001282 0.02 0.03 | | renumber_nodes_and_elem() 4 0.0001 0.000026 0.0001 0.000026 0.00 0.00 | | | | MeshCommunication | | assign_global_indices() 4 0.0661 0.016517 0.0680 0.017001 0.67 0.69 | | compute_hilbert_indices() 1 0.0037 0.003682 0.0037 0.003682 0.04 0.04 | | find_global_indices() 1 0.0004 0.000369 0.0047 0.004721 0.00 0.05 | | parallel_sort() 1 0.0004 0.000436 0.0005 0.000489 0.00 0.00 | | | | MeshOutput | | write_equation_systems() 3 0.0001 0.000018 0.0412 0.013724 0.00 0.42 | | | | MetisPartitioner | | partition() 1 0.0034 0.003437 0.0082 0.008174 0.03 0.08 | | | | NewtonSolver | | solve() 14 0.3192 0.022801 8.5984 0.614175 3.23 86.96 | | | | Parallel | | allgather() 156 0.0082 0.000052 0.0083 0.000053 0.08 0.08 | | broadcast() 1279 0.0060 0.000005 0.0059 0.000005 0.06 0.06 | | gather() 698 0.0074 0.000011 0.0074 0.000011 0.08 0.08 | | max(bool) 59 0.6703 0.011362 0.6703 0.011362 6.78 6.78 | | max(scalar) 48 0.0009 0.000019 0.0009 0.000019 0.01 0.01 | | max(vector) 5 0.0000 0.000007 0.0000 0.000007 0.00 0.00 | | min(vector) 5 0.0003 0.000070 0.0003 0.000070 0.00 0.00 | | probe() 280 0.0009 0.000003 0.0009 0.000003 0.01 0.01 | | receive() 5848 0.0067 0.000001 0.0076 0.000001 0.07 0.08 | | send() 976 0.0016 0.000002 0.0016 0.000002 0.02 0.02 | | send_receive() 298 0.0008 0.000003 0.0035 0.000012 0.01 0.04 | | sum() 76 0.0684 0.000901 0.0684 0.000901 0.69 0.69 | | wait() 280 0.0001 0.000000 0.0001 0.000000 0.00 0.00 | | | | Partitioner | | set_node_processor_ids() 3 0.0006 0.000186 0.0011 0.000359 0.01 0.01 | | set_parent_processor_ids() 1 0.0001 0.000141 0.0001 0.000141 0.00 0.00 | | | | PetscLinearSolver | | solve() 17 0.6399 0.037643 0.6399 0.037643 6.47 6.47 | | | | Residual Assembly | | boundary residual 28616 1.8643 0.000065 2.2205 0.000078 18.85 22.46 | | compute stabilization nodal data 14 0.0480 0.003425 0.5735 0.040961 0.49 5.80 | | init DCO 3825 0.5599 0.000146 0.5599 0.000146 5.66 5.66 | | init nodal data 14 0.0286 0.002043 0.0286 0.002043 0.29 0.29 | | | | Residual+Jacobian Assembly | | compute AiAj 50400 0.2914 0.000006 0.2914 0.000006 2.95 2.95 | | compute AiVj 50400 0.1227 0.000002 0.1227 0.000002 1.24 1.24 | | compute system 3150 1.9787 0.000628 2.4040 0.000763 20.01 24.31 | | init data 3150 1.6888 0.000536 1.6888 0.000536 17.08 17.08 | | | | System | | assemble() 3 0.0006 0.000213 0.0006 0.000213 0.01 0.01 | | solve() 14 0.0002 0.000012 8.5986 0.614188 0.00 86.96 | | | | XdrIO | | read() 1 0.0092 0.009218 0.0095 0.009496 0.09 0.10 | ---------------------------------------------------------------------------------------------------------------- | Totals: 388099 9.8875 100.00 | ---------------------------------------------------------------------------------------------------------------- Successful exit Run complete ------------------------------------------------------------------------------------- (rtest): Running Solution Diff ------------------------------------------------------------------------------------- Solution vectors are equivalent to within tolerance = 5.0000000000000001e-09 Solution relative difference = 2.5288230700145018e-24 Successful exit Run complete ------------------------------------------------------------------------------------- (rtest): PASSED: Test 9 (2D, viscous, thermal nonequilibrium flow over a cylinder with adiabatic noncatalytic wall) -------------------------------------------------------------------------------------