@sales2.1input / do not alter or remove this line! /*========================================================================== /* /* SALES_2 INPUT FILE /* =================== /* /* * Improved input file for the 2-D hydrocode SALE /* /* * Comments may be freely entered in this file following /* the string '/*'. Each entry should be clearly explained /* in this format, greatly easing the task of data entry. /* /* * All data are read in free format, so entries only need /* to be separated by spaces or tabs. /* /*-------------------------------------------------------------------------- /* /* Problem Title Information /* ========================= /* /* Enter a string (up to 80 characters long) that /* describes this computation. Enclose the string in quotes: /* /* jname /*-------------------------------------------------------------------------- 'Simplified lightning simulation - horizontal, no g' /*-------------------------------------------------------------------------- /* /* Start information /* =================== /* /* The following two integers provide information on /* whether the problem is: /* - a completely new problem, /* - a continuation of a previous problem from a /* dump file, /* - or a new problem with initial conditions /* defined in a dump file /* /* ibegcont has values 0, 1, or 2 for the three options /* described above: /* ibegcont = 0, for new problem /* ibegcont = 1, for a continuation /* ibegcont = 2, for restart /* /* If ibegcont is either 1, or 2, ndump is the number of /* the dump file from which the problem data will be /* restored. /* /* The name of the dumpfile must be of the form: sale_dump##.dat /* where ## is "ndump" (a two digit number, e.g. 01 for ndump=1) /* /* ibegcont ndump /*-------------------------------------------------------------------------- 0 0 /*-------------------------------------------------------------------------- /* /* IDL PLOTTING INPUT /* ==================== /* /* The following input values concern the IDL output plots /* which are generated by creating output data files of /* particular quantities for a selection af cells (or vertices), /* at regular time intervals (dtidlp). /* /* These output files may then read by an IDL routine. /* /* ------------------------------------- /* TO DISABLE SET "iidlp" TO ZERO BELOW /* ------------------------------------- /* /* Frequency information - iidlp, dtidlp /* ------------------------------------- /* Please switch the iidlp flag on to generate the data files. /* Please enter the time interval between each data point. /* That is, the run time interval between writing to the /* data files. /* /* Also please enter the number of cells to be examined, ncells /* and the number of plot types to be generated per cell, nplots /* /* ------------------------------------------------------------- /* NOTE: - It is important that for long run times "dtidlp" /* is large enough such that the output files do not /* get too large. /* - A useful guide to the value of dtidlp is: /* dtidlp > twfin/1000. /* - ncells must be less than or equal to 10 /* - nplots must be less than or equal to 6 /* ------------------------------------------------------------- /* /* iidlp dtidlp ncells nplots /*-------------------------------------------------------------------------- 0 0.1 5 3 /*-------------------------------------------------------------------------- /* /* Cell information /* ------------------ /* Please enter the value of i and j corresponding to each of /* the cells/vertices for which the output files will be /* generated. For cells, the i and j required corresponds to /* the vertex coordinates for the bottom left corner of the cell. /* /* ------------------------------------------------------------- /* NOTE: - There *must* be "ncells" sets of coordinates /* - Indexing begins in the bottom left corner of the mesh. /* Therefore, the first cell would have coordinates /* (i_cell,j_cell) = (1,1) /* - For vertex variables enter the vertex coordinates /* (i.e. bottom left coordinate of cell) /* ------------------------------------------------------------- /* /* i_cell j_cell /*-------------------------------------------------------------------------- 10 1 10 5 10 10 10 15 10 20 /*-------------------------------------------------------------------------- /* /* Plot information /* --------------------- /* Please enter a choice of "nplots" quantities to plot for each /* of the cells vs. time. Enter the integer corresponding to /* the plot of choice. /* /* The choices are as follows: /* /* 1 - Pressure /* 2 - Velocity in x-direction /* 3 - Velocity in y direction /* 4 - Principal stress, max /* 5 - Principal stress, int /* 6 - Principal stress, min /* 7 - Density /* 8 - Mass /* 9 - Volume /* 10 - x-coordinate /* 11 - y-coordinate /* 12 - Damage /* /* plot_type(n) - there must be "nplots" number of plot types /*-------------------------------------------------------------------------- 1 2 3 /*-------------------------------------------------------------------------- /* /* Enter data on time step control: /* ================================= /* dt0 = initial time step, recalculated during run /* dtmax = maximum allowed time step /* twfin = problem time to terminate the computation /* dtplot = time interval between plots /* /* Units are in seconds! /* /* inital step max step done! plot time /* dt0 dtmax twfin dtplot /*-------------------------------------------------------------------------- 2.00e-8 2.00e-8 1.0d-7 1.0d-8 /*-------------------------------------------------------------------------- /* /* Enter data on optional outputs: /* ================================== /* Long prints: Complete printout of all important cell /* and vertex quantities. (file = longprint.dat) /* Plus, printout of total energy, internal /* energy, etc. (file = global.dat) /* Tape dumps: Dump to unformatted file of all arrays and /* parameters required for restarting the /* problem from the dump. (file = sale_dump##.dat) /* /* ilprt = 1 forces long printout at print time /* = 0 suppresses long printouts /* dtlprt = problem time interval between long printouts /* idump = 1 forces tape dump at dump time /* = 0 suppresses tape dumps /* dtdump = problem time interval between tape dumps /* /* ------------------------------------------------------------ /* NOTE: - maximum no. of dumps is 24. /* Therefore, dtdump must be greater than /* dtdump > twfin/24 /* ------------------------------------------------------------ /* /* Units of time are in seconds! /* /* Long Printouts? print time Tape Dumps? dump time /* ilprt dtlprt idump dtdump /*------------------------------------------------------------------------- 0 5.0 0 10.0 /*------------------------------------------------------------------------- /* /* ============================================================= /* If continuation or restart the file is read to here only... /* ============================================================= /* /*-------------------------------------------------------------------------- /* /* Enter mesh properties /* ========================= /* The mesh is generated in subroutine celset, however, /* certain parameters can be defined here... /* /* /* Enter the acceleration of gravity: /* ---------------------------------- /* /* units of m/sec**2 /* /* **warning** make sure you get the sign right! conventional /* gravity has gx = 0.0, gy *negative* /* /* gx gy /*-------------------------------------------------------------------------- 0.0 0.0 /*-------------------------------------------------------------------------- /* /* Enter the number of cells in the mesh: /* -------------------------------------- /* nx = the number of cells in the x-direction /* ny = the number of cells in the y-direction /* nproj = the number of cells in the projectile stub /* (choose nproj=0 if mesh has no extension) /* /* ----------------------------------------------------------- /* Note: The value "node" in the common block determines /* the number of vertices allowed in the mesh. /* (nx+1)*(ny+1) should, therefore, NOT exceed this /* value. Current value of node = 4000 /* ----------------------------------------------------------- /* /* nx ny nproj /*-------------------------------------------------------------------------- 25 19 0 /*-------------------------------------------------------------------------- /* /* Enter data on mesh geometry: /* ------------------------------ /* cyl = 1.0 for cylindrical geometry, /* = 0.0 for plane strain geometry /* /* also, enter the flag sphgrd = 1 if you wish the rectangular /* mesh to be mapped to a cylindrical area. If cyl = 1.0 /* this will represent a spherical volume. /* (if used, be sure that ny=2*nx, else mapping is incorrect) /* /* cyl sphgrd /*------------------------------------------------------------------------- 0.0 0 /*------------------------------------------------------------------------- /* /* Enter the dimensions of each rectangular cell: /* ----------------------------------------------- /* All cells are assumed to be the same size before mapping. /* /* Units are in meters! /* /* x-direction y-direction /* dx dy /*------------------------------------------------------------------------- 1.00e-3 1.00e-3 /*------------------------------------------------------------------------- /* /* Enter flags to describe the boundary conditions: /* --------------------------------------------------- /* wb, wl, wr and wt correspond to the boundary conditions /* along the bottom, left, right and top edges of the mesh /* /* 0 = lagrangian surface /* 1 = simple free slip /* 2 = general free slip for curvilinear boundaries /* 3 = no slip /* 4 = continuative outflow /* 5 = specified inflow or outflow /* 6 = applied pressure /* /* bottom left right top /* wb wl wr wt /*------------------------------------------------------------------------- 0 1 1 0 /*------------------------------------------------------------------------- /* /* Enter computational information /* ================================== /* /* Enter flags for the computation type: /* ---------------------------------------- /* Define whether compressible or incompressible flow: /* /* imp = 1 for implicit pressure calculation, /* = 0 for purely explicit calculation /* inc = 1 for incompressible limit variant /* of imp=1 calculation /* irez = rezone flag: /* 0 = eulerian /* 1 = lagrangian /* 2 or greater for a user specified /* continuous rezone /* /* ----------------------------------------------------------- /* Note: rezone operations *never* affect boundaries /* between different materials! /* ----------------------------------------------------------- /* /* imp inc irez /*------------------------------------------------------------------------- 0 0 1 /*------------------------------------------------------------------------- /* /* The following coefficients control the iteration /* process for incompressible computations *only* /* (see SALE manual for further explanation): /* /* om = relaxation coefficient used in pressure iteration /* peps = pressure fraction scaling the relaxation factor rdsdp /* eps = allowed relative error in the pressure iteration /* rf = relaxation factor for continuous grid rezoning /* /* om peps eps rf /*------------------------------------------------------------------------- 1.0 1.0e-4 1.0e-4 0.05 /*------------------------------------------------------------------------- /* /* Enter data on mesh stability control: /* -------------------------------------- /* artvis = artificial (bulk) viscosity coefficient /* anc = alternate node-coupler coefficient /* xi = 1. for 4th-order node-coupling /* = 0. for 2nd-order node-coupling /* /* --------------------------------------------------------- /* Note: anc is reset to zero when cells are elastic /* --------------------------------------------------------- /* /* artvis anc xi /*------------------------------------------------------------------------- 0.1 0 .3 /*------------------------------------------------------------------------- /* /* Enter factors controlling convective fluxing /* ---------------------------------------------- /* (used for Eulerian mode computations only!) /* /* a0, b0 limiting cases- /* 0 0 centered; unstable /* 1 0 full donor cell; stable; diffusive /* 0 1 interpolated donor cell; linearly stable; non-diffusive /* 1 1 stable; diffusive /* /* a0 b0 /*------------------------------------------------------------------------- 1.0 0.0 /*------------------------------------------------------------------------- /* /* Enter the Material information /* ==================================== /* /* Enter the number of different materials employed in /* this computation. /* --------------------------------------------------- /* Actual assignment of material number to different /* portions of the mesh is performed by subroutine celset /* in the code. Edit this routine for custom mesh assignments. /* /* ------------------------------------------------------------ /* Note: Currently set up to consider up to 4 layers of equal /* thickness. /* /* ------------------------ /* LAYER 4 /* ------------------------ /* LAYER 3 /* ------------------------ /* LAYER 2 /* ------------------------ /* LAYER 1 /* ------------------------ /* ------------------------------------------------------------ /* /* MAXIMUM NUMBER OF LAYERS = 4 /* /* numat /*------------------------------------------------------------------------- 3 /*------------------------------------------------------------------------- /* /* Enter a string in quotes describing each material /* employed in this computation (up to 32 characters) /* --------------------------------------------------- /* /* material description /* mat# xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx /*------------------------------------------------------------------------- 1 'Granite' 2 'high energy granite' 3 'Granite' /*------------------------------------------------------------------------- /* /* Enter the parameters describing the equation of state. /* in this case, the Tillotson (1962) equation is used, which /* has 10 input parameters. These express the pressure P as a /* function of the density r and specific internal energy E /* (note that the energy in the reference state is assumed to /* be zero! Negative internal energies are thus possible) /* /* P = A*e + B*e**2 + [a+b/(1.0+E/(E0*h**2))]*r*E /* /* where h=r/r0 and e = h-1.0 is the volume strain. /* This expression is used for compressed states, h>1.0 /* /* For expanded states, h<1.0 and E>Ecv the pressure /* is given by /* /* P = a*r*E + [b*r*E/(1.0+E/(E0*h**2))+A*e*exp(-beta*g)]* /* exp(-alpha*g**2) /* /* where g = r0/r-1.0 /* /* See Appendix II in Melosh, Impact Cratering (1989) for /* more details. /* /* The input parameters are described as: /* /* mat# = material number to which this eos belongs /* r0 = density of the reference state, kg/m**3 /* a = dimensionless parameter /* b = dimensionless paramter: (a + b) = Gruneisen gamma /* A = linear bulk modulus, Pa /* B = coefficient of nonlinear elastic term, Pa /* ezero = parameter categorizing energy dependence, J/kg /* alfa = constant controlling convergence to perfect gas eos /* beta = constant controlling convergence to perfect gas eos /* eiv = specific energy of incipient vaporization, J/kg /* ecv = specific energy of complete vaporizaion, J/kg /* /* To turn this into a stiffened gas equation of state, set /* /* a = gruneisen gamma = ratio of specific heats - 1.0 /* A = bulk modulus, Pa /* b = B = alfa = beta = ezero = 0.0 /* eiv = a very big number /* ecv = a very big number not equal to eiv /* /* To make this a linear elastic solid, set a = 0.0 /* /* To make this a perfect gas, set A = 0.0 /* /* /* mat# r0 a b A B ezero alfa beta eiv ecv /*------------------------------------------------------------------------- 1 2680 0.5 1.3 4e10 1.8e10 1.6e7 5.0 5.0 3.5e6 1.8e7 2 2680 0.5 1.3 4e10 1.8e10 1.6e7 5.0 5.0 3.5e6 1.8e7 3 2680 0.5 1.3 4e10 1.8e10 1.6e7 5.0 5.0 3.5e6 1.8e7 /*------------------------------------------------------------------------- /* /* Enter parameters describing the deviatoric stresses: /* ---------------------------------------------------- /* mat# = material number to which these parameters belong /* visc = linear shear viscosity of material, Pa-sec /* lambda = linear bulk viscosity of material, Pa-s /* mu = elastic shear modulus, Pa /* y0 = hugoniot elastic limit at zero temperature, Pa /* siem = specific internal energy at melting, J/kg /* (the hel is degraded according to the law /* y= y0*(1.-sie/siem)**2) /* /* --------------------------------------------------------- /* Note: - To disable the use of elastic/plastic deviatoric /* stresses, set mu = 0.0 /* --------------------------------------------------------- /* /* mat# visc lambda mu y0 siem csound /*-------------------------------------------------------------------------- /* Material: Westerly granite 1 5.e8 1.e7 2.4e10 0.0 3.3e6 5.183d3 2 5.e8 1.e7 2.4e10 0.0 3.3e6 5.183d3 3 5.e8 1.e7 2.4e10 0.0 3.3e6 5.183d3 /*-------------------------------------------------------------------------- /* /* Enter parameters describing the Grady-Kipp fragmentation model: /* --------------------------------------------------------------- /* mat# = material number to which these parameters belong /* cg = crack velocity, m/sec /* pweib = weibull fracture exponent, dimensionless /* cweib = weibull coefficient k, flaws/m**3 /* /* --------------------------------------------------------- /* Note: - To disable the use of the fragmentation model set /* cweib = 0.0. /* - If mu = 0.0 above, this model will not be applied /* --------------------------------------------------------- /* /* mat# cg pweib cweib /*-------------------------------------------------------------------------- /* Material: Westerly granite (Grady and Lipkin 1980) cweib=4.15e23 1 2.4e3 6.2 0.0 2 2.4e3 6.2 0.0 3 2.4e3 6.2 0.0 /*-------------------------------------------------------------------------- /* /* Enter the initial density and internal energy of each material /* ------------------------------------------------------------- /* /* density, kg/m**3 internal energy, J/kg /* mat# roi siei /*-------------------------------------------------------------------------- 1 2.68e3 0.0 2 2.68e3 3.7e12 3 2.68e3 0.0 /*-------------------------------------------------------------------------- /* /* Enter the specified inflow conditions /* ======================================== /* /* Enter the density and type of any material flowing in /* along a boundary (note that these numbers are ignored /* unless the boundary has been set to type 5, specified /* inflow, above). Note that the material used is that /* already adjacent to the boundary. In units of kg/m**3 /* /* roin /* bottom left right top /*------------------------------------------------------------------------- 0.0 0.0 0.0 0.0 /*------------------------------------------------------------------------- /* /* Enter the internal energy and type of any material /* flowing in along a boundary (note that these numbers are /* ignored unless the boundary has been set to type 5, specified /* inflow, above). Note that material used is that /* already adjacent to the boundary. In units of J/kg /* /* siein /* bottom left right top /*------------------------------------------------------------------------- 0.0 0.0 0.0 0.0 /*------------------------------------------------------------------------- /* /* Enter the pressure applied to a boundary /* (note that these numbers are ignored unless the boundary /* has been set to type 6, applied pressure, above.) /* In units of Pa. /* /* pap /* bottom left right top /*------------------------------------------------------------------------- 0.0 0.0 0.0 0.0 /*------------------------------------------------------------------------- /* /* Enter the x-velocity uin specified at a boundary /* (note that these numbers are ignored unless the boundary /* has been set to type 5, specified outflow, above.) /* In units of m/sec. /* /* Note that if nproj .ne. 0 above, then velocities specified /* for the bottom boundary are applied to the projectile *only*: /* adjacent cells are given zero velocity. Interfacial nodes are /* assigned 1/2 the specified velocity. /* /* uin /* bottom left right top /*------------------------------------------------------------------------- 0.0 0.0 0.0 0.0 /*------------------------------------------------------------------------- /* /* Enter the y-velocity vin specified at a boundary /* (note that these numbers are ignored unless the boundary /* has been set to type 5, specified outflow, above.) /* In units of m/sec. /* /* Note that if nproj .ne. 0 above, then velocities specified /* for the bottom boundary are applied to the projectile *only*: /* adjacent cells are given zero velocity. Interfacial nodes are /* assigned 1/2 the specified velocity. /* /* vin /* bottom left right top /*------------------------------------------------------------------------ 0.0 0.0 0.0 0.0 /*=======================================================================