[Next]:  Wissenschaftlich-technische Dienste  
 [Up]:  Projektbeschreibungen  
 [Previous]:  Bandstrukturberechnungen für Quantum Wells in Halbleiter-Heterostrukturen  
 [Contents]   [Index] 


  
Integrationsalgorithmen für die mehrdimensionalen Eulerschen Gleichungen

  Bearbeiter: W. Dreyer (FG 7), M. Kunik (FG 7), K. Sabelfeld (FG 6) 

Kooperation: N. Simonov (Sibirische Sektion der Russischen Akademie der Wissenschaften, Novosibirsk)

Beschreibung der Forschungsarbeit: Die von Dreyer & Kunik entwickelte Lösungstheorie hyperbolischer Systeme mit darunterliegender Mikrotheorie ([1], [2], [3]) erfordert die Auswertung von ein-, zwei- und dreidimensionalen Integralen. Hierzu sind insbesondere im Mehrdimensionalen effiziente Algorithmen notwendig, die wir gemeinsam mit Nikolai Simonov vom Instititut für Numerische Mathematik und Mathematische Geophysik in Novosibirsk erarbeiten. Begonnen haben wir diese Zusammenarbeit im Berichtszeitraum mit einer Studie über Integrale, die im Zusammenhang mit mehrdimensionalen Eulerschen Gleichungen auftreten. Die hier auftretenden Integranden sind ausnahmslos fast Gaußsche Integrationskerne, die wir durch geschickte Erweiterung auf exakt Gaußsche Kerne umformen konnten. Sodann haben wir zur Auswertung dieser Integrale die Methoden von Stroud/Secrest ([5]) sowie eigens hierzu entwickelte Monte-Carlo-Methoden angewendet.

In dem Spezialfall der zweidimensionalen Eulergleichungen ist zwar die Stroud/Secrest-Methode wesentlich effizienter als die Monte-Carlo-Methode, da sie mit ganz wenigen Gaußschen Integrationspunkten auskommt. Allerdings ist die Monte-Carlo-Methode im Gegensatz hierzu einer Übertragung auf höherdimensionale Probleme mit komplizierter Gebietsgeometrie unmittelbar zugänglich. Zusätzlich kann die Monte-Carlo-Methode durch eine Regularisierung der zufälligen Integrationsknotenpunkte bei bestimmten Anwendungen in ihrer Effizienz erheblich gesteigert werden. Diese Erweiterung wird zur Zeit gemeinsam mit den oben genannten Autoren erarbeitet und auf weitere Anwendungsbeispiele ausgedehnt.

Projektliteratur:

  1. W. DREYER, M. KUNIK, Maximum entropy principle revisited, WIAS-Preprint No. 367, 1997, Contin. Mech. Thermodyn., 10 (1998), No. 6, pp. 331-347.
  2. \dito 
, Reflections of Eulerian shock waves at moving adiabatic boundaries, WIAS-Preprint No. 383, 1997, Monte Carlo Methods and Appl., 4 (1998), No. 3, pp. 231-252.
  3. \dito 
, Initial and boundary value problems of hyperbolic heat conduction, WIAS-Preprint No. 438, 1998, erscheint in: Contin. Mech. Thermodyn.
  4. W. DREYER, M. KUNIK, K. SABELFELD, N. SIMONOV, K. WILMANSKI, Iterative procedure for multidimensional Euler equations, WIAS-Preprint No. 445, 1998, Monte Carlo Methods and Appl., 4 (1998), No. 3, pp. 253-271.
  5. A. H. STROUD, D. SECREST, Approximate integration formulas for certain spherical symmetric region, Math. Comp., 17 (1963), pp. 105-135.


 [Next]:  Wissenschaftlich-technische Dienste  
 [Up]:  Projektbeschreibungen  
 [Previous]:  Bandstrukturberechnungen für Quantum Wells in Halbleiter-Heterostrukturen  
 [Contents]   [Index] 

LaTeX typesetting by I. Bremer
7/30/1999