Researcher profile

R. E. Caflisch

R. E. Caflisch contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.plasm-ph physics.comp-ph physics.flu-dyn math-ph math.DS math.MP math.NA nlin.SI physics.class-ph q-fin.PR

Trust snapshot

Quick read

Trust 19 - Baseline

5works

0followers

10topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

preprint2015arXiv

Complex singularities and PDEs

In this paper we give a review on the computational methods used to characterize the complex singularities developed by some relevant PDEs. We begin by reviewing the singularity tracking method based on the analysis of the Fourier spectrum. We then introduce other methods generally used to detect the hidden singularities. In particular we show some applications of the Padé approximation, of the Kida method, and of Borel-Polya method. We apply these techniques to the study of the singularity formation of some nonlinear dispersive and dissipative one dimensional PDE of the 2D Prandtl equation, of the 2D KP equation, and to Navier-Stokes equation for high Reynolds number incompressible flows in the case of interaction with rigid boundaries.

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.

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.

preprint2012arXiv

European Option Pricing with Transaction Costs and Stochastic Volatility: an Asymptotic Analysis

In this paper the valuation problem of a European call option in presence of both stochastic volatility and transaction costs is considered. In the limit of small transaction costs and fast mean reversion, an asymptotic expression for the option price is obtained. While the dominant term in the expansion it is shown to be the classical Black and Scholes solution, the correction terms appear at $O(\varepsilon^{1/2})$ and $O(\varepsilon)$. The optimal hedging strategy is then explicitly obtained for the Scott's model.

preprint2011arXiv

Beyond the Child-Langmuir Limit

This paper describes a new solution formulation for fully nonlinear and unsteady planar flow of an electron beam in a diode. Using characteristic variables - i.e., variables that follow particle paths - the solution is expressed through an exact analytic, but implicit, formula for any choice of incoming velocity $v_0$, electric field $E_0$ and current $J_0$. For steady solutions, this approach clarifies the origin of the maximal current $J_max$, derived by Child and Langmuir for $v_0=0$ and by Jaffe for $v_0>0$. The implicit formulation is used to find (1) unsteady solutions having constant incoming flux $J_0>J_max$, which leads formation of a virtual cathode, and (2) time-periodic solutions whose average flux exceeds the adiabatic average of $J_max$.