107,603 research outputs found
A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations
Recent work by McClarren & Hauck [29] suggests that the filtered spherical
harmonics method represents an efficient, robust, and accurate method for
radiation transport, at least in the two-dimensional (2D) case. We extend their
work to the three-dimensional (3D) case and find that all of the advantages of
the filtering approach identified in 2D are present also in the 3D case. We
reformulate the filter operation in a way that is independent of the timestep
and of the spatial discretization. We also explore different second- and
fourth-order filters and find that the second-order ones yield significantly
better results. Overall, our findings suggest that the filtered spherical
harmonics approach represents a very promising method for 3D radiation
transport calculations.Comment: 29 pages, 13 figures. Version matching the one in Journal of
Computational Physic
Fourier spectral methods for fractional-in-space reaction-diffusion equations
Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is computationally demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reactiondiffusion equations. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is show-cased by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models,together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator
Immersed boundary methods for numerical simulation of confined fluid and plasma turbulence in complex geometries: a review
Immersed boundary methods for computing confined fluid and plasma flows in
complex geometries are reviewed. The mathematical principle of the volume
penalization technique is described and simple examples for imposing Dirichlet
and Neumann boundary conditions in one dimension are given. Applications for
fluid and plasma turbulence in two and three space dimensions illustrate the
applicability and the efficiency of the method in computing flows in complex
geometries, for example in toroidal geometries with asymmetric poloidal
cross-sections.Comment: in Journal of Plasma Physics, 201
A direct method for the Boltzmann equation based on a pseudo-spectral velocity space discretization
A deterministic method is proposed for solving the Boltzmann equation. The
method employs a Galerkin discretization of the velocity space and adopts, as
trial and test functions, the collocation basis functions based on weights and
roots of a Gauss-Hermite quadrature. This is defined by means of half- and/or
full-range Hermite polynomials depending whether or not the distribution
function presents a discontinuity in the velocity space. The resulting
semi-discrete Boltzmann equation is in the form of a system of hyperbolic
partial differential equations whose solution can be obtained by standard
numerical approaches. The spectral rate of convergence of the results in the
velocity space is shown by solving the spatially uniform homogeneous relaxation
to equilibrium of Maxwell molecules. As an application, the two-dimensional
cavity flow of a gas composed by hard-sphere molecules is studied for different
Knudsen and Mach numbers. Although computationally demanding, the proposed
method turns out to be an effective tool for studying low-speed slightly
rarefied gas flows
- …