search

firemodels / fds

Fire Dynamics Simulator
https://pages.nist.gov/fds-smv/
Other
642 stars 616 forks source link

POROUS_VOLUME_FRACTION #10714

Closed Pethav-7 closed 1 year ago

Pethav-7 commented 2 years ago

Hi.

I just completed two simulations that were exactly the same, exact that the POROUS_VOLUME_FRACTION differd from 0% to a 100%, yet i didnt notice any change in the thermocouples above the material. When the graphs were ploted in excel.

How come the POROUS_VOLUME_FRACTION dosn't effect the movement of hot fire gases through the material?

The following input was used for the material.

&MATL ID='PB', SPECIFIC_HEAT=0.7, CONDUCTIVITY=0.075, DENSITY=450.0, EMISSIVITY=0.93/

&SURF ID='Pyrobubbles', RGB=255,102,0, MATL_ID(1,1)='PB', MATL_MASS_FRACTION(1,1)=1.0, THICKNESS(1)=5.0E-4, GEOMETRY='SPHERICAL', INNER_RADIUS=5.0E-3/

&PART ID='skum', DRAG_LAW='SPHERE', SURF_ID='Pyrobubbles', POROUS_VOLUME_FRACTION=100.0, STATIC=T, PERMEABILITY=1.0E-7,1.E-7,1.E-7 /

The container with the material in it is shown below. Screenshot_1

drjfloyd commented 2 years ago

You have set DRAG_LAW='SPHERE'. This does not invoke the porous media model; DRAG_LAW='POROUS MEDIA' does. The POROUS_VOLUME_FRACTION units are not percent, this should be a fraction from 0 to 1. See the discussion and example in Section 17.4.8.

Since you have not attached a complete input file, I cannot tell if you have properly defined the &INIT for the particles. That could also be incorrect.

Try fixing the issues I noted and report back if you are still having problems (in which case attaching a full simple input file would be helpful).

Pethav-7 commented 2 years ago

There was a problem with uploading the file, but here is the code:

&HEAD CHID='porous_0,35_por_ist_for_0,85'/ &TIME T_END=300.0/ &DUMP DT_RESTART=300.0, DT_SL3D=0.25, NFRAMES=60/ &MISC POROUS_FLOOR=.FALSE./

&MESH ID='Mesh01', IJK=50,50,27, XB=0.0,2.0,0.0,2.0,-0.059231,1.0/

&PART ID='PB_partikel', STATIC=.TRUE., SURF_ID='Pyrobubbles', COLOR='AZURE', AGE=600.0/

&REAC ID='POLYURETHANE', FYI='NFPA Babrauskas', FUEL='REAC_FUEL', C=6.3, H=7.1, O=2.1, N=1.0, SOOT_YIELD=0.1/

