Error when defining new liquid species

220 views
Skip to first unread message

Vishwajeeth Achutha

unread,
Jan 22, 2024, 5:10:33 AM1/22/24
to FDS and Smokeview Discussions
I specified the following lines of code to define a new liquid species to model droplets. This includes all the required parameters required as mentioned in FDS User Guide:

&SPEC ID='New', DENSITY_LIQUID=700, SPECIFIC_HEAT_LIQUID=1.67, VAPORIZATION_TEMPERATURE = 360,MELTING_TEMPERATURE = -24,
HEAT_OF_VAPORIZATION = 950,
ENTHALPY_OF_FORMATION = -400,
H_V_REFERENCE_TEMPERATURE = 316,
VISCOSITY_LIQUID = 0.027,
CONDUCTIVTY_LIQUID = 0.136,
BETA_LIQUID = 0.0098,
MASS_EXTINCTION_COEFFICIENT=3000. /

However.. I have received this error...
-----------------------------------------------------------------
 Starting FDS ...

 MPI Process      0 started on BRLL902043607

 Reading FDS input file ...

WARNING: MPI_PROCESS set for MESH 2 and only one MPI process exists
WARNING: MPI_PROCESS set for MESH 3 and only one MPI process exists
WARNING: MPI_PROCESS set for MESH 4 and only one MPI process exists
WARNING: MPI_PROCESS set for MESH 5 and only one MPI process exists

ERROR: Problem with SPECies number 1, line number 0 (CHID: Droplet_3V)

ERROR: FDS was improperly set-up - FDS stopped (CHID: Droplet_3V)
---------------------------------------------------------------

Kindly do let me know about this.

Thank you
BR,
Vishwa

dr_jfloyd

unread,
Jan 22, 2024, 5:53:17 AM1/22/24
to FDS and Smokeview Discussions
Set T_END=0 so FDS just processes the input.  Then take the / and move it to the midpoint of the input.  For this line that looks to be after HEAT_OF_VAPORIZATION.  Run FDS. If you still get the same error you know it is in the first half of the input. If you don't get the same error, the error is in the second half. Move the / to mid point of the block of input containing the error. Run FDS again. Repeat the process until you have indentified where the error is. Then check that input to make sure that it is in the table of inputs for SPEC, that it is spelled correctly, and that it has the right type in input (i.e., a logical has T or F, an integer does not have a real, etc.).

Vishwajeeth Achutha

unread,
Jan 22, 2024, 6:09:41 AM1/22/24
to FDS and Smokeview Discussions
Thank you. The issue was in 'CONDUCTIVITY_LIQUID'. This parameter is not been recognized by FDS6.7.9. The rest of the parameters seem to work.

dr_jfloyd

unread,
Jan 22, 2024, 6:18:40 AM1/22/24
to FDS and Smokeview Discussions
CONDUCTIVTY_LIQUID
         ^ no I


Vishwajeeth Achutha

unread,
Jan 22, 2024, 6:27:41 AM1/22/24
to FDS and Smokeview Discussions
I corrected the spelling mistake and re-ran it..  'CONDUCTIVITY_LIQUID = 0.136'
It showed the same error.

dr_jfloyd

unread,
Jan 22, 2024, 8:55:19 AM1/22/24
to FDS and Smokeview Discussions
Did you repeat the process to see if you can find another error?

Vishwajeeth Achutha

unread,
Jan 23, 2024, 3:13:34 AM1/23/24
to FDS and Smokeview Discussions
It is working now.. Thanks a lot!!

Vishwajeeth Achutha

unread,
Jan 29, 2024, 2:56:01 AM1/29/24
to FDS and Smokeview Discussions
I have been trying to simulate species droplet evaporation. Droplet dia is app. 0.3micrometer, I have initialized with MASS_PER_VOLUME=0.0001, N_PARTICLES_PER_CELL=3000 droplets in a 3mm volumetric cube. There is also a heat source below this which helps in droplet/species transport. I am unable to see any of these droplets/species in smokeview. I also tried capturing the mass fraction contour of the species, but nothing visible/existing.
Unable to figure this out.

dr_jfloyd

unread,
Jan 29, 2024, 5:46:34 AM1/29/24
to FDS and Smokeview Discussions
Need to know what your inputs are.

Vishwajeeth Achutha

unread,
Jan 29, 2024, 6:16:01 AM1/29/24
to FDS and Smokeview Discussions
I have attached the code snippet here.. I played around with 
 MASS_PER_VOLUME, N_PARTICLES_PER_CELL and MASS_PER_TIME too... But was not able to release droplets or visualize them. Also, does the RAMP_PART work for 'MASS_PER_TIME' too..?


&TIME T_END=10.0, DT=0.01 /

&MESH ID='M1' IJK=176,184,72, XB=0.0,6.4,0.0,4.6,0.0,1.805, MPI_PROCESS=0 /
&MESH ID='M2' IJK=176,92,36, XB=6.4,14.4,0.0,4.6,0.0,1.805, MPI_PROCESS=1 /

&SPEC ID='PARAFFIN', DENSITY_LIQUID=850, SPECIFIC_HEAT_LIQUID=1.67, VAPORIZATION_TEMPERATURE = 400,MELTING_TEMPERATURE = -24,
HEAT_OF_VAPORIZATION = 977  ,
ENTHALPY_OF_FORMATION = -400 ,

