Evaluation of an efficient etack-RLE clustering concept for dynamically adaptive grids

This is the author accepted manuscript. The final version is available from the Society for Industrial and Applied Mathematics via the DOI in this record.Abstract. One approach to tackle the challenge of efficient implementations for parallel PDE simulations on dynamically changing grids is the usage of space-filling curves (SFC). While SFC algorithms possess advantageous properties such as low memory requirements and close-to-optimal partitioning approaches with linear complexity, they require efficient communication strategies for keeping and utilizing the connectivity information, in particular for dynamically changing grids. Our approach is to use a sparse communication graph to store the connectivity information and to transfer data block-wise. This permits efficient generation of multiple partitions per memory context (denoted by clustering) which - in combination with a run-length encoding (RLE) - directly leads to elegant solutions for shared, distributed and hybrid parallelization and allows cluster-based optimizations. While previous work focused on specific aspects, we present in this paper an overall compact summary of the stack-RLE clustering approach completed by aspects on the vertex-based communication that ease up understanding the approach. The central contribution of this work is the proof of suitability of the stack-RLE clustering approach for an efficient realization of different, relevant building blocks of Scientific Computing methodology and real-life CSE applications: We show 95% strong scalability for small-scale scalability benchmarks on 512 cores and weak scalability of over 90% on 8192 cores for finite-volume solvers and changing grid structure in every time step; optimizations of simulation data backends by writer tasks; comparisons of analytical benchmarks to analyze the adaptivity criteria; and a Tsunami simulation as a representative real-world showcase of a wave propagation for our approach which reduces the overall workload by 95% for parallel fully-adaptive mesh refinement and, based on a comparison with SFC-ordered regular grid cells, reduces the computation time by a factor of 7.6 with improved results and a factor of 62.2 with results of similar accuracy of buoy station dataThis work was partly supported by the German Research Foundation (DFG) as part of the Transregional Collaborative Research Centre “Invasive Computing” (SFB/TR 89)

A Theoretical Light-Curve Model for the Recurrent Nova V394 Coronae Austrinae

A theoretical light curve for the 1987 outburst of V394 Coronae Austrinae (V394 CrA) is modeled to obtain various physical parameters of this recurrent nova. We then apply the same set of parametersto a quiescent phase and confirm that these parameters give a unified picture of the binary. The early visual light curve (1-10 days after the optical maximum) is well reproduced by a thermonuclear runaway model on a very massive WD close to the Chandrasekhar limit (1.37 +- 0.01 M_sun). The ensuing plateau phase (10-30 days) is also reproduced by the combination of a slightly irradiated MS and a fully irradiated flaring-up disk with a radius ~1.4 times the Roche lobe size. The best fit parameters are the WD mass 1.37 M_sun, the companion mass 1.5 M_sun (0.8-2.0 M_sun is acceptable), the inclination angle of the orbit i~65-68 degree, and the flaring-up rim ~0.30 times the disk radius. The envelope mass at the optical peak is estimated to be ~6 x 10^{-6} M_sun, which indicates an average mass accretion rate of 1.5 x 10^{-7} M_sun yr^{-1} during the quiescent phase between the 1949 and 1987 outbursts. In the quiescent phase, the observed light curve can be reproduced with a disk size of 0.7 times the Roche lobe size and a rather slim thickness of 0.05 times the accretion disk size at the rim. About 0.5 mag sinusoidal variation of the light curve requires the mass accretion rate higher than ~1.0 x 10^{-7} M_sun yr^{-1}, which is consistent with the above estimation from the 1987 outburst. These newly obtained quantities are exactly the same as those predicted in a new progenitor model of Type Ia supernovae.Comment: 9 pages including 4 figures, to appear in the Astrophysical Journal, Part

Invasive compute balancing for applications with shared and hybrid parallelization

This is the author manuscript. The final version is available from the publisher via the DOI in this record.Achieving high scalability with dynamically adaptive algorithms in high-performance computing (HPC) is a non-trivial task. The invasive paradigm using compute migration represents an efficient alternative to classical data migration approaches for such algorithms in HPC. We present a core-distribution scheduler which realizes the migration of computational power by distributing the cores depending on the requirements specified by one or more parallel program instances. We validate our approach with different benchmark suites for simulations with artificial workload as well as applications based on dynamically adaptive shallow water simulations, and investigate concurrently executed adaptivity parameter studies on realistic Tsunami simulations. The invasive approach results in significantly faster overall execution times and higher hardware utilization than alternative approaches. A dynamic resource management is therefore mandatory for a more efficient execution of scenarios similar to our simulations, e.g. several Tsunami simulations in urgent computing, to overcome strong scalability challenges in the area of HPC. The optimizations obtained by invasive migration of cores can be generalized to similar classes of algorithms with dynamic resource requirements.This work was supported by the German Research Foundation (DFG) as part of the Transregional Collaborative Research Centre ”Invasive Computing” (SFB/TR 89)

