Related papers: A coupled neutronics shock physics solver: impleme…
We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement (AMR)-based multiphysics…
Most massive stars end their lives with core collapse. However, it is not clear which explode as a Core-collapse Supernova (CCSN), leaving behind a neutron star and which collapse to black hole, aborting the explosion. One path to predict…
The newly proposed microstructural constitutive model for polycrystal viscoplasticity in cold and warm regimes (Motaman and Prahl, 2019), is implemented as a microstructural solver via user-defined material subroutine in a finite element…
Fourier solvers have become efficient tools to establish structure-property relations in heterogeneous materials. Introduced as an alternative to the Finite Element (FE) method, they are based on fixed-point solutions of the…
We present mesoscale numerical simulations based on the coupling of the fluctuating lattice Boltzmann method (FLBM) for multicomponent systems with a wetted finite-size particle model. This newly coupled methodologies are used to study the…
In this paper, we present a unified framework to simulate non-Newtonian behaviors. We combine viscous and elasto-plastic stress into a unified particle solver to achieve various non-Newtonian behaviors ranging from fluid-like to solid-like.…
The objective of this article is to address the discretisation of fractured/faulted poromechanical models using 3D polyhedral meshes in order to cope with the geometrical complexity of faulted geological models. A polytopal scheme is…
The recently developed auxiliary field diffusion Monte Carlo method is applied to compute the equation of state and the compressibility of neutron matter. By combining diffusion Monte Carlo for the spatial degrees of freedom and auxiliary…
Coupling of rupture processes in solids with waves also propagating in fluids is a prominent phenomenon arising during tectonic earthquakes. It is executed here in a single `monolithic' model which can asymptotically capture both damageable…
A hybrid Maxwell solver for fully relativistic and electromagnetic (EM) particle-in-cell (PIC) codes is described. In this solver, the EM fields are solved in $k$ space by performing an FFT in one direction, while using finite difference…
The no-core shell model (NCSM) is an \textit{ab initio} method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the…
Field emission coupled with molecular dynamics simulation (FEcMD) software package is a computational tool for studying atomic structure evolution, structural deformation, phase transitions, recrystallization as well as electron emission…
This paper provides an initial description of the Method of Simultaneous Solutions, a Monte Carlo approach that simultaneously solves multiple Boltzmann-transport-like phenomena. Here, it is used to simultaneously solve the neutron…
In this paper, we present a splitting algorithm to solve multicomponent transport models. These models are related to plasma simulations, in which we consider the local thermodynamic equilibrium and weakly ionised plasma-mixture models that…
We present the first two-dimensional general relativistic (GR) simulations of stellar core collapse and explosion with the CoCoNuT hydrodynamics code in combination with the VERTEX solver for energy-dependent, three-flavor neutrino…
We couple two-dimensional hydrodynamics to detailed one-dimensional multigroup flux-limited diffusion neutrino transport to investigate prompt convection in core collapse supernovae. Our initial conditions, time-dependent boundary…
We propose a spectral viscosity method (SVM) to approximate the incompressible Euler equations driven by a multiplicative noise. We show that SVM solution converges to a dissipative measure-valued martingale solution. These solutions are…
Unmodeled beyond Standard Model (BSM) physics in neutrino propagation can masquerade as parameter degeneracies in future precision measurements. Because the upcoming DUNE and T2HK experiments will operate at substantially different…
We propose an efficient algorithm for the evaluation of the potential and its gradient of gravitational/electrostatic $N$-body systems, which we call particle mesh multipole method (PMMM or PM$^3$). PMMM can be understood both as an…
A simplified method to calculate the critical mass of a fissile material sphere is presented. This is a purely pedagogical study, in part to elucidate the historical evolution of criticality calculations. This method employs only elementary…