MATLAB Question: please write Tmat and b matrix for me. Background: Transmissibility Equation: T(-1/2) = d/dx [1/(K(i-1)-1/K(014-1 T (i+1/2) = d/dx [1/K(1) + 1K(i+111-1 Boundary condition: Assumption: in the term of TN+1/2(hn+hn+1), Ki+1 goes to infinity and hn+2 is equal to bbar. Thus, at the outlet boundary condition will be TN-1/2(hn-hn-2) - TN+1/2(hN+hN+1) TN 1/2*(hn-hn-2) -d/dx[1/K(l) + {i+2)14-1 -IN-1/2 * he 2+ IN-1/2+Tr N1+1/2 *hs = TEN+12hbac. MATLAB CODE: Domain size: L = 1000_%domain size [m] N = 100: _%number of cells dx = 4/N:_%cell size [m] x = linspace(dx/2.L-dx/2,N)'; %x-coordinates of cell centers % Boundary conditions Ubac = 1/24/60/60: %prescribed flux boundary (inlet boundary) [m/s] bbar = 100 %prescribed head boundary (outlet boundary) [m] % Fluid condition mu = 10^-3; % viscosity [kg/m/s] g=9.8; % gravitational acceleration (m/s^2) rho = 1000; % fluid density [kg/m^3] % Porous media k = 1e-13; % permeability [m^2] % --- Hydraulic conductivity on cell centers K = rho*g*k/mu; Kvec = K*ones(N, 1); Question 1: please write a matrix (4*4*14*1)=(4*1) for me Let's say (Tmat](hmat]=[b] - 1/2 Matrix [Tmat] | Ti+1/2 (T-1/2+Ti+1/2) -T23/2 0 T-1/2 -T23/2 T32-IN-1/2 0 -TN 1/2 TN-1+TrN+1/2 -TN-1/2 Matrix [hmat) hi IN- hos Matrix [b] Ubar 0 TrN1/2hbac Question 2: What if we have more grid cell (more than N=4) as stated in the Matlab.code. (N=100). Please write Tmat and b MATLAB Question: please write Tmat and b matrix for me. Background: Transmissibility Equation: T(-1/2) = d/dx [1/(K(i-1)-1/K(014-1 T (i+1/2) = d/dx [1/K(1) + 1K(i+111-1 Boundary condition: Assumption: in the term of TN+1/2(hn+hn+1), Ki+1 goes to infinity and hn+2 is equal to bbar. Thus, at the outlet boundary condition will be TN-1/2(hn-hn-2) - TN+1/2(hN+hN+1) TN 1/2*(hn-hn-2) -d/dx[1/K(l) + {i+2)14-1 -IN-1/2 * he 2+ IN-1/2+Tr N1+1/2 *hs = TEN+12hbac. MATLAB CODE: Domain size: L = 1000_%domain size [m] N = 100: _%number of cells dx = 4/N:_%cell size [m] x = linspace(dx/2.L-dx/2,N)'; %x-coordinates of cell centers % Boundary conditions Ubac = 1/24/60/60: %prescribed flux boundary (inlet boundary) [m/s] bbar = 100 %prescribed head boundary (outlet boundary) [m] % Fluid condition mu = 10^-3; % viscosity [kg/m/s] g=9.8; % gravitational acceleration (m/s^2) rho = 1000; % fluid density [kg/m^3] % Porous media k = 1e-13; % permeability [m^2] % --- Hydraulic conductivity on cell centers K = rho*g*k/mu; Kvec = K*ones(N, 1); Question 1: please write a matrix (4*4*14*1)=(4*1) for me Let's say (Tmat](hmat]=[b] - 1/2 Matrix [Tmat] | Ti+1/2 (T-1/2+Ti+1/2) -T23/2 0 T-1/2 -T23/2 T32-IN-1/2 0 -TN 1/2 TN-1+TrN+1/2 -TN-1/2 Matrix [hmat) hi IN- hos Matrix [b] Ubar 0 TrN1/2hbac Question 2: What if we have more grid cell (more than N=4) as stated in the Matlab.code. (N=100). Please write Tmat and b