TMVOC
hello,everyone.I am currently working on a simple pollutant leakage simulation model, and the initial flow field has been simulated. When proceeding to the next step of simulating point source pollutant leakage, an unknown error occurred.As shown in the image below.What might be the reason for this?
4 replies
-
In TMVOC T is always the last primary variable (see table 3 of the user's guide), but in your INCON the last primary variable is 0., which cannot be handled as T shall be higher than the triple point of water.
I guess T=20.°C is assigned as the 3rd primary variable, as it was in the input file of your previous post.
Are you trying to use a steady state flow field with no VOC to initialize the VOC spill by adding one more mass component (the VOC) to the INCON block?
In this case you should assign in MULTI:
NK=3 (H2O, VOC#1, NCG#1)
NEQ=4 (for non isothermal, 3 for isothermal)
NPH= 3
NB=6
NKIN=2 to instruct TMVOC to read an INCON block for only 2 mass components. Look at NKIN in the TMVOC user's guide, page 154.
-
It seems to me you have to assign NKIN=2, while you have not assigned a value to NKIN and then by default NKIN=NK=3. from your output file:
OPTIONS SELECTED ARE: (NK,NEQ,NPH,NB) = ( 3, 3, 3, 8)
NK = 3 - NUMBER OF FLUID COMPONENTS: (1) WATER, NCG NON-CONDENSIBLE GASES, FOLLOWED BY NHC=NK-1-NCG VOCS
NEQ = 3 - NUMBER OF EQUATIONS PER GRID BLOCK
NPH = 3 - NUMBER OF PHASES THAT CAN BE PRESENT
NB = 8 - NUMBER OF SECONDARY PARAMETERS (OTHER THAN COMPONENT MASS FRACTIONS)NKIN= 3 - NUMBER OF MASS COMPONENTS PROVIDED IN SPECIFICATIONS OF INITIAL CONDITIONS (DEFAULT IS NK)
NCG = 1 - NUMBER OF NON-CONDENSIBLE GASES
NHC = 1 - NUMBER OF VOLATILE ORGANIC COMPOUNDSYou need to assign NKIN=2 in block MULTI:
MULTI----1----*----2----*----3----*----4----*----5----*----6----*----7----*----8
3 3 3 8 2Your INCON file is for 2 mass components (H2O, AIR), but you assigned NK=3 for 3 mass components (H2O,DNAPL,AIR). TMVOC shall be informed that it is going to read an INCON file for just 2 components by assigning NKIN=2. Then TMVOC assumes the additional component is a VOC and will reorder the primary variables with the right format.
Your INCON:
A11 1 0.99900000E+00 1
0.1010000000000E+06 0.9900000000000E+00 0.2000000000000E+02P XGair T
The wrong primary variables by using NKIN=3 (default NKIN=NK):
+++++++++ CANNOT FIND PARAMETERS AT ELEMENT *A11 1* IND = 1
+++++++++ XX(M) = 0.101000E+06 0.990000E+00 0.200000E+02 0.000000E+00P XGvoc XGair T
"Cannot find parameters" means TMVOC cannot find a thermodynamic state corresponding to given primary variables. A single gas phase (IND=1) cannot exist with 0.99 mole fraction of DNAPL at P= 1.01 bara and T=0.°C. As it occurs at the very first iteration, the run simply stops.
The correct primary variables assigning NKIN=2 will be
XX(M) = 0.101000E+06 0.000000E+00 0.990000E+00 0.200000E+02
P XGvoc XGair T
Read the detailed explanation at pag. 154 of TMVOC user's guide (2002). I assume the TMVOC-REACT version you are using has not changed NKIN input option.
-
Dear,Alfredo Battistelli . thank you for your assistance. Adding the keyword you mentioned has already yielded results.