&DEVC ID='THCP', QUANTITY='THERMOCOUPLE', XYZ=1.055657,1.304148,0.71/ &DEVC ID='THCP01', QUANTITY='THERMOCOUPLE', XYZ=1.041773,0.69362,0.710088/ &DEVC ID='THCP02', QUANTITY='THERMOCOUPLE', XYZ=1.041773,1.003315,0.71/ &DEVC ID='THCP03', QUANTITY='THERMOCOUPLE', XYZ=1.041773,0.702215,0.694522/ &DEVC ID='THCP04', QUANTITY='THERMOCOUPLE', XYZ=1.041773,1.003315,0.709/ &DEVC ID='THCP05', QUANTITY='THERMOCOUPLE', XYZ=1.055657,1.304148,0.694522/ &DEVC ID='THCP06', QUANTITY='THERMOCOUPLE', XYZ=1.36,0.995988,0.419053, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP07', QUANTITY='THERMOCOUPLE', XYZ=0.671031,0.995988,0.419053, ORIENTATION=1.0,0.0,0.0/ &DEVC ID='THCP08', QUANTITY='THERMOCOUPLE', XYZ=1.04161,1.555658,0.419053, ORIENTATION=0.0,-1.0,0.0/ &DEVC ID='THCP09', QUANTITY='THERMOCOUPLE', XYZ=1.04161,0.41587,0.419053, ORIENTATION=0.0,1.0,-1.0/ &DEVC ID='THCP10', QUANTITY='THERMOCOUPLE', XYZ=1.04161,0.41587,0.626681, ORIENTATION=0.0,1.0,-1.0/ &DEVC ID='THCP11', QUANTITY='THERMOCOUPLE', XYZ=1.36,0.993984,0.626681, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP12', QUANTITY='THERMOCOUPLE', XYZ=0.671031,0.99704,0.626681, ORIENTATION=1.0,0.0,0.0/ &DEVC ID='THCP13', QUANTITY='THERMOCOUPLE', XYZ=1.041483,1.555658,0.626681, ORIENTATION=0.0,-1.0,0.0/ &DEVC ID='THCP14', QUANTITY='THERMOCOUPLE', XYZ=1.384204,0.995988,0.419053, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP15', QUANTITY='THERMOCOUPLE', XYZ=1.384204,0.993984,0.626681, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP16', QUANTITY='THERMOCOUPLE', XYZ=0.671031,0.995988,0.419053, ORIENTATION=1.0,0.0,0.0/

&MATL ID='STEEL', FYI='Drysdale, Intro to Fire Dynamics - ATF NIST Multi-Floor Validation', SPECIFIC_HEAT=0.46, CONDUCTIVITY=45.8, DENSITY=7850.0, EMISSIVITY=0.95/ &MATL ID='PB', SPECIFIC_HEAT=0.7, CONDUCTIVITY=0.075, DENSITY=450.0, EMISSIVITY=0.93/

&SURF ID='Brand', COLOR='RED', HRRPUA=250.0, RAMP_Q='Brand_RAMP_Q', TMP_FRONT=100.0/ &RAMP ID='Brand_RAMP_Q', T=1.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=5.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=6.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=17.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=21.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=22.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=35.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=40.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=45.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=56.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=60.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=61.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=74.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=79.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=84.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=95.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=96.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=97.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=100.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=126.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=137.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=141.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=142.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=155.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=200.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=300.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=315.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=316.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=317.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=320.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=346.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=347.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=351.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=352.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=365.0, F=0.37/ &SURF ID='Stål', RGB=255,52,141, BACKING='VOID', MATL_ID(1,1)='STEEL', MATL_MASS_FRACTION(1,1)=1.0, THICKNESS(1)=1.0E-3/ &SURF ID='Pyrobubbles', RGB=255,102,0, MATL_ID(1,1)='PB', MATL_MASS_FRACTION(1,1)=1.0, THICKNESS(1)=5.0E-4, GEOMETRY='SPHERICAL', INNER_RADIUS=5.0E-3/

&PART ID='skum', DRAG_LAW='POROUS MEDIA', SURF_ID='Pyrobubbles', POROUS_VOLUME_FRACTION=0.35, STATIC=T, DRAG_COEFFICIENT=0.1,0.1,0.1 PERMEABILITY=1.0E-7,1.E-7,1.E-7 /

&INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.661034,1.38413,0.41587,1.55587,0.320765,0.639705/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=1.222775,1.38413,0.41587,1.55587,1.110223E-16,0.320765/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=1.166992,1.222773,0.778304,1.228273,0.205,0.320765/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.860773,0.914425,0.778522,1.228561,0.205,0.320765/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.661034,0.860714,0.41587,1.55587,0.0,0.320765/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.860714,1.222586,1.22856,1.55587,0.0,0.320005/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.860714,1.222586,0.416388,0.777676,0.0,0.320005/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.91456,1.167396,1.08,1.228487,0.204483,0.320005/ &INIT ID='Particle Cloud', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=10, MASS_PER_VOLUME=450.0, XB=0.91456,1.16,0.778043,0.926531,0.204343,0.320005/

