Boundary condition problem for diffusion equation

Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
3 replies · 2K views
brambram
Messages
7
Reaction score
0
BOUNDARY CONDITION PROBLEM
I have came up with matrix for numerical solution for a problem where chemical is introduced to channel domain, concentration equation:

δc/δt=D*((δ^2c)/(δx^2))-kc

assuming boundary conditions for c(x,t) as : c(0,t)=1, c(a,t)=0. Where a is channel's length.

What if domain is not-infinite (as it is above) and we cannot assume that the chemical is not dispersed at the end of the channel- so the boundary condition c(a,t)=0 is no longer valid?

What boundary condition can we use when we cannot assume that the chemical is not dispersed at the outlet of the channel? Can you state it or not? How about computational fluid dynamics boundary conditions options?
 
Physics news on Phys.org
At the two ends of the sample, you can have a wide variety of different combinations of boundary conditions:

Specified concentration c = C(t)
Zero flux: δc/δx = 0
Specified flux: -D ∂c/∂x = [itex]\phi[/itex](t)
Convective mass transfer: -D ∂c/∂x = k (c - c)

There are an infinite number of others also, but these are the ones you run into most often in practice.
 
If there is no flux crossing a boundary, then the gradient of the flux normal to the boundary is 0.

For a 1-D problem like yours, that means ##\partial c/\partial x = 0##