Register to reply

Alternative boundary conditions - Thomas-algorithm

Share this thread:
Huibert
#1
Dec1-10, 10:26 AM
P: 6
Hello,

I have to solve a diffusion equation:
MatrixL * Csim(:,i+1) = MatrixR * Csim(:,i) + BoundaryConditions
where Csim = concentration, j = location, i = time.

Boundary conditions are of type Dirichlet (Csim = 5 at j = 1, Csim = 0 at j = end). So I used:

MatrixL(1,:) = (1 0 0 . . 0) 
MatrixR(1,:) = (0 0 0 . . 0) 
BoundaryConditions(1) = (5)
so
Matrices, first line: (1 0 . 0) * Csim(:,i+1) = (0 0 . 0) * Csim(:,i) + 5
Matrices, last line: (. . 0 1) * Csim(:,i+1) = (. . 0 1) * Csim(:,i) + 0

--------------------------------------------------------------------------------

To solve this problem with the Thomas Algorithm, I have to write the equation as

MatrixL2*Csim(:,i+1)  = Csim(:,i)
So MatrixL2 = inv(MatrixR)*MatrixL

However, it is not possible to calculate inv(MatrixR) when MatrixR(1,:)= MatrixR(end,:)= 0

--------------------------------------------------------------------------------

So I tried to describe the boundary conditions in the folloing way :

Matrices, first line: (1 0 0 . .) * Csim(:,i+1) = (1 0 0 . .) * Csim(:,i)
Matrices, last line: (. . 0 0 1) * Csim(:,i+1) = (. . 0 0 1) * Csim(:,i)

But that cat won't jump. So could you please help me to find what's wrong with this?

--------------------------------------------------------------------------------

Because the problem is a littlebit complicated to explain, a longer description is in the attachment.
Attached Files
File Type: pdf Question.pdf (119.9 KB, 1 views)
Phys.Org News Partner Science news on Phys.org
Experts defend operational earthquake forecasting, counter critiques
EU urged to convert TV frequencies to mobile broadband
Sierra Nevada freshwater runoff could drop 26 percent by 2100

Register to reply

Related Discussions
Boundary conditions Aerospace Engineering 2
Find Flux Density On One Side of Dielectric Boundary Given Boundary Conditions Advanced Physics Homework 0
Pde with boundary conditions Differential Equations 1
What's the difference between initial conditions and boundary conditions? Differential Equations 9
Boundary conditions. Calculus & Beyond Homework 4