Équations différentielles non linéaires (2nd ordre)

Un problème, une question, un nouveau théorème ?

Modérateurs : JeanN, Michel Quercia

Répondre
AchilleTalon
Messages : 31
Enregistré le : lun. juil. 30, 2018 8:37 am
Classe : Psi

Équations différentielles non linéaires (2nd ordre)

Message par AchilleTalon » sam. mai 18, 2019 9:43 pm

Bonjour,

je travaille actuellement sur la mission IKAROS : Un satellite utilisant la pression photonique (pression exercée par les photons du soleil) pour garder son orbite autour du soleil.

Physiquement, j'arrive cependant à un jeu d'équations que je n'arrive pas à résoudre :

https://image.noelshack.com/fichiers/20 ... ations.jpg

Il s'agit d'équations non linéaires couplées (je crois), si quelqu'un saurait les résoudre, ou encore mieux, pourrait me donner des indications pour tracer cela sur Python (Si j'arrive à obtenir le profil de l'angle théta et de la distance satellite-soleil, r, en fonction du temps, je pourrais sûrement me débrouiller pour tracer la trajectoire du satellite)
Je précise que les seules inconnues sont la distance r (satellite-soleil) et l'angle théta.

Tous les commentaires constructifs me seraient sûrement d'une grande aide alors merci d'avance aux matheux qui pourront m'aider !
Clem

matmeca_mcf1
Messages : 1414
Enregistré le : mar. févr. 13, 2018 10:22 am

Re: Équations différentielles non linéaires (2nd ordre)

Message par matmeca_mcf1 » sam. mai 18, 2019 10:00 pm

Avez-vous essayer des schémas numériques pour résoudre numériquement le système d'EDO?
Ancien ENS Cachan (maths) 1999--2003
Enseignant-Chercheur à l'Enseirb-Matmeca (Bordeaux INP)
Les opinions exprimées ci-dessus sont miennes et ne reflètent pas la position officielle de l'école dans laquelle j'enseigne.

AchilleTalon
Messages : 31
Enregistré le : lun. juil. 30, 2018 8:37 am
Classe : Psi

Re: Équations différentielles non linéaires (2nd ordre)

Message par AchilleTalon » sam. mai 18, 2019 10:13 pm

matmeca_mcf1 a écrit :
sam. mai 18, 2019 10:00 pm
Avez-vous essayer des schémas numériques pour résoudre numériquement le système d'EDO?
Je ne connais pas le système EDO, cependant, je connais les algorithmes d'Euler pour la résolution des équations différentielles linéaires d'ordre 1/2.
Auriez vous une idée d'algorithme Python pouvant venir à bout de ces équations ?

matmeca_mcf1
Messages : 1414
Enregistré le : mar. févr. 13, 2018 10:22 am

Re: Équations différentielles non linéaires (2nd ordre)

Message par matmeca_mcf1 » sam. mai 18, 2019 11:32 pm

EDO = équations différentielles ordinaires. Pour de la mécanique céleste, Euler va donner des résultats non satisfaisants. Avec Euler explicite, toutes les planètes quittent le système solaire. Avec Euler implicite, toutes les planètes s'écrasent dans le Soleil.

Essayez RK4: ce ne sera pas très satisfaisant pour des temps longs (plusieurs milliers de révolutions) mais cela restera raisonnablement précis aux temps courts. Pour de la mécanique céleste, il faudrait utiliser des méthodes symplectiques. Vous pouvez essayer Störmer-Verlet. Ce sera moins précis au temps court mais plus stable sur les temps longs.
Ancien ENS Cachan (maths) 1999--2003
Enseignant-Chercheur à l'Enseirb-Matmeca (Bordeaux INP)
Les opinions exprimées ci-dessus sont miennes et ne reflètent pas la position officielle de l'école dans laquelle j'enseigne.

AchilleTalon
Messages : 31
Enregistré le : lun. juil. 30, 2018 8:37 am
Classe : Psi

Re: Équations différentielles non linéaires (2nd ordre)

Message par AchilleTalon » dim. mai 19, 2019 11:31 am

matmeca_mcf1 a écrit :
sam. mai 18, 2019 11:32 pm
EDO = équations différentielles ordinaires. Pour de la mécanique céleste, Euler va donner des résultats non satisfaisants. Avec Euler explicite, toutes les planètes quittent le système solaire. Avec Euler implicite, toutes les planètes s'écrasent dans le Soleil.

Essayez RK4: ce ne sera pas très satisfaisant pour des temps longs (plusieurs milliers de révolutions) mais cela restera raisonnablement précis aux temps courts. Pour de la mécanique céleste, il faudrait utiliser des méthodes symplectiques. Vous pouvez essayer Störmer-Verlet. Ce sera moins précis au temps court mais plus stable sur les temps longs.
Merci beaucoup, il me semble effectivement que RK4 est adaptée, je cherche activement comment l'implémenter sur Maple si quelqu'un à des conseils (c'est hors programme CPGE) !
Edit : j'ai trouvé effectivement des choses intéressantes sur Stormer-Verlet, cependant, le fait que mes équations soient couplées me pose problème...

matmeca_mcf1
Messages : 1414
Enregistré le : mar. févr. 13, 2018 10:22 am

Re: Équations différentielles non linéaires (2nd ordre)

Message par matmeca_mcf1 » dim. mai 19, 2019 2:31 pm

Ne faites pas du numérique avec Maple. Maple est un logiciel de calcul symbolique. Faites le en python avec numpy, ou en octave ou en C/fortran.
Ancien ENS Cachan (maths) 1999--2003
Enseignant-Chercheur à l'Enseirb-Matmeca (Bordeaux INP)
Les opinions exprimées ci-dessus sont miennes et ne reflètent pas la position officielle de l'école dans laquelle j'enseigne.

Répondre

Qui est en ligne

Utilisateurs parcourant ce forum : Aucun utilisateur enregistré et 7 invités