RnD/reactor
3
0
Fork 0
Code Issues 1 Pull Requests Actions Packages Projects Releases Wiki Activity

Fix boundary conditions in gradient and second derivative operators for Reactor_PFR

Browse Source
This commit is contained in:
p.vanderwilt
2025-06-24 16:38:07 +02:00
parent 2e76f733a8
commit 9f13229785
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

8
dependencies/reactor_class.js vendored
View File

@@ -179,20 +179,20 @@ class Reactor_PFR {
}
makeDoperator() { // create the upwind scheme gradient operator
const I = math.diag(Array(this.n_x).fill(-1), 0);
const A = math.resize(math.diag(Array(this.n_x).fill(1), 1), [this.n_x, this.n_x]);
const I = math.diag(Array(this.n_x).fill(1), 0);
const A = math.resize(math.diag(Array(this.n_x).fill(-1), -1), [this.n_x, this.n_x]);
I[0][0] = 0;
I[0][1] = -1;
I[this.n_x-1][this.n_x-1] = 0; // Neumann boundary condition at x=L
return math.add(I, A);
}
makeD2operator() { // create the upwind scheme second derivative operator
makeD2operator() { // create the central second derivative operator
const I = math.diag(Array(this.n_x).fill(-2), 0);
const A = math.resize(math.diag(Array(this.n_x).fill(1), 1), [this.n_x, this.n_x]);
const B = math.resize(math.diag(Array(this.n_x).fill(1), -1), [this.n_x, this.n_x]);
I[0][0] = 0;
I[0][1] = -1;
I[0][0] = -1; // Dichelet boundary condition at outlet
return math.add(I, A, B);
}
}