&OBST ID='Låda nedre del', XB=0.661031,0.671031,0.41587,1.54587,0.0,0.665, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Batteri', XB=0.924273,1.159273,0.928563,1.078563,0.205,0.315, RGB=12,0,0, SURF_ID='Brand'/ &OBST ID='Låda nedre del', XB=0.68,1.4,0.4,0.44,0.019231,0.646923, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Låda nedre del', XB=1.36,1.4,0.44,1.56,0.019231,0.646923, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Låda nedre del', XB=0.68,1.36,1.56,1.56,0.019231,0.646923, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=0.64,0.68,0.44,1.56,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=0.64,1.4,0.4,0.44,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=1.36,1.4,0.44,1.56,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=0.68,1.36,1.56,1.56,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda golv', XB=0.68,1.4,0.4,1.56,-0.02,0.019231, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,0.92,0.84,1.16,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,0.44,0.84,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,1.16,1.56,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.44,0.36,0.4,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=1.16,1.4,0.84,1.16,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=1.4,1.44,0.4,1.6,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,0.4,0.44,0.686154,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,1.56,1.6,0.686154,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.92,1.16,0.84,1.16,0.725385,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,1.6,1.6,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,0.64,0.44,1.6,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,0.64,0.4,0.44,0.646923,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,1.24,0.76,1.24,0.215385,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,1.24,0.76,0.76,0.019231,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,1.24,1.24,1.24,0.019231,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,0.88,0.76,1.24,0.019231,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=1.24,1.24,0.76,1.24,0.019231,0.215385, SURF_ID='Stål'/

&HOLE ID='Hole', XB=0.64067,0.670671,0.42677,1.546534,0.665,0.67/ &HOLE ID='Hole', XB=1.374204,1.404205,0.42593,1.545694,0.665,0.67/ &HOLE ID='Hole', XB=0.660796,0.671136,0.415652,1.556099,0.647197,0.665/ &HOLE ID='Hole', XB=1.374047,1.384293,0.415525,1.55613,0.647108,0.665/ &HOLE ID='Hole', XB=0.660796,1.38533,1.54556,1.556227,0.647197,0.665/ &HOLE ID='Hole', XB=0.660796,1.384846,0.415285,0.425952,0.647197,0.665/ &HOLE ID='Hole', XB=0.64,1.4,0.4,0.44,0.646923,0.686154/ &HOLE ID='Hole', XB=0.64,1.4,1.56,1.6,0.646923,0.686154/

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

&SLCF QUANTITY='VISIBILITY', PBY=1.0/ &SLCF QUANTITY='VISCOSITY', PBY=1.0/ &SLCF QUANTITY='TEMPERATURE', PBX=1.38/ &SLCF QUANTITY='TEMPERATURE', PBX=0.68/ &SLCF QUANTITY='TEMPERATURE', PBZ=0.662692/ &SLCF QUANTITY='TEMPERATURE', VECTOR=.TRUE., XB=0.661032,1.38413,0.41587,1.55587,1.110223E-16,0.96441, FYI='Temperatur'/

&TAIL /

drjfloyd commented 2 years ago

What is the issue you have with this input file? Your original issue was a problem with porous media; however, you have defined the particle class 'skum' as being porous media in this input but you have not placed any of these particles in the domain. Are you reporting a new issue with the particle class PB_partikel?

Pethav-7 commented 2 years ago

Hello, Sorry for the mistake. I have now corrected the code and rerun the simulations. To repeat, the scenario I’m trying to simulate is smoke and heat transport from a battery fire (red arrow below), through a fixed granular media between the battery and the top of the steel box (area indicated with yellow arrow below, this area is also bolded in the code below), to for example a TC on the outside of a steel box. Of course the smoke doesn’t penetrate the steel box but it is the temperature on the outside of the steel box, measured with TC3, which is important. As can be seen in the code I have filled the steel box with several zones with granular media. But the main point I think is the yellow zone that should act as a porous media between the battery and the top of the steel box. image

