Source author record

Yajuan Sun

Yajuan Sun appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

Hamiltonian Particle-in-Cell methods for Vlasov-Poisson equations

In this paper, Particle-in-Cell algorithms for the Vlasov-Poisson system are presented based on its Poisson bracket structure. The Poisson equation is solved by finite element methods, in which the appropriate finite element spaces are taken to guarantee that the semi-discretized system possesses a well defined discrete Poisson bracket structure. Then, splitting methods are applied to the semi-discretized system by decomposing the Hamiltonian function. The resulting discretizations are proved to be Poisson bracket preserving. Moreover, the conservative quantities of the system are also well preserved. In numerical experiments, we use the presented numerical methods to simulate various physical phenomena. Due to the huge computational effort of the practical computations, we employ the strategy of parallel computing. The numerical results verify the efficiency of the new derived numerical discretizations.

preprint2016arXiv

Explicit high-order non-canonical symplectic particle-in-cell algorithms for Vlasov-Maxwell systems

Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure derived from the Lagrangian of the particle-field system, which is naturally discrete in particles. The electromagnetic field is spatially-discretized using the method of discrete exterior calculus with high-order interpolating differential forms for a cubic grid. The resulting time-domain Lagrangian assumes a non-canonical symplectic structure. It is also gauge invariant and conserves charge. The system is then solved using a splitting method discovered by He et al., which produces five exactly-soluable sub-systems, and high-order structure- preserving algorithms follow by combinations. The explicit, high-order, and conservative nature of the algorithms is especially suitable for long-term simulations of particle-field systems with extremely large number of degrees of freedom on massively parallel supercomputers. The algorithms have been tested and verified by the two physics problems, i.e., the nonlinear Landau damping and the electron Bernstein wave.

preprint2016arXiv

Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations

In this paper, we develop Hamiltonian particle-in-cell methods for Vlasov-Maxwell equations by applying conforming finite element methods in space and splitting methods in time. For the spatial discretisation, the criteria for choosing finite element spaces are presented such that the semi-discrete system possesses a discrete non-canonical Poisson structure. We apply a Hamiltonian splitting method to the semi-discrete system in time, then the resulting algorithm is Poisson preserving and explicit. The conservative properties of the algorithm guarantee the efficient and accurate numerical simulation of the Vlasov-Maxwell equations over long-time.

preprint2016arXiv

High order volume-preserving algorithms for relativistic charged particles in general electromagnetic fields

We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general construction of volume-preserving methods that can be applied to systems with time-dependent electromagnetic fields. The newly derived methods provide numerical solutions with good accuracy and conservative properties over long time of simulation. Furthermore, because of the use of processing technique the high order methods are explicit, and cost less than the methods derived from standard compositions, thus are more efficient. The results are verified by the numerical experiments. Linear stability analysis of the methods show that the high order processed method allows larger time step size during integration.

preprint2015arXiv

Canonical symplectic particle-in-cell method for long-term large-scale simulations of the Vlasov-Maxwell system

Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretizing its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g., $10^{9}$, degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves.

preprint2015arXiv

Explicit non-canonical symplectic algorithms for charged particle dynamics

We study the non-canonical symplectic structure, or K-symplectic structure inherited by the charged particle dynamics. Based on the splitting technique, we construct non-canonical symplectic methods which is explicit and stable for the long-term simulation. The key point of splitting is to decompose the Hamiltonian as four parts, so that the resulting four subsystems have the same structure and can be solved exactly. This guarantees the K-symplectic preservation of the numerical methods constructed by composing the exact solutions of the subsystems. The error convergency of numerical solutions is analyzed by means of the Darboux transformation. The numerical experiment display the long-term stability and efficiency for these methods.

preprint2015arXiv

Hamiltonian integration methods for Vlasov-Maxwell equations

Hamiltonian integration methods for the Vlasov-Maxwell equations are developed by a Hamiltonian splitting technique. The Hamiltonian functional is split into five parts, i.e., the electrical energy, the magnetic energy, and the kinetic energy in three Cartesian components. Each of the subsystems is a Hamiltonian system with respect to the Morrison-Marsden-Weinstein Poisson bracket and can be solved exactly. Compositions of the exact solutions yield Poisson structure preserving, or Hamiltonian, integration methods for the Vlasov-Maxwell equations, which have superior long-term fidelity and accuracy.

preprint2013arXiv

Enhanced HBVMs for the numerical solution of Hamiltonian problems with multiple invariants

Recently, the class of energy-conserving Runge-Kutta methods named Hamiltonian Boundary Value Methods (HBVMs), has been proposed for the efficient solution of Hamiltonian problems, as well as for other types of conservative problems. In this paper, we report further advances concerning such methods, resulting in their enhanced version (Enhanced HBVMs, or EHBVMs). The basic theoretical results are sketched, along with a few numerical tests on a Hamiltonian problem, taken from the literature, possessing multiple invariants.

