Roman Samulyak

2papers

2 Papers

NAMar 13, 2016
AP-Cloud: Adaptive Particle-in-Cloud Method for Optimal Solutions to Vlasov-Poisson Equation

Xingyu Wang, Roman Samulyak, Xiangmin Jiao et al.

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.

NAMar 17, 2016
Lagrangian Particle Method for Compressible Fluid Dynamics

Hsin-Chiang Chen, Roman Samulyak, Wei Li

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.