Examples |
Table of contents |
Two interacting blast waves Radiative compact supernova remnant Hawley-Zabusky angled shock problem Supernova shock instability Galactic cooling flows Disk-planet interaction |
Two interacting blast waves
References: Woodward, P.R., Colella, P., 1984, JCP, 54, 115
|
Problem | |
---|---|
All zones are plotted in the case of AMR ; original postscript files
are also available
(ppm 6.45M,
amr 1.65M).
The movie presents temporal evolution of density and velocity on 5 levels of grids with frames taken each 5 time steps.
|
Movies | |||||
---|---|---|---|---|---|
truncation error/resolution | |||||
10-1 | 10-2 | 10-3 | 10-4 | 6400 | |
2.7 M | 2.7 M | 2.8 M | 3.0 M | 2.6 M |
Performance | ||
---|---|---|
IBM RS/6000 590 | ||
SGI Indigo2 R4400 | ||
SGI Origin200 R10000 |
Radiative compact supernova remnant
References: Plewa, T., 1995, MNRAS, 275, 143
|
Problem | |
---|---|
The evolution of compact supernova remnant (cSNR) is one of the most
challenging problems of astrophysical hydrodynamics due to fact that
the cooling time scale becomes much shorter than the dynamical time
scale as the remnant enters its radiative phase of evolution. The
cooled material condensates very rapidly and the thin
shell forms behind the forward shock. At the end of this
simulation density of
the thin shell is nearly 4 orders of magnitude higher than the
density of ambient gas. Another consequence of the catastrophic
cooling is the formation of secondary shocks which modify the flow
pattern near the shell. The cSNRs are one of the most luminuous objects in the Universe with peak luminosities reaching 10+43 erg/s. |
Movies | ||||
---|---|---|---|---|
4 levels | 5 levels | |||
whole domain | 4.2 M | 4.0 M | ||
shocked gas | 3.9 M | 3.8 M |
Performance | ||
---|---|---|
IBM RS/6000 590 | ||
SGI Indigo2 R4400 | ||
SGI Origin200 R10000 |
Hawley-Zabusky angled shock problem
References: Hawley, J.F., Zabusky, N.J., 1989, Phys. Rev. Let., 63, 1241
|
Problem | |
---|---|
Initially Mach 1.2 shock is positioned at x=5 and contact discontinuity touches y=0 at x=15. The density jump across discontinuity is equal to 3 and its inclination is equal to 30 degrees. The gas is ideal and gamma=1.4. The evolution starts at t=0 and is followed up to t=620 inside the computational domain of size 62x372. Lower and upper boundaries are reflecting; at the left-hand boundary we imposed fixed inflow with the state set equal to post-shock values and the right-hand boundary is transmitting. |
Movies | ||||||
---|---|---|---|---|---|---|
truncation error/resolution | ||||||
10-1 | 10-2 | 10-3 | 960x160 | |||
density | 0.6 M | 0.6 M | 0.6 M | 0.6 M | ||
fluid tracer | 0.5 M | 0.5 M | 0.5 M | 0.5 M | ||
schlieren-type | 0.8 M | 0.7 M | 0.7 M | 0.8 M | ||
grid levels | 0.8 M | 0.7 M | 0.7 M |
Performance | ||
---|---|---|
SGI Origin200 R10000 |
Supernova shock instability
References: Müller, E., Fryxell, B., Arnett, D., 1991, A&A, 251, 505
|
Problem | |
---|---|
We consider evolution of strongly perturbed supernova in two dimensions assuming spherical symmetry. We start 300 s after the core bounce and continue calculations to t=13000 s. Development of Rayleigh-Taylor instabilities and subsequent large-scale mixing is observed and believed as being responsible for some spectral features of SN 1987A and most probably common for other Type-II supernovae. |
Movies | ||||
---|---|---|---|---|
768x200 | 3072x400 | |||
density | 0.6 M | 0.6 M | ||
Lagrangean coordinate | 0.5 M | 0.5 M | ||
schlieren-type | 0.8 M | 0.7 M | ||
grid levels | 0.8 M | 0.7 M |
Performance | ||
---|---|---|
SGI Indigo2 R4400 | ||
SGI Origin200 R10000 |
Galactic cooling flows
References: Kritsuk, A., Plewa, T., Müller, E., 2000, astro-ph/0008017
|
Problem | |
---|---|
We use hydrodynamic simulations with adaptive grid refinement to study the dependence of hot gas flows in X-ray luminous giant elliptical galaxies on the efficiency of heat supply to the gas. We consider a number of potential heating mechanisms including Type Ia supernovae and sporadic nuclear activity of a central supermassive black hole. As a starting point for this research we use an equilibrium hydrostatic recycling model (Kritsuk 1996). We show that a compact cooling inflow develops, if the heating is slightly insufficient to counterbalance radiative cooling of the hot gas in the central few kiloparsecs. An excessive heating in the centre, instead, drives a convectively unstable outflow. We model the onset of the instability and a quasi-steady convective regime in the core of the galaxy in two-dimensions assuming axial symmetry. Provided the power of net energy supply in the core is not too high, the convection remains subsonic. The convective pattern is dominated by buoyancy driven large-scale mushroom-like structures. Unlike in the case of a cooling inflow, the X-ray surface brightness of an (on average) isentropic convective core does not display a sharp maximum at the centre. A hybrid model, which combines a subsonic peripheral cooling inflow with an inner convective core, appears to be stable. We also discuss observational implications of these results. |
Results | |
---|---|
Hydrodynamics of galactic cooling flows web page with lots of computer generated movies and still pictures. |
Disk-planet interaction
References: Ciecielag, P., Plewa, T., Rózyczka, M., 2000, Astron. Nachr., 321, 171
|
Problem | |
---|---|
A problem of mass flow in the immediate vicinity of a planet embedded in a protoplanetary disk is studied numerically in two dimensions. Large differences in temporal and spatial scales involved suggest that a specialized discretization method for solution of hydrodynamical equations may offer great savings in computational resources, and can make extensive parameter studies feasible. Preliminary results obtained with help of Adaptive Mesh Refinement technique and high-order explicit Eulerian solver are presented. This combination of numerical techniques appears to be an excellent tool which allows for direct simulations of mass flow in vicinity of the accretor at moderate computational cost. The present communication is focused on the structure of the outer part of the circumplanetary disk. We employ the multifluid option to the hydrodynamical solver to prove that the circumplanetary disk is composed of the gas transfered into it from the edges of the gap. |
Results | |
---|---|
AVS visualization |