Optical Detection of Two Intermediate Mass Binary Pulsar Companions

We report the detection of probable optical counterparts for two Intermediate Mass Binary Pulsar (IMBP) systems, PSR J1528-3146 and PSR J1757-5322. Recent radio pulsar surveys have uncovered a handful of these systems with putative massive white dwarf companions, thought to have an evolutionary history different from that of the more numerous class of Low Mass Binary Pulsars (LMBPs) with He white dwarf companions. The study of IMBP companions via optical observations offers us several new diagnostics: the evolution of main sequence stars near the white-dwarf-neutron star boundary, the physics of white dwarfs close to the Chandrasekhar limit, and insights into the recycling process by which old pulsars are spun up to high rotation frequencies. We were unsuccessful in our attempt to detect optical counterparts of PSR J1141-6545, PSR J1157-5112, PSR J1435-6100, and PSR J1454-5846.Comment: 9 pages, 2 figures, accepted for publication in ApJ

A Precise Distance to IRAS 00420+5530 via H2O Maser Parallax with the VLBA

We have used the VLBA to measure the annual parallax of the H2O masers in the star-forming region IRAS 00420+5530. This measurement yields a direct distance estimate of 2.17 +/- 0.05 kpc (<3%), which disagrees substantially with the standard kinematic distance estimate of ~4.6 kpc (according to the rotation curve of Brand and Blitz 1993), as well as most of the broad range of distances (1.7-7.7 kpc) used in various astrophysical analyses in the literature. The 3-dimensional space velocity of IRAS 00420+5530 at this new, more accurate distance implies a substantial non-circular and anomalously slow Galactic orbit, consistent with similar observations of W3(OH) (Xu et al., 2006; Hachisuka et al. 2006), as well as line-of-sight velocity residuals in the rotation curve analysis of Brand and Blitz (1993). The Perseus spiral arm of the Galaxy is thus more than a factor of two closer than previously presumed, and exhibits motions substantially at odds with axisymmetric models of the rotating Galaxy.Comment: 33 pages, 12 figures; Accepted by ApJ (to appear March 2009

Mid-Infrared Imaging of NGC 6334 I

We present high-resolution (<0.5") mid-infrared Keck II images of individual sources in the central region of NGC 6334 I. We compare these images to images at a variety of other wavelengths from the near infrared to cm radio continuum and speculate on the nature of the NGC 6334 I sources. We assert that the cometary shape of the UCHII region here, NGC 6334 F, is due to a champagne-like flow from a source on the edge of a molecular clump and not a due to a bow shock caused by the supersonic motion of the UCHII region through the interstellar medium. The mid-infrared emission in concentrated into an arc of dust that define the boundary between the UCHII region and the molecular clump. This dust arc contains a majority of the masers in the region. We discuss the nature of the four near-infrared sources associated with IRS-I 1, and suggest that one of the sources, IRS1E, is responsible for the heating and ionizing of the UCHII region and the mid-infrared dust arc. Infrared source IRS-I 2, which has been thought to be a circumstellar disk associated with a linear distribution of methanol masers, is found not to be directly coincident with the masers and elongated at a much different position angle. IRS-I 3 is found to be a extended source of mid-infrared emission coming from a cluster of young dusty sources seen in the near-infrared.Comment: Accepted for publication by the Astrophysical Journal, 27 pages, 9 figure

An empirical temperature calibration for the Delta a photometric system. II. The A-type and mid F-type star

With the Delta a photometric system, it is possible to study very distant galactic and even extragalactic clusters with a high level of accuracy. This can be done with a classical color-magnitude diagram and appropriate isochrones. The new calibration presented in this paper is a powerful extension. For open clusters, the reddening is straightforward for an estimation via Isochrone fitting and is needed in order to calculate the reddening-free, temperature sensitive, index (g1-y)0. As a last step, the calibration can be applied to individual stars. Because no a-priori reddening-free photometric parameters are available for the investigated spectral range, we have applied the dereddening calibrations of the Stromgren uvbybeta system and compared them with extinction models for the Milky Way. As expected from the sample of bright stars, the extinction is negligible for almost all objects. As a next step, already established calibrations within the Stromgren uvbybeta, Geneva 7-color, and Johnson UBV systems were applied to a sample of 282 normal stars to derive a polynomial fit of the third degree for the averaged effective temperatures to the individual (g1-y)0 values with a mean of the error for the whole sample of Delta T(eff) is 134K, which is lower than the value in Paper I for hotter stars. No statistically significant effect of the rotational velocity on the precision of the calibration was found.Comment: 5 pages, 2 figures, accepted by A&

Polarimetry of an Intermediate-age Open Cluster: NGC 5617

We present polarimetric observations in the UBVRI bands of 72 stars located in the direction of the medium age open cluster NGC 5617. Our intention is to use polarimetry as a tool membership identification, by building on previous investigations intended mainly to determine the cluster's general characteristics rather than provide membership suitable for studies such as stellar content and metallicity, as well as study the characteristics of the dust lying between the Sun and the cluster. The obsevations were carried out using the five-channel photopolarimeter of the Torino Astronomical Observatory attached to the 2.15m telescope at the Complejo Astron\'omico El Leoncito (CASLEO; Argentina. We are able to add 32 stars to the list of members of NGC 5617, and review the situation for others listed in the literature. In particular, we find that five blue straggler stars in the region of the cluster are located behind the same dust as the member stars are and we confirm the membership of two red giants. The proposed polarimetric memberships are compared with those derived by photometric and kinematical methods, with excellent results. Among the observed stars, we identify 10 with intrinsic polarization in their light. NGC 5617 can be polarimetrically characterized with $P_{max}= 4.40$ and $\theta_{v}= 73^\circ.1$. The spread in polarization values for the stars observed in the direction of the cluster seems to be caused by the uneven distribution of dust in front of the cluster's face. Finally, we find that in the direction of the cluster, the interstellar medium is apparently free of dust, from the Sun's position up to the Carina-Sagittarius arm, where NGC 5617 seems to be located at its farthest border

A Near-Infrared (JHK) Survey of the Vicinity of the HII region NGC 7538: Evidence for a Young Embedded Cluster

We describe the results of two near infrared (K-band) imaging surveys and a three color (JHK) survey of the vicinity of NGC 7538. The limiting magnitudes are K ~ 16.5 and K ~ 17.5 mag for the K-band surveys and K ~ 15 mag for the JHK survey. We identify more than 2000 and 9000 near-infrared (NIR) sources on the images of the two K-band surveys and 786 NIR sources in the JHK survey. From color-color diagrams, we derive a reddening law for background stars and identify 238 stars with NIR excesses. Contour maps indicate a high density peak coincident with a concentration of stars with NIR excesses. We identify this peak as a young, embedded cluster and confirm this result with the K-band luminosity function, color histograms, and color-magnitude diagrams. The center of the cluster is at RA = 23:13:39.34, DEC = 61:29:18.9. The cluster radius is $\sim$ 3' ~ 2.5 pc for an adopted distance, d ~ 2.8 kpc. For d = 2.8 kpc, and reddening, E_{J-K} = 0.55 mag, the slope of the logarithmic K-band luminosity function (KLF) of the cluster, s ~ 0.32 +- 0.03, agrees well with previous results for L1630 (s = 0.34) and M17 (s = 0.26).Comment: 26 pages with 11 figures. Accepted by Astronomical Journa

The Temperature and Cooling Age of the White-Dwarf Companion to the Millisecond Pulsar PSR B1855+09

We report on Keck and {\em Hubble Space Telescope} observations of the binary millisecond pulsar PSR B1855+09. We detect its white-dwarf companion and measure \mv=25.90\pm0.12 and \mi=24.19\pm0.11 (Vega system). From the reddening-corrected color, (\mv-\mi)_0=1.06\pm0.21, we infer a temperature \Teff=4800\pm800 K. The white-dwarf mass is known accurately from measurements of the Shapiro delay of the pulsar signal, \Mcomp=0.258^{+0.028}_{-0.016} \Msun. Hence, given a cooling model, one can use the measured temperature to determine the cooling age. The main uncertainty in the cooling models for such low-mass white dwarfs is the amount of residual nuclear burning, which is set by the thickness of the hydrogen layer surrounding the helium core. From the properties of similar systems, it has been inferred that helium white dwarfs form with thick hydrogen layers, with mass \simgt3\times10^{-3} \Msun, which leads to significant additional heating. This is consistent with expectations from simple evolutionary models of the preceding binary evolution. For PSR B1855+09, though, such models lead to a cooling age of $\sim10$Gyr, which is twice the spin-down age of the pulsar. It could be that the spin-down age were incorrect, which would call the standard vacuum dipole braking model into question. For two other pulsar companions, however, ages well over 10 Gyr are inferred, indicating that the problem may lie with the cooling models. There is no age discrepancy for models in which the white dwarfs are formed with thinner hydrogen layers (\simlt3\times10^{-4} \Msun).Comment: 7 pages, 1 figure, aas4pp2.sty. Accepted for publication in ApJ