[Date Index]
[Thread Index]
[Author Index]
Re: Help me : Solve a simple PDE in Mathematica
*To*: mathgroup at smc.vnet.net
*Subject*: [mg110532] Re: Help me : Solve a simple PDE in Mathematica
*From*: schochet123 <schochet123 at gmail.com>
*Date*: Wed, 23 Jun 2010 01:54:17 -0400 (EDT)
*References*: <hvmvoo$pp1$1@smc.vnet.net> <hvn9ms$7iv$1@smc.vnet.net>
Please note the recent message with Subject: Help me : Solve a simple
PDE in Mathematica
by Oliver Ruebenkoenig from Wolfram
who points out his Package available at
http://library.wolfram.com/infocenter/Conferences/7549/
It uses a finite element method to discretize multi-dimensional pure
boundary value problems and then solves the resulting system of linear
equations
This will probably solve your problem.
It can also solve time-dependent problems by using NDSolve on the
discretized system.
Very nice job, Oliver
Steve
On Jun 22, 2:02 pm, thaihang le <thaihang... at gmail.com> wrote:
> Dear Steve,
>
> Thank you for your great help.
>
> With your answer, no i understand. Solving PDE is not as simple
> although we know the equation and boundary. The most important is how
> to imply it in Mathematica.
>
> Mathematica uses Method of lines and solve very well simple PDE. But
> as my problem, the boundary condition is split in two domains :
>
> At z= 0 :
> When r < 1 : The dirichlet boundary is used : u[r,z] =0
> When r > 1 : the Neumann is used : D[u[r,z],{r,1}] = 0
>
> That's why i dont know how to imply this in Mathematica. For example,
> if the boundary is the same like this :
>
> r < 1 : u[r,z]=0
> r> 1 : u[r,z] =1
>
> I can imply in Mathematica as : z=0 : u[r,0] == If[r > 1, 1, 0].
>
> Thanks again Steve.
> I try to go by another approach.
Prev by Date:
**Re: Difficulty with saving Package .m files**
Next by Date:
**Re: Book?**
Previous by thread:
**Re: Help me : Solve a simple PDE in Mathematica**
Next by thread:
**Re: Floating point exception when taking a reciprocal**
| |