Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity

Two finite difference discretization schemes for approximating the spatial derivatives in the diffusion equation in spherical coordinates with variable diffusivity are presented and analyzed. The numerical solutions obtained by the discretization schemes are compared for five cases of the functional...

Full description

Bibliographic Details
Main Authors: Ford-Versypt, Ashlee Nicole, Braatz, Richard D
Other Authors: Massachusetts Institute of Technology. Department of Chemical Engineering
Format: Article
Language:en_US
Published: Elsevier 2017
Online Access:http://hdl.handle.net/1721.1/106943
https://orcid.org/0000-0003-4304-3484
_version_ 1811080980924989440
author Ford-Versypt, Ashlee Nicole
Braatz, Richard D
author2 Massachusetts Institute of Technology. Department of Chemical Engineering
author_facet Massachusetts Institute of Technology. Department of Chemical Engineering
Ford-Versypt, Ashlee Nicole
Braatz, Richard D
author_sort Ford-Versypt, Ashlee Nicole
collection MIT
description Two finite difference discretization schemes for approximating the spatial derivatives in the diffusion equation in spherical coordinates with variable diffusivity are presented and analyzed. The numerical solutions obtained by the discretization schemes are compared for five cases of the functional form for the variable diffusivity: (I) constant diffusivity, (II) temporally dependent diffusivity, (III) spatially dependent diffusivity, (IV) concentration-dependent diffusivity, and (V) implicitly defined, temporally and spatially dependent diffusivity. Although the schemes have similar agreement to known analytical or semi-analytical solutions in the first four cases, in the fifth case for the variable diffusivity, one scheme produces a stable, physically reasonable solution, while the other diverges. We recommend the adoption of the more accurate and stable of these finite difference discretization schemes to numerically approximate the spatial derivatives of the diffusion equation in spherical coordinates for any functional form of variable diffusivity, especially cases where the diffusivity is a function of position.
first_indexed 2024-09-23T11:39:54Z
format Article
id mit-1721.1/106943
institution Massachusetts Institute of Technology
language en_US
last_indexed 2024-09-23T11:39:54Z
publishDate 2017
publisher Elsevier
record_format dspace
spelling mit-1721.1/1069432022-10-01T05:08:04Z Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity Ford-Versypt, Ashlee Nicole Braatz, Richard D Massachusetts Institute of Technology. Department of Chemical Engineering Ford-Versypt, Ashlee Nicole Braatz, Richard D Two finite difference discretization schemes for approximating the spatial derivatives in the diffusion equation in spherical coordinates with variable diffusivity are presented and analyzed. The numerical solutions obtained by the discretization schemes are compared for five cases of the functional form for the variable diffusivity: (I) constant diffusivity, (II) temporally dependent diffusivity, (III) spatially dependent diffusivity, (IV) concentration-dependent diffusivity, and (V) implicitly defined, temporally and spatially dependent diffusivity. Although the schemes have similar agreement to known analytical or semi-analytical solutions in the first four cases, in the fifth case for the variable diffusivity, one scheme produces a stable, physically reasonable solution, while the other diverges. We recommend the adoption of the more accurate and stable of these finite difference discretization schemes to numerically approximate the spatial derivatives of the diffusion equation in spherical coordinates for any functional form of variable diffusivity, especially cases where the diffusivity is a function of position. 2017-02-15T16:49:16Z 2017-02-15T16:49:16Z 2014-08 2014-04 Article http://purl.org/eprint/type/JournalArticle 0098-1354 http://hdl.handle.net/1721.1/106943 Ford Versypt, Ashlee N., and Richard D. Braatz. “Analysis of Finite Difference Discretization Schemes for Diffusion in Spheres with Variable Diffusivity.” Computers & Chemical Engineering 71 (2014): 241–252. https://orcid.org/0000-0003-4304-3484 en_US http://dx.doi.org/10.1016/j.compchemeng.2014.05.022 Computers & Chemical Engineering Creative Commons Attribution-NonCommercial-NoDerivs License http://creativecommons.org/licenses/by-nc-nd/4.0/ application/pdf Elsevier PMC
spellingShingle Ford-Versypt, Ashlee Nicole
Braatz, Richard D
Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
title Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
title_full Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
title_fullStr Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
title_full_unstemmed Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
title_short Analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
title_sort analysis of finite difference discretization schemes for diffusion in spheres with variable diffusivity
url http://hdl.handle.net/1721.1/106943
https://orcid.org/0000-0003-4304-3484
work_keys_str_mv AT fordversyptashleenicole analysisoffinitedifferencediscretizationschemesfordiffusioninsphereswithvariablediffusivity
AT braatzrichardd analysisoffinitedifferencediscretizationschemesfordiffusioninsphereswithvariablediffusivity