ff3d-users
[Top][All Lists]

## Re: [ff3d-users] Code to solve Fick's law, what's wrong??!!

 From: Stephane Del Pino Subject: Re: [ff3d-users] Code to solve Fick's law, what's wrong??!! Date: Sat, 15 May 2004 18:48:06 +0200 User-agent: KMail/1.6.2

```-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

Hello Arun.

Le samedi 15 Mai 2004 17:24, Arun Viswanathan a écrit :
> Hey Stephane,
>
> I am so sorry to be bothering you so often, the whole
> thing is new to me, that's why I have too many
> doubts....here's yet another...:)...
>
> In the example of the Euler's equation that you sent
> me, 'dt' is simply declared as any other variable. In
> my equation of Fick's law, I don't think I can do
> that. The equation I am trying to solve goes this way:
No. This example does not deal with Euler's equations but is an implicit Euler
*scheme* wrote to solve your problem.

As it is the case in every freefem version, functions do not depend on time,
only on space. But this is not a problem.

> d/dt [C(x,t)] = D * d^2/dx [C(x,t)]
In ff3d, to solve this problem the spacial discretization is done
independently. So, to solve
an implicit Euler scheme consists in writing
where n is the time step.
For convenience, this can be written:

This is exactly what solves the programm that I sent to you.

I hope that it is now more clear.
Best regards,
Stephane.

PS : please do not cross post your requests. This is not related with the -dev
mailing list and I read more often the -users mailing list that this
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.2.4 (GNU/Linux)

iD8DBQFApknHbZyhKhZwdc0RAsreAJ4qxgRJXQsfyCDaI6dCP05csRJaHQCgjQLG
TGMxaWbcVG4pi8q0iOEVySg=
=bOXW
-----END PGP SIGNATURE-----

```