https://github.com/ashleefv/fd_spheres_variable_diffusivity
Finite difference discretization scheme for diffusion in a sphere with variable diffusivity
Science Score: 23.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
○CITATION.cff file
-
○codemeta.json file
-
○.zenodo.json file
-
✓DOI references
Found 2 DOI reference(s) in README -
✓Academic publication links
Links to: zenodo.org -
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (5.3%) to scientific vocabulary
Repository
Finite difference discretization scheme for diffusion in a sphere with variable diffusivity
Basic Info
- Host: GitHub
- Owner: ashleefv
- License: bsd-3-clause
- Language: MATLAB
- Default Branch: master
- Size: 56.6 KB
Statistics
- Stars: 2
- Watchers: 3
- Forks: 1
- Open Issues: 0
- Releases: 2
Metadata Files
README.md
FDspheresvariable_diffusivity
Finite difference discretization scheme for diffusion in a sphere with variable diffusivity
Overview
This code solves the PDE for Fickian diffusion with variable diffusivity within a radially symmetric sphere. The PDE is solved numerically using method of lines, the spherical finite difference discretization scheme defined and analyzed in the reference, and any built-in MATLAB ODE solver. Examples are provided for variable diffusivity cases I-IV presented in the reference listed below. The use of ode45 solver is shown in this code. Any ODE solver could be selected.
Finite Difference Discretization Scheme for Spherical Diffusion with Variable Diffusivity
Code Author
Ashlee N. Ford Versypt, School of Chemical Engineering, Oklahoma State University ashleefv@okstate.edu
Reference
When using this software, please cite:
A. N. Ford Versypt & R. D. Braatz, Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity, Computers & Chemical Engineering, 71 (2014) 241-252, https://doi.org/10.1016/j.compchemeng.2014.05.022
Diffusivity cases (defined computationally in FDspheresvariable_diffusivity.m)
DCASE = 1; Case I: constant diffusivity with alpha = alpha0
DCASE = 2; Case II: time-dependent diffusivity alpha = alpha0 for t =< tau and alpha = alpha0exp(k(t - tau)) for t > tau
DCASE = 3; Case III: spatially dependent diffusivity alpha = alpha0*r^2
DCASE = 4; CASE IV: concentration-dependent diffusivity alpha = alpha0exp(k-kc)
Code
solveFDspheresvariablediffusivity.m solves the diffusion PDE in spherical coordinates for c(r,t) with variable diffusivity defined by DCASE
FDspheresvariable_diffusivity uses the finite difference discretization scheme defined in the reference as Scheme 1 and the variable diffusivity Cases I-IV analyzed in that work.
Owner
- Name: Ashlee Ford Versypt
- Login: ashleefv
- Kind: user
- Location: Buffalo, NY
- Company: University at Buffalo, The State University of New York
- Website: https://sites.google.com/site/ashleefordversypt/
- Twitter: FordVersyptLab
- Repositories: 1
- Profile: https://github.com/ashleefv