BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

L. F. Ricketson

L. F. Ricketson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NAphysics.plasm-phphysics.comp-ph

Trust snapshot

Quick read

Trust 19 - Baseline
5works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2015arXiv

A multilevel Monte Carlo method for a class of McKean-Vlasov processes

We generalize the multilevel Monte Carlo (MLMC) method of Giles to the simulation of systems of particles that interact via a mean field. When the number of particles is large, these systems are described by a McKean-Vlasov process - a stochastic differential equation (SDE) whose coefficients depend on expectations of the solution as well as pathwise data. In contrast to standard MLMC, the new method uses mean field estimates at coarse levels to inform the fine level computations. Using techniques from the theory of propagation of chaos, we prove convergence and complexity results for the algorithm in a special case. We find that the new method achieves $L^1$ error of size $\varepsilon$ with $O(\varepsilon^{-2} (\log \varepsilon)^5)$ complexity, in contrast to the $O(\varepsilon^{-3})$ complexity of standard methods. We also prove a variance scaling result that strongly suggests similar performance of the algorithm in a more general context. We present numerical examples from applications and observe the expected behavior in each case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate Derivative Evaluation for any Grad-Shafranov Solver

We present a numerical scheme that can be combined with any fixed boundary finite element based Poisson or Grad-Shafranov solver to compute the first and second partial derivatives of the solution to these equations with the same order of convergence as the solution itself. At the heart of our scheme is an efficient and accurate computation of the Dirichlet to Neumann map through the evaluation of a singular volume integral and the solution to a Fredholm integral equation of the second kind. Our numerical method is particularly useful for magnetic confinement fusion simulations, since it allows the evaluation of quantities such as the magnetic field, the parallel current density and the magnetic curvature with much higher accuracy than has been previously feasible on the affordable coarse grids that are usually implemented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An entropy based thermalization scheme for hybrid simulations of Coulomb collisions

We formulate and test a hybrid fluid-Monte Carlo scheme for the treatment of elastic collisions in gases and plasmas. While our primary focus and demonstrations of applicability are for moderately collisional plasmas, as described by the Landau-Fokker-Planck equation, the method is expected to be applicable also to collision processes described by the Boltzmann equation. This scheme is similar to the previously discussed velocity-based scheme [R. Caflisch, et. al, Multiscale Modeling & Simulation 7, 865, (2008)] and the scattering-angle-based scheme [A.M. Dimits, et. al, Bull. APS 55, no. 15 (2010, Abstract: XP9.00006)], but with a firmer theoretical basis and without the inherent limitation to the Landau-Fokker-Planck case. It gives a significant performance improvement (e.g., error for a given computational effort) over the velocity-based scheme. These features are achieved by assigning passive scalars to each simulated particle and tracking their evolution through collisions. The method permits a detailed error analysis that is confirmed by numerical results. The tests performed are for the evolution from anisotropic Maxwellian and a bump-on-tail distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Multilevel Monte Carlo simulation of Coulomb collisions

We present a new, for plasma physics, highly efficient multilevel Monte Carlo numerical method for simulating Coulomb collisions. The method separates and optimally minimizes the finite-timestep and finite-sampling errors inherent in the Langevin representation of the Landau-Fokker-Planck equation. It does so by combining multiple solutions to the underlying equations with varying numbers of timesteps. For a desired level of accuracy epsilon, the computational cost of the method is order(epsilon^{-2}) or order(epsilon^{-2} (\ln epsilon)^2), depending on the underlying discretization, Milstein or Euler-Maruyama respectively. This is to be contrasted with a cost of order(epsilon^{-3}) for direct simulation Monte Carlo or binary collision methods. We successfully demonstrate the method with a classic beam diffusion test case in 2D, making use of the Levy area approximation for the correlated Milstein cross terms, and generating a computational saving of a factor of 100 for epsilon = 10^{-5}. We discuss the importance of the method for problems in which collisions constitute the computational rate limiting step, and its limitations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Three Improvements to Multi-Level Monte Carlo Simulation of SDE Systems

We introduce three related but distinct improvements to multilevel Monte Carlo (MLMC) methods for the solution of systems of stochastic differential equations (SDEs). Firstly, we show that when the payoff function is twice continuously differentiable, the computational cost of the scheme can be dramatically reduced using a technique we call `Ito linearization'. Secondly, by again using Ito's lemma, we introduce an alternative to the antithetic method of Giles et. al [M.B. Giles, L. Szpruch. arXiv preprint arXiv:1202.6283, 2012] that uses an approximate version of the Milstein discretization requiring no Levy area simulation to obtain the theoretically optimal cost-to-error scaling. Thirdly, we generalize the antithetic method of Giles to arbitrary refinement factors. We present numerical results and compare the relative strengths of various MLMC-type methods, including each of those presented here.

0 citations0 reviewsNo code linkedOpen access