Some results are shown in the figure below. image

The simulation matrix is as follows: image

As can be seen in the graphs a change in the conductivity affects the results, as expected. But it can also be seen that a change in porosity does not affect the results, which is very unexpected, I think. Have I misunderstood something about the porosity?

The code is as follows:

Simulering.fds Generated by PyroSim - Version 2022.1.0422 2022-jun-16 22:24:51

&HEAD CHID='kon_0,1365'/ &TIME T_END=300/ &DUMP DT_RESTART=300.0, DT_SL3D=0.25, NFRAMES=60/ &MISC POROUS_FLOOR=.FALSE./

&MESH ID='Mesh01', IJK=50,50,27, XB=0.0,2.0,0.0,2.0,-0.059231,1.0/

&PART ID='PB_partikel', STATIC=.TRUE., SURF_ID='Pyrobubbles', COLOR='AZURE', AGE=600.0/

&REAC ID='POLYURETHANE', FYI='NFPA Babrauskas', FUEL='REAC_FUEL', C=6.3, H=7.1, O=2.1, N=1.0, SOOT_YIELD=0.1/

&DEVC ID='THCP', QUANTITY='THERMOCOUPLE', XYZ=1.055657,1.304148,0.71/ &DEVC ID='THCP01', QUANTITY='THERMOCOUPLE', XYZ=1.041773,0.69362,0.710088/ &DEVC ID='THCP02', QUANTITY='THERMOCOUPLE', XYZ=1.041773,1.003315,0.71/ &DEVC ID='THCP03', QUANTITY='THERMOCOUPLE', XYZ=1.041773,0.702215,0.694522/ &DEVC ID='THCP04', QUANTITY='THERMOCOUPLE', XYZ=1.041773,1.003315,0.709/ &DEVC ID='THCP05', QUANTITY='THERMOCOUPLE', XYZ=1.055657,1.304148,0.694522/ &DEVC ID='THCP06', QUANTITY='THERMOCOUPLE', XYZ=1.36,0.995988,0.419053, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP07', QUANTITY='THERMOCOUPLE', XYZ=0.671031,0.995988,0.419053, ORIENTATION=1.0,0.0,0.0/ &DEVC ID='THCP08', QUANTITY='THERMOCOUPLE', XYZ=1.04161,1.555658,0.419053, ORIENTATION=0.0,-1.0,0.0/ &DEVC ID='THCP09', QUANTITY='THERMOCOUPLE', XYZ=1.04161,0.41587,0.419053, ORIENTATION=0.0,1.0,-1.0/ &DEVC ID='THCP10', QUANTITY='THERMOCOUPLE', XYZ=1.04161,0.41587,0.626681, ORIENTATION=0.0,1.0,-1.0/ &DEVC ID='THCP11', QUANTITY='THERMOCOUPLE', XYZ=1.36,0.993984,0.626681, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP12', QUANTITY='THERMOCOUPLE', XYZ=0.671031,0.99704,0.626681, ORIENTATION=1.0,0.0,0.0/ &DEVC ID='THCP13', QUANTITY='THERMOCOUPLE', XYZ=1.041483,1.555658,0.626681, ORIENTATION=0.0,-1.0,0.0/ &DEVC ID='THCP14', QUANTITY='THERMOCOUPLE', XYZ=1.384204,0.995988,0.419053, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP15', QUANTITY='THERMOCOUPLE', XYZ=1.384204,0.993984,0.626681, ORIENTATION=-1.0,0.0,0.0/ &DEVC ID='THCP16', QUANTITY='THERMOCOUPLE', XYZ=0.671031,0.995988,0.419053, ORIENTATION=1.0,0.0,0.0/

