DC conduction equation in Matlab

  • Context: MATLAB 
  • Thread starter Thread starter Mattw
  • Start date Start date
  • Tags Tags
    Conduction Dc Matlab
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
5 replies · 2K views
Mattw
Messages
8
Reaction score
0
Hello

I am solving DC conduction equation in Matlab. My equation is 1/r * ∂/∂r [r*σ* ∂V/∂r]=∂ρ/∂t
My equation is 1D(only varies with radius) and time dependent. i am using pdepe solver.The space charge ( ρ) is the time dependent variable. But the form of MATLAB pdepe solver the time dependent is the unknown term V. Can alyone help me how to solve this equation with pdepe or other solvers.

Thanks
 
Physics news on Phys.org
Good Morning Mattw,

So, a few questions:

Your unknown is [itex]V\left(t,r\right )[/itex]?
[itex]\rho \left (t,r\right)[/itex] is known for all [itex]t[/itex] and [itex]r[/itex]?

J.
 
Hi jfgobin

the space charge distribution ρ(t,r) is unknown.
 
Mattw, so you are trying to solve for both [itex]V[/itex] and [itex]\rho[/itex].

What is your second equation?
 
Or is is just Maxwell's?
 
Hi jfgobin

ρ(t,r) is just the continuti equation from Maxwell's.

I am solving this equation with matlab, and i am using pdepe solver. at the beginning if i assume my problem is in steady state then i can make ∂ρ/∂t zero. But the problem is i can't make "C"( coefficient of the time dependent term) zero in pdepe solve. Do you have any idea if there is MATLAB PDE solver for 1/r * ∂/∂r [r*σ* ∂V/∂r]=0 OR how to make "c" zero in pdepe solver?

Thanks