Open AccessSimulation of stellar instabilities with vastly different time‐scales using domain decomposition
Author(s) -
Grott M.,
Chernigovski S.,
Glatzel W.
Publication year - 2003
Publication title -
monthly notices of the royal astronomical society
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 2.058
H-Index - 383
eISSN - 1365-2966
pISSN - 0035-8711
DOI - 10.1046/j.1365-8711.2003.06896.x
Subject(s) - physics , instability , shock wave , envelope (radar) , domain decomposition methods , shock (circulatory) , stars , time domain , term (time) , stability (learning theory) , decomposition , domain (mathematical analysis) , mechanics , statistical physics , scale (ratio) , astrophysics , mathematical analysis , aerospace engineering , medicine , ecology , radar , mathematics , quantum mechanics , finite element method , machine learning , biology , computer science , engineering , computer vision , thermodynamics
Strange mode instabilities in the envelopes of massive stars lead to shock waves, which can oscillate on a much shorter time‐scale than that associated with the primary instability. The phenomenon is studied by direct numerical simulation using a, with respect to time, implicit Lagrangian scheme, which allows for a variation of the dependent variables by several orders of magnitude. The time‐step for the simulation of the system is reduced appreciably by the shock oscillations and prevents its long‐term study. A procedure based on domain decomposition is proposed to surmount the difficulty of vastly different time‐scales in various regions of the stellar envelope and thus to enable the desired long‐term simulations. Criteria for domain decomposition are derived and the proper treatment of the resulting inner boundaries is discussed. Tests of the approach are presented and its viability is demonstrated by application to a model for the star P Cygni. In this investigation primarily the feasibility of domain decomposition for the problem considered is studied. We intend to use the results as the basis of an extension to two‐dimensional simulations.