Network model of time-dependent solute transport in the brain

This folder contains code needed to run the time dependent solute transport network model. Similar code for the fluid transport hydraulic resistance network model can be found at https://zenodo.org/records/6540366 for the steady state solute transport model can be found at https://zenodo.org/records/11385333 Scripts: Run_LaplaceSolver_Trail.m: A script that allows you to load the parameters from "param_1g.mat" and run the hydraulic network and the time-dependent solute transport code using Laplace transform. As written, this script will output a file containing the solution, make a plot of the average concentration in each segment at different time, and will print the median concentration in each vessel type (pial, penetrating, capillary, or parenchymal) Revised_Figure1.m: Comparison of the results from the numerical solver and the semi-analytical Laplace transform solver Toolboxes: The Laplace_Solver_Refined_Toolbox contains the solver for the semi-analytical solution of the network solute transport using Laplace Transform The Numerical_Solver_Toolbox contains the numerical solver for the network solute transport using spectral method %Parameter sets: "param_1g.mat": brain network hydraulic resistance parameters generated using the code for the fluid transport hydraulic resistance network model can be found at https://zenodo.org/records/6540366 %Data SteadyInjection_1_Dispersion_Laplace_vectorized.mat: Result from the Laplace Transform solver SteadyInjection_1_increasedflow_dt1.mat, SteadyInjection_1_increasedflow_dt5, SteadyInjection_1_increasedflow_dt10: Results from the numerical solver using different time step (dt)