Related papers: Computational geometric methods for preferential c…
A new approach for integration of motion in many-body systems of interacting polyatomic molecules is proposed. It is based on splitting time propagation of pseudo-variables in a modified phase space, while the real translational and…
We study suspensions of oblate rigid particles in a viscous fluid for different values of the particle volume fractions. Direct numerical simulations have been performed using a direct-forcing immersed boundary method to account for the…
Pinched flow fractionation is shown to be an efficient and selective way to quickly separate particles by size in a very polydisperse semi-concentrated suspension. In an effort to optimize the method, we discuss the quantitative influence…
Particle-In-Cell codes are widely used for plasma physics simulations. It is often the case that particles within a computational cell need to be split to improve the statistics or, in the case of non-uniform meshes, to avoid the…
In this paper we demonstrate a computational method to solve the inverse scattering problem for a star-shaped, smooth, penetrable obstacle in 2D. Our method is based on classical ideas from computational geometry. First, we approximate the…
In this work it is shown how the immersed boundary method of (Peskin2002) for modeling flexible structures immersed in a fluid can be extended to include thermal fluctuations. A stochastic numerical method is proposed which deals with…
We consider a composite medium, which consists of a homogeneous matrix containing a statistically homogeneous set of multimodal spherical inclusions. This model is used to represent the morphology of heterogeneous solid propellants (HSP)…
We discuss a spectral method for the numerical solution of the Vlasov-Poisson system where the velocity space is decomposed by means of an Hermite basis. We describe a semi-implicit time discretization that extends the range of numerical…
A class of variational schemes for the hydrodynamic-electrodynamic model of lossless free-electron gas in a quasineutral background is developed for high-quality simulations of surface plasmon polaritons. The Lagrangian density of lossless…
We propose a one-way coupled model that tracks individual primary particles in a conceptually simple cellular flow setup to predict flocculation in turbulence. This computationally efficient model accounts for Stokes drag, lubrication,…
We consider the geometric numerical integration of Hamiltonian systems subject to both equality and "hard" inequality constraints. As in the standard geometric integration setting, we target long-term structure preservation. We…
By numerical integration, we study the relaxation dynamics of degenerate harmonic oscillator modes dispersively coupled to particle positions. Depending on whether the effective inertial potential induced by the oscillators keep the…
A numerical method is presented for first-principle simulations of charged colloidal dispersions in electrolyte solutions. Utilizing a smoothed profile for colloid-solvent boundaries, efficient mesoscopic simulations are enabled for…
We investigate the motion of a suspended non-Brownian sphere past a fixed cylindrical or spherical obstacle in the limit of zero Reynolds number for arbitrary particle-obstacle aspect ratios. We consider both a suspended sphere moving in a…
Computational modeling and simulation of fluid-structure interactions constitute a fundamental cornerstone for advancing aerospace engineering endeavors. This paper addresses the notion and implementation of the immersed boundary method for…
Nonlocal diffusion model provides an appropriate description of the diffusion process of solute in the complex medium, which cannot be described properly by classical theory of PDE. However, the operators in the nonlocal diffusion models…
We present a new particle-based (discrete element) numerical method for the simulation of granular dynamics, with application to motions of particles on small solar system body and planetary surfaces. The method employs the parallel N-body…
A kinetic flux-splitting procedure used in conjunction with local thermodynamic equilibrium in a finite volume allows us to investigate numerically discrete-velocity gas flows. The procedure, outlined for a general discrete-velocity gas, is…
The most difficult aspect of the realistic modeling of granular materials is how to capture the real shape of the particles. Here we present a method to simulate granular materials with complex-shaped particles. The particle shape is…
Calculating cost-effective solutions to particle dynamics in viscous flows is an important problem in many areas of industry and nature. We implement a second-order symmetric splitting method on the governing equations for a rigid…