Input Data of 2-D Adaptive Mesh Generator (for interaction with the MODEL error estimator code.) Reference: H.C. Huang, A.S. Usmani, Finite Element Analysis for Heat Transfer, London, Springer-Verlag, 1994. This program requires data files: h_geom.dat to define the geometry h_mesh.dat to define the types and values of boundary conditions densit.dat to control the element sizes to be created A. Homogeneous regions without holes File h_geom contains: <<----- <<----- <<----- <<----- <<----- 0 1 0.00001 3 ! holes, regions, tolerance, smoothing iterations 1 ! region number 4 ! number of line segments bounding the region 1 ! 1st segment type (1 line_to, 2 arc_to) 0.,0. 100.,0. ! segment end coordinates 1 ! 2nd segment type 100.,3. ! end coordinates, from last point 1 ! 3rd segment type 0.,3. ! end coordinates, from last point 1 ! last segment type, closes to 1st point ! other remarks can be stored here. ! Donald, E_P_S_V ! ! v ^ 3 ------------(seg 3)------------ ! | | ! (seg 4) Region 1 (seg 2) ! | | ! |- 0 ------------(seg 1)------------ 100 ----> s File h_mesh.dat contains: <<----- <<----- <<----- <<----- <<----- "Title: Donald, E_P_S_V, 1/31/03" 1. 1. ! fake material properties for first region (then 2nd, 3rd if any) 2 ! Number of different Dirichlet boundary condition values 10. ! 1st Dirichlet value (left edge) 0. ! 2nd Dirichlet value (right edge) 2 ! Number of segments with Dirichlet boundary conditions 2 +2 ! 1st segment number, DBC number (- for 1st point only) 4 +1 ! 1st segment number, DBC number (- for 1st point only) 0 ! Number of segements with non-zero Neumann conditions ! other remarks can be stored here. ! Neumann = 0 ! v ^ 3 ------------(seg 3)------------ ! | | ! DBC=10 (seg 4) Region 1 (seg 2) DBC = 0 ! | | ! |- 0 ------------(seg 1)------------ 100 ----> s ! Neumann = 0 File densit.dat contains (to also be read by Matlab): <<----- <<----- NUM_D 1 1 = Number of density conntrol points, two dummy unity terms. X, Y, El_size = Point where element size is given, and size. (Manual input for generation fo first mesh, but for later runs rename msh_new_el_size.tmp, from MODEL error estimator, to densit.dat.) Run with Mesh_adapt (or /net/screech/akin/h_mesh/new_mesh) Program will issue prompt: please input element type. 1 -- triangle with three nodes; 2 -- triangle with six nodes; 3 -- quadrilateral with four nodes; 1 Program responds with for subdomain 1 enter 0 for no smoothing enter x for x smoothing iterations (Ignore, and) NOTE: NODE FLAG & BC_VALUE SAVED TO msh_bc_values.dat NOTE: NODE FLAG & COORDINATES SAVED TO msh_bc_xyz.dat NOTE: CONNECTIVITY SAVED TO msh_typ_nodes.dat These are used by MODEL to run and Matlab to plot. ================================================================ Below are an example of the 3 input files. See ref for figures. ================================================================ (h_geom.dat) 0 1 0.00001 3 ! nbnds natrib tolrn nitsm [holes regions tol smooth] 1 ! region 4 ! noseg [number of sides for this hole] 1 ! ityp [first type 1 line 2 arc+- 3 user 4 other side] 0. 0. 160. 0. ! first line x,y x,y 1 ! ityp [type 1 line_to 2 arc_to+- 3 user 4 other side] 110. 120. ! ending point, radius 1 ! ityp [type 1 line_to 2 arc_to+- 3 user 4 other side] 50. 120. ! ending point, radius 1 ! ityp [type 1 line_to 2 arc_to+- 3 user 4 other side] Allaire fin (h_mesh.dat) Allaire fin convection 1.0 1. ! (CONDY (IM), CAPCY (IM), IM = 1, NMATR) 1 ! NOBCD [DIRICHLET BOUNDARY CONDITIONS] 111. ! bc 1 flag 1 ! sides that have a bc 1 1 ! boundary, bc number applied (- if first pt only) 0 ! NOBCN [FLUX, CONVECTION, RADIATION] (densit.dat) 4 1 1 0. 0. 5. 160. 0. 5. 110. 120. 10. 50. 120. 7.5 ========================== end =================================