老师,建了一个树木模型,但树就是着不起来,向您请教一下啥原因呢!
&MESHID='Mesh01', IJK=50,40,25, XB=0.0,50.0,0.0,40.0,0.0,25.0/
&SPEC ID='WATERVAPOR', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID='CARBONMONOXIDE', LUMPED_COMPONENT_ONLY=.TRUE./
&SPECID='NITROGEN', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID='CARBONDIOXIDE', LUMPED_COMPONENT_ONLY=.TRUE./
&SPECID='OXYGEN', LUMPED_COMPONENT_ONLY=.TRUE./
&SPECID='mucai', FORMULA='CH1.7O0.74N0.002'/
&SPEC ID='SOOT',FORMULA='C', LUMPED_COMPONENT_ONLY=.TRUE./
&SPEC ID='AIR',BACKGROUND=.TRUE.,
SPEC_ID(1)='NITROGEN',
SPEC_ID(2)='OXYGEN',
VOLUME_FRACTION(1)=3.76,
VOLUME_FRACTION(2)=1.0/
&SPECID='PRODUCT',
SPEC_ID(1)='CARBON DIOXIDE',
SPEC_ID(2)='CARBON MONOXIDE',
SPEC_ID(3)='NITROGEN',
SPEC_ID(4)='SOOT',
SPEC_ID(5)='WATER VAPOR',
VOLUME_FRACTION(1)=0.964384,
VOLUME_FRACTION(2)=3.655E-3,
VOLUME_FRACTION(3)=3.84076,
VOLUME_FRACTION(4)=0.031961,
VOLUME_FRACTION(5)=0.85/
&REACID='MUCAI',
HEAT_OF_COMBUSTION=1.8E4,
FUEL='mucai',
SPEC_ID_NU='AIR','mucai','PRODUCT',
NU=-1.02121,-1.0,1.0/
&MATLID='shugan',
SPECIFIC_HEAT=3.0,
CONDUCTIVITY=0.28,
DENSITY=900.0,
HEAT_OF_COMBUSTION=8000.0,
N_REACTIONS=1,
HEAT_OF_REACTION=1500.0,
SPEC_ID(1,1)='PRODUCT',
NU_SPEC(1,1)=1.0,
A=5.6666E4,
E=75.6/
&MATLID='shuguan',
SPECIFIC_HEAT=3.2,
CONDUCTIVITY=0.05,
DENSITY=15.0,
HEAT_OF_COMBUSTION=1.3E4,
N_REACTIONS=1,
HEAT_OF_REACTION=1500.0,
SPEC_ID(1,1)='PRODUCT',
NU_SPEC(1,1)=1.0,
A=5.6666E4,
E=75.79/
&SURF ID='wind',
RGB=26,101,217,
VEL=-4.0,
PROFILE='PARABOLIC'/
&SURFID='shuagan',
RGB=24,34,19,
BURN_AWAY=.TRUE.,
BACKING='VOID',
MATL_ID(1,1)='shugan',
MATL_MASS_FRACTION(1,1)=1.0,
THICKNESS(1)=0.15/
&SURFID='shuguan',
RGB=146,237,103,
BACKING='VOID',
MATL_ID(1,1)='shuguan',
MATL_MASS_FRACTION(1,1)=1.0,
THICKNESS(1)=0.15/
&SURFID='burner',
COLOR='RED',
HRRPUA=4.6E4,
RAMP_Q='burner_RAMP_Q',
TMP_FRONT=300.0/
&RAMPID='burner_RAMP_Q', T=0.0, F=0.0/
&RAMPID='burner_RAMP_Q', T=5.0, F=1.0/
&RAMPID='burner_RAMP_Q', T=10.0, F=1.0/
&RAMPID='burner_RAMP_Q', T=15.0, F=1.0/
&RAMPID='burner_RAMP_Q', T=20.0, F=0.0/
&OBSTID='Obstruction', XB=0.5,12.5,2.5,3.0,0.0,2.0, SURF_ID='INERT'/
&OBST ID='Slab',XB=1.0,2.0,3.0,4.0,1.0,2.0, THICKEN=.TRUE., SURF_ID='shuguan'/
&OBST ID='Slab',XB=2.0,3.0,4.0,5.0,1.0,2.0, THICKEN=.TRUE., SURF_ID='shuguan'/
&OBST ID='Slab',XB=2.0,4.0,1.0,2.0,1.0,2.0, THICKEN=.TRUE., SURF_ID='shuguan'/
&OBST ID='Slab',XB=2.0,4.0,2.0,4.0,1.0,4.0, THICKEN=.TRUE., SURF_ID='shuguan'/
&OBST ID='Slab',XB=3.0,4.0,4.0,5.0,1.0,3.0, THICKEN=.TRUE., SURF_ID='shuguan'/
&OBST ID='Slab',XB=4.0,5.0,2.0,4.0,1.0,2.0, THICKEN=.TRUE., SURF_ID='shuguan'/
&VENT ID='MeshVent: Mesh01 [XMAX]', SURF_ID='OPEN', XB=50.0,50.0,0.0,40.0,0.0,25.0/
&VENT ID='MeshVent: Mesh01 [XMIN]', SURF_ID='OPEN', XB=0.0,0.0,0.0,40.0,0.0,25.0/
&VENT ID='MeshVent: Mesh01 [YMAX]', SURF_ID='OPEN', XB=0.0,50.0,40.0,40.0,0.0,25.0/
&VENT ID='MeshVent: Mesh01 [YMIN]', SURF_ID='wind', XB=0.0,50.0,0.0,0.0,0.0,25.0/
&VENT ID='MeshVent: Mesh01 [ZMAX]', SURF_ID='OPEN', XB=0.0,50.0,0.0,40.0,25.0,25.0/
&VENTID='burner', SURF_ID='burner', XB=2.5,12.5,2.0,3.0,2.0,2.0, SPREAD_RATE=2.0,XYZ=7.5,2.5,2.0/