A diffusion Monte Carlo method for charge density on a conducting surface at non-constant potentials

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Unjong Yu ◽  
Hoseung Jang ◽  
Chi-Ok Hwang

Abstract We develop a last-passage Monte Carlo algorithm on a conducting surface at non-constant potentials. In the previous researches, last-passage Monte Carlo algorithms on conducting surfaces with a constant potential have been developed for charge density at a specific point or on a finite region and a hybrid BIE-WOS algorithm for charge density on a conducting surface at non-constant potentials. In the hybrid BIE-WOS algorithm, they used a deterministic method for the contribution from the lower non-constant potential surface. In this paper, we modify the hybrid BIE-WOS algorithm to a last-passage Monte Carlo algorithm on a conducting surface at non-constant potentials, where we can avoid the singularities on the non-constant potential surface very naturally. We demonstrate the last-passage Monte Carlo algorithm for charge densities on a circular disk and the four rectangle plates with a simple voltage distribution, and update the corner singularities on the unit square plate and cube.

2020 ◽  
Vol 3 (8) ◽  
pp. 2000075
Author(s):  
Hoseung Jang ◽  
Unjong Yu ◽  
Youngjoo Chung ◽  
Chi‐Ok Hwang

2005 ◽  
Vol 08 (05) ◽  
pp. 553-574 ◽  
Author(s):  
SEBASTIAN E. FERRANDO ◽  
ARIEL J. BERNAL

A new simulation based algorithm to approximate prices of path dependent European options is introduced. The algorithm is defined for tree-like approximations to the underlying process and makes extensive use of structural properties of the discrete approximation. We indicate the advantages of the new algorithm in comparison to standard Monte Carlo algorithms. In particular, we prove a probabilistic error bound that compares the quality of both approximations. The algorithm is of general applicability and, for a large class of options, it has the same computational complexity as Monte Carlo.


Author(s):  
Andrey Yu. Ambos ◽  
Galiya Lotova ◽  
Guennady Mikhailov

AbstractA Monte Carlo algorithm admitting parallelization is constructed for estimation of probability moments of the spectral radius of the operator of the integral equation describing transfer of particles with multiplication in a random medium. A randomized homogenization method is developed with the same aim on the base of the theory of small perturbations and diffusive approximation. Test calculations performed for a one-group spherically symmetric model system have shown a satisfactory concordance of results obtained from two models.


1961 ◽  
Vol 57 (3) ◽  
pp. 623-627 ◽  
Author(s):  
W. D. Collins

Copson(1) has shown that the problem of determining the surface density of electric charge induced on a thin circular disk maintained at a constant potential in an external electrostatic field can be solved by two applications of the known solution of Abel's integral equation. This note shows that Copson's method can be extended to determine the surface charge density induced on the disk when situated inside an earthed coaxial infinitely long hollow cylinder. Whilst this problem is not solved in closed form, it is shown to be governed by a Fredholm integral equation of the second kind, which can be solved approximately by iteration when the radius of the cylinder is large compared with that of the disk.


Biometrika ◽  
2020 ◽  
Vol 107 (4) ◽  
pp. 997-1004
Author(s):  
Qifan Song ◽  
Yan Sun ◽  
Mao Ye ◽  
Faming Liang

Summary Stochastic gradient Markov chain Monte Carlo algorithms have received much attention in Bayesian computing for big data problems, but they are only applicable to a small class of problems for which the parameter space has a fixed dimension and the log-posterior density is differentiable with respect to the parameters. This paper proposes an extended stochastic gradient Markov chain Monte Carlo algorithm which, by introducing appropriate latent variables, can be applied to more general large-scale Bayesian computing problems, such as those involving dimension jumping and missing data. Numerical studies show that the proposed algorithm is highly scalable and much more efficient than traditional Markov chain Monte Carlo algorithms.


2005 ◽  
Vol 8 ◽  
pp. 217-250
Author(s):  
Sophie Ambrose ◽  
Max Neunhöffer ◽  
Cheryl E. Praeger ◽  
Csaba Schneider

AbstractThis paper presents a generalisation of the sifting procedure introduced originally by Sims for computation with finite permutation groups, and now used for many computational procedures for groups, such as membership testing and finding group orders. The new procedure is a Monte Carlo algorithm, and it is presented and analysed in the context of black-box groups. It is based on a chain of subsets instead of a subgroup chain. Two general versions of the procedure are worked out in detail, and applications are given for membership tests for several of the sporadic simple groups. The authors' major objective was that the procedures could be proved to be Monte Carlo algorithms, and the costs computed. In addition, they explicitly determined suitable subset chains for six of the sporadic groups, and then implemented the algorithms involving these chains in the GAP computational algebra system. It turns out that sample imple-mentations perform well in practice. The implementations will be made available publicly in the form of a GAP package.


2002 ◽  
Vol 13 (03) ◽  
pp. 343-365 ◽  
Author(s):  
TETSUYA TAKAISHI ◽  
PHILIPPE de FORCRAND

We discuss hybrid Monte Carlo algorithms for odd-flavor lattice QCD simulations. The algorithms include a polynomial approximation, which enables us to simulate odd-flavor QCD in the framework of the hybrid Monte Carlo algorithm. In order to make the algorithms exact, the correction factor to the polynomial approximation is also included in an economical, stochastic way. We test the algorithms for nf = 1, 1 + 1 and 2 + 1 flavors and compare results with other algorithms.


Sign in / Sign up

Export Citation Format

Share Document