&MATL ID='STEEL', FYI='Drysdale, Intro to Fire Dynamics - ATF NIST Multi-Floor Validation', SPECIFIC_HEAT=0.46, CONDUCTIVITY=45.8, DENSITY=7850.0, EMISSIVITY=0.95/ &MATL ID='PB', SPECIFIC_HEAT=0.7, CONDUCTIVITY=0.1365, DENSITY=450.0, EMISSIVITY=0.93/

&SURF ID='Brand', COLOR='RED', HRRPUA=250.0, RAMP_Q='Brand_RAMP_Q', TMP_FRONT=100.0/ &RAMP ID='Brand_RAMP_Q', T=1.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=5.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=6.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=17.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=21.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=22.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=35.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=40.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=45.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=56.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=60.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=61.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=74.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=79.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=84.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=95.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=96.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=97.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=100.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=126.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=137.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=141.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=142.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=155.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=200.0, F=0.37/ &RAMP ID='Brand_RAMP_Q', T=300.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=315.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=316.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=317.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=320.0, F=0.0/ &RAMP ID='Brand_RAMP_Q', T=346.0, F=0.48/ &RAMP ID='Brand_RAMP_Q', T=347.0, F=0.51/ &RAMP ID='Brand_RAMP_Q', T=351.0, F=1.0/ &RAMP ID='Brand_RAMP_Q', T=352.0, F=0.92/ &RAMP ID='Brand_RAMP_Q', T=365.0, F=0.37/

&SURF ID='Stål', RGB=255,52,141, BACKING='VOID', MATL_ID(1,1)='STEEL', MATL_MASS_FRACTION(1,1)=1.0, THICKNESS(1)=1.0E-3/ &SURF ID='Pyrobubbles', RGB=255,102,0, MATL_ID(1,1)='PB', MATL_MASS_FRACTION(1,1)=1.0, THICKNESS(1)=5.0E-4, GEOMETRY='SPHERICAL', INNER_RADIUS=5.0E-3/

&PART ID='skum', DRAG_LAW='POROUS MEDIA', SURF_ID='Pyrobubbles', POROUS_VOLUME_FRACTION=0.85, STATIC=T,DRAG_COEFFICIENT=0.1,0.1,0.1, PERMEABILITY=1.0E-7,1.E-7,1.E-7 /

&INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.661034,1.38413,0.41587,1.55587,0.320765,0.639705/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=1.222775,1.38413,0.41587,1.55587,1.110223E-16,0.320765/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=1.166992,1.222773,0.778304,1.228273,0.205,0.320765/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.860773,0.914425,0.778522,1.228561,0.205,0.320765/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.661034,0.860714,0.41587,1.55587,0.0,0.320765/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.860714,1.222586,1.22856,1.55587,0.0,0.320005/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.860714,1.222586,0.416388,0.777676,0.0,0.320005/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.91456,1.167396,1.08,1.228487,0.204483,0.320005/ &INIT ID='skum', PART_ID='PB_partikel', N_PARTICLES_PER_CELL=1, MASS_PER_VOLUME=450.0, CELL_CENTERED=T, XB=0.91456,1.16,0.778043,0.926531,0.204343,0.320005/

