별지제 호서식 연구결과보고서 과제명 소속소방산업기술연구소연구책임자권성필 연구기간 연구목표 연구배경

별지제 호서식 연구결과보고서 과제명 소속소방산업기술연구소연구책임자권성필 연구기간 연구목표 연구배경

연구내용

β

β

연구결과 참고문헌

붙임 ρ ρ α ρ α ρ α ρ α α

α α ρ α α α ρ α α α

β α αβ α α αβ αβ β ρ exp max αβ αβ α α αβ α exp

κ π max P r

Δ

ρ δ max P r

δ ρ δ

붙임 #include <math.h> #include <stdlib.h> #include <stdio.h> void mesh(file *fp) { int i, j, k; float x_ini, x_end, y_ini, y_end, z_ini, z_end; float dx=0.25, dy=0.25, dz=0.50; for(i=0; i<4; i++){ z_ini=i*dz; z_end=(i+1)*dz; for(j=0; j<4; j++){ y_ini=j*dy; y_end=(j+1)*dy; } } for(k=0; k<4; k++){ x_ini=k*dx; x_end=(k+1)*dx; fprintf(fp, "&MESH IJK = 25, 25, 50,"); fprintf(fp, "XB=%3.2f, %3.2f, %3.2f, %3.2f, %3.2f, %3.2f,", x_ini, x_end, y_ini, y_end, z_ini, z_end); fprintf(fp, " / \n"); } fprintf(fp, " \n"); } void matl(file *fp) { fprintf(fp, "&MATL ID = 'CONCRETE' \n"); fprintf(fp, " DENSITY = 2200. \n"); fprintf(fp, " CONDUCTIVITY = 1.2 \n"); fprintf(fp, " SPECIFIC_HEAT = 0.88 /&MATL ID = 'CONCRETE' \n"); fprintf(fp, " DENSITY = 2200. \n"); fprintf(fp, " CONDUCTIVITY = 1.2 \n"); fprintf(fp, " SPECIFIC_HEAT = 0.88 / \n\n"); fprintf(fp, "&MISC SURF_DEFAULT fprintf(fp, "&SURF ID = 'WALL' / \n\n"); = 'BURNER', HRRPUA=1000., PART_ID='smoke' / \n\n"); fprintf(fp, "&REAC ID = 'ETHANOL VAPOR' \n"); fprintf(fp, " FYI = 'Ethanol, C_2 H_6 O' \n"); fprintf(fp, " HEAT_OF_COMBUSTION = 26780. \n"); fprintf(fp, " O = 1. \n"); fprintf(fp, " C = 2. \n");

fprintf(fp, " H = 6. \n"); fprintf(fp, " SOOT_YIELD = 0.008 / \n\n"); fprintf(fp, "&SURF ID = 'BURNER', HRRPUA=1000. / \n\n"); fprintf(fp, "&MATL ID = 'CELLULOSE' \n"); fprintf(fp, " FYI = 'Properties completely fabricated' \n"); fprintf(fp, " SPECIFIC_HEAT = 2.3 \n"); fprintf(fp, " CONDUCTIVITY = 0.29 \n"); fprintf(fp, " DENSITY = 40.0 \n"); fprintf(fp, " N_REACTIONS = 1 \n"); fprintf(fp, " NU_FUEL = 1. \n"); fprintf(fp, " REFERENCE_TEMPERATURE = 500. \n"); fprintf(fp, " HEAT_OF_REACTION = 418. \n"); fprintf(fp, " HEAT_OF_COMBUSTION = 8000. / \n\n"); fprintf(fp, "&MATL ID = 'LIGNIN' \n"); fprintf(fp, " FYI = 'Properties completely fabricated' \n"); fprintf(fp, " SPECIFIC_HEAT = 1.1 \n"); fprintf(fp, " CONDUCTIVITY = 0.1 \n"); fprintf(fp, " DENSITY = 550. \n"); fprintf(fp, " N_REACTIONS = 1 \n"); fprintf(fp, " NU_FUEL = 1. \n"); fprintf(fp, " REFERENCE_TEMPERATURE = 500. \n"); fprintf(fp, " HEAT_OF_REACTION = 418. \n"); fprintf(fp, " HEAT_OF_COMBUSTION = 8000. / \n\n"); fprintf(fp, "&MATL ID = 'STEEL' \n"); fprintf(fp, " EMISSIVITY = 1.0 \n"); fprintf(fp, " DENSITY = 7850. \n"); fprintf(fp, " CONDUCTIVITY = 45.8 \n"); fprintf(fp, " SPECIFIC_HEAT = 0.46 / \n\n"); fprintf(fp, "&MATL ID = 'ETHANOL LIQUID' \n"); fprintf(fp, " EMISSIVITY = 1.0 \n"); fprintf(fp, " NU_FUEL = 0.97 \n"); fprintf(fp, " HEAT_OF_REACTION= 880. \n"); fprintf(fp, " CONDUCTIVITY = 0.17 \n"); fprintf(fp, " SPECIFIC_HEAT = 2.45 \n"); fprintf(fp, " DENSITY = 787. \n"); fprintf(fp, " ABSORPTION_COEFFICIENT = 40. \n"); fprintf(fp, " BOILING_TEMPERATURE = 76. / \n\n"); fprintf(fp, "&MATL ID = 'WATER' \n"); fprintf(fp, " EMISSIVITY = 1.0 \n"); fprintf(fp, " DENSITY = 1000. \n"); fprintf(fp, " CONDUCTIVITY = 0.6 \n"); fprintf(fp, " SPECIFIC_HEAT = 4.19 \n"); fprintf(fp, " N_REACTIONS = 1 \n"); fprintf(fp, " A = 1E20 \n"); fprintf(fp, " E = 1.62E+05 \n");

