Fan not passing extra species

[GoogleCodeExporter]

Application Version: 5.3.1
SVN Revision Number: 3729
Compile Date: Wed, 08 Apr 2009
Operating System: Windows Vista Ultimate SP2 32-bit

Describe details of the issue below:

We had a customer modeling a ventilation system with a contaminant gas. 
The volume fraction was not being calculated as expected.  I traced the
problem to a fan (porous SURF with specified velocity applied to a thin
OBST) that was between two meshes that made up the duct system.  It appears
that if the fan is at the boundary of a mesh, then the composition of extra
species is not passed to the other side.  It works as expected if the
obstruction is moved slightly to one side or the other.  I did not test
with combustion products, so I do not know if other parts of the gas phase
composition are being passed properly.

In the test case below, the concentration of extra species 1 is zero in
mesh_2 when the obstruction is at x=0.  If the obstruction is at moved
slightly to x=1.8, then the concentration in mesh_2 is 5% as expected.

====================================

&TIME T_END=60.00/
&DUMP DT_PL3D=5.00, PLOT3D_QUANTITY='VOLUME
FRACTION','U-VELOCITY','V-VELOCITY','W-VELOCITY','PRESSURE',
PLOT3D_SPEC_ID(1)='ExSPEC'/

&MESH ID='MESH_1', IJK=10,5,5, XB=0.00,2.00,0.00,1.00,0.00,1.00/
&MESH ID='MESH_2', IJK=10,5,5, XB=2.00,4.00,0.00,1.00,0.00,1.00/

&SPEC ID='ExSPEC'/

&SURF ID='supply',
      RGB=26,204,26,
      VEL=-1.00,
      MASS_FRACTION(1)=0.0500,
      TAU_MF(1)=1.00/
&SURF ID='fan',
      RGB=51,51,204,
      VEL=1.00,
      POROUS=.TRUE./

&DEVC ID='C1', QUANTITY='VOLUME FRACTION', SPEC_ID='ExSPEC',
XYZ=1.50,0.50,0.50/
&DEVC ID='C2', QUANTITY='VOLUME FRACTION', SPEC_ID='ExSPEC',
XYZ=2.50,0.50,0.50/
&OBST XB=2.00,2.00,0.00,1.00,0.00,1.00, SURF_ID='fan'/ Obstruction

&VENT SURF_ID='supply', XB=0.00,0.00,0.00,1.00,0.00,1.00/ Vent
&VENT SURF_ID='OPEN', XB=4.00,4.00,0.00,1.00,0.00,1.00/ Vent[1]

&SLCF QUANTITY='VELOCITY', VECTOR=.TRUE., PBY=0.3000/
&SLCF QUANTITY='VOLUME FRACTION', SPEC_ID='ExSPEC', PBY=0.3000/

Original issue reported on code.google.com by hardeman...@gmail.com on 20 Jul 2009 at 9:27

[GoogleCodeExporter]

Thanks, I'll take a look. In future, it's best to post the input file. Cutting 
and 
pasting can lead to end of line errors or confusion.

Original comment by mcgra...@gmail.com on 20 Jul 2009 at 9:59

• Changed state: Accepted

[GoogleCodeExporter]

Well, you know us Thunderhead folks -- I just copied it from the FDS view 
inside of
PyroSim.  I can save it to a text file next time.

Original comment by hardeman...@gmail.com on 20 Jul 2009 at 10:17

[GoogleCodeExporter]

Tricky one -- the logic for these porous obstructions and obstructions that 
straddle 
two meshes was mutually exclusive. I added some logic to allow both conditions 
to be 
applied.

Original comment by mcgra...@gmail.com on 24 Jul 2009 at 7:01

• Changed state: Fixed

[GoogleCodeExporter]

Assume verified.

Original comment by mcgra...@gmail.com on 23 Nov 2009 at 9:46

• Changed state: Verified