comp.lang.idl-pvwave archive: archive » How to implement periodic boundaries in IMSL_PDE

Home » Public Forums » archive » How to implement periodic boundaries in IMSL_PDE_MOL?

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

How to implement periodic boundaries in IMSL_PDE_MOL? [message #75437]

Wed, 09 March 2011 01:32

Tatcher
Messages: 7
Registered: September 2008

Junior Member

Hi!

I try to use IMSL_PDE_MOL to solve the diffusion equation Ut=Uxx to
become acquainted with the routine.

Initial values:
t=0
0<=x<=1
U(x,0)=gauss function

If I am using zero flux boundaries all is working smoothly und the
distribution decays nicely.
u(0,t)=u(1,t)=0
the values for 'f_bc' are then: alpha=1, beta=0,gammap=0

If I want to use periodic boundaries, I am lost. I want to do
something like this U[0,n]=U[X-2,n] ; U[N-1,n]=U[1,n] How do I get the
"ghost variables" of U in 'f_bc' and how to update alpha, beta and
gammap?

Cheers,

Christian

Report message to a moderator

[Message index]

How to implement periodic boundaries in IMSL_PDE_MOL?

By: Tatcher on Wed, 09 March 2011 01:32

Previous Topic:	Re: IDL 8.0 image() scaling
Next Topic:	Extract Array positions for a set of Values

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 03:33:19 PDT 2025

Total time taken to generate the page: 0.72046 seconds