An alle Mathematica-Kundigen, ich als Mathematica-Neuling versuche gerade, eine numerische Lösung zu einem System partieller Differentialgleichungen zu erhalten, bislang aber ohne Erfolg. Es handelt sich um die Beschreibung eines "unendlich langen" magnetischen Zylinders, dessen Längsachse sich vertikal zu einem sinusförmigen magnetischen Wechselfeld befindet. Der Zylinder wird je nach Parameter um seine Längsachse eine Bewegung ausführen, wobei dies Stillstand, Oszillation, Rotation und chaotische Bewegung beinhaltet. Die Funktion Phi[tau] beschreibt die Position des Zylinders im Bogenmaß, v[r,tau] das Strömungsfeld. Die beiden Gleichungen sind über eine Randbedingung miteinander verknüpft. Mit der Funktion NDSolve funktioniert das leider nicht. Das entsprechende notebook hänge ich an. Woran liegt das? Als Anfänger bin ich für jeden Hinweis dankbar! Noch einen Hinweise: Die Funktion v(r,tau) kann man "von Hand" mittels einer Fourier-Reihe lösen. Damit kommt man aber in der zweiten Gleichung auch nicht weiter, da der Ausdruck Sin[tau] Cos[Phi[tau]] nicht integrierbar ist. Mit freundlichen Grüßen, Christian Reichert
Zylinder.nb
Description: Mathematica Notebook document
********************************************************************** Dipl.-Ing. Christian Reichert Institut für Technische Chemie - Bereich Wasser- und Geotechnologie (ITC-WGT) Forschungszentrum Karlsruhe GmbH Postfach 3640 D-76021 Karlsruhe Tel: 07247/82-2984 Fax: 07247/82-3478 http://www.fzk.de/itc-wgt **********************************************************************
