################################################################################ TEST FOR LONGITUDIAL MECHANICAL DISPERSION OF TEMPERATURE The problem is not described in any paper. See a similar problem for solute transport in the following paper: Afanasyev, A.A. 2018. Numerical modelling of solute flow dispersion in porous media using simulator MUFITS // IOP Conf. Series: Journal of Physics: Conf. Series. 1129: 012002 DOI:10.1088/1742-6596/1129/1/012002. -------------------------------------------------------------------------------- Last update on 17.07.2019 Program Version: 2019.B01 ################################################################################ RUNSPEC ################### RUNSPEC section begins here ###################### GASSTORE ^^^^^^ We switch the program to the GASSTORE EOS module ^^^^^^^^^^^^^^ The fluid consists of WATER - water METRIC ^^^^^^^^^^^^^^^^^^^^^ We use METRIC units ^^^^^^^^^^^^^^^^^^^^^^^^^^^^ MECHDISP We enable the mechanical dispersion modelling option. FAST We use FAST option. AUXARRAY We create auxiliary array #XCOORD which #XCOORD / will be used for reporting x coordinate / fro GRID ##################### GRID section begins here ####################### The grid is specified within brackets MAKE-ENDMAKE MAKE <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< -- cartesian We select Cartesian gridding -- grid nx ny nz option and specify the number of CART 12000 1 1 / along every axis. XYZBOUND -- xmin-xmax ymin-ymax zmin-zmax We specify the domain extent in metres. -20 100 0 1 0 1 / The domain length is 120 metres. BOUNDARY -- ---- ----------- ------------ -------- -------- 102 6* 'I+' 5* INFTHIN 4* 1 2 / -- ---- ----------- ------------ -------- -------- / We create infinetely thin gost grid block at the right boundary. The parameters in this block are fixed (ACTNUM=2), will be used as open boundary. Further we can refer to this block with FLUXNUM=102. SRCSPECG -- name i j k MYSRC 1 1 1 / We create a point source with the name / MYSRC and place it in the grid block i=1, j=1 and k=1. ENDMAKE >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> EQUALS We equate PORO 0.1 / - porosity to 0.1; PERMX 100 / - X permeability to 100 mD; PERMY 100 / - Y permeability to 100 mD. / RPTGRID We reporrt these properties from the DEPTH PORO PERMX PERMZ / GRID section. PROPS ####################### PROPS section begins here #################### Rock properties are specified within brackets ROCK-ENDROCK ROCK <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / ROCKDH We specify that rock density and heat 2896.5 1.0 / capacity are 2896.5 kg/m3 and 1 kJ/kg/K. MDISPL This keyword specifies the longitudial -- phase -- sat -- dispersivity for the LIQ phase. The LIQ 0 1 / dispersivity is equal to the true fluid -- vel ------------- velocity and it does not depend on the 0 0 0 / phase saturation. 1 1 1 / / ENDROCK >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> The relative permeabilities are specified within brackets SAT-ENDSAT SAT <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< / ENDSAT >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> INIT ####################### INIT section begins here ##################### REGALL ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ The following OPERAREG keyword will be applied to all grid blocks. OPERAREG -- -------- -------- -------- -------- -------- -------- -------- -------- PRES 1* ROCKNUM 1 EQUALS 50 / TEMPC 1* ROCKNUM 1 EQUALS 30 / TEMPC XCOORD FLUXNUM 1 SETINT 31 -3000 0 / TEMPC 1* FLUXNUM 101 EQUALS 31 / TEMPC 1* FLUXNUM 102 EQUALS 30 / -- -------- -------- -------- -------- -------- -------- -------- -------- / With the above keyword we specfy that at t=0: - pressure = 50 bar in the whole domain; - temperature is equal to 31 C at x<0; - temperature is equal to 30 C at x>=0. EQUALNAM The temperature of the fluid injected throuh PRES 50 MYSRC / the point source is 31 C at 50 bar. TEMPC 31 / / COPY We copy x coordinate to the user created array XCOORD #XCOORD / #XCOORD. / RPTSUM PRES TEMPC PHST DENT DEN$1 ETM$1 #XCOORD / We report these parameters at every report time. SCHEDULE ##################################################################### REPORTS NOTHING / No reports to LOG-file. Linear solver setup: ILUTFILL - Ilut-fill parameter is 5; 5 / ILUTDROP - Ilut-drop parameter is 1e-4. 1E-4 / TUNING The maximum time step is 50 days. 1* 50 / SRCINJE Through the point source we inject MYSRC RATEIN 1* 100 1E-4 / fluid at a rate of 1e-4 m3/day. / TSTEP We advance simulation to 365699.4489 100*3656.994489 / days reporting results every 3656.994489 days. POST ######################## POST section begins here #################### CONVERT We convert the output to ParaView compatible format. MNEMPARA The change the output format of the TEMPC 1PE12.5 / TEMPC mnemonic. / RPTPOST The report list consits of two columns, #XCOORD TEMPC / x coordinate and temperature. POSTSPEC With these kewords we report the 1* TIME25.CSV 25 / temperature distribution at the 25th POSTBLK report time. 1750 1 1 / 2000 1 1 / 2250 1 1 / 2500 1 1 / 2750 1 1 / 3000 1 1 / 3250 1 1 / 3500 1 1 / 3750 1 1 / 4000 1 1 / 4250 1 1 / 4500 1 1 / 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / / POSTSPEC With these kewords we report the 1* TIME50.CSV 50 / temperature distribution at the 50th POSTBLK report time. 2750 1 1 / 3000 1 1 / 3250 1 1 / 3500 1 1 / 3750 1 1 / 4000 1 1 / 4250 1 1 / 4500 1 1 / 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / 6500 1 1 / 6750 1 1 / 7000 1 1 / 7250 1 1 / / POSTSPEC With these kewords we report the 1* TIME75.CSV 75 / temperature distribution at the 75th POSTBLK report time. 3750 1 1 / 4000 1 1 / 4250 1 1 / 4500 1 1 / 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / 6500 1 1 / 6750 1 1 / 7000 1 1 / 7250 1 1 / 7500 1 1 / 7750 1 1 / 8000 1 1 / 8250 1 1 / / POSTSPEC With these kewords we report the 1* TIME100.CSV 100 / temperature distribution at the 100th POSTBLK report time. 4750 1 1 / 5000 1 1 / 5250 1 1 / 5500 1 1 / 5750 1 1 / 6000 1 1 / 6250 1 1 / 6500 1 1 / 6750 1 1 / 7000 1 1 / 7250 1 1 / 7500 1 1 / 7750 1 1 / 8000 1 1 / 8250 1 1 / 8500 1 1 / 8750 1 1 / 9000 1 1 / 9250 1 1 / / END #######################################################################