Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userConservative Integrators for Piecewise Smooth Systems with Transversal Dynamics
65
0
0.0
(
0
)
Added by
Anil Hirani
Publication date
2021
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We introduce conservative integrators for long term integration of piecewise smooth systems with transversal dynamics and piecewise smooth conserved quantities. In essence, for a piecewise dynamical system with piecewise defined conserved quantities such that its trajectories cross transversally to its interface, we combine Mannshardts transition scheme and the Discrete Multiplier Method to obtain conservative integrators capable of preserving conserved quantities up to machine precision and accuracy order. We prove that the order of accuracy of the integrators is preserved after crossing the discontinuity in the case of codimension one number of conserved quantities. Numerical examples illustrate the preservation of accuracy order.
rate research
Read More
The paper deals with numerical discretizations of separable nonlinear Hamiltonian systems with additive noise. For such problems, the expected value of the total energy, along the exact solution, drifts linearly with time. We present and analyze a time integrator having the same property for all times. Furthermore, strong and weak convergence of the numerical scheme along with efficient multilevel Monte Carlo estimators are studied. Finally, extensive numerical experiments illustrate the performance of the proposed numerical scheme.
We perform a numerical analysis of a class of randomly perturbed {H}amiltonian systems and {P}oisson systems. For the considered additive noise perturbation of such systems, we show the long time behavior of the energy and quadratic Casimirs for the exact solution. We then propose and analyze a drift-preserving splitting scheme for such problems with the following properties: exact drift preservation of energy and quadratic Casimirs, mean-square order of convergence one, weak order of convergence two. These properties are illustrated with numerical experiments.
An explicit numerical strategy that practically preserves invariants is derived for conservative systems by combining an explicit high-order Runge-Kutta (RK) scheme with a simple modification of the standard projection approach, which is named the explicit invariants-preserving (EIP) method. The proposed approach is shown to have the same order as the underlying RK method, while the error of invariants is analyzed in the order of $mathcal{O}left(h^{2(p+1)}right),$ where $h$ is the time step and $p$ represents the order of the method. When $p$ is appropriately large, the EIP method is practically invariants-conserving because the error of invariants can reach the machine accuracy. The method is illustrated for the cases of single and multiple invariants, with regard to both ODEs and high-dimensional PDEs. Extensive numerical experiments are presented to verify our theoretical results and demonstrate the superior behaviors of the proposed method in a long time numerical simulation. Numerical results suggest that the fourth-order EIP method preserves much better the qualitative properties of the flow than the standard fourth-order RK method and it is more efficient in practice than the fully implicit integrators.
This article presents and analyses an exponential integrator for the stochastic Manakov equation, a system arising in the study of pulse propagation in randomly birefringent optical fibers. We first prove that the strong order of the numerical approximation is $1/2$ if the nonlinear term in the system is globally Lipschitz-continuous. Then, we use this fact to prove that the exponential integrator has convergence order $1/2$ in probability and almost sure order $1/2$, in the case of the cubic nonlinear coupling which is relevant in optical fibers. Finally, we present several numerical experiments in order to support our theoretical findings and to illustrate the efficiency of the exponential integrator as well as a modified version of it.
We develop Random Batch Methods for interacting particle systems with large number of particles. These methods use small but random batches for particle interactions, thus the computational cost is reduced from $O(N^2)$ per time step to $O(N)$, for a system with $N$ particles with binary interactions. On one hand, these methods are efficient Asymptotic-Preserving schemes for the underlying particle systems, allowing $N$-independent time steps and also capture, in the $N to infty$ limit, the solution of the mean field limit which are nonlinear Fokker-Planck equations; on the other hand, the stochastic processes generated by the algorithms can also be regarded as new models for the underlying problems. For one of the methods, we give a particle number independent error estimate under some special interactions. Then, we apply these methods to some representative problems in mathematics, physics, social and data sciences, including the Dyson Brownian motion from random matrix theory, Thomsons problem, distribution of wealth, opinion dynamics and clustering. Numerical results show that the methods can capture both the transient solutions and the global equilibrium in these problems.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
