BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shay I. Heizler

Shay I. Heizler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.stat-mechphysics.comp-phphysics.plasm-phastro-ph.IMcond-mat.mtrl-scicond-mat.softphysics.flu-dyn

Catalog footprint

What is connected

8works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

8 published item(s)

preprint2022arXiv

A New Discrete Implicit Monte Carlo Scheme for Simulating Radiative Transfer Problems

We present a new algorithm for radiative transfer-based on a statistical Monte Carlo approach-that does not suffer from teleportation effects, on the one hand, and yields smooth results, on the other hand. Implicit Monte Carlo (IMC) techniques for modeling radiative transfer have existed from the 1970s. When they are used for optically thick problems, however, the basic algorithm suffers from "teleportation" errors, where the photons propagate faster than the exact physical behavior, due to the absorption-blackbody emission processes. One possible solution is to use semianalog Monte Carlo, in its new implicit form (ISMC), which uses two kinds of particles, photons and discrete material particles. This algorithm yields excellent teleportation-free results, but it also produces noisier solutions (relative to classic IMC), due to its discrete nature. Here, we derive a new Monte Carlo algorithm, Discrete Implicit Monte Carlo (DIMC), which also uses the idea of two kinds of discrete particles, and thus does not suffer from teleportation errors. DIMC implements the IMC discretization and creates new radiation photons for each time step, unlike ISMC. Using the continuous absorption technique, DIMC yields smooth results like classic IMC. One of the main elements of the algorithm is the avoidance of the explosion of the particle population, by using particle merging. We test the new algorithm on 1D and 2D cylindrical problems, and show that it yields smooth, teleportation-free results. We finish by demonstrating the power of the new algorithm on a classic radiative hydrodynamic problem-an opaque radiative shock wave. This demonstrates the power of the new algorithm for astrophysical scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multi-Frequency Implicit Semi-analog Monte-Carlo (ISMC) Radiative Transfer Solver in Two-Dimensions (without Teleportation)

We study the multi-dimensional radiative transfer phenomena using the ISMC scheme, in both gray and multi-frequency problems. Implicit Monte-Carlo (IMC) schemes have been in use for five decades. The basic algorithm yields teleportation errors, where photons propagate faster than the correct heat front velocity. Recently [Poëtte and Valentin, J. Comp. Phys., 412, 109405 (2020)], a new implicit scheme based on the semi-analog scheme was presented and tested in several one-dimensional gray problems. In this scheme, the material energy of the cell is carried by material-particles, and the photons are produced only from existing material particles. As a result, the teleportation errors vanish, due to the infinite discrete spatial accuracy of the scheme. We examine the validity of the new scheme in two-dimensional problems, both in Cartesian and Cylindrical geometries. Additionally, we introduce an expansion of the new scheme for multi-frequency problems. We show that the ISMC scheme presents excellent results without teleportation errors in a large number of benchmarks, especially against the slow classic IMC convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Supersonic-subsonic transition region in radiative heat flow via self-similar solutions

We study the radiative hydrodynamics flow of radiation-driven heat waves in hot dense plasmas, using approximate self-similar solutions. Specifically, we have focused on the intermediate regime between pure radiative supersonic flow and the pure subsonic regime. These two regimes were investigated both using exact self-similar solutions and numerical simulations, however, most of the study used numerical simulations, mainly because the radiative heat wave and the shock regions are not self-similar altogether. In a milestone work [J. Garnier et al., Phys. Plas., 13, 092703 (2006)], it was found that for a specific power law dependency temperature profile, a unique exact self-similar solution exists, that is valid for all physical regimes. In this work we approximate Garnier's exact solution for a general power-law temperature-dependency, using simple analytical considerations. This approximate solution yields a good agreement compared to numerical simulations for the different thermodynamic profiles within the expected range of validity. In addition, we offer an approximate solution for the energies absorbed in the matter, again, for a general power-law temperature profile. Our approximate self-similar solution for the energy yields very good results comparing to exact numerical simulations for both gold and $\mathrm{Ta_2O_5}$. We also set a comparison of our self-similar solutions with the results of an experiment for radiation temperature measurement in a hohlraum in low-density foams that is addressed directly, to the intermediate regime, yielding a good agreement and similar trends. The different models as well as the numerical simulations are powerful tools to analyze the supersonic-subsonic transition region.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Key to understanding supersonic radiative Marshak waves using simple models and advanced simulations

