Source author record

Roman Samulyak

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

ResearcherUnclaimed source record

math.NA physics.comp-ph physics.plasm-ph math-ph math.MP nlin.CD

Catalog footprint

What is connected

11works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Computational Model for Granular Flows based on the Lagrangian Particle Method

A numerical model and parallel software for 3D simulations of granular flows have been developed based on the Lagrangian particle (LP) method [R.Samulyak, X. Wang, H.-C. Chen, Lagrangian particle method for compressible fluid dynamics, J. Comput. Phys. 362 (2018) 1-19], originally developed for compressible hydrodynamic flows, including free surface / multiphase problems. LP uses local polynomial least square fittings on particle-based stencils that ensure numerical convergence to the prescribed order. The granular flow model implements continuum equations with a $μ(I)$-rheology closure that is capable of describing two-directional transitions of the flow between various regimes characterized by solid-, liquid-, and gas-like features. The granular flow code has been parallelized for distributed memory supercomputers and validated by comparing 3D simulations to experimental data on the collapse of granular columns. Numerical simulations showing the granular flow transition to a gas-like regime have also been presented.

preprint2022arXiv

Massively-Parallel Lagrangian Particle Code and Applications

Massively-parallel, distributed-memory algorithms for the Lagrangian particle hydrodynamic method [R. Samulyak, X. Wang, H.-C. Chen, Lagrangian particle method for compressible fluid dynamics, J. Comput. Phys., 362 (2018), 1-19] have been developed, verified, and implemented. The key component of parallel algorithms is a particle management module that includes a parallel construction of octree databases, dynamic adaptation and refinement of octrees, and particle migration between parallel subdomains. The particle management module is based on the p4est (parallel forest of k-trees) library. The massively-parallel Lagrangian particle code has been applied to a variety of fundamental science and applied problems. A summary of Lagrangian particle code applications to the injection of impurities into thermonuclear fusion devices and to the simulation of supersonic hydrogen jets in support of laser-plasma wakefield acceleration projects has also been presented.

preprint2021arXiv

SPACE: 3D Parallel Solvers for Vlasov-Maxwell and Vlasov-Poisson Equations for Relativistic Plasmas with Atomic Transformations

A parallel, relativistic, three-dimensional particle-in-cell code SPACE has been developed for the simulation of electromagnetic fields, relativistic particle beams, and plasmas. In addition to the standard second-order Particle-in-Cell (PIC) algorithm, SPACE includes efficient novel algorithms to resolve atomic physics processes such as multi-level ionization of plasma atoms, recombination, and electron attachment to dopants in dense neutral gases. SPACE also contains a highly adaptive particle-based method, called Adaptive Particle-in-Cloud (AP-Cloud), for solving the Vlasov-Poisson problems. It eliminates the traditional Cartesian mesh of PIC and replaces it with an adaptive octree data structure. The code's algorithms, structure, capabilities, parallelization strategy and performances have been discussed. Typical examples of SPACE applications to accelerator science and engineering problems are also presented.

preprint2016arXiv

AP-Cloud: Adaptive Particle-in-Cloud Method for Optimal Solutions to Vlasov-Poisson Equation

We propose a new adaptive Particle-in-Cloud (AP-Cloud) method for obtaining optimal numerical solutions to the Vlasov-Poisson equation. Unlike the traditional particle-in-cell (PIC) method, which is commonly used for solving this problem, the AP-Cloud adaptively selects computational nodes or particles to deliver higher accuracy and efficiency when the particle distribution is highly non-uniform. Unlike other adaptive techniques for PIC, our method balances the errors in PDE discretization and Monte Carlo integration, and discretizes the differential operators using a generalized finite difference (GFD) method based on a weighted least square formulation. As a result, AP-Cloud is independent of the geometric shapes of computational domains and is free of artificial parameters. Efficient and robust implementation is achieved through an octree data structure with 2:1 balance. We analyze the accuracy and convergence order of AP-Cloud theoretically, and verify the method using an electrostatic problem of a particle beam with halo. Simulation results show that the AP-Cloud method is substantially more accurate and faster than the traditional PIC, and it is free of artificial forces that are typical for some adaptive PIC techniques.