preprint2013arXiv

Multiple invariants conserving Runge-Kutta type methods for Hamiltonian problems

In a recent series of papers, the class of energy-conserving Runge-Kutta methods named Hamiltonian BVMs (HBVMs) has been defined and studied. Such methods have been further generalized for the efficient solution of general conservative problems, thus providing the class of Line Integral Methods (LIMs). In this paper we derive a further extension, which we name Enhanced Line Integral Methods (ELIMs), more tailored for Hamiltonian problems, allowing for the conservation of multiple invariants of the continuous dynamical system. The analysis of the methods is fully carried out and some numerical tests are reported, in order to confirm the theoretical achievements.

preprint2012arXiv

An Input-Output Simulation Approach to Controlling Multi-AffineSystems for Linear Temporal Logic Specifications

This paper presents an input-output simulation approach to controlling multi-affine systems for linear temporal logic (LTL) specifications, which consists of the following steps. First, we partition the state space into rectangles, each of which satisfies atomic LTL propositions. Then, we study the control of multi-affine systems on rectangles including the control of driving all trajectories starting from a rectangle to exit through a facet and the control of stabilizing the system towards a desired point. With the proposed controllers, a finitely abstracted transition system is constructed which is shown to be input-output simulated by the rectangular transition system of the multi-affine system. Since input-output simulation preserves LTL properties, the controller synthesis of the multi-affine system for LTL specifications is achieved by designing a nonblocking supervisor for the abstracted transition system and by continuously implementing the resulting supervisor for the original multi-affine system.

preprint2012arXiv

Bisimilarity Enforcing Supervisory Control for Deterministic Specifications

This paper investigates the supervisory control of nondeterministic discrete event systems to enforce bisimilarity with respect to deterministic specifications. A notion of synchronous simulation-based controllability is introduced as a necessary and sufficient condition for the existence of a bisimilarity enforcing supervisor, and a polynomial algorithm is developed to verify such a condition. When the existence condition holds, a supervisor achieving bisimulation equivalence is constructed. Furthermore, when the existence condition does not hold, two different methods are provided for synthesizing maximal permissive sub-specifications.

preprint2012arXiv

The minimal stage, energy preserving Runge-Kutta method for polynomial Hamiltonian systems is the Averaged Vector Field method

No Runge-Kutta method can be energy preserving for all Hamiltonian systems. But for problems in which the Hamiltonian is a polynomial, the Averaged Vector Field (AVF) method can be interpreted as a Runge-Kutta method whose weights $b_i$ and abscissae $c_i$ represent a quadrature rule of degree at least that of the Hamiltonian. We prove that when the number of stages is minimal, the Runge-Kutta scheme must in fact be identical to the AVF scheme.

preprint2011arXiv

Computation for Supremal Simulation-Based Controllable and Strong Observable Subautomata

Bisimulation relation has been successfully applied to computer science and control theory. In our previous work, simulation-based controllability and simulation-based observability are proposed, under which the existence of bisimilarity supervisor is guaranteed. However, a given specification automaton may not satisfy these conditions, and a natural question is how to compute a maximum permissive subspecification. This paper aims to answer this question and investigate the computation of the supremal simulation-based controllable and strong observable subautomata with respect to given specifications by the lattice theory. In order to achieve the supremal solution, three monotone operators, namely simulation operator, controllable operator and strong observable operator, are proposed upon the established complete lattice. Then, inequalities based on these operators are formulated, whose solution is the simulation-based controllable and strong observable set. In particular, a sufficient condition is presented to guarantee the existence of the supremal simulation-based controllable and strong observable subautomata. Furthermore, an algorithm is proposed to compute such subautomata.

preprint2011arXiv

Decentralized Supervisory Control of Discrete Event Systems for Bisimulation Equivalence

In decentralized systems, branching behaviors naturally arise due to communication, unmodeled dynamics and system abstraction, which can not be adequately captured by the traditional sequencing-based language equivalence. As a finer behavior equivalence than language equivalence, bisimulation not only allows the full set of branching behaviors but also explicitly specifies the properties in terms of temporal logic such as CTL* and mu-calculus. This observation motivates us to consider the decentralized control of discrete event systems (DESs) for bisimulation equivalence in this paper, where the plant and the specification are taken to be nondeterministic and the supervisor is taken to be deterministic. An automata-based control framework is formalized, upon which we develop three architectures with respect to different decision fusion rules for the decentralized bisimilarity control, named a conjunctive architecture, a disjunctive architecture and a general architecture. Under theses three architectures, necessary and sufficient conditions for the existence of decentralized bisimilarity supervisors are derived respectively, which extend the traditional results of supervisory control from language equivalence to bisimulation equivalence. It is shown that these conditions can be verified with exponential complexity. Furthermore, the synthesis of bisimilarity supervisors is presented when the existence condition holds.