Quasisteady and steady states in global gyrokinetic particle-in-cell simulations

Jolliet, S.; McMillan, B. F.; Vernay, T.; Villard, L.; Bottino, A.; Angelino, P.
May 2009
Physics of Plasmas;May2009, Vol. 16 Issue 5, p052307
Academic Journal
Collisionless delta-f gyrokinetic particle-in-cell simulations suffer from the entropy paradox, in which the entropy grows linearly in time while low-order moments are saturated. As a consequence, these simulations do not reach a steady state and are unsuited to make quantitative predictions. A solution to this issue is the introduction of artificial dissipation. The notion of steady state in gyrokinetic simulations is studied by deriving an evolution equation for the fluctuation entropy and applying it to the global collisionless particle-in-cell code ORB5 [S. Jolliet et al., Comput. Phys. Commun. 177, 409 (2007)]. It is shown that a recently implemented noise-control algorithm [B. F. McMillan et al., Phys. Plasmas 15, 052308 (2008)] based on a W-stat provides the necessary dissipation to reach a steady state. The two interesting situations of decaying and driven turbulence are considered. In addition, it is shown that a separate heating algorithm, not based on a W-stat, does not lead to a statistical steady state.


Related Articles

  • Can Water Clusters Ionize HCl under the Conditions of Arctic Stratosphere? 1. Intermolecular Interactions. Shevkunov, S. V. // Colloid Journal;Mar/Apr2004, Vol. 66 Issue 2, p216 

    Pseudopotential model was constructed to simulate the H3O+(H2O)nCl– clusters at room and stratosphere temperatures using the Monte Carlo method. Numerical values of interaction parameters were restored from the experimental data on the free energy and entropy of vapor nucleation on ions...

  • Amorphous polymeric nitrogen—toward an equation of state. Yakub, L. // Low Temperature Physics;Apr2009, Vol. 35 Issue 4, p326 

    The thermodynamic properties of amorphous polymeric nytrogen are analyzed theoretically on the basis of a free-energy expansion. An approach allowing estimation of the Helmholtz free energy of the high-pressure disordered condensed phases of nitrogen is proposed. The approach is based on...

  • Numerical entropy and phason elastic constants of plane random tilings with any 2D-fold symmetry. Destainville, N. // European Physical Journal B -- Condensed Matter;Jul2006, Vol. 52 Issue 1, p119 

    We perform Transition matrix Monte Carlo simulations to evaluate the entropy of rhombus tilings with fixed polygonal boundaries and 2D-fold rotational symmetry. We estimate the large-size limit of this entropy for D=4 to 10. We confirm analytic predictions of [N. Destainville et al., J. Stat....

  • Measuring the chemical potential of polymer solutions and melts in computer simulations. Müller, M.; Paul, W. // Journal of Chemical Physics;1/1/1994, Vol. 100 Issue 1, p719 

    We are presenting a new method to measure the chemical potential of polymer chains in computer simulations which has virtually no limitations in chain length. The method uses a thermodynamic integration over the excluded volume interaction of a ghost chain to determine its excess chemical...

  • One dimensional model for water and aqueous solutions. Part V. Monte Carlo simulation of dilute solutions of hard rod in waterlike particles. Lyubartsev, Alexander; Ben-Naim, Arieh // Journal of Chemical Physics;11/28/2009, Vol. 131 Issue 20, p204507 

    We have carried out Monte Carlo simulation on the primitive one dimensional model for water described earlier [A. Ben-Naim, J. Chem. Phys. 128, 024506 (2008)]. We show that by taking into account second nearest neighbor interactions, one can obtain the characteristic anomalous solvation...

  • Excess entropy and structural transitions in a two-dimensional square-shoulder fluid. Singh, Murari; Hongjun Liu; Kumar, Sanat K.; Ganguly, Abir; Chakravarty, Charusita // Journal of Chemical Physics;2/21/2010, Vol. 132 Issue 7, p074503 

    Metropolis Monte Carlo simulations on the square-shoulder fluid of Malescio and Pellicane are used to trace the temperature dependent excess entropy, the heat capacity, and configurational energy along several isochores, including those for which mechanically stable zero-temperature structures...

  • Polymer distribution in connected spherical domains. Sousa, A. F.; Pais, A. A. C. C.; Linse, P. // Journal of Chemical Physics;6/1/2005, Vol. 122 Issue 21, p214902 

    The distribution of neutral and charged polymers with different flexibilities between two spheres of varying volume connected by a short and narrow cylinder has been investigated by Monte Carlo simulations. The uncharged chain displayed mostly a single-sphere occupancy due to the high...

  • The configurational entropy of mixing of interstitials solid solutions. Garcés, Jorge // Applied Physics Letters;4/19/2010, Vol. 96 Issue 16, p161904 

    A simple and general expression for the configurational entropy of mixing of multicomponent interstitials solid solutions is deduced based on a probabilistic approach to compute the configurational entropy. The good agreement obtained with previous experimental and theoretical works based on...

  • Configurational entropy of hydrogen-disordered ice polymorphs. Herrero, Carlos P.; Ramírez, Rafael // Journal of Chemical Physics;6/21/2014, Vol. 140 Issue 23, p234502-1 

    The configurational entropy of several H-disordered ice polymorphs is calculated by means of a thermodynamic integration along a path between a totally H-disordered state and one fulfilling the Bernal-Fowler ice rules. A Monte Carlo procedure based on a simple energy model is used, so that the...


Read the Article


Sorry, but this item is not currently available from your library.

Try another library?
Sign out of this library

Other Topics