4,800 research outputs found
High Performance Computing With A Conservative Spectral Boltzmann Solver
We present new results building on the conservative deterministic spectral method for the space inhomogeneous Boltzmann equation developed by Gamba and Tharkabhushaman. This approach is a two-step process that acts on the weak form of the Boltzmann equation, and uses the machinery of the Fourier transform to reformulate the collisional integral into a weighted convolution in Fourier space. A constrained optimization problem is solved to preserve the mass, momentum, and energy of the resulting distribution. We extend this method to second order accuracy in space and time, and explore how to leverage the structure of the collisional formulation for high performance computing environments. The locality in space of the collisional term provides a straightforward memory decomposition, and we perform some initial scaling tests on high performance computing resources. We also use the improved computational power of this method to investigate a boundary-layer generated shock problem that cannot be described by classical hydrodynamics.Mathematic
Conservative Deterministic Spectral Boltzmann Solver Near The Grazing Collisions Limit
We present new results building on the conservative deterministic spectral method for the space homogeneous Boltzmann equation developed by Gamba and Tharkabhushaman. This approach is a two-step process that acts on the weak form of the Boltzmann equation, and uses the machinery of the Fourier transform to reformulate the collisional integral into a weighted convolution in Fourier space. A constrained optimization problem is solved to preserve the mass, momentum, and energy of the resulting distribution. Within this framework we have extended the formulation to the case of more general case of collision operators with anisotropic scattering mechanisms, which requires a new formulation of the convolution weights. We also derive the grazing collisions limit for the method, and show that it is consistent with the Fokker-Planck-Landau equations as the grazing collisions parameter goes to zero.Mathematic
On the stability of homogeneous solutions to some aggregation models
Vasculogenesis, i.e. self-assembly of endothelial cells leading to capillary network formation, has been the object of many experimental investigations in recent years, due to its relevance both in physiological and in pathological conditions. We performed a detailed linear stability analysis of two models of in vitro vasculogenesis, with the aim of checking their potential for structure formation starting from initial data representing a continuum cell monolayer. The first model turns out to be unstable at low cell densities, while pressure stabilizes it at high densities. The second model is instead stable at low cell densities. Detailed information about the instability regions and the structure of the critical wave numbers are obtained in several interesting limiting cases. We expect that altogether, this information will be useful for further comparisons of the two models with experiments
Test of a simple and flexible molecule model for alpha-, beta- and gamma-S8 crystals
S8 is the most stable compound of elemental sulfur in solid and liquid
phases, at ambient pressure and below 400K. Three crystalline phases of S8 have
been clearly identified in this range of thermodynamic parameters, although no
calculation of its phase diagram has been performed yet. alpha- and gamma-S8
are orientationally ordered crystals while beta-S8 is measured as
orientationally disordered. In this paper we analyze the phase diagram of S8
crystals, as given by a simple and flexible molecule model, via a series of
molecular dynamics (MD) simulations.
The calculations are performed in the constant pressure- constant temperature
ensemble, using an algorithm that is able to reproduce structural phase
transitions.Comment: RevTex,7 pages, 5 figures,to appear in J. Chem. Phy
A Fast Conservative Spectral Solver For The Nonlinear Boltzmann Collision Operator
We present a conservative spectral method for the fully nonlinear Boltzmann collision operator based on the weighted convolution structure in Fourier space developed by Gamba and Tharkabhushnanam.. This method can simulate a broad class of collisions, including both elastic and inelastic collisions as well as angularly dependent cross sections in which grazing collisions play a major role. The extension presented in this paper consists of factorizing the convolution weight on quadrature points by exploiting the symmetric nature of the particle interaction law, which reduces the computational cost and memory requirements of the method to O(M(2)N(4)logN) from the O(N-6) complexity of the original spectral method, where N is the number of velocity grid points in each velocity dimension and M is the number of quadrature points in the factorization, which can be taken to be much smaller than N. We present preliminary numerical results.Mathematic
A discontinuous Galerkin method for the Vlasov-Poisson system
A discontinuous Galerkin method for approximating the Vlasov-Poisson system
of equations describing the time evolution of a collisionless plasma is
proposed. The method is mass conservative and, in the case that piecewise
constant functions are used as a basis, the method preserves the positivity of
the electron distribution function and weakly enforces continuity of the
electric field through mesh interfaces and boundary conditions. The performance
of the method is investigated by computing several examples and error estimates
associated system's approximation are stated. In particular, computed results
are benchmarked against established theoretical results for linear advection
and the phenomenon of linear Landau damping for both the Maxwell and Lorentz
distributions. Moreover, two nonlinear problems are considered: nonlinear
Landau damping and a version of the two-stream instability are computed. For
the latter, fine scale details of the resulting long-time BGK-like state are
presented. Conservation laws are examined and various comparisons to theory are
made. The results obtained demonstrate that the discontinuous Galerkin method
is a viable option for integrating the Vlasov-Poisson system.Comment: To appear in Journal for Computational Physics, 2011. 63 pages, 86
figure
Linear Theory of Electron-Plasma Waves at Arbitrary Collisionality
The dynamics of electron-plasma waves are described at arbitrary
collisionality by considering the full Coulomb collision operator. The
description is based on a Hermite-Laguerre decomposition of the velocity
dependence of the electron distribution function. The damping rate, frequency,
and eigenmode spectrum of electron-plasma waves are found as functions of the
collision frequency and wavelength. A comparison is made between the
collisionless Landau damping limit, the Lenard-Bernstein and Dougherty
collision operators, and the electron-ion collision operator, finding large
deviations in the damping rates and eigenmode spectra. A purely damped entropy
mode, characteristic of a plasma where pitch-angle scattering effects are
dominant with respect to collisionless effects, is shown to emerge numerically,
and its dispersion relation is analytically derived. It is shown that such a
mode is absent when simplified collision operators are used, and that
like-particle collisions strongly influence the damping rate of the entropy
mode.Comment: 23 pages, 10 figures, accepted for publication on Journal of Plasma
Physic
- …