@@@@@ @@ @ @ @@@ @ @ @@ @@@ @ @ @ @ @ @ @@ @@@@@ @ @@ @ @ @@@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @@ @@ @ @ @ @ @ @ @ @ @ @@ @ @ @ @ @ @@ @ @ @ @ @ @ @ @@ @@@@ @ @@ @ @ @ @ @ @ @ @@@@ @ @ @ @ @ @ @ @@@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @@ @ @ @ @ @ @@ @ @@ @@ @@@ @ @ @@@@ @@@ @ @ @ @@ @@@@ @ @ @ @ @@ @ @ @ @ @@ @ @ TOUGH2 IS A PROGRAM FOR MULTIPHASE MULTICOMPONENT FLOW IN PERMEABLE MEDIA, INCLUDING HEAT FLOW. IT IS A MEMBER OF THE MULKOM FAMILY OF CODES, DEVELOPED AT LAWRENCE BERKELEY NATIONAL LABORATORY. ******************************************************************************** ******************* ******************** ******************* TOUGH2 - VERSION 2.0 (OCTOBER 1999) ******************** ******************* T2CG2 Solver Package ******************** ******************* ******************** ******************************************************************************** Copyright 1999 by The Regents of the University of California (subject to approval by the U.S. Department of Energy). NOTICE: This software was developed under funding from the U.S. Department of Energy and the U.S. Government consequently retains certain rights as follows: the U.S. Government has been granted for itself and others acting on its behalf a paid-up, nonexclusive, irrevocable, worldwide license in the software to reproduce, prepare derivative works, and perform publicly and display publicly. Beginning five (5) years after the date permission to assert copyright is obtained from the U.S. Department of Energy, and subject to any subsequent five (5) year renewals, the U.S. Government is granted for itself and others acting on its behalf a paid-up, nonexclusive, irrevocable, worldwide license in the software to reproduce, prepare derivative works, distribute copies to the public, perform publicly and display publicly, and to permit others to do so. PARAMETERS FOR FLEXIBLE DIMENSIONING OF MAJOR ARRAYS (MAIN PROGRAM) ARE AS FOLLOWS MNEL = 60000 MNCON = 170000 MNEQ = 6 MNK = 5 MNPH = 3 MNB = 8 MNOGN = 2000 MGTAB = 2000 =================================================================================================================================== MAXIMUM NUMBER OF VOLUME ELEMENTS (GRID BLOCKS): MNEL = 60000 MAXIMUM NUMBER OF CONNECTIONS (INTERFACES): MNCON = 170000 MAXIMUM LENGTH OF PRIMARY VARIABLE ARRAYS: MPRIM = 360000 MAXIMUM NUMBER OF GENERATION ITEMS (SINKS/SOURCES): MNOGN = 2000 MAXIMUM NUMBER OF TABULAR (TIME-DEPENDENT) GENERATION DATA: MGTAB = 2000 LENGTH OF SECONDARY PARAMETER ARRAY: MSEC =17220000 MAXIMUM NUMBER OF JACOBIAN MATRIX ELEMENTS: MNZ =14400000 LARGE LINEAR EQUATION ARRAYS: LENGTH OF IRN IS LIRN =14400000 LENGTH OF ICN AND CO IS LICN =14400000 =================================================================================================================================== array dimensioning is made according to the needs of the conjugate gradient solvers when using MA28 or LUBAND, only a smaller-size problem can be accommodated restriction with MA28 is: {number of elements} + 2 * {number of connections} < {MNEL + 2* MNCON}/4 =================================================================================================================================== SUMMARY OF DISK FILES FILE *VERS* DOES NOT EXIST --- OPEN AS A NEW FILE FILE *MESH* DOES NOT EXIST --- OPEN AS A NEW FILE FILE *INCON* DOES NOT EXIST --- OPEN AS A NEW FILE FILE *GENER* DOES NOT EXIST --- OPEN AS A NEW FILE FILE *SAVE* DOES NOT EXIST --- OPEN AS A NEW FILE FILE *LINEQ* DOES NOT EXIST --- OPEN AS A NEW FILE FILE *TABLE* DOES NOT EXIST --- OPEN AS A NEW FILE =================================================================================================================================== PROBLEM TITLE: *SAM7C3* ... one-d radial injection into saturated region. DOMAIN NO. 1 MATERIAL NAME -- domen DOMAIN NO. 2 MATERIAL NAME -- bndry WRITE FILE *GENER* FROM INPUT DATA WRITE FILE *INCON* FROM INPUT DATA WRITE FILE *MESH* FROM INPUT DATA ************************************************************************************** * * * Comments in Input Data (have encountered 10 records with unknown keywords) * * (print up to 50 of them) * * * ************************************************************************************** * * * ----*----1 MOP: 123456789*123456789*1234 ---*----5----*----6----*----7----*----8 * * * * IE(14) = 1, use SUPST and GASEOS for enthalpy. IE(14) = 0, use GASEOS only. * * IE(15) = 1, use Peng-Robinson eqn. of state. IE(15) = 2,3 use RK,SRK, resp. * * IE(16) = 1, NCG is CO2. IE(16) = 2, NCG is N2. * * * * xOLVR ---1----*----2----*----3----*----4----*----5----*----6----*----7----*----8 * * 3 Z2 O2 * * * * * * * ************************************************************************************** ************************************************************************************ * EVALUATE FLOATING POINT ARITHMETIC * ************************************************************************************ * * * FLOATING POINT PROCESSOR HAS APPROXIMATELY 77 SIGNIFICANT DIGITS * * * * DEFAULT VALUE OF INCREMENT FACTOR FOR NUMERICAL DERIVATIVES IS DFAC = 0.1036E-38 * * DEFAULT VALUE FOR DFAC WILL BE USED * * * ************************************************************************************ * The numbers of connections in the X-, Y- and Z-directions are 11, 0 and 0 respectively * * * ********************************************************************************************************************************** all NCON = 11 connections read from file *MESH* reference known elements, and have been initialized to the data arrays RFILE DEBUG: Finished reading MESH file. NEL = 12 RFILE DEBUG: ELEM array sample from MESH (first 15 or NEL): RFILE DBG: ELEM( 1) = >ELE11< RFILE DBG: ELEM( 2) = >ELE12< RFILE DBG: ELEM( 3) = >ELE13< RFILE DBG: ELEM( 4) = >ELE14< RFILE DBG: ELEM( 5) = >ELE15< RFILE DBG: ELEM( 6) = >ELE16< RFILE DBG: ELEM( 7) = >ELE17< RFILE DBG: ELEM( 8) = >ELE18< RFILE DBG: ELEM( 9) = >ELE19< RFILE DBG: ELEM( 10) = >ELE20< RFILE DBG: ELEM( 11) = >ELE21< RFILE DBG: ELEM( 12) = >ELE22< RFILE DEBUG: --- End of ELEM array sample from MESH --- &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& Summary of capabilities for random permeability modification &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& Modification of absolute permeability on a grid block-by-grid block basis will be made when a domain "SEED " is present in data block "ROCKS", as follows. k ---> k' = k*m Here, k is the absolute permeability specified for the reservoir domain to which the grid block belongs. Parameter m is a "permeability modifier" which can be internally generated or externally prescribed by the user on a block-by-block basis. When permeability modification is in effect, the strength of capillary pressure will, following Leverett (1941), automatically be scaled as Pcap ---> Pcap' = Pcap/SQRT(m). User-supplied permeability modifiers have to be entered as parameter "PMX" in columns 41-50 of an ELEMEnt record. Permeability modification options are selected through parameters in data block "ROCKS". &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& Summary of available permeability modification options (with s - random number between 0 and 1; PMX - user-supplied modifiers in data block "ELEME"): (1) externally supplied: m = PMX - PER(2) (2) "linear" (DROK.ne.0): m = PER(1) * s - PER(2) (3) "logarithmic" (DROK.eq.0): m = exp(- PER(1) * s) - PER(2) &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&& if a domain "SEED " is present, permeability modification will be made &&&& if no domain "SEED " is present, no permeability modification will be made >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< >>>>>>>>>>>>>>>>>>>>>>>>>>> domain = "SEED " is not present, no permeability modification will be made <<<<<<<<<<<<<<<<<<<<<<<<<< >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< Data provided in domain "SEED " are used to select the following options. DROK = *** random number seed for internal generation of "linear" permeability modifiers. = 0: (default) no internal generation of "linear" permeability modifiers. > 0: perform "linear" permeability modification; random modifiers are generated internally with DROK as seed. POR = *** random number seed for internal generation of "logarithmic" permeability modifiers, = 0: (default) no internal generation of "logarithmic" permeability modifiers. > 0: perform "logarithmic" permeability modification; random modifiers are generated internally with POR as seed. &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&& note: if both DROK and POR are specified as non-zero, DROK takes precedence &&&&& &&&&& if both DROK and POR are zero, permeability modifiers as supplied through "ELEME" data will be used &&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& PER(1) = *** scale factor (optional) for internally generated permeability modifiers. = 0: (defaults to PER(1) = 1): permeability modifiers are generated as random numbers in the interval (0, 1). > 0: permeability modifiers are generated as random numbers in the interval (0, PER(1)). PER(2) = *** shift (optional) for internal or external permeability modifiers. = 0: (default) no shift is applied to permeability modifiers. > 0: permeability modifiers are shifted according to m' = m - PER(2). All m' < 0 are set equal to zero. REFERENCE TO UNKNOWN ELEMENT *A1 1* AT SOURCE # 1 --- WILL IGNORE SOURCE REFERENCE TO UNKNOWN ELEMENT *A1 1* AT SOURCE # 2 --- WILL IGNORE SOURCE RFILE DEBUG: Finished reading GENER file. NOGN = 43176 RFILE DEBUG: START path. Initializing X with DEP/INDOM. RFILE DEBUG: DEP array values (first nkin1): 0.00000E+00 0.65000E+02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 1)="ELE11", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 2)="ELE12", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 3)="ELE13", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 4)="ELE14", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 5)="ELE15", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 6)="ELE16", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 7)="ELE17", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 8)="ELE18", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 9)="ELE19", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 10)="ELE20", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 11)="ELE21", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: DEP Loop - ELEM( 12)="ELE22", DEP_Temp_C= 0.00000E+00, Set X_Temp_C= 0.00000E+00 RFILE DEBUG: Finished pre-initializing X with DEP/INDOM. RFILE DEBUG: After DEP/INDOM, ELEM( 1)="ELE11", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 2)="ELE12", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 3)="ELE13", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 4)="ELE14", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 5)="ELE15", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 6)="ELE16", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 7)="ELE17", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 8)="ELE18", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 9)="ELE19", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 10)="ELE20", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 11)="ELE21", T_C= 0.00000E+00 RFILE DEBUG: After DEP/INDOM, ELEM( 12)="ELE22", T_C= 0.00000E+00 RFILE DEBUG:Now reading specific INCON from RFILE DEBUG: Read header from temp INCON file RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> RFILE DEBUG: For EL=>ELE**< DEPU Temp_C= 0.65000E+02 WILL IGNORE INITIAL CONDITION AT UNKNOWN ELEMENT ELE** RFILE DEBUG: NO MATCH for INCON EL=>ELE** RFILE DEBUG: Read from unit 1: EL=> WWWWW WARNING WWWWW: INCON data at 12 unknown elements have been ignored Blank EL found on unit 1,INCON part returning. MESH HAS 12 ELEMENTS ( 12 ACTIVE) AND 11 CONNECTIONS (INTERFACES) BETWEEN THEM GENER HAS 2 SINKS/SOURCES ********************************************************************************************************************************** * * * M A T R I X S O L V E R A N D R E L E V A N T I N F O R M A T I O N * * * ********************************************************************************************************************************** The solver is determined from MOP(21) The solution method indicator MATSLV = 0 - reset internally to the default, MATSLV = 3 The solution method indicator MATSLV = 3 MATSLV = 1: MA28 NOT AVAILABLE MATSLV = 2: SUBROUTINE DSLUBC - BI-CONJUGATE GRADIENT SOLVER Incomplete LU factorization preconditioning MATSLV = 3: SUBROUTINE DSLUCS (DEFAULT) - Lanczos-type Conjugate Gradient Squared solver Incomplete LU factorization preconditioning MATSLV = 4: SUBROUTINE DSLUGM - Generalized Minimum Residual Conjugate Gradient solver Incomplete LU factorization preconditioning MATSLV = 5: SUBROUTINE DLUSTB - STABILIZED BI-CONJUGATE GRADIENT SOLVER Incomplete LU factorization preconditioning MATSLV = 6: SUBROUTINE LUBAND - Direct solver using LU decomposition RITMAX: Maximum # of CG iterations as fraction of the total number of equations = 1.00000E-01 (0.0 < RITMAX <= 1.0, Default = 0.1) CLOSUR: Convergence criterion for the CG iterations = 1.00000E-06 (1.0e-12 <= CLOSUR <= 1.0e-6, Default = 1.0e-6) NMAXIT: Maximum # of CG iterations - not to exceed the total number of equations NELA*NEQ = 20 (20 < NMAXIT <= NREDM) The matrix Z-preprocessing system is ZPROCS = Z1 ZPROCS = Z0: No Z-preprocessing; default for NEQ = 1 and for MATSLV = 1, 6 ZPROCS = Z1: Replacement of zeros on the main-diagonal by a small number; default for NEQ > 1 and for 1 < MATSLV < 6 ZPROCS = Z2: Linear combination of equations in each element to produce non-zero main diagonal entries ZPROCS = Z3: Normalization of equations, followed by Z2 ZPROCS = Z4: Same as in OPROCS = O4 The matrix O-preprocessing system is OPROCS = O0 OPROCS = O0: No O-preprocessing; default for NEQ = 1 and for MATSLV = 1, 6 OPROCS = O1: Elimination of lower half of the main-diagonal submatrix with center pivoting OPROCS = O2: O1+Elimination of upper half of the main-diagonal submatrix with center pivoting OPROCS = O3: O2+Normalization - Results in unit main-diagonal submatrices OPROCS = O4: pre-processing which results in unit main-diagonal submatrices without center pivoting END OF TOUGH2 INPUT JOB --- ELAPSED TIME = 0.0000 SECONDS *********************************************************************************************************************************** * ARRAY *MOP* ALLOWS TO GENERATE MORE PRINTOUT IN VARIOUS SUBROUTINES, AND TO MAKE SOME CALCULATIONAL CHOICES. * *********************************************************************************************************************************** MOP(1) = 0 *** ALLOWS TO GENERATE A SHORT PRINTOUT FOR EACH NEWTON-RAPHSON ITERATION = 0, 1, OR 2: GENERATE 0, 1, OR 2 LINES OF PRINTOUT MORE PRINTOUT IS GENERATED FOR MOP(I) > 0 IN THE FOLLOWING SUBROUTINES (THE LARGER MOP IS, THE MORE WILL BE PRINTED). MOP(2) = 3 *** CYCIT MOP(3) = 0 *** MULTI MOP(4) = 1 *** QU MOP(5) = 0 *** EOS MOP(6) = 0 *** LINEQ MOP(8) = 2 *** DISF (T2DM ONLY) MOP(7) = 0 *** IF UNEQUAL ZERO, WILL GENERATE A PRINTOUT OF INPUT DATA CALCULATIONAL CHOICES OFFERED BY MOP ARE AS FOLLOWS: MOP(9) = 0 *** CHOOSES FLUID COMPOSITION ON WITHDRAWAL (PRODUCTION). = 0: ACCORDING TO RELATIVE MOBILITIES. = 1: ACCORDING TO COMPOSITION IN PRODUCING ELEMENT. MOP(10) = 0 *** CHOOSES INTERPOLATION FORMULA FOR DEPENDENCE OF THERMAL CONDUCTIVITY ON LIQUID SATURATION (SL). = 0: K = KDRY + SQRT(SL)*(KWET-KDRY) = 1: K = KDRY + SL*(KWET-KDRY) MOP(11) = 0 *** CHOOSES EVALUATION OF MOBILITY AND ABSOLUTE PERMEABILITY AT INTERFACES. = 0: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS UPSTREAM WEIGHTED. = 1: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS UPSTREAM WEIGHTED. = 2: MOBILITIES ARE UPSTREAM WEIGHTED WITH WUP. (DEFAULT IS WUP = 1.0). PERMEABILITY IS HARMONIC WEIGHTED. = 3: MOBILITIES ARE AVERAGED BETWEEN ADJACENT ELEMENTS. PERMEABILITY IS HARMONIC WEIGHTED. = 4: MOBILITY * PERMEABILITY PRODUCT IS HARMONIC WEIGHTED. MOP(12) = 0 *** CHOOSES PROCEDURE FOR INTERPOLATING GENERATION RATES FROM A TIME TABLE. = 0: TRIPLE LINEAR INTERPOLATION. = 1: "STEP FUNCTION" OPTION. MOP(13) = 0 *** T2DM ONLY. SPECIFIES ASSIGNMENT OF COMPONENTS OF BOUNDARY VELOCITY AND CONCENTRATION GRADIENT VECTORS. AFFECTS ONLY THE INTERPOLATED COMPONENTS. DIRECT COMPONENTS ARE USED WHERE AVAILABLE. = 0: VELOCITY AND GRADIENT AT BOUNDARY ARE ZERO. = 1: VELOCITY IS ZERO; GRADIENT AT BOUNDARY IS NEAREST-NEIGHBOR. = 2: VELOCITY IS NEAREST NEIGHBOR; GRADIENT AT BOUNDARY IS ZERO. = 3: VELOCITY AND GRADIENT AT BOUNDARY ARE NEAREST-NEIGHBOR. MOP(14) = 0 *** SPECIFIES THE HANDLING OF PIVOT FAILURES IN THE LINEAR EQUATION SOLUTION (MA28 only) = 0: PERFORM NEW MATRIX DECOMPOSITION AFTER PIVOT FAILURE > 0: IGNORE PIVOT FAILURE AND PROCEED MOP(15) = 0 *** ALLOWS TO SELECT A SEMI-ANALYTICAL HEAT EXCHANGE CALCULATION WITH CONFINING BEDS. = 0: NO SEMI-ANALYTICAL HEAT EXCHANGE > 0: SEMI-ANALYTICAL HEAT EXCHANGE ENGAGED (WHEN A SPECIAL SUBROUTINE MODULE *QLOSS* IS PRESENT) MOP(16) = 1 *** PERMITS TO CHOOSE TIME STEP SELECTION OPTION = 0: USE TIME STEPS EXPLICITLY PROVIDED AS INPUT. > 0: INCREASE TIME STEP BY AT LEAST A FACTOR 2, IF CONVERGENCE OCCURS IN .LE. MOP(16) ITERATIONS. MOP(17) = 0 *** HANDLES SCALING OPTIONS. = 0: NO SCALING. = 7: SCALING. MOP(18) = 0 *** ALLOWS TO SELECT HANDLING OF INTERFACE DENSITY. = 0: PERFORM UPSTREAM WEIGHTING FOR INTERFACE DENSITY. > 0: COMPUTE INTERFACE DENSITY AS AVERAGE OF THE TWO GRID BLOCK DENSITIES. HOWEVER, WHEN ONE OF THE TWO PHASE SATURATIONS IS ZERO, DO UPSTREAM WEIGHTING. MOP(19) = 0 *** is used in some EOS-modules for selecting different options MOP(20) = 0 *** is used in some EOS-modules for selecting different options MOP(21) = 0 *** PERMITS TO SELECT LINEAR EQUATION SOLVER = 0: DEFAULTS TO MOP(21) = 3 = 1: DIRECT SOLVER MA28 IS NOT AVAILABLE IN THIS VERSION - INTERNALLY RESET TO 6 (LUBAND DIRECT SOLVER) = 2: SUBROUTINE DSLUBC: BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 3: SUBROUTINE DSLUCS: BI-CONJUGATE GRADIENT SOLVER - LANCZOS TYPE; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 4: SUBROUTINE DSLUGM: GENERALIZED MINIMUM RESIDUAL CONJUGATE GRADIENTS; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 5: SUBROUTINE DLUSTB: Stabilized bi-conjugate gradient solver; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 6: SUBROUTINE LUBAND: Direct solver using LU decomposition MOP(22) = 0 *** T2DM ONLY. SPECIFIES METHOD OF SUMMATION OF DUPLICATE ELEMENTS IN THE JACOBIAN. *** = 0: USE SUMDUP2. = 1: USE SUMDUP. = 2: IF MATSLV = 1, USE MA28 INTERNAL SUMMATION. MOP(23) = 0 *** T2DM ONLY. HANDLES EFFECTS OF NON-CONNECTED NEIGHBOR GRID BLOCKS ON DISPERSIVE FLUX IN DISF. *** = 0: INCLUDE INFLUENCE OF NEIGHBOR GRID BLOCKS. (MORE ACCURATE JACOBIAN, SLOWER LINEAR EQUATION SOLUTION). = 1: NEGLECT INFLUENCE OF NEIGHBOR GRID BLOCKS. (LESS ACCURATE JACOBIAN, FASTER LINEAR EQUATION SOLUTION). MOP(24) = 0 *** PERMITS TO SELECT HANDLING OF MULTIPHASE DIFFUSIVE FLUXES AT INTERFACES = 0: HARMONIC WEIGHTING OF FULLY-COUPLED EFFECTIVE MULTIPHASE DIFFUSIVITY = 1: SEPARATE HARMONIC WEIGHTING OF GAS AND LIQUID PHASE DIFFUSIVITIES. *********************************************************************************************************************************** Use diffusion coefficients from data block DIFFU. DEBUG T0 VALUE BEFORE ASSIGNMENT TO TBASE: T0 = 25.0000000000000 DEBUG: He diff. by GASDIF. DEBUG EOS ICALL1: TBASE_FOR_DIFFUSIVITY = 1.709061403020980E-312 MOP(20) = 0 *** NO LONGER USED IN EOS7C. *** *********************************************************************************************************************************** * EOS7C: EQUATION OF STATE FOR MIXTURES OF WATER/BRINE/NON-CONDENSIBLE GAS/TRACER/HELIUM * *********************************************************************************************************************************** ENTHALPY OPTION: SUPST (IE(14) = 1) EQUATION OF STATE: PR (IE(15) = 1) NON-CONDENSIBLE GAS: CO2 (IE(16) = 1) OPTIONS SELECTED ARE: (NK,NEQ,NPH,NB,NKIN) = (5,6,2,8,5) NK = 5 - NUMBER OF COMPONENTS NEQ = 6 - NUMBER OF EQUATIONS PER GRID BLOCK NPH = 2 - NUMBER OF PHASES THAT CAN BE PRESENT NB = 8 - NUMBER OF SECONDARY PARAMETERS (OTHER THAN' COMPONENT MASS FRACTIONS) NKIN = 5 - number of components for initializing thermodynamic conditions (default is NKIN = NK) For NB = 6, diffusion is "off", for NB = 8, diffusion is "on" AVAILABLE OPTIONS for (NK,NEQ,NPH,NB): (5,5,2,6 or 8) - WATER, BRINE, NCG, TRC, HE; ISOTHERMAL; VARIABLES (P, XB, XNCG, XTRC, X OR S+10, T) (5,6,2,6 or 8) - WATER, BRINE, NCG, TRC, CH4; NON-ISOTHERMAL; VARIABLES (P, XB, XNCG, XTRC, X OR S+10, T) (4,4,2,6 or 8) - WATER, BRINE, NCG, TRC, NO CH4; ISOTHERMAL; VARIABLES (P, XB, XNCG, XTRC, T) (4,5,2,6 or 8) - WATER, BRINE, NCG, TRC, NO CH4; NON-ISOTHERMAL; VARIABLES (P, XB, XNCG, XTRC, T) NKIN = NK or NKIN = NK-2. Default options are (5,5,2,8) - isothermal, diffusion "on", NKIN=NK THE NK = 4 ("NO CH4") OPTIONS MAY ONLY BE USED FOR PROBLEMS WITH SINGLE-PHASE LIQUID CONDITIONS THROUGHOUT. THE NORMAL NUMBER OF SECONDARY PARAMETERS OTHER THAN MASS FRACTIONS IS 6 PER PHASE. IN EOS7C, WE OPTIONALLY ADD TO THIS A SATURATION-DEPENDENT TORTUOSITY FOR EACH PHASE, AS WELL AS TEMPERATURE AND PRESSURE DEPENDENCE OF THE DIFFUSION COEFFICIENT. *********************************************************************************************************************************** NKIN = 5 *** ALLOWS INITIALIZATION WITH DIFFERENT SETS OF PRIMARY VARIABLES. *** *** THIS IS USEFUL FOR STARTING EOS7C SIMULATIONS FROM EOS7 INITIAL CONDITIONS. *** = NK (default): (P,XB,XNCG,XTRC,XCH4,T) FOR SINGLE PHASE, (P,XB,XNCG,XTRC,S+10,T) FOR TWO-PHASE. (EOS7C FORMAT). = NK-2: (P,XB,XCH4,T) FOR SINGLE PHASE, (P,XB,S+10,T) FOR TWO-PHASE. (EOS7 FORMAT). WILL INITIALIZE XNCG = XTRC = 0. *********************************************************************************************************************************** THE PRIMARY VARIABLES ARE P - PRESSURE T - TEMPERATURE XB - BRINE MASS FRACTION XNCG - MASS FRACTION OF NON-CONDENSIBLE GAS XTRC - MASS FRACTION OF VOLATILE TRACER S+10. - (GAS PHASE SATURATION + 10.) X - HELIUM MASS FRACTION IN GAS T - TEMPERATURE ****************************** ***************************************************************** * COMPONENTS * * FLUID PHASE CONDITION PRIMARY VARIABLES * ****************************** ***************************************************************** * * * * * # 1 - WATER * * SINGLE-PHASE GAS (#) P, XB, XNCG, XTRC, X, T * * * * * * # 2 - BRINE * * SINGLE-PHASE LIQUID (*) P, XB, XNCG, XTRC, X, T * * * * * * # 3 - NCG * * TWO-PHASE (*) P, XB, XNCG, XTRC, S+10., T * * * * * * # 4 - TRC * ***************************************************************** * * (#) XB, XNCG, XTRC AND X ARE MASS FRACTIONS IN THE GAS PHASE. * # 5 - CH4 * (*) XB, XNCG, XTRC AND X ARE MASS FRACTIONS IN THE AQUEOUS PHASE. * * * # 6 - HEAT * * * *********************************************************************************************************************************** NEGATIVE REFERENCE PRESSURE OF -.100000E+01 PA WAS SPECIFIED, THUS BRINE PROPERTIES ARE IDENTICAL TO WATER FOR ALL SALINITIES. PROPERTIES OF NCG AND TRACER: DOMAIN NCG TRC HALF-LIFE (SECONDS): -ALL- NO DECAY 0.6000E+01 MOLECULAR WEIGHT (GM/MOLE): -ALL- 0.4401E+02 0.2000E+02 INVERSE HENRY CONST. (1/PA): -ALL- VARIABLE 0.6000E+01 HEAT CAPACITY (J/(KG K)): -ALL- VARIABLE NEGLIGIBLE GAS PHASE DIFFUSIVITY (M**2/S): -ALL- 0.1000E-04 0.1000E-04 AQ. PHASE DIFFUSIVITY (M**2/S): -ALL- 0.1000E-09 0.1000E-09 DISTRIBUTION COEFF. (M**3/KG): domen 0.0000E+00 0.0000E+00 DISTRIBUTION COEFF. (M**3/KG): bndry 0.0000E+00 0.0000E+00 MOLECULAR DIFFUSIVITY OF WATER, BRINE, NCG, TRC, AND HELIUM THROUGH THE GASEOUS AND AQUEOUS PHASES, (FDDIAG(PHASE,COMP)) [M**2/S]: PHASE 1 = GAS; PHASE 2 = AQUEOUS PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP PHASE COMP -1- -1- -1- -2- -1- -3- -1- -4- -1- -5- -2- -1- -2- -2- -2- -3- -2- -4- -2- -5- 0.10000E-04 0.00000E+00 0.10000E-04 0.10000E-04 0.85536E-03 0.10000E-09 0.10000E-09 0.10000E-09 0.10000E-09 0.10000E-09 TEMPERATURE = 0.000000E+00 OUT OF RANGE IN SAT +++++++++ CANNOT FIND PARAMETERS AT ELEMENT *ELE11* XX(M) = 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 0.000000E+00 !!!!!! ERRONEOUS DATA INITIALIZATION !!!!!! CHECK INCON FORMAT AND ASSIGNMENTS OF NK AND NKIN. ---------- STOP EXECUTION ----------