Researcher profile

Tomasz Plewa

Tomasz Plewa contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.SR physics.flu-dyn physics.plasm-ph astro-ph physics.comp-ph astro-ph.HE astro-ph.EP astro-ph.GA physics.chem-ph physics.data-an

Trust snapshot

Quick read

Trust 21 - Emerging

9works

0followers

10topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

A hybrid adaptive multiresolution approach for the efficient simulation of reactive flows

Computational studies that use block-structured adaptive mesh refinement (AMR) approaches suffer from unnecessarily high mesh resolution in regions adjacent to important solution features. This deficiency limits the performance of AMR codes. In this work a novel hybrid adaptive multiresolution (HAMR) approach to AMR-based calculations is introduced to address this issue. The multiresolution (MR) smoothness indicators are used to identify regions of smoothness on the mesh where the computational cost of individual physics solvers may be decreased by replacing direct calculations with interpolation. We suggest an approach to balance the errors due to the adaptive discretization and the interpolation of physics quantities such that the overall accuracy of the HAMR solution is consistent with that of the MR-driven AMR solution. The performance of the HAMR scheme is evaluated for a range of test problems, from pure hydrodynamics to turbulent combustion.

preprint2020arXiv

SN Ia DDT Explosions Powered by the Zeldovich Reactivity Gradient Mechanism

The aim of this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyze weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for deflagration-to-detonation transition (DDT) to occur. We observe formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel'dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated either aided by reactivity gradients or by collisions of carbon detonations. The observed evolutionary timescales are found sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to DDT might be of consequence for the chemical composition of the outer ejecta regions of SN Ia from the single degenerate channel, and offer potential for validation of the proposed model.

preprint2017arXiv

Detonability of white dwarf plasma: turbulence models at low densities

We study the conditions required to produce self-sustained detonations in turbulent, carbon-oxygen degenerate plasma at low densities. We perform a series of three-dimensional hydrodynamic simulations of turbulence driven with various degrees of compressibility. The average conditions in the simulations are representative of models of merging binary white dwarfs. We find that material with very short ignition times is abundant in the case that turbulence is driven compressively. This material forms contiguous structures that persist over many ignition times, and that we identify as prospective detonation kernels. Detailed analysis of prospective kernels reveals that these objects are centrally-condensed and their shape is characterized by low curvature, supportive of self-sustained detonations. The key characteristic of the newly proposed detonation mechanism is thus high degree of compressibility of turbulent drive. The simulated detonation kernels have sizes notably smaller than the spatial resolution of any white dwarf merger simulation performed to date. The resolution required to resolve kernels is 0.1 km. Our results indicate a high probability of detonations in such well-resolved simulations of carbon-oxygen white dwarf mergers. These simulations will likely produce detonations in systems of lower total mass, thus broadening the population of white dwarf binaries capable of producing Type Ia supernovae. Consequently, we expect a downward revision of the lower limit of the total merger mass that is capable of producing a prompt detonation. We review application of the new detonation mechanism to various explosion scenarios of single, Chandrasekhar-mass white dwarfs.

preprint2014arXiv

Toward connecting core-collapse supernova theory with observations: I. Shock revival in a 15 Msun blue supergiant progenitor with SN 1987A energetics

We study the evolution of the collapsing core of a 15 Msun blue supergiant supernova progenitor from the core bounce until 1.5 seconds later. We present a sample of hydrodynamic models parameterized to match the explosion energetics of SN 1987A. We find the spatial model dimensionality to be an important contributing factor in the explosion process. Compared to two-dimensional simulations, our three-dimensional models require lower neutrino luminosities to produce equally energetic explosions. We estimate that the convective engine in our models is 4% more efficient in three dimensions than in two dimensions. We propose that the greater efficiency of the convective engine found in three-dimensional simulations might be due to the larger surface-to-volume ratio of convective plumes, which aids in distributing energy deposited by neutrinos. We do not find evidence of the standing accretion shock instability nor turbulence being a key factor in powering the explosion in our models. Instead, the analysis of the energy transport in the post-shock region reveals characteristics of penetrative convection. The explosion energy decreases dramatically once the resolution is inadequate to capture the morphology of convection on large scales. This shows that the role of dimensionality is secondary to correctly accounting for the basic physics of the explosion. We also analyze information provided by particle tracers embedded in the flow, and find that the unbound material has relatively long residency times in two-dimensional models, while in three dimensions a significant fraction of the explosion energy is carried by particles with relatively short residency times.

