DEMOCRITUS: An adaptive particle in cell (PIC) code for object-plasma interactions

Lapenta, Giovanni

doi:10.1016/j.jcp.2011.02.041

Cited by 24 publications

(21 citation statements)

References 38 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The physical to logical space transformation uses the tri-linear interpolation method described in [46]. The concept of performing particle assignment and interpolation in PIC codes has been successfully demonstrated by the CPIC [18] and DEMOCRITUS [37] codes.…”

Section: Composite Linear Volume Methodsmentioning

confidence: 99%

“…In the context of near-Earth plasma-body interactions, the interaction may be considered electrostatic and unmagnetized [37,50,55,56]. Under these assumptions, Maxwell's equations reduce to Poisson's equation for the electric potential φ,…”

Section: Governing Equationsmentioning

confidence: 99%

“…A majority of PIC codes developed for spacebased applications (PicUp3D [22], CPIC [18], DEMOCRITUS [37], and DiP2/3D [42]) have been developed to study the charging of satellites and instrument calibration. As a result, they do not include gas-surface interaction models found in DSMC codes.…”

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

pdFOAM: A PIC-DSMC code for near-Earth plasma-body interactions

et al. 2017

View full text Add to dashboard Cite

Understanding the interaction of the near-Earth space environment with orbiting bodies is critical, both from a design and scientific perspective. In Low Earth Orbit (LEO), the interaction between the ionosphere and orbiting objects is well studied from a charging perspective. Not well understood is the effect of the ionosphere on the motion of LEO objects i.e. charged aerodynamics. This paper presents the implementation, validation, and verification of the hybrid electrostatic Particle-in-Cell (PIC) -Direct Simulation Monte Carlo (DSMC) code, pdFOAM, to study both the neutral and charged particle aerodynamics of LEO objects. The 2D aerodynamic interaction of a cylinder with a fixed uniform surface potential of −50 V and mesothermal O + and H + plasmas representative of ionospheric conditions is investigated. New insights into the role of bounded ion jets and their effect on surface forces are presented. O + bounded ion jets are observed to cause a 4.4% increase ion direct Charged Particle Drag (dCPD), while H + ion jets produce a net reduction in H + drag by 23.7% i.e. they cause a thrust force. As a result, we conclude that past work, primarily based on Orbital Motion Limited theory, does not adequately capture the physics of LEO charged aerodynamics. Hence, we recommend a revisit of conclusions regarding the significance of CPD to LEO objects -pdFOAM being an appropriate tool for this purpose.

show abstract

Section: Composite Linear Volume Methodsmentioning

confidence: 99%

Section: Governing Equationsmentioning

confidence: 99%

See 1 more Smart Citation

pdFOAM: A PIC-DSMC code for near-Earth plasma-body interactions

et al. 2017

View full text Add to dashboard Cite

show abstract

“…Means the higher the numerical approaches of each individual element of the cycle, the higher approach is also the underlying full error of the cycle, see further ideas of adaptation and their underlying errors [56].…”

Section: )mentioning

confidence: 99%

Engineering Applications

Geiser

2015

Multicomponent and Multiscale Systems

View full text Add to dashboard Cite

In this chapter, we discuss the different engineering applications related to multicomponent and multiscale models, that occur in different categories (microscopic, mesoscopic and macroscopic). As we have described in the Introduction on p. xxv. That such models can be used as basic model to couple to more complicate models, describing materials, interfaces, etc., see Rosso and de Baas (Review of materials modelling: what makes a material function? Let me compute the ways, 2014, [1]).We deal with the following engineering applications with the two classified multiscaling approaches, see also the Introduction on p. xxv and in Fig. 5.1:• Different time-or spatial scales of same model (mono model or basic model).• Linking of different models (different models or multimodel).One of the main contributions to link the different scales and models together are the coupling techniques. In our motivation, such coupling of different time-and spatial scales or coupling of different models, need the suggested methods, e.g. multiscale methods, multicomponent methods, that allow a data transfer between the different scales and models. Such methods, we have introduced in the previous sections and now, we will close the gap between the theoretical discussions of methods and their applications to engineering problems. In such a stage, we have to adapt the numerical schemes with respect to the real-life properties and we obtain truly working multiscale approaches, that we solve the engineering complexities, see [1].Based on the real-life applications, we could study such helpful coupling of different scales or different models. Therefore, we overcome the gap of the recent problem in linking different scales of complicated engineering models in the industrial applications. Multiscale Methods for Langevin-Like EquationsAbstract In this section, we discuss multiscale methods, that solve Langevin-like equations, see [2]. The underlying ideas are to split into a deterministic and stochastic part of the Langevin equation, see [3]. The splitting methods are based on additive and iterative schemes, which are discussed with respect of their benefits and drawbacks, see [4]. We are motivated to reduce computational time for Coulomb collisions in plasma done in particle simulations. We extend splitting schemes, which are well-known in deterministic applications, to stochastic applications and modify the methods with respect to the stochastic terms. Such an idea allows to solve the multiscale behaviour of the coarse deterministic and fine stochastic timescales in an adequate computational time. Introduction of the ProblemIn the underlying problem, we are motivated to develop fast algorithms to solve the Coulomb collisions in plasma simulations, see [5].Recently in the literature, we can find two main ideas to deal with the Coulomb collisions in particle simulations. Here we have the following two ideas:• Binary algorithm: Particles in a finite cell selected into binary pairs. The collision is computed by the scattered velocities by an ...

show abstract

“…For example, the PEGASUS code [9] simulates astrophysical plasma using the hybrid-PIC approach. Object-plasma interactions are investigated in DEMOCRITUS [10] also using the PIC approach. QUICKPIC [11] has been developed to study the wakefield acceleration in plasmas.…”

Section: Introductionmentioning

confidence: 99%