preprint2016arXiv

Lagrangian Particle Method for Compressible Fluid Dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) an upwinding second-order particle-based algorithm with limiter, providing accuracy and long term stability, (c) elimination of the dependence on artificial parameters such as the smoothening length in SPH, causing difficulties in the case of large density changes, and (d) accurate resolution of states at free interfaces. Numerical verification test demonstrating the convergence order are presented as well as examples of complex free surface flows.

preprint2015arXiv

On Stability of Targets for Plasma Jet Induced Magnetoinertial Fusion

The compression and stability of plasma targets for the plasma jet-induced magneto-inertial fusion (PJMIF) have been investigated via large scale simulations using the FronTier code capable of explicit tracking of material interfaces. In the PJMIF concept, a plasma liner, formed by the merger of a large number of radial, highly supersonic plasma jets, implodes on a magnetized plasma target and compresses it to conditions of the fusion ignition. A multi-stage computational approach for simulations of the liner-target interaction and the compression of plasma targets has been developed to minimize computing time. Simulations revealed important features of the target compression process, including instability and disintegration of targets. The non-uniformity of the leading edge of the liner, caused by plasma jets as well as oblique shock waves between them, leads to instabilities during the target compression. By using front tracking, the evolution of targets has been studied in 3-dimensional simulations. Optimization studies of target compression with different number of jets have also been performed.

preprint2015arXiv

Transonic ablation flow regimes of high-Z pellets

In this letter, we report results of numerical studies of the ablation of argon and neon pellets in tokamaks and compare them with theoretical predictions and studies of deuterium pellets. Results demonstrate the influence of atomic physics processes on the pellet ablation process.

preprint2013arXiv

A Comparison Study of Two Methods for Elliptic Boundary Value Problems

In this paper, we perform a comparison study of two methods (the embedded boundary method and several versions of the mixed finite element method) to solve an elliptic boundary value problem.

preprint2013arXiv

An Embedded Boundary Method for Two Phase Incompressible Flow

We develop an embedded boundary method (EBM) to solve the two-phase incompressible flow with piecewise constant density. The front tracking method is used to track the interface. The fractional step methods are used to solve the incompressible Navier-Stokes equations while the EBM is used in the projection step to solve an elliptic interface problem for the pressure with a jump equal to the surface tension force across the interface. Several examples are used to verify the accuracy of the method.

preprint1999arXiv

New mathematical models for particle flow dynamics

A new class of integro-partial differential equation models is derived for the prediction of granular flow dynamics. These models are obtained using a novel limiting averaging method (inspired by techniques employed in the derivation of infinite-dimensional dynamical systems models) on the Newtonian equations of motion of a many-particle system incorporating widely used inelastic particle-particle force formulas. By using Taylor series expansions, these models can be approximated by a system of partial differential equations of the Navier-Stokes type. The exact or approximate governing equations obtained are far from simple, but they are less complicated than most of the continuum models now being used to predict particle flow behavior. Solutions of the new models for granular flows down inclined planes and in vibrating beds are compared with known experimental and analytical results and good agreement is obtained.

preprint1998arXiv

The integrability of Lie-invariant geometric objects generated by ideals in the Grassmann algebra

We investigate closed ideals in the Grassmann algebra serving as bases of Lie-invariant geometric objects studied before by E. Cartan. Especially, the E. Cartan theory is enlarged for Lax integrable nonlinear dynamical systems to be treated in the frame work of the Wahlquist Estabrook prolongation structures on jet-manifolds and Cartan-Ehresmann connection theory on fibered spaces. General structure of integrable one-forms augmenting the two-forms associated with a closed ideal in the Grassmann algebra is studied in great detail. An effective Maurer-Cartan one-forms construction is suggested that is very useful for applications. As an example of application the developed Lie-invariant geometric object theory for the Burgers nonlinear dynamical system is considered having given rise to finding an explicit form of the associated Lax type representation.