Plot Page: Cout_t_Gisf_0
This is a plot of the outflow concentration as a function
of time. Gisf, consumption of the metabolite, is always zero.
As PSg (the exchange between the plasma and the isf) increases,
the peak decreases an appears at a later time. This happens
because as PSg increases, the model becomes effectively a
btex10 model (single capillary).
***********************************************************
Nested Plot: the NestedPlotPars should be loaded (reduces
t.max to 10.)
Nested Plot Cp(t,x,PSg,Gisf) is delineated with 25 plots.
Each plot sows the concentration in the capillary (time is
horizontal axis, down stream distance is the vertical axis.
The contour lines are for increasing the concentration in
the plasma as a function of time and distance in the capillary.
The axes for t and x are labeled X1 and Y1 respectively. The
axes for X2 and Y2 are PSg and Gisf respectively.
The first column has PSg=0 (exchange between plasma and isf).
Therefor changes in Gisf have no effect.
Nested Plot: Cout_t_PSg_Gisf are line plots of the outflow
concentration where PSG ranges from
PSG=0, black
PSg=3, red
PSg=6, orange
PSg=9, yellow
PSg=15, green.
An alternative way to do these plot would be to leave PSg and
Gisf as parameterss (not realDomains) and loop over them in the
Loops Graphical User interface.
// MODEL NUMBER:
// MODEL NUMBER: 0366
// MODEL NAME: NestedPlots
// SHORT DESCRIPTION: Models a tissue cylinder consisting of two regions: plasma,
// and interstitial fluid. Contains nested plots.
import nsrunit; unit conversion on;
math NestedPlots { // Based on BTEX20
/* PARAMETERS AND KEY TO NAMES
p = PLASMA
isf = INTERSTITIAL FLUID REGION*/
// INDEPENDENT VARIABLES
realDomain t sec ; t.min=0; t.max=30; t.delta=0.1;
realDomain x cm; real L=0.1 cm, Ngrid=61; x.min=0; x.max=L; x.ct=Ngrid;
private x.min, x.max, x.ct;
// ADDITIONAL DOMAINS
// PS means Permeability-surface area product between two regions
// between p and isf
realDomain PSg ml/(g*min); PSg.min=1; PSg.max=5; PSg.delta=1;
// Gisf is the consumption rate coefficient in each the isf(Gulosity)
realDomain Gisf ml/(g*min); Gisf.min=0; Gisf.max=10; Gisf.delta=2.5;
// OTHER PARAMETERS
real Fp = 1 ml/(g*min), // Plasma flow
// VOLUMES (a p after the region name implies "prime", a virtual volume)
Vp = 0.05 ml/g, // p
Visfp = 0.15 ml/g, // isf (Visf')
// G is the consumption rate coefficient in each region (Gulosity)
Gp = 0 ml/(g*min), // p
// D is the axial diffusion coefficient
Dp = 1.0e-5 cm^2/sec, // p
Disf = 1.0e-6 cm^2/sec; // isf
// The hVolumes protect against zero divides
private real hVp =if(Vp>0) Vp else (1e-6 ml/g);
private real hVisfp =if(Visfp>0) Visfp else (1e-6 ml/g);
// INFLOWING CONCENTRATION
extern real Cin(t) mM;
// CONCENTRATION VARIABLES
real Cp(t,x,PSg,Gisf) mM, // p
Cisf(t,x,PSg,Gisf) mM, // isf
Cout(t,PSg,Gisf) mM; // Outflow Concentration from plasma region
// BOUNDARY CONDITIONS (Note total flux BC for inflowing region.)
when (x=x.min) { (-Fp*L/hVp)*(Cp-Cin)+Dp*Cp:x = 0; Cisf:x = 0; }
when (x=x.max) { Cp:x = 0; Cisf:x = 0; Cout = Cp;}
// INITIAL CONDITIONS
when (t=t.min) { Cp =if(x=x.min) Cin else 0; Cisf = 0; }
// PARTIAL DIFFERENTIAL EQUATIONS
Cp:t = -Fp*L/hVp*Cp:x -Gp/hVp*Cp+ Dp*Cp:x:x
+ PSg/hVp*(Cisf-Cp);
Cisf:t = -Gisf/hVisfp*Cisf + Disf*Cisf:x:x
+ PSg/hVisfp*(Cp-Cisf) ;
}
/*
DETAILED DESCRIPTION:
This program is a version of BTEX20, the two region
partial differential model. Two of the model parameters,
PSg (exchange between plasma and isf) and Gisf
(consumption in the isf) have been modified to be
realDomains instead. This necessitates making
Cp(t,x) into Cp(t,x,PSg,Gisf),
Cisf(t,x) into Cisf(t,x,PSg,Gisf) and
Cout(t) into Cout(t,PSg,Gisf).
The details of making nested plots can be found at
http://www.physiome.org/jsim/docs/User_Nested.html
The "Nesting" menu is on the second line of the nesting plot
page, i.e.
Message | Plotname | ...
PlotName: | File | Nesting | View
KEY WORDS:
nested, nested plots, world within worlds, BTEX20,PDE,convection,diffusion,
permeation, reaction,distributed,capillary,plasma,isf,interstitial fluid
REFERENCES:
Bosan, Sorel and Harris, Thomas R.. A Visualizetion-Based Analysis Method
for Multiparameter Models of Capillary Tissue-Exchange. Annals of Biomedical
Engineering, Vol. 24, pp.124-138, 1996.
W.C. Sangren and C.W. Sheppard. A mathematical derivation of the
exchange of a labelled substance between a liquid flowing in a
vessel and an external compartment. Bull Math BioPhys, 15, 387-394,
1953.
C.A. Goresky, W.H. Ziegler, and G.G. Bach. Capillary exchange modeling:
Barrier-limited and flow-limited distribution. Circ Res 27: 739-764, 1970.
J.B. Bassingthwaighte. A concurrent flow model for extraction
during transcapillary passage. Circ Res 35:483-503, 1974.
B. Guller, T. Yipintsoi, A.L. Orvis, and J.B. Bassingthwaighte. Myocardial
sodium extraction at varied coronary flows in the dog: Estimation of
capillary permeability by residue and outflow detection. Circ Res 37: 359-378, 1975.
C.P. Rose, C.A. Goresky, and G.G. Bach. The capillary and
sarcolemmal barriers in the heart--an exploration of labelled water
permeability. Circ Res 41: 515, 1977.
J.B. Bassingthwaighte, C.Y. Wang, and I.S. Chan. Blood-tissue
exchange via transport and transformation by endothelial cells.
Circ. Res. 65:997-1020, 1989.
Poulain CA, Finlayson BA, Bassingthwaighte JB.,Efficient numerical methods
for nonlinear-facilitated transport and exchange in a blood-tissue exchange
unit, Ann Biomed Eng. 1997 May-Jun;25(3):547-64.
REVISION HISTORY:
Revised by BEJ 04/14/09
Boundary Conditions
Revised by GR 09/22/09
Added Statistics and Reformatted.
JSim SOFTWARE COPYRIGHT AND REQUEST FOR ACKNOWLEDGMENT OF USE:
JSim software was developed with support from NIH grants HL088516,
and HL073598. Please cite these grants in any publication for which
this software is used and send one reprint of published abstracts or
articles to the address given below. Academic use is unrestricted.
Software may be copied so long as this copyright notice is included.
Copyright (C) 1999-2009 University of Washington.
Contact Information:
The National Simulation Resource,
Director J. B. Bassingthwaighte,
Department of Bioengineering,
University of Washington, Seattle, WA
98195-5061
*/
NESTED PLOTS:
It is highly recommended that the user read
http://www.physiome.org/jsim/docs/User_Nested.html
before starting.
This program is a version of BTEX20, the two region
partial differential model. Two of the model parameters,
PSg (exchange between plasma and isf) and Gisf
(consumption in the isf) have been modified to be
realDomains instead. This necessitates making
Cp(t,x) into Cp(t,x,PSg,Gisf),
Cisf(t,x) into Cisf(t,x,PSg,Gisf) and
Cout(t) into Cout(t,PSg,Gisf).
The details of making nested plots can be found at
http://www.physiome.org/jsim/docs/User_Nested.html
The "Nesting" menu is on the second line of the nesting plot
page, i.e.
Message ]]| Plotname | ...
PlotName: | File | Nesting | View