Kinetic Solvers with Adaptive Mesh in Phase Space

An Adaptive Mesh in Phase Space (AMPS) methodology has been developed for solving multi-dimensional kinetic equations by the discrete velocity method. A Cartesian mesh for both configuration (r) and velocity (v) spaces is produced using a tree of trees data structure. The mesh in r-space is automatically generated around embedded boundaries and dynamically adapted to local solution properties. The mesh in v-space is created on-the-fly for each cell in r-space. Mappings between neighboring v-space trees implemented for the advection operator in configuration space. We have developed new algorithms for solving the full Boltzmann and linear Boltzmann equations with AMPS. Several recent innovations were used to calculate the discrete Boltzmann collision integral with dynamically adaptive mesh in velocity space: importance sampling, multi-point projection method, and the variance reduction method. We have developed an efficient algorithm for calculating the linear Boltzmann collision integral for elastic and inelastic collisions in a Lorentz gas. New AMPS technique has been demonstrated for simulations of hypersonic rarefied gas flows, ion and electron kinetics in weakly ionized plasma, radiation and light particle transport through thin films, and electron streaming in semiconductors. We have shown that AMPS allows minimizing the number of cells in phase space to reduce computational cost and memory usage for solving challenging kinetic problems

Robust event-driven particle tracking in complex geometries

The authors gratefully acknowledge the support of the German Research Foundation (DFG) through Grants PO 472/20 and SFB-814. FPS and ZISC are thanked for support. We would like to thank Prapanch Nair for helpful discussions.Peer reviewedPostprin

Cosmological Radiative Transfer Codes Comparison Project I: The Static Density Field Tests

Radiative transfer simulations are now at the forefront of numerical astrophysics. They are becoming crucial for an increasing number of astrophysical and cosmological problems; at the same time their computational cost has come to the reach of currently available computational power. Further progress is retarded by the considerable number of different algorithms (including various flavours of ray-tracing and moment schemes) developed, which makes the selection of the most suitable technique for a given problem a non-trivial task. Assessing the validity ranges, accuracy and performances of these schemes is the main aim of this paper, for which we have compared 11 independent RT codes on 5 test problems: (0) basic physics, (1) isothermal H II region expansion and (2) H II region expansion with evolving temperature, (3) I-front trapping and shadowing by a dense clump, (4) multiple sources in a cosmological density field. The outputs of these tests have been compared and differences analyzed. The agreement between the various codes is satisfactory although not perfect. The main source of discrepancy appears to reside in the multi-frequency treatment approach, resulting in different thicknesses of the ionized-neutral transition regions and different temperature structure. The present results and tests represent the most complete benchmark available for the development of new codes and improvement of existing ones. To this aim all test inputs and outputs are made publicly available in digital form.Comment: 32 pages, 39 figures (all color), comments welcom

A Three-dimensional Particle-in-Cell Methodology on Unstructured Voronoi Grids with Applications to Plasma Microdevices

The development and numerical implementation of a three-dimensional Particle-In-Cell (PIC) methodology on unstructured Voronoi-Delauney tetrahedral grids is presented. Charge assignment and field interpolation weighting schemes of zero- and first-order are formulated based on the theory of long-range constraints for three-dimensional unstructured grids. The algorithms for particle motion, particle tracing, particle injection, and loading are discussed. Solution to Poisson\u27s equation is based on a finite-volume formulation that takes advantage of the Voronoi-Delauney dual. The PIC methodology and code are validated by application to the problem of current collection by cylindrical Langmuir probes in stationary and moving collisionless plasmas. Numerical results are compared favorably with previous numerical and analytical solutions for a wide range of probe radius to Debye length ratios, probe potentials, and electron to ion temperature ratios. A methodology for evaluation of the heating, slowing-down and deflection times in 3D PIC simulations is presented. An extensive parametric evaluation is performed and the effects of the number of computational particles per cell, the ratio of cell-edge to Debye length, and timestep are investigated. The unstructured PIC code is applied to the simulation of Field Emission Array (FEA) cathodes. Electron injection conditions are obtained from a Field Emission microtip model and the simulation domain includes the FEA cathode and anode. Currents collected by the electrodes are compared to theoretical values. Simulations show the formation of the virtual cathode and three-dimensional effects under certain injection conditions. The unstructured PIC code is also applied to the simulation of a micro-Retarding Potential Analyzer. For simple cases the current at the collector plate is compared favorably with theoretical predictions. The simulations show the complex structure of the potential inside the segmented microchannel, the phase space of plasma species and the space-charge effects not captured by the theory

