This Article 
 Bibliographic References 
 Add to: 
A Fast Method for Estimating Discrete Field Values in Early Engineering Design
March 1996 (vol. 2 no. 1)
pp. 35-43

Abstract—Much of the analysis done in engineering design involves the solution of partial differential equations (PDEs) that are subject to initial-value or boundary-value conditions; generically these are called "field problems." Finite-element and finite-difference methods (FEM, FDM) are the predominant solution techniques, but these are often too expensive or too tedious to use in the early phases of design. What's needed is a fast method to compute estimates of field values at a few critical points that uses simple and robust geometric tools. This paper describes such a method. It is based on an old technique—integrating PDEs through stochastic (Monte Carlo) sampling—that is accelerated through the use of ray representations. In the first (pre-processing) stage, the domain (generally a mechanical part) is coherently sampled to produce a ray-rep. The second stage involves the usual stochastic sampling of the field, which is now enhanced by exploiting the semi-discrete character of ray-reps. The method is relatively insensitive to the complexity of the shape being analyzed, and it has adjustable precision. Its mechanics and advantages are illustrated by using Laplace's equation as an example.

[1] R.B. Tilove,"Set Membership Classification: A Unified Approach to Geometric Intersection Problems," IEEE Trans. on Computers, vol. 29, no. 20, pp. 874-883, Oct. 1980.
[2] T.J. Hoffman and N.E. Banks,"Monte Carlo Solution to the Dirichlet Problem with the Double-Layer Potential Density," Trans. Am. Nucl. Soc., vol. 18, pp. 136-137, 1974.
[3] T.J. Hoffman and N.E. Banks,"Extension of the Surface Density Monte Carlo Technique to Heat Transfer in Nonconvex Regions," Trans. Am. Nucl. Soc., vol. 19, p. 164, 1974.
[4] O.D. Kellogg,Foundations of Potential Theory. Dover Publications, Inc., 1929.
[5] J. Arvo and D. Kirk,"A Survey of Ray Tracing Acceleration Techniques," ACM SIGGRAPHE '88 Tutorial Notes, vol. 7, New York, 1988.
[6] A. Fujimoto, T. Takayu, and K. Iwata, "ARTS: Accelerated Ray-Tracing System," IEEE Computer Graphics and Applications, Vol. 6, No. 4, April 1986, pp. 16-26.
[7] J.G. Cleary and G. Wyvill,"Analysis of an Algorithm for Fast Ray-Tracing Using Uniform Space Subdivision," Visual Computer, vol. 4, pp. 65-83, 1988.
[8] J.L. Ellis,G. Kedem,T.C. Lyerly,D.G. Thielman,R.J. Marisa,J.P. Menon, and H.B. Voelcker,"The RayCasting Engine and Ray Representations: A Technical Summary," Int'l J. Computational Geometry and Applications, vol. 1, no. 4, pp. 347-380, Dec. 1991.
[9] J.P. Menon,R.J. Marisa, and J. Zagajac,"More Powerful Solid Modeling Through Ray Representations," IEEE Computer Graphics and Applications, vol. 14, no. 3, pp. 22-35, May 1994.
[10] J. Zagajac,"Distributed Raycasting on the SP2," Tech. Report CPA95-5, Sibley School of Mechanical Eng., Cornell University, Mar. 1996.
[11] J.E. Bresenham,"Algorithm for Computer Control of a Digital Plotter," IBM Syst. J., vol. 4, no. 1, pp. 25-30, 1965.

Index Terms:
Ray tracing, Monte Carlo, random walk, engineering analysis.
Jovan Zagajac, "A Fast Method for Estimating Discrete Field Values in Early Engineering Design," IEEE Transactions on Visualization and Computer Graphics, vol. 2, no. 1, pp. 35-43, March 1996, doi:10.1109/2945.489385
Usage of this product signifies your acceptance of the Terms of Use.