Solve Reaction-Diffusion Problem: This Problem Is Already Done, But If You Had To Start From Scratch You Would

Solve Reaction-Diffusion problem
solved using collocation at Gauss points (Legendre points) and Moments method
This problem is already done, but if you had to start from scratch you would:
1 Open visual basic editor and import OCCdll.bas and LUFactor.bas
2 Make sure Excel can find the dll by:
(a) putting it in the "current directory" (click on file open to determine the current directory)
(b) using an explicit path in the basic code
(c) putting it somewhere on the path
3 use OCC_Setup with supplied values for symmetry, type and geometry
4 The use supplied functions to get coefficients and solve the linear system
! Symmetry: NonSymmetric = 0, Symmetric = 1 0 Summary of Flux Calculation

! Type: Gauss=1, Lobatto=2 1 left right
! Geometry: Planar=0, Cylindrical=1, Spherical=2 0 0.05062 0.13368
OCC_Setup #VALUE! #VALUE! #VALUE!
Number Interior Points 4 #VALUE! #VALUE!
Number of Points Total 6 #VALUE! #VALUE!
0.04903 0.12654
set a to 1 for no variation, 0.2 for variation in paper #VALUE! #VALUE!
a error: #VALUE! #VALUE!
Spatial variation 0.2 3.2 #VALUE! #VALUE!
Thiele parameters 5 100 -3.15% -5.34%
#VALUE! #VALUE!
Conventional formulation
Xc Wc R(x) B array
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
Symmetric formulation
C array Eq. (1.10) (note: symmetric matrix)
A array Mmult(A,y) D mass matrix

#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
Moment method with interpolated rates
Transpose of rate coefficients from above
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! < Rx transposed
Dc*rate - interpolated rate terms (note: not symmetric) Solution matrix with interpolated rates
Moment method with exact integration of rate terms

Full exact Moments Mass matrix from below; Solution matrix
Calculation of Mass Matrix for Moments Method using higher Gaussian quadrature:
8 use 2 degrees higher to get exact integration
X W R W*R
#VALUE! #VALUE! #VALUE! #VALUE! L* - interp. polynomials thru interior pts L- interp. polynomials thru all pts
#VALUE! #VALUE! #VALUE! #VALUE!
W L* (WL*) transpose
WRL* (WRL*) transpose
x & y exact solution from large n for comparison

0 0.020032 0.0660995 0.1355657 0.2246803 0.328638 0.441834 0.558166 0.671362
0 0.097512 0.294848 0.5310374 0.742405 0.886061 0.958471 0.983554 0.981023
Calculation of interpolated solution Collocation Moments

x L1 L2 L3 L4 L5 L6 Ycol Yint
0.00 #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
press ctrl-alt-F9 if Excel fails to update when needed
1.20000
unfortunately, graph does not like to update, normal Microsoft
Note: to execute an array valued function
1. mark the area for the function result, NxN
where N is the total number of points
2. specify the function in the usual way 1.00000
he current directory) 3. use ctrl-shift-enter to execute the function,
otherwise only one box will be filled in
4. if only one box is filled in, mark the area for
the result, press F2 and then ctrl-shift-enter
0.80000
ry of Flux Calculation
total
0.18429 exact (from large n)
0.60000
#VALUE! Orthogonal Collocation first derivative
#VALUE! Orthogonal Collocation Eq (1.24)
#VALUE! Moments interpolated rates
0.17557 Moments, approx. quadrature (5 pt Gauss)
#VALUE! Moments, exact rate integration (6 pt Gauss) 0.40000
#VALUE! collocation
#VALUE! Moments interpolated rates
-4.73% Moments, approx. quadrature (5 pt Gauss)
#VALUE! Moments, exact rate integration
0.20000
Eq. (1.6) RHS Yc

#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! 0.00000
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! 0.0 0.1 0.2 0.3 0.
Exact
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
Moments Int
Moments
(note: symmetric matrix) RHS Yc Mmult Residual

#VALUE! #VALUE! #VALUE! #VALUE! 0.00000 #VALUE! #VALUE!
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
D mass matrix bug in this function, no first and last columns (calculation below)
#VALUE! #VALUE! #VALUE! #VALUE! #VALUE!
interpolated rates RHS Solution Mmult Residual
solution
RHS 0.0000 1st Deriv.
#VALUE! #VALUE! #VALUE!
#VALUE! #VALUE! #VALUE!
0.0000
L- interp. polynomials thru all pts

Pure mass matrix: (no reaction term)

Mass matrix exact integration row sum

0.77532 0.864434 0.933901 0.979968 1

0.944804 0.832326 0.584843 0.234696 0
Ymom
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
#VALUE!
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
Exact Collocation Collocation
Moments Int Moments Int Moments
Moments

Solve Reaction-Diffusion Problem: This Problem Is Already Done, But If You Had To Start From Scratch You Would

Uploaded by

Copyright:

Available Formats

Solve Reaction-Diffusion Problem: This Problem Is Already Done, But If You Had To Start From Scratch You Would

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Solve Reaction-Diffusion Problem: This Problem Is Already Done, But If You Had To Start From Scratch You Would

Uploaded by

Copyright:

Available Formats

Solve Reaction-Diffusion problem

! Symmetry: NonSymmetric = 0, Symmetric = 1 0 Summary of Flux Calculation

A array Mmult(A,y) D mass matrix

Moment method with exact integration of rate terms

x & y exact solution from large n for comparison

Calculation of interpolated solution Collocation Moments

Eq. (1.6) RHS Yc

(note: symmetric matrix) RHS Yc Mmult Residual

L- interp. polynomials thru all pts

Pure mass matrix: (no reaction term)

Mass matrix exact integration row sum

0.77532 0.864434 0.933901 0.979968 1

You might also like