@@@@@ @@ @ @ @@@ @ @ @@ @@@ @ @ @ @ @ @ @@ @@@@@ @ @@ @ @ @@@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @@ @@ @ @ @ @ @ @ @ @ @ @@ @ @ @ @ @ @@ @ @ @ @ @ @ @ @@ @@@@ @ @@ @ @ @ @ @ @ @ @@@@ @ @ @ @ @ @ @ @@@ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @ @@ @ @ @ @ @ @@ @ @@ @@ @@@ @ @ @@ @@@ @ @ @ @@ @@@@ @ @ @ @ @@ @ @ @ @ @@ @ @ TOUGH3 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 LABORATORY BY KARSTEN PRUESS ET AL. ******************************************************************************** ********************** TOUGH3 VER. 1.0.0 (January 2018) ********************* ******************************************************************************** SUMMARY OF DISK FILES 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 *TABLE* DOES NOT EXIST --- OPEN AS A NEW FILE =================================================================================================================================== PROBLEM TITLE: PINNs for CCS in 2D model (initial condition generator) (Reading ROCKS block) DOMAIN NO. 1 MATERIAL NAME -- Brick (Reading MULTI block) (Reading SELEC block) (Reading START block) HAVE READ UNKNOWN BLOCK LABEL "----*" --- IGNORE THIS, AND CONTINUE READING INPUT DATA (Reading PARAM block) HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA (Reading ELEME block) WRITE FILE *MESH* FROM INPUT DATA (Reading CONNE block) (Reading INCON block) WRITE FILE *INCON* FROM INPUT DATA (Reading OUTPU block) HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA (Reading TIMES block) HAVE READ UNKNOWN BLOCK LABEL " " --- IGNORE THIS, AND CONTINUE READING INPUT DATA (Reading ENDCY block) ************************************************************************************ * EVALUATE FLOATING POINT ARITHMETIC * ************************************************************************************ * * * FLOATING POINT PROCESSOR HAS APPROXIMATELY 15 SIGNIFICANT DIGITS * * * * DEFAULT VALUE OF INCREMENT FACTOR FOR NUMERICAL DERIVATIVES IS DFAC = 0.1490E-07 * * DEFAULT VALUE FOR DFAC WILL BE USED * * * ************************************************************************************ ALL NCON = 1 CONNECTIONS READ FROM FILE *MESH* REFERENCE KNOWN ELEMENTS PARAMETERS FOR FLEXIBLE DIMENSIONING OF MAJOR ARRAYS (MAIN PROGRAM) ARE AS FOLLOWS MNEL = 2 MNCON = 1 MNEQ = 3 MNK = 3 MNPH = 3 MNB = 6 MNOGN = 0 MGTAB = 0 =================================================================================================================================== MAXIMUM NUMBER OF VOLUME ELEMENTS (GRID BLOCKS): MNEL = 2 MAXIMUM NUMBER OF CONNECTIONS (INTERFACES): MNCON = 1 MAXIMUM LENGTH OF PRIMARY VARIABLE ARRAYS: MPRIM = 8 MAXIMUM NUMBER OF GENERATION ITEMS (SINKS/SOURCES): MNOGN = 0 MAXIMUM NUMBER OF TABULAR (TIME-DEPENDENT) GENERATION DATA: MGTAB = 0 LENGTH OF SECONDARY PARAMETER ARRAY: MNSEC = 232 MESH HAS 2 ELEMENTS AND 1 CONNECTIONS (INTERFACES) BETWEEN THEM GENER HAS 0 SINKS/SOURCES END OF PART 1 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) = 1 *** 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) = 0 *** CYCIT MOP(3) = 0 *** MULTI MOP(4) = 0 *** QU MOP(5) = 0 *** EOS MOP(6) = 0 *** LINEQ MOP(7) = 0 *** IF UNEQUAL ZERO, WILL GENERATE A PRINTOUT OF INPUT DATA CALCULATIONAL CHOICES OFFERED BY MOP ARE AS FOLLOWS: MOP(8) = 0 *** IF ISOT IS NEGATIVE, CHOOSES OPTION FOR REDUCING FRACTURE-MATRIX INTERFACE AREA. 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) = 2: K = C0+C1*T+C2*Sw+C3*POR MOP(11) = 4 *** 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) = 2 *** CHOOSES PROCEDURE FOR INTERPOLATING GENERATION RATES FROM A TIME TABLE. = 0: TRIPLE LINEAR INTERPOLATION. = 1: "STEP FUNCTION" OPTION. = 2: RIGOROUS STEP RATE OPTION. MOP(13) = 1 *** DEFINES CONTENT OF INCON AND SAVE FILE. = 0: STANDARD CONTENT. = 2: READS PARAMETERS OF HYSTERESIS MODEL FROM FILE INCON. MOP(15) = 0 *** ALLOWS TO SELECT A SEMI-ANALYTICAL HEAT EXCHANGE CALCULATION WITH CONFINING BEDS. = 0: NO SEMI-ANALYTICAL HEAT EXCHANGE. = 1: SEMI-ANALYTICAL LINEAR HEAT EXCHANGE ENGAGED. INITIAL TEMPERATURE OF CONFINING BEDS IS UNIFORM. = 2: SEMI-ANALYTICAL LINEAR HEAT EXCHANGE ENGAGED. INITIAL TEMPERATURE OF CONFINING LAYERS IS NOT UNIFORM. = 5: SEMI-ANALYTICAL RADIAL HEAT EXCHANGE WITH PROPERTIES GIVEN IN MATERIAL QLOSS = 6: SEMI-ANALYTICAL RADIAL HEAT EXCHANGE WITH DEPTH-DEPENDENT PROPERTIES (DEPTH, RADIUS, TEMPERATURE, CONDUCTIVITY, DENSITY, CAPACITY) PROVIDED ON FILE radqloss.dat MOP(16) = 2 *** PERMITS TO CHOOSE TIME STEP SELECTION OPTION = 0: AUTOMATIC TIME STEPPING BASED ON MAXIMUM CHANGE IN SATURATION. = 1: AUTOMATIC TIME STEPPING BASED ON NUMBER OF ITERATIONS NEEDED FOR CONVERGENCE. > 1: INCREASE TIME STEP BY AT LEAST A FACTOR 2, IF CONVERGENCE OCCURS IN .LE. MOP(16) ITERATIONS. MOP(17) = 0 *** HANDLES TIME STEPPING AFTER LINEAR EQUATION SOLVER FAILURE. = 0: NO TIME STEP REDUCTION DESPITE LINEAR EQUATION SOLUTION FAILURE. = 9: REDUCE TIME STEP AFTER LINEAR EQUATION SOLUTION FAILURE. 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(21) = 4 *** PERMITS TO SELECT LINEAR EQUATION SOLVER FROM PACKAGE < 2: DEFAULTS TO MOP(21) = 3 = 2: DSLUBC: BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 3: DSLUCS: BI-CONJUGATE GRADIENT SOLVER - LANCZOS TYPE; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 4: DSLUGM: GENERALIZED MINIMUM RESIDUAL CONJUGATE GRADIENTS; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 5: DLUSTB: STABILIZED BI-CONJUGATE GRADIENT SOLVER; PRECONDITIONER: INCOMPLETE LU FACTORIZATION = 6: LUBAND: DIRECT SOLVER USING LU DECOMPOSITION = 7: AZTEC: PARALLEL ITERATIVE SOLVER = 8: PETSc: PARALLEL ITERATIVE SOLVER 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 FOR EACH PHASE FLUX *********************************************************************************************************************************** *********************************************************************************************************************************** &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& &&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 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. FILE *CO2TAB* EXISTS --- OPEN AS AN OLD FILE ************************************************* * CO2 PROPERTIES READ FROM FILE CO2TAB * * NPK=127; NTK= 51 * * P(1) = 0.100000E+06; P(127) = 0.600000E+08 * * T(1) = 0.304000E+01; T( 51) = 0.103040E+03 * ************************************************* *********************************************************************************************************************************** * ECO2: EQUATION OF STATE FOR MIXTURES OF WATER/SODIUM CHLORIDE/CO2 * *********************************************************************************************************************************** OPTIONS SELECTED ARE: (NK,NEQ,NPH,NB) = (3,3,3, 6) NK = 3 - NUMBER OF FLUID COMPONENTS NEQ = 3 - NUMBER OF EQUATIONS PER GRID BLOCK NPH = 3 - NUMBER OF PHASES THAT CAN BE PRESENT NB = 6 - NUMBER OF SECONDARY PARAMETERS (OTHER THAN COMPONENT MASS FRACTIONS) FOR NB = 6, DIFFUSION IS "OFF", FOR NB = 8, DIFFUSION IS "ON" AVAILABLE OPTIONS: (NK,NEQ,NPH,NB)=(3,4,3,6 OR 8) -WATER, NACL, CO2; NON-ISOTHERMAL; VAR. (P, XS OR SS+10, X3 OR SG+10, T) (3,3,3,6 OR 8) -WATER, NACL, CO2; ISOTHERMAL; VAR. (P, XS OR SS+10, X3 OR SG+10, T) DEFAULT OPTIONS ARE (3,4,3,6) - NON-ISOTHERMAL, DIFFUSION "OFF" *********************************************************************************************************************************** THE PRIMARY VARIABLES ARE P - PRESSURE XSM - SALT MASS FRACTION IN TWO-COMPONENT H2O-NaCl (XS) OR SOLID SATURATION (SS+10.) X3 - CO2 MASS FRACTION OR SG+10. - GAS PHASE SATURATION + 10. T - TEMPERATURE ****************************** ********************************************************* * COMPONENTS * * FLUID PHASE CONDITION PRIMARY VARIABLES * ****************************** ********************************************************* * * * * * # 1 - WATER * * SINGLE-PHASE GAS P, XSM, X3, T * * * * * * # 2 - NACL * * SINGLE-PHASE AQUEOUS P, XSM, X3, T * * * * * * # 3 - CO2 * * TWO-PHASE P, XSM, SG+10., T * * * * * * # 4 - HEAT * ********************************************************* * * ****************************** *********************************************************************************************************************************** * ARRAY *IE* ALLOWS TO MAKE CHOICES AMONG DIFFERENT EOS OPTIONS * *********************************************************************************************************************************** IE(1) = 0: *** NUMBER OF ADDITIONAL RECORDS READ IN DATA BLOCK SELEC. IE(3) = 0: *** ALLOWS CHOICE OF OPTIONS FOR CALCULATION OF BRINE VISCOSITY. = 0: PHILLIPS ET AL. (1981) (DEFAULT). = 1: PALLISER AND MCKIBBIN (1998). = 2: MAO AND SUN (2006) = 3: POTTER (1976). IE(4) = 1: *** ALLOWS TO CHOOSE AMONG DIFFERENT CORRELATIONS FOR THE DENSITY OF COMPRESSED BRINE. = 1: ANDERSEN ET AL. (1992). (DEFAULT=0) = 2: PRITCHETT (1993). = 3: BRINE COMPRESSIBILITY EQUAL TO WATER COMPRESSIBILITY AT THE SAME REDUCED TEMPERATURE. = 4: BRINE COMPRESSIBILITY EQUAL TO WATER COMPRESSIBILITY AT THE SAME TEMPERATURE = 5: BATZLE AND WANG (1992). = 6: DRIESNER (2007). IE(11) = 0: *** ALLOWS CHOICE OF DEPENDENCE OF PERMEABILITY ON ACTIVE FLOW POROSITY. = 0: PERMEABILITY DOES NOT VARY WITH FLOW POROSITY. = 1: PERMEABILITY VARIES AS PHIF**FE(1). = 2: SERIES FRACTURE MODEL. = 3: SERIES TUBE MODEL. IE(12) = 0: *** ALLOWS CHOICE OF MODEL FOR WATER SOLUBILITY IN CO2. = 0: AFTER SPYCHER AND PRUESS (2005). = 1: EVAPORATION MODEL (WATER DENSITY IN GAS BASED ON DENSITY OF VAPOR IN EQUILIBRIUM WITH BRINE). IE(13) = 0: *** ALLOWS CHOICE OF DEPENDENCE OF BRINE DENSITY ON DISSOLVED CO2. = 0: BRINE DENSITY VARIES WITH DISSOLVED CO2, ACCORDING TO THE CORRELATION FOR CO2 MOLAR VOLUME IN SOLUTION FROM J. GARCIA (LBNL-49023). = 1: BRINE DENSITY IS INDEPENDENT OF CO2 CONCENTRATION. IE(14) = 3: *** ALLOWS CHOICE OF TREATMENT OF THERMOPHYSICAL PROPERTIES DEPENDENCE ON SALT CONTENT. = 0: FULL DEPENDENCE (DEFAULT). > 0: SUPPRESSION OF ALL THERMOPHYSICAL PROPERTIES DEPENDENCE OTHER THAN SALT SOLUBILITY. IE(15) = 4: *** ALLOWS CHOICE OF CORRELATION FOR VAPOR SATURATED BRINE ENTHALPY. = 1: MICHAELIDES (1981). = 2: MILLER (1978). (OBSOLETE) = 3: PHILLIPS ET AL. (1981). = 4: LORENZ, MARIC AND RIRSHL (2000). (DEFAULT=0) = 5: DRIESNER (2007). *********************************************************************************************************************************** ********** VOLUME- AND MASS-BALANCES ******************************************************************************************** ********** [KCYC,ITER] = [ 0, 0] ***** THE TIME IS 0.00000E+00 SECONDS, OR 0.00000E+00 DAYS Balances for entire system PHASES PRESENT COMPONENT MASS IN PLACE (kg) ********************************************************* ************************************************************************ PHASES * GAS AQUEOUS SOLID COMPONENTS * WATER SALT CO2 HEAT ********************************************************* ************************************************************************ * PHASES * VOL. (m^3) * 0.10000000E+03 0.30000000E+03 0.00000000E+00 GAS PHASE * 0.21981018E+03 0.00000000E+00 0.83971172E+05 0.47152579E+11 MASS (kg) * 0.84190982E+05 0.30107436E+06 0.00000000E+00 AQUEOUS * 0.29511747E+06 0.00000000E+00 0.59568939E+04 0.55805279E+11 ********************************************************* SOLID SALT * 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 ROCK * 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.18000000E+12 TOTAL * 0.29533728E+06 0.00000000E+00 0.89928066E+05 0.28295786E+12 ************************************************************************ ********************************************************************************************************************************** A11 1( 1, 1) ST = 0.100000E-01 DT = 0.100000E-01 A11 1( 2, 1) ST = 0.300000E-01 DT = 0.200000E-01 A11 1( 3, 1) ST = 0.700000E-01 DT = 0.400000E-01 A11 1( 4, 1) ST = 0.150000E+00 DT = 0.800000E-01 A11 1( 5, 1) ST = 0.310000E+00 DT = 0.160000E+00 A11 1( 6, 1) ST = 0.630000E+00 DT = 0.320000E+00 A11 1( 7, 1) ST = 0.127000E+01 DT = 0.640000E+00 ...ITERATING... AT [ 8, 1] --- DELTEX = 0.128000E+01 MAX. RES. = 0.135612E-04 AT ELEMENT A11 1 EQUATION 3 A11 1( 8, 2) ST = 0.255000E+01 DT = 0.128000E+01 ...ITERATING... AT [ 9, 1] --- DELTEX = 0.256000E+01 MAX. RES. = 0.270819E-04 AT ELEMENT A11 1 EQUATION 3 A11 1( 9, 2) ST = 0.511000E+01 DT = 0.256000E+01 ...ITERATING... AT [ 10, 1] --- DELTEX = 0.512000E+01 MAX. RES. = 0.540023E-04 AT ELEMENT A11 1 EQUATION 3 A11 1( 10, 2) ST = 0.102300E+02 DT = 0.512000E+01 ...ITERATING... AT [ 11, 1] --- DELTEX = 0.102400E+02 MAX. RES. = 0.107362E-03 AT ELEMENT A11 1 EQUATION 3 A11 1( 11, 2) ST = 0.204700E+02 DT = 0.102400E+02 ...ITERATING... AT [ 12, 1] --- DELTEX = 0.204800E+02 MAX. RES. = 0.212187E-03 AT ELEMENT A11 1 EQUATION 3 A11 1( 12, 2) ST = 0.409500E+02 DT = 0.204800E+02 ...ITERATING... AT [ 13, 1] --- DELTEX = 0.409600E+02 MAX. RES. = 0.414464E-03 AT ELEMENT A11 1 EQUATION 3 A11 1( 13, 2) ST = 0.819100E+02 DT = 0.409600E+02 ...ITERATING... AT [ 14, 1] --- DELTEX = 0.819200E+02 MAX. RES. = 0.791132E-03 AT ELEMENT A11 1 EQUATION 3 A11 1( 14, 2) ST = 0.163830E+03 DT = 0.819200E+02 ...ITERATING... AT [ 15, 1] --- DELTEX = 0.163840E+03 MAX. RES. = 0.144449E-02 AT ELEMENT A11 1 EQUATION 3 A11 1( 15, 2) ST = 0.327670E+03 DT = 0.163840E+03 ...ITERATING... AT [ 16, 1] --- DELTEX = 0.327680E+03 MAX. RES. = 0.242742E-02 AT ELEMENT A11 1 EQUATION 3 A11 1( 16, 2) ST = 0.655350E+03 DT = 0.327680E+03 ...ITERATING... AT [ 17, 1] --- DELTEX = 0.655360E+03 MAX. RES. = 0.352230E-02 AT ELEMENT A11 1 EQUATION 3 A11 1( 17, 2) ST = 0.131071E+04 DT = 0.655360E+03 ...ITERATING... AT [ 18, 1] --- DELTEX = 0.131072E+04 MAX. RES. = 0.402245E-02 AT ELEMENT A11 1 EQUATION 3 A11 1( 18, 2) ST = 0.262143E+04 DT = 0.131072E+04 ...ITERATING... AT [ 19, 1] --- DELTEX = 0.262144E+04 MAX. RES. = 0.322772E-02 AT ELEMENT A11 1 EQUATION 3 A11 1( 19, 2) ST = 0.524287E+04 DT = 0.262144E+04 ...ITERATING... AT [ 20, 1] --- DELTEX = 0.524288E+04 MAX. RES. = 0.162424E-02 AT ELEMENT A11 1 EQUATION 3 A11 1( 20, 2) ST = 0.104858E+05 DT = 0.524288E+04 ...ITERATING... AT [ 21, 1] --- DELTEX = 0.104858E+05 MAX. RES. = 0.467393E-03 AT ELEMENT A11 1 EQUATION 3 A11 1( 21, 2) ST = 0.209715E+05 DT = 0.104858E+05 ...ITERATING... AT [ 22, 1] --- DELTEX = 0.209715E+05 MAX. RES. = 0.723808E-04 AT ELEMENT A11 1 EQUATION 3 A11 1( 22, 2) ST = 0.419430E+05 DT = 0.209715E+05 A11 1( 23, 1) ST = 0.838861E+05 DT = 0.419430E+05 ...ITERATING... AT [ 24, 1] --- DELTEX = 0.838861E+05 MAX. RES. = 0.116543E-04 AT ELEMENT A11 1 EQUATION 3 A11 1( 24, 2) ST = 0.167772E+06 DT = 0.838861E+05 A11 1( 25, 1) ST = 0.335544E+06 DT = 0.167772E+06 A11 1( 26, 1) ST = 0.671089E+06 DT = 0.335544E+06 A11 1( 27, 1) ST = 0.864000E+06 DT = 0.192911E+06 PINNs for CCS in 2D model (initial condition generator) OUTPUT DATA AFTER ( 27, 1)-2-TIME STEPS THE TIME IS 0.100000E+02 DAYS 0@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ TOTAL TIME KCYC ITER ITERC KON DX1M DX2M DX3M MAX. RES. NER KER DELTEX 0.864000E+06 27 1 43 2 0.00000E+00 0.00000E+00 0.00000E+00 0.27808E-06 2 3 0.19291E+06 0@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ELEM. INDEX PRES_G PRES_L SAT_G SAT_L TEMP PCAP_GL REL_G REL_L X_WATER_G X_CO2_G X_WATER_L X_CO2_L VIS_G VIS_L DEN_G DEN_L (PA) (PA) (-) (-) (DEC-C) (PA) (-) (-) (-) (-) (-) (-) (PA*S) (PA*S) (KG/M**3) (KG/M**3) A11 0 1 0.22855E+08 0.2286E+08 0.5008E+00 0.4992E+00 0.4500E+02 0.0000E+00 0.4357E+00 0.8827E-02 0.2608E-02 0.9974E+00 0.9412E+00 0.5879E-01 0.7968E-04 0.5994E-03 0.8407E+03 0.1012E+04 A11 1 2 0.22855E+08 0.2286E+08 0.0000E+00 0.1000E+01 0.4500E+02 0.0000E+00 0.0000E+00 0.1000E+01 0.1000E+01 0.0000E+00 0.1000E+01 0.4596E-04 0.1000E+01 0.5994E-03 0.0000E+00 0.1000E+04 0@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ ********** VOLUME- AND MASS-BALANCES ******************************************************************************************** ********** [KCYC,ITER] = [ 27, 1] ***** THE TIME IS 0.86400E+06 SECONDS, OR 0.10000E+02 DAYS Balances for entire system PHASES PRESENT COMPONENT MASS IN PLACE (kg) ********************************************************* ************************************************************************ PHASES * GAS AQUEOUS SOLID COMPONENTS * WATER SALT CO2 HEAT ********************************************************* ************************************************************************ * PHASES * VOL. (m^3) * 0.10015277E+03 0.29984723E+03 0.00000000E+00 GAS PHASE * 0.21958333E+03 0.00000000E+00 0.83976590E+05 0.47178755E+11 MASS (kg) * 0.84196174E+05 0.30106910E+06 0.00000000E+00 AQUEOUS * 0.29511764E+06 0.00000000E+00 0.59514569E+04 0.55759108E+11 ********************************************************* SOLID SALT * 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 ROCK * 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.18000000E+12 TOTAL * 0.29533722E+06 0.00000000E+00 0.89928047E+05 0.28293786E+12 ************************************************************************ ********************************************************************************************************************************** TOUGH STATUS: Write SAVE file for time = 864000.00000000000 EEE EEE Number of processors = 1 EEE Time perform model computation = 0.0000000000000000 EEE of which spent in lin. solv. = 0.0000000000000000 EEE and spent on other = 0.0000000000000000 EEE EEE Total number of time steps = 27 EEE Average time in solver per time step = 0.0000000000000000 EEE Average time spent on other per time step = 0.0000000000000000 EEE EEE Total number Newton steps = 16 EEE Average number of Newton steps per time step 0.59259259259259256 EEE Average time per Newton step = 0.0000000000000000 EEE Average time spent on other per Newton st = 0.0000000000000000 EEE EEE Total number of iter = 16 EEE Average number of iter per call 1.0000000000000000 EEE Average time per iter = 0.0000000000000000 EEE EEE ============================================= EEE EEE END OF TOUGH3 SIMULATION RUN ELAPSED TIME = 0.000 SEC DISTRIBUTION + CALCULATION TIME = 0.000 SEC DATA INPUT TIME = 0.000 SEC