Langevin PDF simulation of particle deposition in a turbulent pipe flow

The paper deals with the description of particle deposition on walls from a turbulent flow over a large range of particle diameter, using a Langevin PDF model. The first aim of the work is to test how the present Langevin model is able to describe this phenomenon and to outline the physical as- pects which play a major role in particle deposition. The general features and characteristics of the present stochastic model are first recalled. Then, results obtained with the standard form of the model are presented along with an analysis which has been carried out to check the sensitivity of the predictions on different mean fluid quantities. These results show that the physical repre- sentation of the near-wall physics has to be improved and that, in particular, one possible route is to introduce specific features related to the near-wall coherent structures. In the following, we propose a simple phenomenological model that introduces some of the effects due to the presence of turbulent coherent structures on particles in a thin layer close to the wall. The results obtained with this phenomenological model are in good agreement with experimental evidence and this suggests to pursue in that direction, towards the development of more general and rigorous stochastic models that provide a link between a geometrical description of turbulent flow and a statistical one.Comment: 40 pages, 8 figure

Development of an Unstructured 3-D Direct Simulation Monte Carlo/Particle-in-Cell Code and the Simulation of Microthruster Flows

This work is part of an effort to develop an unstructured, three-dimensional, direct simulation Monte Carlo/particle-in-cell (DSMC/PIC) code for the simulation of non-ionized, fully ionized and partially-ionized flows in micropropulsion devices. Flows in microthrusters are often in the transitional to rarefied regimes, requiring numerical techniques based on the kinetic description of the gaseous or plasma propellants. The code is implemented on unstructured tetrahedral grids to allow discretization of arbitrary surface geometries and includes an adaptation capability. In this study, an existing 3D DSMC code for rarefied gasdynamics is improved with the addition of the variable hard sphere model for elastic collisions and a vibrational relaxation model based on discrete harmonic oscillators. In addition the existing unstructured grid generation module of the code is enhanced with grid-quality algorithms. The unstructured DSMC code is validated with simulation of several gaseous micronozzles and comparisons with previous experimental and numerical results. Rothe s 5-mm diameter micronozzle operating at 80 Pa is simulated and results are compared favorably with the experiments. The Gravity Probe-B micronozzle is simulated in a domain that includes the injection chamber and plume region. Stagnation conditions include a pressure of 7 Pa and mass flow rate of 0.012 mg/s. The simulation examines the role of injection conditions in micronozzle simulations and results are compared with previous Monte Carlo simulations. The code is also applied to the simulation of a parabolic planar micronozzle with a 15.4-micron throat and results are compared with previous 2D Monte Carlo simulations. Finally, the code is applied to the simulation of a 34-micron throat MEMS-fabricated micronozzle. The micronozzle is planar in profile with sidewalls binding the upper and lower surfaces. The stagnation pressure is set at 3.447 kPa and represents an order of magnitude lower pressure than used in previous experiments. The simulation demonstrates the formation of large viscous boundary layers in the sidewalls. A particle-in-cell model for the simulation of electrostatic plasmas is added to the DSMC code. Solution to Poisson\u27s equation on unstructured grids is obtained with a finite volume implementation. The Poisson solver is validated by comparing results with analytic solutions. The integration of the ionized particle equations of motion is performed via the leapfrog method. Particle gather and scatter operations use volume weighting with linear Lagrange polynomial to obtain an acceptable level of accuracy. Several methods are investigated and implemented to calculate the electric field on unstructured meshes. Boundary conditions are discussed and include a formulation of plasma in bounded domains with external circuits. The unstructured PIC code is validated with the simulation of a high voltage sheath formation