&HEAD CHID='onetree_surf_1mesh' TITLE='Single processor One tree over surface fuel, uses fuel element model' / - Specify number of grid cells (IJK) and domain size (XB, meters) in x,y,z directons &MESH IJK=160,30,60, XB=-4,12,-3,0,0,6 / - Stretch computation grid in the z direction &TRNZ IDERIV=1,CC=0,PC=0.5 / &TRNZ IDERIV=2,CC=0,PC=0.0 / - Duration of simulation in seconds &TIME T_END=120 / - Miscellaneous parameters &MISC BAROCLINIC=.TRUE. / - Gaseous combustion parameters &REAC ID='WOOD' FYI='Ritchie, et al., 5th IAFSS, C_3.4 H_6.2 O_2.5, dHc = 15MW/kg' SOOT_YIELD = 0.02 O = 2.5 C = 3.4 H = 6.2 HEAT_OF_COMBUSTION = 17700. / - Account for water vapor released during drying &SPEC ID='WATER VAPOR' / -Tree crown vegetation based on NIST Douglas fir tree burn experiments. The bulk density is higher than normally found in most trees in the field &PART ID='TREE VEG',TREE=.TRUE.,QUANTITIES='VEG_TEMPERATURE', VEG_INITIAL_TEMPERATURE=20., VEG_SV=3900.,VEG_MOISTURE=1.00,VEG_CHAR_FRACTION=0.25, VEG_DRAG_COEFFICIENT=0.375,VEG_DENSITY=514.,VEG_BULK_DENSITY=2.76, VEG_BURNING_RATE_MAX=0.4,VEG_DEHYDRATION_RATE_MAX=0.4, VEG_REMOVE_CHARRED=.TRUE. / &TREE PART_ID='TREE VEG',XYZ=4,0,0,FUEL_GEOM='CONE', CROWN_WIDTH=3,CROWN_BASE_HEIGHT=1,TREE_HEIGHT=6, OUTPUT_TREE=.TRUE.,LABEL='tree_veg' / -Surface vegetation. Pinus Ponderosa pine needles. Uses Catchepole CST, 131, pp. 1-37, 1998 values of case PPMC59 &PART ID='GROUND VEG',TREE=.TRUE.,QUANTITIES='VEG_TEMPERATURE', VEG_INITIAL_TEMPERATURE=20., VEG_SV=5710.,VEG_MOISTURE=0.07,VEG_CHAR_FRACTION=0.25, VEG_DRAG_COEFFICIENT=0.375,VEG_DENSITY=510.,VEG_BULK_DENSITY=10, VEG_BURNING_RATE_MAX=0.4,VEG_DEHYDRATION_RATE_MAX=0.4, VEG_REMOVE_CHARRED=.TRUE. / &TREE PART_ID='GROUND VEG',FUEL_GEOM='RECTANGLE', XB = 0,8,-2,0,0,0.051,OUTPUT_TREE=.TRUE.,LABEL='ground-veg' / - Ignitor fire location and heat release rate magnitude and history 260 kW/m^2 fire on bottom of domain and at 0 <= x <= 0.2 meters; -2 <= y <= 0 meters &SURF ID='IGN FIRE',HRRPUA=260,RAMP_Q='RAMPFIRE' / &RAMP ID='RAMPFIRE',T= 0.0,F=0.0 / &RAMP ID='RAMPFIRE',T= 1.0,F=0.0 / &RAMP ID='RAMPFIRE',T= 1.5,F=1.0 / &RAMP ID='RAMPFIRE',T=12.5,F=1.0 / &RAMP ID='RAMPFIRE',T=13.0,F=0.0 / &VENT XB=0,0.2,-2,0,0,0,SURF_ID='IGN FIRE' / - Set inflow velocity characteristics, 2 m/s wind into domain &SURF ID='INFLOW',VEL=-2,RAMP_V='RAMPVEL' / &RAMP ID='RAMPVEL',T=0.0,F=0.0 / &RAMP ID='RAMPVEL',T=0.5,F=1.0 / - Boundary conditions on computational domain - Domain-Boundary conditions &VENT XB=-4,-4,-3, 0,0,6,SURF_ID='INFLOW' / inflow through x=-4 (left side) &VENT XB=12,12,-3, 0,0,6,SURF_ID='OPEN' / &VENT XB=-4,12,-3,-3,0,6,SURF_ID='OPEN' / &VENT XB=-4,12, 0, 0,0,6,SURF_ID='MIRROR' / symmetry along y=0 plane &VENT XB=-4,12,-3, 0,6,6,SURF_ID='OPEN' / - Output data to be viewed using Smokeview -- time intervals, in seconds, at which various quantities are written &DUMP DT_SLCF=0.1,DT_PART=0.1,DT_BNDF=0.1,DT_PL3D=20. / -- two-dimensional slice files &SLCF PBY= 0, QUANTITY='TEMPERATURE',VECTOR=.TRUE. / &SLCF PBY= 0, QUANTITY='soot density' / &SLCF PBY= 0, QUANTITY='WATER VAPOR' / H2O from veg drying -- boundary files &BNDF QUANTITY='NET HEAT FLUX'/ &BNDF QUANTITY='RADIATIVE HEAT FLUX'/ &BNDF QUANTITY='CONVECTIVE HEAT FLUX'/ - Declare end of input file &TAIL /