Post shut-in arrest and recession solutions for a deflating hydraulic fracture in a permeable elastic medium

The file names embed the two dimensionless parameters: omega=Ts/tmm~=Ts (since tmm~=1) and phiV. To avoid decimal points in the file name I have multiplied the phiV value by 100. So the above data file is for the case omega=10^{-3} and phiV=0.5.

Input = struct('Ep',Ep,... % Pa plane strain modulus 'mup',mup,... % Pa*s, alternate fluid viscosity 'Cp',Cp,... % m/s^1/2, alternate Carter's leak-off coefficient 'Kp',Kp,... % Pa*m^1/2, alternate fracture toughness 'Q0',Q0,... % m/s, injection rate 'Ts',Ts,... % shut-in time 'omega',omega,... % dimensionless shut-in time 'phiV',phiV,... % arrest regime parameter 'Nr',Nr,... % number of grid points in r direction 'Nt',itcol); % number of time steps till collapse Results=struct('pt',P(1,1:itcol),... % wellbore pressure versus t 'Rt',R(1:itcol),... % fracture radius versus t 'wt',W(1,1:itcol),... % wellbore aperture versus t 'eta',eta(1:itcol),... % efficiency versus t 'pr',P(:,1:itcol),... % fluid pressure versus r at all times Nt 'wr',W(:,1:itcol),... % fracture width versus r at all times Nt (because of the moving mesh to plot in real space use plot(rho*R(it),wr(:,it))) 'rho',rho,... % lateral spatial coordinate 't',time(1:itcol),... % time 'keyindx',[its ita itd itcol],... % key indices keyindx(1)=its (shut-in index), keyindx(2)=ita (arrest), keyindx(3)=itr (recession), keyindx(4)=itc (collapse) 'Input',Input); % Input Structure