A New Software and Hardware Parallelized Floating Random-Walk Algorithm for the Modified Helmholtz Equation Subject to Neumann and Mixed Boundary Conditions
Keywords:
A New Software and Hardware Parallelized Floating Random-Walk Algorithm for the Modified Helmholtz Equation Subject to Neumann and Mixed Boundary ConditionsAbstract
A new floating random-walk algorithm for the one-dimensional modified Helmholtz equation subject to Neumann and mixed boundary conditions problems is developed in this paper. Traditional floating random-walk algorithms for Neumann and mixed boundary condition problems have involved “reflecting boundaries” resulting in relatively large computational times. In a recent paper, we proposed the elimination of the use of reflecting boundaries through the use of novel Green’s functions that mimic the boundary conditions of the problem of interest. The methodology was validated by a solution of the one-dimensional Laplace’s equation. In this paper, we extend the methodology to the floating random-walk solution of the onedimensional modified Helmholtz equation, and excellent agreement has been obtained between an analytical solution and floating random-walk results. The algorithm has been parallelized and a near linear rate of parallelization has been obtained with as many as thirty-two processors. These results have previously been published in [1]. In addition, a GPU implementation employing 4096 simultaneous threads displayed a similar near-linear parallelization gain and a one to two orders of magnitude improvement over the CPU implementation. An immediate application of this research is in the numerical solution of the electromagnetic diffusion equation in magnetically permeable and electrically conducting objects with applications in dielectrometry and magnetometry sensors that have the ability to detect sub-surface objects such as landmines. The ultimate goal of this research is the application of this methodology to the solution of aerodynamical flow problems.
Downloads
References
K. Chatterjee, M. Sandora, C. W. Yu, S.
Srinivasan, J. Poggie, “A New Parallelized
Floating Random-Walk Algorithm for the
Modified Helmholtz Equation Subject to
Neumann and Mixed Boundary Conditions:
Validation with a 1D Benchmark Problem,”
The 25th International Review of Progress in
Applied Computational Electromagnetics,
March 2009.
Y. L. Le Coz and R. B. Iverson, “A
Stochastic Algorithm for High Speed
Capacitance Extraction in Integrated
Circuits,” Solid-State Electronics, Vol. 35,
pp. 1005-1012, 1992.
M. Sobol, A Primer for the Monte Carlo
Method, CRC Press: Boca Raton, 1994.
R. Haberman, Elementary Applied Partial
Differential Equations, 3rd ed., Prentice-
Hall: New Jersey, 1998.
A. Haji-Sheikh, Application of Monte Carlo
Methods to Thermal Conduction Problems,
Ph.D. dissertation, University of Minnesota,
, pp. 106-108.
K. Chatterjee, C. Yu, S. Srinivasan and J.
Poggie, “A New Floating Random Walk
Methodology for Neumann and Mixed
Boundary Condition Problems Without
Reflections at Boundaries: Validation with
Laplace’s Equation in One Dimension,”
Far East Journal of Applied Mathematics,
Vol. 26, No. 3, pp. 705-713, 2007.
NVIDIA, 2007, Compute Unified Device
Architecture:
http://developer.download.nvidia.com/comp
ute/cuda/1_0/NVIDIA_CUDA_Programmin
g_Guide_1.0.pdf.
V. Podlozhnyuk, Parallel Mersenne Twister:
http://developer.download.nvidia.com/comp
ute/cuda/sdk/website/projects/MersenneTwi
ster/doc/MersenneTwister.pdf, 2007.