Field_2D_Types_Ave_Lib(1) (Echo of msh_remarks.txt) =================== Begin Application Remarks =========================== Patch test with eight T6 elements Numerical check T6 isoparametric Poisson element via patch test with only interior nodes 7, 10-12, 15, 22-25 unknown 5--8-13-16-21 (4,4) Mesh shown to left. : / : / : Exact solution u = 1 + 3x - 4y 4 24 12 25 20 :/ :/ : du/dx = 3, du/dy = -4 3--7-11-15-19 are constant everywhere : / : / : 2 22 10 23 18 :/ :/ : 1--6--9-14--17 ->X (0,0) Nodal dof: u-temperature Element type = 1 (a Poisson, everywhere) Element connection: six nodes per element Element properties 4 (columns in msh_properties.tmp): conductivity (k), volumetric source (Q), density (Rho), thickness (Thick) ==================== End Application Remarks ============================ Begin Phase 2, input mesh, properties, sources and EBC Read 25 nodes. (Echo node, file msh_bc_xyz.txt) node, bc_flags, 2 coordinates 1, 1 0 0 2, 1 0 1 3, 1 0 2 4, 1 0 3 5, 1 0 4 6, 1 1 0 7, 0 1 2 8, 1 1 4 9, 1 2 0 10, 0 2 1 11, 0 2 2 12, 0 2 3 13, 1 2 4 14, 1 3 0 15, 0 3 2 16, 1 3 4 17, 1 4 0 18, 1 4 1 19, 1 4 2 20, 1 4 3 21, 1 4 4 22, 0 1 1 23, 0 3 1 24, 0 1 3 25, 0 3 3 Read 8 elements with its type & 6 nodes each. Element number, type, connection list 1, 1 1 9 11 6 10 22 2, 1 1 11 3 22 7 2 3, 1 9 17 19 14 18 23 4, 1 9 19 11 23 15 10 5, 1 3 11 13 7 12 24 6, 1 3 13 5 24 8 4 7, 1 11 19 21 15 20 25 8, 1 11 21 13 25 16 12 Read 1 materials with 5 properties 5 homogeneous material values per elem type (Echoing text file msh_properties.txt) 1 1 0 0 0 Above properties are: TYPE 1 or 2: Kx, Ky, G, Q (Q2,...,Qn) For -[d/dx (Kx du/dx) + d/dy (Ky du/dy)] + G u - Q = 0 TYPE 3: 1=q_in, 2==0, 3==0, 4==0 For -Kn du/dn = q_in (postive when entering) TYPE 4: 1=g_in, 2==h_in, 3==0, 4==0 For -Kn du/dn = g_in + h_in * u, where for convection: h_in = h the convection coeff), g_in = h * u_infinity Note: expecting 16 essential BC values. Read 16 EBC data sets (Echo of file text file msh_ebc.txt) Node, DOF, Value. 1 1 1 2 1 -3 3 1 -7 4 1 -11 5 1 -15 6 1 4 8 1 -12 9 1 7 13 1 -9 14 1 10 16 1 -6 17 1 13 18 1 9 19 1 5 20 1 1 21 1 -3 CReated file plot_input_2d_mesh.png NOTE: no debug outputs requested. Begin Phase 3, assemble (scatter) element arrays Displaying the first 2 sets of element matrices Element matrix S_e = 0.5000 0.1667 0.0000 -0.6667 -0.0000 -0.0000 0.1667 1.0000 0.1667 -0.6667 -0.6667 -0.0000 0.0000 0.1667 0.5000 -0.0000 -0.6667 -0.0000 -0.6667 -0.6667 -0.0000 2.6667 0.0000 -1.3333 -0.0000 -0.6667 -0.6667 0.0000 2.6667 -1.3333 -0.0000 -0.0000 -0.0000 -1.3333 -1.3333 2.6667 Element matrix C_e = 0 0 0 0 0 0 Element matrix S_e = 0.5000 -0.0000 0.1667 0.0000 -0.0000 -0.6667 -0.0000 0.5000 0.1667 0.0000 -0.6667 -0.0000 0.1667 0.1667 1.0000 0 -0.6667 -0.6667 0.0000 0.0000 0 2.6667 -1.3333 -1.3333 -0.0000 -0.6667 -0.6667 -1.3333 2.6667 0.0000 -0.6667 -0.0000 -0.6667 -1.3333 0.0000 2.6667 Element matrix C_e = 0 0 0 0 0 0 Total input source terms = 0 Begin Phase 4, modify matrices for EBC Begin Phase 5, solve modified system Computed Solution: Node, 1 results per node 1, 1.000e+00 2, -3.000e+00 3, -7.000e+00 4, -1.100e+01 5, -1.500e+01 6, 4.000e+00 7, -4.000e+00 8, -1.200e+01 9, 7.000e+00 10, 3.000e+00 11, -1.000e+00 12, -5.000e+00 13, -9.000e+00 14, 1.000e+01 15, 2.000e+00 16, -6.000e+00 17, 1.300e+01 18, 9.000e+00 19, 5.000e+00 20, 1.000e+00 21, -3.000e+00 22, -1.323e-15 23, 6.000e+00 24, -8.000e+00 25, -2.000e+00 Created file color_scalar_result.png Created file result_surface_38.png Reactions at essential BCs Node, DOF, Reaction Value 1 1 3.3333e-01 2 1 -4.0000e+00 3 1 -2.0000e+00 4 1 -4.0000e+00 5 1 -2.3333e+00 6 1 5.3333e+00 8 1 -5.3333e+00 9 1 2.6667e+00 13 1 -2.6667e+00 14 1 5.3333e+00 16 1 -5.3333e+00 17 1 2.3333e+00 18 1 4.0000e+00 19 1 2.0000e+00 20 1 4.0000e+00 21 1 -3.3333e-01 Totals = 3.7970e-14 Begin Phase 6, use the answers Element Post-processing: El, Q_Pt, Coordinates 1, 1, 1.33333 0.666667 El, Q_Pt, Gradient Vector 1, 1, 3 -4 El, Q_Pt, Flux Vector 1, 1, -3 4 El, Q_Pt, Coordinates 1, 2, 0.8 0.4 El, Q_Pt, Gradient Vector 1, 2, 3 -4 El, Q_Pt, Flux Vector 1, 2, -3 4 El, Q_Pt, Coordinates 1, 3, 1.6 1.2 El, Q_Pt, Gradient Vector 1, 3, 3 -4 El, Q_Pt, Flux Vector 1, 3, -3 4 El, Q_Pt, Coordinates 1, 4, 1.6 0.4 El, Q_Pt, Gradient Vector 1, 4, 3 -4 El, Q_Pt, Flux Vector 1, 4, -3 4 El, Q_Pt, Coordinates 2, 1, 0.666667 1.33333 El, Q_Pt, Gradient Vector 2, 1, 3 -4 El, Q_Pt, Flux Vector 2, 1, -3 4 El, Q_Pt, Coordinates 2, 2, 0.4 0.8 El, Q_Pt, Gradient Vector 2, 2, 3 -4 El, Q_Pt, Flux Vector 2, 2, -3 4 El, Q_Pt, Coordinates 2, 3, 0.4 1.6 El, Q_Pt, Gradient Vector 2, 3, 3 -4 El, Q_Pt, Flux Vector 2, 3, -3 4 El, Q_Pt, Coordinates 2, 4, 1.2 1.6 El, Q_Pt, Gradient Vector 2, 4, 3 -4 El, Q_Pt, Flux Vector 2, 4, -3 4 El, Q_Pt, Coordinates 3, 1, 3.33333 0.666667 El, Q_Pt, Gradient Vector 3, 1, 3 -4 El, Q_Pt, Flux Vector 3, 1, -3 4 El, Q_Pt, Coordinates 3, 2, 2.8 0.4 El, Q_Pt, Gradient Vector 3, 2, 3 -4 El, Q_Pt, Flux Vector 3, 2, -3 4 El, Q_Pt, Coordinates 3, 3, 3.6 1.2 El, Q_Pt, Gradient Vector 3, 3, 3 -4 El, Q_Pt, Flux Vector 3, 3, -3 4 El, Q_Pt, Coordinates 3, 4, 3.6 0.4 El, Q_Pt, Gradient Vector 3, 4, 3 -4 El, Q_Pt, Flux Vector 3, 4, -3 4 El, Q_Pt, Coordinates 4, 1, 2.66667 1.33333 El, Q_Pt, Gradient Vector 4, 1, 3 -4 El, Q_Pt, Flux Vector 4, 1, -3 4 El, Q_Pt, Coordinates 4, 2, 2.4 0.8 El, Q_Pt, Gradient Vector 4, 2, 3 -4 El, Q_Pt, Flux Vector 4, 2, -3 4 El, Q_Pt, Coordinates 4, 3, 2.4 1.6 El, Q_Pt, Gradient Vector 4, 3, 3 -4 El, Q_Pt, Flux Vector 4, 3, -3 4 El, Q_Pt, Coordinates 4, 4, 3.2 1.6 El, Q_Pt, Gradient Vector 4, 4, 3 -4 El, Q_Pt, Flux Vector 4, 4, -3 4 El, Q_Pt, Coordinates 5, 1, 1.33333 2.66667 El, Q_Pt, Gradient Vector 5, 1, 3 -4 El, Q_Pt, Flux Vector 5, 1, -3 4 El, Q_Pt, Coordinates 5, 2, 0.8 2.4 El, Q_Pt, Gradient Vector 5, 2, 3 -4 El, Q_Pt, Flux Vector 5, 2, -3 4 El, Q_Pt, Coordinates 5, 3, 1.6 3.2 El, Q_Pt, Gradient Vector 5, 3, 3 -4 El, Q_Pt, Flux Vector 5, 3, -3 4 El, Q_Pt, Coordinates 5, 4, 1.6 2.4 El, Q_Pt, Gradient Vector 5, 4, 3 -4 El, Q_Pt, Flux Vector 5, 4, -3 4 El, Q_Pt, Coordinates 6, 1, 0.666667 3.33333 El, Q_Pt, Gradient Vector 6, 1, 3 -4 El, Q_Pt, Flux Vector 6, 1, -3 4 El, Q_Pt, Coordinates 6, 2, 0.4 2.8 El, Q_Pt, Gradient Vector 6, 2, 3 -4 El, Q_Pt, Flux Vector 6, 2, -3 4 El, Q_Pt, Coordinates 6, 3, 0.4 3.6 El, Q_Pt, Gradient Vector 6, 3, 3 -4 El, Q_Pt, Flux Vector 6, 3, -3 4 El, Q_Pt, Coordinates 6, 4, 1.2 3.6 El, Q_Pt, Gradient Vector 6, 4, 3 -4 El, Q_Pt, Flux Vector 6, 4, -3 4 El, Q_Pt, Coordinates 7, 1, 3.33333 2.66667 El, Q_Pt, Gradient Vector 7, 1, 3 -4 El, Q_Pt, Flux Vector 7, 1, -3 4 El, Q_Pt, Coordinates 7, 2, 2.8 2.4 El, Q_Pt, Gradient Vector 7, 2, 3 -4 El, Q_Pt, Flux Vector 7, 2, -3 4 El, Q_Pt, Coordinates 7, 3, 3.6 3.2 El, Q_Pt, Gradient Vector 7, 3, 3 -4 El, Q_Pt, Flux Vector 7, 3, -3 4 El, Q_Pt, Coordinates 7, 4, 3.6 2.4 El, Q_Pt, Gradient Vector 7, 4, 3 -4 El, Q_Pt, Flux Vector 7, 4, -3 4 El, Q_Pt, Coordinates 8, 1, 2.66667 3.33333 El, Q_Pt, Gradient Vector 8, 1, 3 -4 El, Q_Pt, Flux Vector 8, 1, -3 4 El, Q_Pt, Coordinates 8, 2, 2.4 2.8 El, Q_Pt, Gradient Vector 8, 2, 3 -4 El, Q_Pt, Flux Vector 8, 2, -3 4 El, Q_Pt, Coordinates 8, 3, 2.4 3.6 El, Q_Pt, Gradient Vector 8, 3, 3 -4 El, Q_Pt, Flux Vector 8, 3, -3 4 El, Q_Pt, Coordinates 8, 4, 3.2 3.6 El, Q_Pt, Gradient Vector 8, 4, 3 -4 El, Q_Pt, Flux Vector 8, 4, -3 4 ans = 0 Total domain volume = 16 Begin SCP averaging flux at the nodes NOTE: maximum number of element neighbors = 6 NOTE: All elements have neighbors Range of averaged flux values at mesh nodes FLUX #: MAXIMUM, NODE MINIMUM, NODE 1: -3.0000e+00, 6 -3.0000e+00, 17 2: 4.0000e+00, 17 4.0000e+00, 10 Patch averaged fluxes at mesh nodes: Node, 2 flux components per node 1 -3.000e+00 4.000e+00 2 -3.000e+00 4.000e+00 3 -3.000e+00 4.000e+00 4 -3.000e+00 4.000e+00 5 -3.000e+00 4.000e+00 6 -3.000e+00 4.000e+00 7 -3.000e+00 4.000e+00 8 -3.000e+00 4.000e+00 9 -3.000e+00 4.000e+00 10 -3.000e+00 4.000e+00 11 -3.000e+00 4.000e+00 12 -3.000e+00 4.000e+00 13 -3.000e+00 4.000e+00 14 -3.000e+00 4.000e+00 15 -3.000e+00 4.000e+00 16 -3.000e+00 4.000e+00 17 -3.000e+00 4.000e+00 18 -3.000e+00 4.000e+00 19 -3.000e+00 4.000e+00 20 -3.000e+00 4.000e+00 21 -3.000e+00 4.000e+00 22 -3.000e+00 4.000e+00 23 -3.000e+00 4.000e+00 24 -3.000e+00 4.000e+00 25 -3.000e+00 4.000e+00 NOTE: ave fluxes saved to file scp_node_ave_fluxes.txt For color_scp_flux_step.m,quiver_scp_ave_pt_flux.m etc Using a scale of 1 and vector increment of 1 Read 25 mesh coordinate pairs Read 8 elements connections Retained the first 8 with type < 3 Read 25 nodal flux sets Created file scp_flux_vec_mesh.png quit