fprintf(fp, " NU_WATER = 1.0 \n"); fprintf(fp, " HEAT_OF_REACTION = 2260. / \n\n"); fprintf(fp, "&SURF ID = 'WOOD_CELL' \n"); fprintf(fp, " COLOR = 'BURNT UMBER' \n"); fprintf(fp, " BURN_AWAY =.TRUE. \n"); fprintf(fp, " MATL_ID(1,1:3) = 'CELLULOSE','WATER','LIGNIN' \n"); fprintf(fp, " MATL_MASS_FRACTION(1,1:3) = 0.70,0.10,0.20 \n"); fprintf(fp, " MATL_ID(2,1) = 'LIGNIN' \n"); fprintf(fp, " THICKNESS(1:2) = 0.1,0.002 \n"); fprintf(fp, " PART_ID = 'smoke' / \n\n"); fprintf(fp, "&SURF ID = 'WALL' \n"); fprintf(fp, " RGB = 200,200,200 \n"); fprintf(fp, " COLOR = 'GRAY' \n"); fprintf(fp, " MATL_ID = 'CONCRETE' \n"); fprintf(fp, " THICKNESS = 0.012 / \n\n"); fprintf(fp, "&SURF ID = 'STEEL SHEET' \n"); fprintf(fp, " COLOR = 'SILVER' \n"); fprintf(fp, " MATL_ID = 'STEEL' \n"); fprintf(fp, " THICKNESS = 0.01 \n"); fprintf(fp, " BACKING = 'EXPOSED'/ \n\n"); fprintf(fp, "&SURF ID = 'ETHANOL POOL' \n"); fprintf(fp, " FYI = '4 kg of ethanol in a 0.7 m x 0.8 m pan' \n"); fprintf(fp, " COLOR = 'STEEL BLUE' \n"); fprintf(fp, " MATL_ID = 'ETHANOL LIQUID' \n"); fprintf(fp, " THICKNESS = 0.01 / \n\n"); } fprintf(fp, "&PART ID = 'smoke', MASSLESS =.TRUE., SAMPLING_FACTOR = 1 / \n\n"); void obst(file *fp) { int i, j, k; float x_ini, x_end, y_ini, y_end, z_ini, z_end; fprintf(fp, "&OBST XB = 0.20, 0.80, 0.20, 0.23, 0.17, 0.20, SURF_ID = 'STEEL SHEET' / \n"); fprintf(fp, "&OBST XB = 0.20, 0.80, 0.77, 0.80, 0.17, 0.20, SURF_ID = 'STEEL SHEET' / \n"); fprintf(fp, "&OBST XB = 0.20, 0.23, 0.20, 0.80, 0.17, 0.20, SURF_ID = 'STEEL SHEET' / \n"); fprintf(fp, "&OBST XB = 0.77, 0.80, 0.20, 0.80, 0.17, 0.20, SURF_ID = 'STEEL SHEET' / \n\n"); fprintf(fp, "&OBST XB = 0.77, 0.80, 0.20, 0.23, 0.00, 0.17, SURF_ID = 'STEEL SHEET' / \n"); fprintf(fp, "&OBST XB = 0.77, 0.80, 0.77, 0.80, 0.00, 0.17, SURF_ID = 'STEEL SHEET' / \n"); fprintf(fp, "&OBST XB = 0.20, 0.23, 0.20, 0.23, 0.00, 0.17, SURF_ID = 'STEEL SHEET' / \n"); fprintf(fp, "&OBST XB = 0.20, 0.23, 0.77, 0.80, 0.00, 0.17, SURF_ID = 'STEEL SHEET' / \n\n"); fprintf(fp, "&OBST XB = 0.20, 0.80, 0.20, 0.80, 0.05, 0.08,"); fprintf(fp, " SURF_IDS = 'ETHANOL POOL', 'STEEL SHEET', 'STEEL SHEET'/ \n"); fprintf(fp, "&OBST XB = 0.20, 0.80, 0.20, 0.80, 0.07, 0.07, SURF_ID = 'BURNER' / \n\n");

z_ini=0.17; for(k=0; k<26; k++){ z_ini+=0.03; z_end=z_ini+0.03; if(k%4==0){ x_ini=0.09; for(i=0; i<6; i++){ if(i%2) x_ini+=0.12; else x_ini+=0.11; x_end=x_ini+0.02; fprintf(fp, "&OBST XB = %3.2f, %3.2f, 0.20, 0.80, %3.2f, %3.2f,", x_ini, x_end, z_ini, z_end); fprintf(fp, " SURF_ID = 'WOOD_CELL' / \n"); } fprintf(fp, " \n"); } if(k%4==1){ y_ini=0.09; for(j=0; j<6; j++){ if(j%2) y_ini+=0.12; else y_ini+=0.11; y_end=y_ini+0.02; fprintf(fp, "&OBST XB = 0.20, 0.80, %3.2f, %3.2f, %3.2f, %3.2f,", y_ini, y_end, z_ini, z_end); fprintf(fp, " SURF_ID = 'WOOD_CELL' / \n"); } fprintf(fp, " \n"); } if(k%4==2){ x_ini=0.15; for(i=0; i<5; i++){ if(i%2) x_ini+=0.12; else x_ini+=0.11; x_end=x_ini+0.02; fprintf(fp, "&OBST XB = %3.2f, %3.2f, 0.20, 0.80, %3.2f, %3.2f,", x_ini, x_end, z_ini, z_end); fprintf(fp, " SURF_ID = 'WOOD_CELL' / \n"); } fprintf(fp, " \n"); } if(k%4==3){ y_ini=0.15; for(j=0; j<5; j++){ if(j%2) y_ini+=0.12; else y_ini+=0.11; y_end=y_ini+0.02; fprintf(fp, "&OBST XB = 0.20, 0.80, %3.2f, %3.2f, %3.2f, %3.2f,",

} } } fprintf(fp, " \n"); } y_ini, y_end, z_ini, z_end); fprintf(fp, " SURF_ID = 'WOOD_CELL' / \n"); void output(file *fp) { fprintf(fp, "Open up the domain on all sides except the bottom (ZMIN) \n\n"); fprintf(fp, "&VENT MB = 'XMIN', SURF_ID = 'OPEN' / \n"); fprintf(fp, "&VENT MB = 'XMAX', SURF_ID = 'OPEN' / \n"); fprintf(fp, "&VENT MB = 'YMIN', SURF_ID = 'OPEN' / \n"); fprintf(fp, "&VENT MB = 'YMAX', SURF_ID = 'OPEN' / \n"); fprintf(fp, "&VENT MB = 'ZMAX', SURF_ID = 'OPEN' / \n\n"); } fprintf(fp, "&SLCF PBX = 0.260, QUANTITY = 'TEMPERATURE',VECTOR=.TRUE. / \n"); fprintf(fp, "&SLCF PBX = 0.260, QUANTITY = 'HRRPUV' / \n"); fprintf(fp, "&SLCF PBX = 0.260, QUANTITY = 'MIXTURE_FRACTION' / \n\n"); main() { FILE *fp; fp = fopen("model_04.fds", "w+"); fprintf(fp, "&HEAD CHID = 'model_04', TITLE = 'Test Case for 64 CPUs, "); fprintf(fp, "SVN $Revision: 935 $' / \n\n"); fprintf(fp, "&TIME TWFIN = 60 / \n\n"); mesh(fp); matl(fp); obst(fp); output(fp); fprintf(fp, "&TAIL \n\n"); } fclose(fp);