Main Page | Report this Page
 
   
Science Forum Index  »  Math - Numerical Analysis Forum  »  Solution to diffusion eqtn in a disk...
Page 1 of 1    
Author Message
Gib Bogle...
Posted: Wed May 14, 2008 2:52 am
Guest
The 2D diffusion equation is to be solved within the problem domain: a
disk of radius a. There are N constant point sources of mass within the
disk, and Dirichlet b.c. at the edge r = a, given by the function
f(theta) = U(a,theta).

Am I correct in thinking that the steady state solution for U(r,theta)
can be computed as the sum of the N Green's functions corresponding to
the sources (with the b.c. U(a,theta) = 0) and the Poisson's integral
formula corresponding to the boundary function f?

Thanks.
Peter Spellucci...
Posted: Wed May 14, 2008 2:52 am
Guest
In article <g0e5md$bko$1 at (no spam) lust.ihug.co.nz>,
Gib Bogle <bogle at (no spam) ihug.too.much.spam.co.nz> writes:
Quote:
The 2D diffusion equation is to be solved within the problem domain: a
disk of radius a. There are N constant point sources of mass within the
disk, and Dirichlet b.c. at the edge r = a, given by the function
f(theta) = U(a,theta).

Am I correct in thinking that the steady state solution for U(r,theta)
can be computed as the sum of the N Green's functions corresponding to
the sources (with the b.c. U(a,theta) = 0) and the Poisson's integral
formula corresponding to the boundary function f?

Thanks.
yes:


Delta u(i) = f(i) on D , u=0 on \partial D i=1,..,N
f(i) representing the source terms
Delta u(N+1) = 0 on D , u=f on \partial D

let
v= sum_{i=1,...,N+1} u(i)
=>
Delta v = sum_{i=1,..,N} f(i) on D, v=f on \partial D


hth
peter
Gib Bogle...
Posted: Wed May 14, 2008 5:01 pm
Guest
Peter Spellucci wrote:
Quote:
In article <g0e5md$bko$1 at (no spam) lust.ihug.co.nz>,
Gib Bogle <bogle at (no spam) ihug.too.much.spam.co.nz> writes:
The 2D diffusion equation is to be solved within the problem domain: a
disk of radius a. There are N constant point sources of mass within the
disk, and Dirichlet b.c. at the edge r = a, given by the function
f(theta) = U(a,theta).

Am I correct in thinking that the steady state solution for U(r,theta)
can be computed as the sum of the N Green's functions corresponding to
the sources (with the b.c. U(a,theta) = 0) and the Poisson's integral
formula corresponding to the boundary function f?

Thanks.
yes:

Delta u(i) = f(i) on D , u=0 on \partial D i=1,..,N
f(i) representing the source terms
Delta u(N+1) = 0 on D , u=f on \partial D

let
v= sum_{i=1,...,N+1} u(i)
=
Delta v = sum_{i=1,..,N} f(i) on D, v=f on \partial D


hth
peter


Thanks Peter, I was fairly sure, just wanted confirmation.
 
Page 1 of 1       All times are GMT - 5 Hours
The time now is Tue Oct 14, 2008 10:11 am