如何实现/测出一个材料的温度阶跃？（热点火？）

Question

各位大佬：

请教一个问题，我想模拟一个材料的热分解及点火，即达到一定温度后，材料的温度会突然升高直到燃尽。

但是我不知如何进行设置，燃烧热设置得过大的话FDS又会报错：数值不稳定。

另外像这样的温度测点，我使用Inside wall temperature测温，数据就会显示在室温而测不出来任何值。

请各位大佬赐教。

Answer 1

在材料的热解反应里，设置材料热解后的产物为燃料 ，这样材料在热解之前不会燃烧，热解后发生燃烧温度自然会上升

Comments

您好，我按照步骤设置了，但是还是测不到温度的陡升，可否请您看一看呢？代码如下：

&HEAD CHID='28'/

&TIME T_END=800.0/

&DUMP DT_RESTART=300.0, DT_SL3D=0.25/

&MESH ID='Mesh01', IJK=40,40,40, XB=-1.0,1.0,-1.0,1.0,0.0,2.0/

&SPEC ID='PROPANE'/

&REAC ID='FDS6 METHANE',

      FYI='FDS6 Predefined',

      FUEL='METHANE'/

&DEVC ID='SOLID101', QUANTITY='ADIABATIC SURFACE TEMPERATURE', XYZ=0.0,0.0,0.7, IOR=3/

&DEVC ID='SOLID10', DEPTH=0.4, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,0.3, IOR=-3/

&DEVC ID='SOLID11', QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.0,0.0,0.3, IOR=-3/

&DEVC ID='SOLID2', DEPTH=0.399, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=-0.199,-0.199,0.3, IOR=-3/

&DEVC ID='SOLID6', DEPTH=1.0E-3, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=-0.199,-0.199,0.3, IOR=-3/

&DEVC ID='SOLID3', DEPTH=0.399, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.199,-0.199,0.3, IOR=-3/

&DEVC ID='SOLID7', DEPTH=1.0E-3, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.199,-0.199,0.3, IOR=-3/

&DEVC ID='SOLID8', DEPTH=1.0E-3, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.199,0.199,0.3, IOR=-3/

&DEVC ID='SOLID4', DEPTH=0.399, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=0.199,0.199,0.3, IOR=-3/

&DEVC ID='SOLID9', DEPTH=1.0E-3, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=-0.199,0.199,0.3, IOR=-3/

&DEVC ID='SOLID5', DEPTH=0.399, QUANTITY='INSIDE WALL TEMPERATURE', XYZ=-0.199,0.199,0.3, IOR=-3/

&DEVC ID='SOLID111', QUANTITY='ADIABATIC SURFACE TEMPERATURE', XYZ=0.0,0.0,0.3, IOR=-3/

&MATL ID='1',

      SPECIFIC_HEAT=1.176,

      CONDUCTIVITY=0.4644,

      DENSITY=1845.0,

      HEAT_OF_COMBUSTION=4.78E+12,

      N_REACTIONS=1,

      SPEC_ID(1,1)='PROPANE',

      NU_SPEC(1,1)=1.0,

      A=4780.0,

      E=1.439E+5/

&SURF ID='1',

      RGB=146,130,148,

      TEXTURE_MAP='psm_carpet_grey.jpg',

      TEXTURE_WIDTH=0.6096,

      TEXTURE_HEIGHT=0.6096,

      BURN_AWAY=.TRUE.,

      MATL_ID(1,1)='1',

      MATL_MASS_FRACTION(1,1)=1.0,

      THICKNESS(1)=0.4/

&SURF ID='burner',

      COLOR='RED',

      HRRPUA=1000.0,

      TMP_FRONT=300.0/

&OBST ID='Obstruction', XB=-0.2,0.2,-0.2,0.2,0.3,0.7, SURF_ID='1'/ 

&VENT ID='Mesh Vent: Mesh01 [XMAX]', SURF_ID='OPEN', XB=1.0,1.0,-1.0,1.0,0.0,2.0/ 

&VENT ID='Mesh Vent: Mesh01 [XMIN]', SURF_ID='OPEN', XB=-1.0,-1.0,-1.0,1.0,0.0,2.0/ 

&VENT ID='Mesh Vent: Mesh01 [YMAX]', SURF_ID='OPEN', XB=-1.0,1.0,1.0,1.0,0.0,2.0/ 

&VENT ID='Mesh Vent: Mesh01 [YMIN]', SURF_ID='OPEN', XB=-1.0,1.0,-1.0,-1.0,0.0,2.0/ 

&VENT ID='Mesh Vent: Mesh01 [ZMAX]', SURF_ID='OPEN', XB=-1.0,1.0,-1.0,1.0,2.0,2.0/ 

&VENT ID='Vent01', SURF_ID='burner', XB=-0.5,0.5,-0.5,0.5,1.0E-3,1.0E-3/ 

&BNDF QUANTITY='ADIABATIC SURFACE TEMPERATURE'/

&BNDF QUANTITY='BACK WALL TEMPERATURE'/

&BNDF QUANTITY='GAS TEMPERATURE'/

&BNDF QUANTITY='NET HEAT FLUX'/

&BNDF QUANTITY='WALL TEMPERATURE'/

&SLCF QUANTITY='TEMPERATURE', ID='1', PBX=0.0/

&SLCF QUANTITY='TEMPERATURE', ID='2', PBY=0.0/

&TAIL /