Observations of Gaia16apd revealed extremely luminous ultraviolet emission among superluminous supernovae (SLSNe). Using radiation hydrodynamics simulations we perform a comparison of UV light curves, color temperatures and photospheric velocities between the most popular SLSN models: pair-instability supernova, magnetar and interaction with circumstellar medium. We find that the interaction model is the most promising to explain the extreme UV luminosity of Gaia16apd. The differences in late-time UV emission and in color evolution found between the models can be used to link an observed SLSN event to the most appropriate model. Observations at UV wavelengths can be used to clarify the nature of SLSNe and more attention should be paid to them in future follow-up observations.

We investigate the possibility of generating and studying turbulence in plasma by means of high-energy density laser-driven experiments. Our focus is to create supersonic, self-magnetized turbulence with characteristics that resemble those found in the interstellar medium (ISM). We consider a target made of a spherical core surrounded by a shell made of denser material. The shell is irradiated by a sequence of laser pulses sending inward-propagating shocks that convert the inner core into plasma and create turbulence. In the context of the evolution of the ISM, the shocks play the role of supernova remnant shocks and the core represents the ionized interstellar medium. We consider the effects of both pre-existing and self-generating magnetic fields and study the evolution of the system by means of two-dimensional numerical simulations. We find that the evolution of the turbulent core is generally, subsonic with rms-Mach number $M_t\approx 0.2$. We observe an isotropic, turbulent velocity field with an inertial range power spectra of $P(k)\propto k^{-2.3}$. We account for the effects of self-magnetization and find that the resulting magnetic field has characteristic strength $\approx 3\times 10^{4}$ G. The corresponding plasma beta is $\approx 1\times 10^{4}$--$1\times 10^{5}$, indicating that the magnetic field does not play an important role in the dynamical evolution of the system. The natural extension of this work is to study the system evolution in three-dimensions, with various laser drive configurations, and targets with shells and cores of different masses. The latter modification may help to increase the turbulent intensity and possibly create transonic turbulence. One of the key challenges is to obtain transonic turbulent conditions in a quasi-steady state environment.

(abridged) There exists a substantial disagreement between computer simulation results and high-energy density laboratory experiments of the Rayleigh-Taylor instability Kuranz et al. (2010). We adopt the Braginskii formulation for transport in hot, dense plasma, implement and verify the additional physics modules, and conduct a computational study of a single-mode RTI in two dimensions with various combinations of the newly implemented modules. We find that magnetic fields reach levels on the order of 11 MG in the absence of thermal conduction. We observe denting of the RT spike tip and generation of additional higher order modes as a result of these fields. Contrary to interpretation presented in earlier work Nishiguchi (2002), the additional mode is not generated due to modified anisotropic heat transport effects but due to dynamical effect of self-generated magnetic fields. The main effects of thermal conduction are a reduction of the RT instability growth rate (by about 20% for conditions considered here) and inhibited mixing on small scales. In this case, the maximum self-generated magnetic fields are weaker (approximately 1.7 MG). These self-generated magnetic fields are of very similar strength compared to magnetic fields observed recently in HED laboratory experiments Manuel et al. (2012). We find that thermal conduction plays the dominant role in the evolution of the model RTI system considered. It smears out small-scale structure and reduces the RTI growth rate. This may account for the relatively featureless RT spikes seen in experiments, but does not explain mass extensions observed in experiments. Resistivity and related heat source terms were not included in the present work, but we estimate their impact on RTI as modest and not affecting our main conclusions. Resistive effects will be discussed in detail in the next paper in the series.

We present a study of numerical behavior of a thickened flame used in Flame Capturing (FC, Khokhlov (1995)) for tracking thin unresolved physical flames in deflagration simulations. We develop a steady-state procedure for calibrating the flame model used, and test it against analytical results. We observe numerical noises generated by original realization of the technique. Alternative artificial burning rates are discussed, which produce acceptably quiet flames. Two new quiet models are calibrated to yield required "flame" speed and width, and further studied in 2D and 3D setting. Landau-Darrieus type instabilities of the flames are observed. One model also shows significantly anisotropic propagation speed on the grid, both effects increasingly pronounced at larger matter expansion as a result of burning; this makes the model unacceptable for use in type Ia supernova simulations. Another model looks promising for use in flame capturing at fuel to ash density ratio of order 3 and below. That "Model B" yields flames completely localized within a region 6 cells wide at any expansion. We study Markstein effect for flame models described, through direct numerical simulations and through quasi-steady technique developed. Comparison demonstrates that Markstein effect dominates instability effects on curved flame speed for Model B in 2D simulations for fuel to ash density ratio of about 2.5 and below. We find critical wavelength of LD instability by direct simulations of perturbed nearly planar flames; these agree with analytical estimates when Markstein number values found in this work are used. We conclude that the behavior of model B is well understood, and optimal for FC applications among all flame models proposed to date.

Various flame tracking techniques are often used in hydrodynamic simulations. Their use is indispensable when resolving actual scale of the flame is impossible. We show that parameters defining "artificial flame" propagation found from model systems may yield flame velocities several times distinct from the required ones, due to matter expansion being ignored in the models. Integral effect of material expansion due to burning is incorporated into flame capturing technique (FCT) [Khokhlov(1995)]. Interpolation formula is proposed for the parameters governing flame propagation yielding 0.2% accurate speed and width for any expansion (and at least 0.01% accurate for expansions typical in type Ia supernova explosions.) Several models with simple burning rates are studied with gas expansion included. Plausible performance of the technique in simulations is discussed. Its modification ensuring finite flame width is found. Implementation suggestions are summarized, main criterion being the scheme performance being insensitive to expansion parameter (thus absence of systematic errors when the burning progresses from inner to outer layers); in this direction promising realizations are found, leading to flame structure not changing while flame evolves through the whole range of densities in the white dwarf.