H_V_REFERENCE_TEMPERATURE = 316,
VISCOSITY_LIQUID = 0.027,
MASS_EXTINCTION_COEFFICIENT=9100,
CONDUCTIVITY_LIQUID = 0.136/

&PART ID='droplets', SPEC_ID='PAR',INITIAL_TEMPERATURE=315, MONODISPERSE= .TRUE.,
DIAMETER= 0.3,
QUANTITIES(1:3)='PARTICLE DIAMETER','PARTICLE TEMPERATURE', 'PARTICLE VELOCITY'/

&INIT PART_ID='droplets', XB=0.100,0.316,2.300,2.516,0.439,0.449, MASS_PER_VOLUME=0.001, N_PARTICLES_PER_CELL=50000 /

RAMP ID='cycle', T= 0, F=0.0 /
RAMP ID='cycle', T= 21.9, F=0.0 /
RAMP ID='cycle', T= 22, F= 1 /
RAMP ID='cycle', T=27.9, F= 1 /
RAMP ID='cycle', T=28, F=0.0 /

RAMP ID='cycle', T= 35.9, F=0.0 /
RAMP ID='cycle', T= 36, F= 1 /
RAMP ID='cycle', T=38.9, F= 1 /
RAMP ID='cycle', T=39, F=0.0 /

RAMP ID='cycle', T= 44.9, F=0.0 /
RAMP ID='cycle', T= 45, F= 1 /
RAMP ID='cycle', T=46.9, F= 1 /
RAMP ID='cycle', T=47, F=0.0 /

RAMP ID='cycle', T= 52.9, F=0.0 /
RAMP ID='cycle', T= 53, F= 1 /
RAMP ID='cycle', T=53.9, F= 1 /
RAMP ID='cycle', T=54, F=0.0 /


--- Geometric name
&OBST ID='Generator_full', XB=0.100,0.316,2.300,2.516,0.02,0.438, COLOR='BLACK' /

&SURF ID='C_top', VEL=-0.75, TMP_FRONT=66. /
&VENT XB=0.100,0.316,2.300,2.516,0.438,0.438, SURF_ID = 'C_top',COLOR='RED' /

&SURF ID='C_bottom', VEL=0.48, TMP_FRONT=28. /
&VENT XB=0.100,0.316,2.300,2.516,0.02,0.02,SURF_ID = 'C_bottom',COLOR='BLUE' /


--- Beam detector
&PROP ID='Smoke', QUANTITY='PATH OBSCURATION', SPEC_ID='PAR', PART_ID='droplets' /
&DEVC XB= 5.26,5.26,1.8,2.8,1.68,1.68, PROP_ID='Smoke', ID='beam_1' /


--- Slice files

&SLCF PBY=2.408, QUANTITY='MASS FRACTION', SPEC_ID='PAR' /

&SLCF XB=0.0,3.0,0.0,4.0,0.0,1.6, QUANTITY='MASS FRACTION', SPEC_ID='PAR' /


DUMP SMOKE3D_QUANTITY='MASS FRACTION', SMOKE3D_SPEC_ID='PAR' /
&ISOF QUANTITY='TEMPERATURE',VALUE(1)=50., VALUE(2)=200., VALUE(3)=500., VALUE(4)=1000., /
&ISOF QUANTITY='VELOCITY',VALUE(1)=0.5, VALUE(2)=1 ,VALUE(3)=10 /

&DUMP DT_PL3D = 10, MASS_FILE=.TRUE. ,  PLOT3D_QUANTITY(1:5)='TEMPERATURE',
'U-VELOCITY','V-VELOCITY','W-VELOCITY','MASS FRACTION', PLOT3D_SPEC_ID(5)='PAR'/

&TAIL /

Thank you

Randy McDermott

unread,
Jan 29, 2024, 9:46:36 AM1/29/24
to fds...@googlegroups.com
50,000 particles per cell?  I don't think you need this.  FDS will weight the particles to get the right mass.  In the attached test case, I've used 500 part per cell.  And I fixed a couple of other problem with your input file.

--
You received this message because you are subscribed to the Google Groups "FDS and Smokeview Discussions" group.
To unsubscribe from this group and stop receiving emails from it, send an email to fds-smv+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/fds-smv/4b7b28b6-295a-4259-a0ca-92774c24cc39n%40googlegroups.com.
test.fds
Screenshot 2024-01-29 at 9.43.54 AM.png

dr_jfloyd

unread,
Jan 29, 2024, 10:40:28 AM1/29/24
to FDS and Smokeview Discussions
You haven't defined any gas phase properties for the paraffiin. Make sure that the assumptions FDS is making for the gas phase species properties are appropriate for your needs.

The units for DIAMETER on PART is micrometers. Do you have the correct units for your input? 0.3 micrometers is a very small droplet.

Vishwajeeth Achutha

unread,
Jan 31, 2024, 4:54:04 AM1/31/24
to FDS and Smokeview Discussions
Thank you for the responses. Yes, I require to define only liquid properties.. and the droplet value also corresponds to 0.3micrometers.
I tried working with the 'test.fds'. Thanks a lot for the input. The issue was the droplets were released only in the initial few seconds, during which even the species was observable and then nothing visible.. However, 
1) I would like to release the droplets continuously from the domain for the entire duration (10s). 
2) And would like to visualize the 3D smoke for the entire duration.

Thank you.

Best regards,
Vishwajeeth

Randy McDermott

unread,
Jan 31, 2024, 6:57:48 AM1/31/24
to fds...@googlegroups.com
Reply all
Reply to author
Forward
0 new messages