How to implemented boundary condition in solving convection diffusion equation in cylindical coordinate
For C_c at r = R_c, you can consider that node as C(length(r1),j)
For C_a at r = R_c, you can consider that node as C(length(r1)+1,j)
Then, via central difference,
dC_c/dr (at r = R_c) can be written as (C(length(r1)+1,j)-C(length(r1)-1,j))/(del_r)
Similarly, dC_a/dr (at r = R_c) can be written as (C(length(r1)+2,j)-C(length(r1),j))/(del_r)
Try to suitably implement the above in the condition mass flux continuity condition at r = R_c, and check what happens. Also, ensure that the other conditions are correctly implemented, and then run the code.
I want to know how to write this in matlab so that I will not get error. I have aslo attached my matlab file hereFor C_c at r = R_c, you can consider that node as C(length(r1),j)
For C_a at r = R_c, you can consider that node as C(length(r1)+1,j)
Then, via central difference,
dC_c/dr (at r = R_c) can be written as (C(length(r1)+1,j)-C(length(r1)-1,j))/(del_r)
Similarly, dC_a/dr (at r = R_c) can be written as (C(length(r1)+2,j)-C(length(r1),j))/(del_r)
Try to suitably implement the above in the condition mass flux continuity condition at r = R_c, and check what happens. Also, ensure that the other conditions are correctly implemented, and then run the code.
I want to know how to write this in matlab so that I will not get error. I have aslo attached my matlab file here For C_c at r = R_c, you can consider that node as C(length(r1),j)
For C_a at r = R_c, you can consider that node as C(length(r1)+1,j)
Then, via central difference,
dC_c/dr (at r = R_c) can be written as (C(length(r1)+1,j)-C(length(r1)-1,j))/(del_r)
Similarly, dC_a/dr (at r = R_c) can be written as (C(length(r1)+2,j)-C(length(r1),j))/(del_r)
Try to suitably implement the above in the condition mass flux continuity condition at r = R_c, and check what happens. Also, ensure that the other conditions are correctly implemented, and then run the code.
I want to know how to write this in matlab so that I will not get error. I have aslo attached my matlab file here boundary condition, pde MATLAB Answers — New Questions
