This model simulates the diffusion of a substance through a region with a constant diffusivity and different solubilities inside and outside the region.
Model number: 0176
|Run Model: ||    Help running a JSim model.|
(JSim model applet may take 10-20 seconds to load.)
where C is the concentration of the diffusing species, D is the rate of diffusion of the species in the axial direction, and x and t are the spatial and time domain respectively. The initial conditions are:
where Clh = 10 mM and Crh = 3 mM are the initial concentration at the left hand and right hand boundaries, respectively. To properly represent the variation of the concentration at the boundaries we must impose a concentration flux boundary condition at the left hand and right hand boundaries. Simply fixing the boundaries at the outside concentration divided by the partition coefficient will represent the solution properly in the steady state but will not accurately describe the concentrations close to the boundary at times significantly less than that to reach the steady state. We have used a central difference approximation to the flux at the boundaries establish the boundary conditions. More formally we have:
where Qlh and Qrh are the concentration fluxes at the left and right hand boundaries respectively and are given by the central difference approximation:
The equations for this model may also be viewed by running the JSim model applet and clicking on the Source tab at the bottom left of JSim's Run Time graphical user interface. The equations are written in JSim's Mathematical Modeling Language (MML). See the Introduction to MML and the MML Reference Manual. Additional documentation for MML can be found by using the search option at the Physiome home page.
Bassingthwaighte JB. Transport in Biological Systems, Springer Verlag, New York, 2007.
- Diffusion Tutorial,
- 1-D Diffusion modeled as a partial differential equation,
- 1-D Diffusion with asymmetrical Consumption modeled as a partial differential equation,
- 1-D diffusion-advection equation with Robin boundary condition
- Random Walks of multiple particles in 1 dimension
- Random Walk of single particle in 2 dimensions
- Fractional Brownian Motion Walk in 2 dimensions
- Diffusion in a uniform slab
- Two Slab diffusion: Different diffusion coeffs in adjacent slabs require special boundary conditions
- Heat equation in two dimensions with Dirichlet boundary conditions
- Safford 1977 Dead end pore model for Calcium diffusion in muscle
- Safford 1978 Water diffusion in heart
- Suenson 1974 Diffusion in heart tissue, sucrose and water
- Facilitated diffusion through 2 regions
- Barrer Diffusion: Diffusion through 1-D slab with recipient chamber on right
We welcome comments and feedback for this model. Please use the button below to send comments:
Model HistoryGet Model history in CVS.
Please cite www.physiome.org in any publication for which this software is used and send one reprint to the address given below:
The National Simulation Resource, Director J. B. Bassingthwaighte, Department of Bioengineering, University of Washington, Seattle WA 98195-5061.
[This page was last modified 14Mar18, 3:17 pm.]
Model development and archiving support at physiome.org provided by the following grants: NIH U01HL122199 Analyzing the Cardiac Power Grid, 09/15/2015 - 05/31/2020, NIH/NIBIB BE08407 Software Integration, JSim and SBW 6/1/09-5/31/13; NIH/NHLBI T15 HL88516-01 Modeling for Heart, Lung and Blood: From Cell to Organ, 4/1/07-3/31/11; NSF BES-0506477 Adaptive Multi-Scale Model Simulation, 8/15/05-7/31/08; NIH/NHLBI R01 HL073598 Core 3: 3D Imaging and Computer Modeling of the Respiratory Tract, 9/1/04-8/31/09; as well as prior support from NIH/NCRR P41 RR01243 Simulation Resource in Circulatory Mass Transport and Exchange, 12/1/1980-11/30/01 and NIH/NIBIB R01 EB001973 JSim: A Simulation Analysis Platform, 3/1/02-2/28/07.