This article studies the propagation of supersonic radiative Marshak waves. These waves are radiation dominated, and play an important role in inertial confinement fusion and in astrophysical and laboratory systems. For that reason, this phenomenon has attracted considerable experimental attention in recent decades in several different facilities. The present study integrates the various experimental results published in the literature, demonstrating a common physical base. A new simple semi-analytic model is derived and presented along with advanced radiative hydrodynamic implicit Monte Carlo direct numerical simulations, which explain the experimental results. This study identifies the main physical effects dominating the experiments, notwithstanding their different apparatuses and different physical regimes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Self-similar solution of the subsonic radiative heat equations using a binary equation of state

Radiative subsonic heat waves, and their radiation driven shock waves, are important hydro-radiative phenomena. The high pressure, causes hot matter in the rear part of the heat wave to ablate backwards. At the front of the heat wave, this ablation pressure generates a shock wave which propagates ahead of the heat front. Although no self-similar solution of both the ablation and shock regions exists, a solution for the full problem was found in a previous work. Here, we use this model in order to investigate the effect of the equation of state (EOS) on the propagation of radiation driven shocks. We find that using a single ideal gas EOS for both regions, as used in previous works, yields large errors in describing the shock wave. We use the fact that the solution is composed of two different self-similar solutions, one for the ablation region and one for the shock, and apply two ideal gas EOS (binary-EOS), one for each region, by fitting a detailed tabulated EOS to power laws at different regimes. By comparing the semi-analytic solution with a numerical simulation using a full EOS, we find that the semi-analytic solution describes both the heat and the shock regions well.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Full self-similar solutions of the subsonic radiative heat equations

We study the phenomenon of diffusive radiative heat waves (Marshak waves) under general boundary conditions. In particular, we derive full analytic solutions for the subsonic case, that include both the ablation and the shock wave regions. Previous works in this regime, based on the work of [R. Pakula and R. Sigel, Phys. Fluids. 443, 28, 232 (1985)], present self-similar solutions for the ablation region alone, since in general, the shock region and the ablation region are not self-similar together. Analytic results for both regions were obtained only for the specific case in which the ratio between the ablation front velocity and the shock velocity is constant. In this work, we derive a full analytic solution for the whole problem in general boundary conditions. Our solution is composed of two different self-similar solutions, one for each region, that are patched at the heat front. The ablative region of the heat wave is solved in a manner similar to previous works. Then, the pressure at the front, which is derived from the ablative region solution, is taken as a boundary condition to the shock region, while the other boundary is described by Hugoniot relations. The solution is compared to full numerical simulations in several representative cases. The numerical and analytic results are found to agree within $1\%$ in the ablation region, and within $2-5\%$ in the shock region. This model allows better prediction of the physical behavior of radiation induced shock waves, and can be applied for high energy density physics experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microbranching in mode-I fracture using large scale simulations of amorphous and perturbed lattice models

We study the high-velocity regime mode-I fracture instability when small microbranches start to appear near the main crack, using large scale simulations. Some of the features of those microbranches have been reproduced qualitatively in smaller scale studies (using ${\cal O}(10^4)$ atoms) on both a model of an amorphous materials (via the continuous random network model) and using perturbed lattice models. In this study, larger scale simulations (${\cal O}(10^6)$ atoms) were performed using multi-threading computing on a GPU device, in order to achieve more physically realistic results. First, we find that the microbranching pattern appears to be converging with the lattice width. Second, the simulations reproduce the growth of the size of a microbranch as a function of the crack velocity, as well as the increase of the amplitude of the derivative of the electrical resistance RMS with respect to the time as a function of the crack velocity. In addition, the simulations yield the correct branching angle of the microbranches, and the power law governing the shape of the microbranches seems to be lower than one, so that the side cracks turn over in the direction of propagation of the main crack as seen in experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Micro-branching in mode-I fracture in a randomly perturbed lattice

We study mode-I fracture in lattices with noisy bonds. In contrast to previous attempts, by using a small parameter that perturbs the force-law between the atoms in perfect lattices and using a 3-body force law, simulations reproduce the qualitative behavior of the beyond steady-state cracks in the high velocity regime, including reasonable micro-branching. As far as the physical properties such as the structure factor $g(r)$, the radial or angular distributions, these lattices share the physical properties of perfect lattices rather than that of an amorphous material (e.g., the continuous random network model). A clear transition can be seen between steady-state cracks, where a single crack propagates in the midline of the sample and the regime of unstable cracks, where micro-branches start to appear near the main crack, in line with previous experimental results. This is seen both in a honeycomb lattice and a fully hexagonal lattice. This model reproduces the main physical features of propagating cracks in brittle materials, including the behavior of velocity as a function of driving displacement and the increasing amplitude of oscillations of the electrical resistance. In addition, preliminary indications of power-law behavior of the micro-branch shapes can be seen, potentially reproducing one of the most intriguing experimental results of brittle fracture.

0 citations0 reviewsNo code linkedOpen access