preprint2013arXiv

Prospects of Turbulence Studies in High-Energy Density Laser-Generated Plasma: Numerical Investigations in Two Dimensions

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.

preprint2013arXiv

The Braginskii model of the Rayleigh-Taylor instability. I. Effects of self-generated magnetic fields and thermal conduction in two dimensions

(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.

preprint2008arXiv

The Physics of Wind-Fed Accretion

We provide a brief review of the physical processes behind the radiative driving of the winds of OB stars and the Bondi-Hoyle-Lyttleton capture and accretion of a fraction of the stellar wind by a compact object, typically a neutron star, in detached high-mass X-ray binaries (HMXBs). In addition, we describe a program to develop global models of the radiatively-driven photoionized winds and accretion flows of HMXBs, with particular attention to the prototypical system Vela X-1. The models combine XSTAR photoionization calculations, HULLAC emission models appropriate to X-ray photoionized plasmas, improved models of the radiative driving of photoionized winds, FLASH time-dependent adaptive-mesh hydrodynamics calculations, and Monte Carlo radiation transport. We present two- and three-dimensional maps of the density, temperature, velocity, ionization parameter, and emissivity distributions of representative X-ray emission lines, as well as synthetic global Monte Carlo X-ray spectra. Such models help to better constrain the properties of the winds of HMXBs, which bear on such fundamental questions as the long-term evolution of these binaries and the chemical enrichment of the interstellar medium.

preprint2006arXiv

Detonating Failed Deflagration Model of Thermonuclear Supernovae II. Comparison to Observations

We develop and demonstrate the methodology of testing multi-dimensional supernova models against observations by studying the properties of one example of the detonation from failed deflagration (DFD) explosion model of thermonuclear supernovae. Using time-dependent multi-dimensional radiative transfer calculations, we generate the synthetic broadband optical light curves, near-infrared light curves, color evolution curves, full spectral time-series, and spectropolarization of the model, as seen from various viewing angles. All model observables are critically evaluated against examples of well-observed, standard Type Ia supernovae (SNe Ia). We explore the consequences of the intrinsic model asphericity by studying the dependence of the model emission on viewing angle, and by quantifying the resulting dispersion in (and internal correlations between) various model observables. These statistical properties of the model are also evaluated against those of the available observational sample of SNe Ia. On the whole, the DFD model shows good agreement with a broad range of SN Ia observations. Certain deficiencies are also apparent, and point to further developments within the basic theoretical framework. We also identify several intriguing orientation effects in the model which suggest ways in which the asphericity of SNe Ia may contribute to their photometric and spectroscopic diversity and, conversely, how the relative homogeneity of SNe Ia constrains the degree of asymmetry allowable in the models. The comprehensive methodology adopted in this work proves an essential component of developing and validating theoretical supernova models, and helps motivate and clearly define future directions in both the modeling and the observation of SNe Ia.

preprint1994arXiv

Fast Radiative Shocks in Dense Media. III. Properties of the Emission

Evolution of fast, radiative shocks in high density medium is presented. Ionizing spectra and approximate broad band light curves of the shocked gas are calculated. Emergent shock spectra, as seen by a distant observer, are obtained from photoionization models. The emergent spectra have a power-law shape $F_ν\proptoν^{-α}$ with mean spectral index $α\sim0.6-1.0$ in the energy range $0.01-10$ keV, and have a high-energy cutoff corresponding to the original shock velocity. It is shown that the models exhibit promising features that may account for some photometric and spectral properties of Active Galactic Nuclei.