&OBST ID='Låda nedre del', XB=0.661031,0.671031,0.41587,1.54587,0.0,0.665, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Batteri', XB=0.924273,1.159273,0.928563,1.078563,0.205,0.315, RGB=12,0,0, SURF_ID='Brand'/ &OBST ID='Låda nedre del', XB=0.68,1.4,0.4,0.41,0.019231,0.646923, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Låda nedre del', XB=1.39,1.4,0.44,1.56,0.019231,0.646923, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Låda nedre del', XB=0.68,1.36,1.56,1.56,0.019231,0.646923, COLOR='INVISIBLE', SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=0.64,0.68,0.44,1.56,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=0.64,1.4,0.4,0.44,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=1.36,1.4,0.44,1.56,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda övre del', XB=0.68,1.36,1.56,1.56,0.646923,0.686154, SURF_ID='Stål'/ &OBST ID='Låda golv', XB=0.68,1.4,0.4,1.56,-0.02,0.019231, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,0.92,0.84,1.16,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,0.44,0.84,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,1.16,1.56,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.44,0.36,0.4,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=1.16,1.4,0.84,1.16,0.686154,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=1.4,1.44,0.4,1.6,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,0.4,0.44,0.686154,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,1.56,1.6,0.686154,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.92,1.16,0.84,1.16,0.725385,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,1.4,1.6,1.6,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,0.64,0.44,1.6,0.646923,0.686154, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Lock', XB=0.64,0.64,0.4,0.44,0.646923,0.725385, RGB=0,153,153, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,1.24,0.76,1.24,0.215385,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,1.24,0.76,0.76,0.019231,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,1.24,1.24,1.24,0.019231,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=0.88,0.88,0.76,1.24,0.019231,0.215385, SURF_ID='Stål'/ &OBST ID='Batterihållare', XB=1.24,1.24,0.76,1.24,0.019231,0.215385, SURF_ID='Stål'/

&HOLE ID='Hole', XB=0.64067,0.670671,0.42677,1.546534,0.665,0.67/ &HOLE ID='Hole', XB=1.374204,1.404205,0.42593,1.545694,0.665,0.67/ &HOLE ID='Hole', XB=0.660796,0.671136,0.415652,1.556099,0.647197,0.665/ &HOLE ID='Hole', XB=1.374047,1.384293,0.415525,1.55613,0.647108,0.665/ &HOLE ID='Hole', XB=0.660796,1.38533,1.54556,1.556227,0.647197,0.665/ &HOLE ID='Hole', XB=0.660796,1.384846,0.415285,0.425952,0.647197,0.665/ &HOLE ID='Hole', XB=0.64,1.4,0.4,0.44,0.646923,0.686154/ &HOLE ID='Hole', XB=0.64,1.4,1.56,1.6,0.646923,0.686154/

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

&SLCF QUANTITY='VISIBILITY', PBY=1.0/ &SLCF QUANTITY='VISCOSITY', PBY=1.0/ &SLCF QUANTITY='TEMPERATURE', PBX=1.38/ &SLCF QUANTITY='TEMPERATURE', PBX=0.68/ &SLCF QUANTITY='TEMPERATURE', PBZ=0.662692/ &SLCF QUANTITY='TEMPERATURE', VECTOR=.TRUE., XB=0.661032,1.38413,0.41587,1.55587,1.110223E-16,0.96441, FYI='Temperatur'/

&TAIL /

drjfloyd commented 2 years ago

Try removing the MASS_PER_VOLUME. This is not in the example in the User's Guide. You are specifying twice the mass to place in the volume: once by defining the porous volume fraction and a second time by giving MASS_PER_VOLUME.

As a note, you defined your particles as a sphere with an 0.5 mm outer layer of PB around a 5 mm radius void space. The MATL PB has a density of 450 kg/m3 making the effective density of a single spherical particle ~112 kg/m3. For a sphere, the maximum close packing volume fraction is ~64 %. Therefore, the maximum mass that physically makes sense to specify as a MASS_PER_VOLUME would be ~72 kg/m3.

Pethav-7 commented 2 years ago

I tried removing the MASS_PER_VOLUME parameter for the &INIT (bolded in previous post) and did a rerun of the simulations. I still get the same temperature measurments for the "base"-simulation and the "lower porosity"-simulation.

Any other idea on why there's no difference between them?

Resaults displayed in picture below. Screenshot_2

mcgratta commented 2 years ago

POROUS MEDIA is rarely used, as far as I can tell. If you believe that there is something wrong with the algorithm, create a very simple test case for which there is a known solution. That will enable us to better diagnose a problem. To just say that a particular property did not make a big difference does not tell us if there is something really wrong.

drjfloyd commented 1 year ago

No further responses from OP. Closing.