Modify capillary pressure unit by unit
Hi, everyone
I have recently read a lot of coupling programs carried out based on the TOUGH series of software, and I noticed that they have an equation for modifying the capillary pressure, which means that they have corrected the capillary pressure for each cell based on porosity and permeability, may I ask where I can assign a cell-by-cell value to the capillary pressure?
(Rutqvist et al., 2002, Environmental Geology)
4 replies
-
Most people coupling flac3d with TOUGH at each time-step. For such a case, you need to update the capillary pressure at the beginning of each time step in time stepping loop. For example, in TOUGH3, you need update the parameters at subroutine CYCIT in file Main_Comp.f90 near the following lines:
do while (.not. stop_cycle)
if (KON .eq. 2) then
KCYC=KCYC+1
!
!-----COME HERE FOR NEW TIME STEP.
!
KC=KC+1! you may insert the parameters updating due to coupling effect here.