Discussion :

[yann_m]

Bien le bonjour,

quelqu'un aurait-il eu vent d'une methode permettant la resolution par elements finis d'une EDP comportant un bilaplacien ?

La seule methode qui me vient en tete pour le moment est de calculer la derrivee troisieme de mon inconnue (ordre 3 en formulation faible) a partir des valeurs discretes de cettes inconnue sur les noeuds voisins.
Cependant j'imagine qu'a cet ordre, cette facon de faire est tres/trop sensible aux erreurs numeriques...

Si quelqu'un a une idee je l'en remercie par avance.

Bien cordialement,
Yann.

[FR119492]

Salut!

Lorsque tu appliques la méthode des éléments finis à une équation de Laplace ou de Poisson avec des mailles triangulaires, tu minimise une fonctionnelle qui fait intervenir les dérivées premières et tu interpoles linéairement dans chaque triangle, ce qui te conduit à résoudre un système linéaire dont les inconnues sont les valeurs de la fonction recherchée sur les sommets des triangles.

Alors, bien que je n'aie jamais eu à traiter ce problème, je te suggère d'essayer de faire une interpolation quadratique sur chaque triangle, les inconnues étant les valeurs de la fonction recherchée sur les sommets des triangles et au milieu de chacun de leurs côtés.

Sans garantie, mais bonne chance quand même.
Jean-Marc Blanc

[yann_m]

Salut,

merci beaucoup de m'avoir répondu.
Oui en effet une interpolation quadratique sera meilleure en terme d'erreur. Cependant, ma question porte plutôt sur la façon de développer le gradient d'un Laplacien en FEM...
Question bête en fait, mais hélas je n'ai rien vu de tel auparavant...

Bien cordialement,
Yann.

[Aleph69]

Bonjour,

tu ne doit pas avoir de dérivée d'ordre 3 dans ta formulation faible : il faut appliquer deux fois la formule de la divergence de Green. Ainsi, lorsque tu testes le gradient du laplacien, tu obtiens l'opposé du laplacien testé par un gradient, sommé avec un terme de bord que tu dois traiter à l'aide des conditions aux limites imposées à ton problème.

[yann_m]

Salut,
je ne comprends pas bien ta réponse, tu veux dire faire deux intégrations par parties successives afin de mettre l'intégrale volumique du bilaplacien sous forme faible ?
Dans ce cas, de la seconde intégration par partie nous obtenons deux intégrales surfaciques et une intégrale de volume, cette dernière composé de deux Laplaciens... L'une sur mon inconnue, l'autre sur la fonction d'interpolation.
Je ne suis pas un familier de cette méthode numérique, je crois que mon problème vient de là. Je ne vois pas bien comment programmer cette intégrale...
Merci beaucoup pour ta réponse en tout cas !
Yann.

[Aleph69]

Bonjour,

pour un problème unidimensionnel, appliquer deux fois la formule de la divergence de Green revient bien à intégrer deux fois par parties. L'approximation par éléments finis de ta formulation faible peut ensuite être effectuée de manière classique; je ne vois pas le point bloquant pour toi : il faudrait que tu développes.

Quelques remarques de vocabulaire si tu es amené à rédiger un document sur ton travail par la suite :
1. mettre une intégrale (quantité globale) sous forme faible ne veut rien dire, cela concerne uniquement des équations différentielles (quantités locales);
2. mieux vaut parler d'approximation d'une intégrale que de programmation d'une intégrale.

[yann_m]

Salut,
merci pour le vocabulaire, en effet c'est important et j'ai tendance a prendre quelques libertes avec, je m'en excuse...

Pour developper un peu voila ce qui me bloque :

Le peu que j'ai fait avec les elements finis concernait des EDPs comportant un laplacien, par exemple l'equation de Poisson. La formulation faible de cette equation (obtenue comme tu l'as fait remarquer grace a la formule de la divergence de Green), donne une integrale surfacique, et une integrale volumique. Cette derniere comportant un produit de deux Laplaciens. Celui de la variable recherchee, et celui de la fonction d'interpolation. Dans ce cas, pas de probleme.

Maintenant si je suis le meme raisonnement avec une EDP comportant un bilaplacien. Je realise la premiere integration par partie (formule de la divergence de Green), le premier terme est une integrale de surface faisant intervenir le produit d'un gradient d'ordre 3 avec une fonction test, et la seconde une integrale volumique faisant intervenir le produit d'un gradient d'ordre 3 avec le gradient de la fonction test.

A ce point pas de probleme, ce que tu m'as suggere si j'ai bien compris, c'est d'utiliser une seconde fois la formule de la divergence de Green sur l'integrale volumique.
Dans ce cas une seconde integrale de surface apparait, ainsi qu'une nouvelle integrale volumique, composee d'un produit de deux Laplaciens. Le Laplacien de mon inconnue et le Laplcien de ma fonction test.

C'est a ce point ou je suis bloque, comment faire pour approximer l'integrale de deux Laplaciens par elements finis ?

Encore merci de prendre de ton temps,
bien cordialement,
Yann.