Répondre à la discussion
Affichage des résultats 1 à 2 sur 2

probleme de condition limite sous mathematica



  1. #1
    azeit4

    Smile probleme de condition limite sous mathematica


    ------

    Bonjour
    Dans le cadre de mon projet de 3ème année d'école d'ingé, je dois effectuer une modélisation de l'échauffement d'une carotte d'argilite lors du carrotage.

    Pour cela, je demande à mathematica de résoudre une fonction de 2 variables (r et t).

    Je définis les conditions aux limites qui me paraissent nécessaires mais mathématica me revoit l'erreur suivante (en rouge).

    Merci d'avance de toute forme d'aide que vous pourrez m'apporter.


    A = D[Y[r, t], t] == k*(D[Y[r, t], {r, 2}] + (1/r)*D[Y[r, t], r])

    NDSolve[{A, Y[60, t] == 300, Y[r, 0] == 11, Y[r, 180] == 300}, Y, {r,
    19, 60}, {t, 0, 180}]


    NDSolve::ivone: Boundary values may only be specified for one \
    independent variable. Initial values may only be specified at one \
    value of the other independent variable

    -----

  2. #2
    Rincevent

    Re : probleme de condition limite sous mathematica

    salut,

    si je comprends bien tu cherches à résoudre une équation type diffusion donc parabolique...

    bien que n'ayant jamais joué avec mathematica pour résoudre des PDE, je crois comprendre un peu ce que tu fais (je connais quand même Nsolve) et voici quelques remarques/questions :

    - tu sembles imposer la condition initiale Y(r,t=0)=11 mais en même temps une condition "finale" Y(r,t=180)=300 : si c'est bien ce que tu fais, c'est effectivement interdit pour une équation parabolique

    - si je comprends bien, tu imposes aussi la condition limite Y(r=60,t)=300. Si c'est bien une condition limite que tu imposes, elle est incompatible avec ta condition initiale, ce qui peut poser problème... car tu demandes à avoir Y(60,0)=11 (car valable pour tout r) mais aussi Y(60,0)=300 (car valable pour tout t).
    Ceux qui manquent de courage ont toujours une philosophie pour le justifier. A.C.

Discussions similaires

  1. condition initiale sous proteus isis
    Par simonsays dans le forum Électronique
    Réponses: 4
    Dernier message: 24/01/2009, 17h29
  2. Vibration de poutres en flexion sous mathematica
    Par stephaninho dans le forum Technologies
    Réponses: 1
    Dernier message: 01/05/2008, 23h36
  3. Petit probleme avec mathematica
    Par Olorin dans le forum Logiciel - Software - Open Source
    Réponses: 0
    Dernier message: 24/01/2008, 15h09
  4. Hydrodynamique, problème de condition aux limites
    Par Alexander dans le forum Physique
    Réponses: 2
    Dernier message: 25/08/2007, 16h09
  5. problème de maximisation avec mathematica
    Par mmarie3 dans le forum Mathématiques du supérieur
    Réponses: 1
    Dernier message: 04/07/2007, 20h36