T h i s paper presents a method for optimum placement of nodes in grid generation for process simulation together with an algorithm for updating the grid after each addition of a new node to ensure that the Delaunay property is satisfied. Placement of nodes is decided on by considering the optimum error in evaluating the integral sv C(z)dz. The best error estimate is obtained when the node coincides with the centroid (the center of mass) of its own Voronoi region and moreover when the Voronoi region is symmetric. After addition of a node the grid is updated to maintain the Delaunay property using the Delaunay-Voronoi algorithm. I t ro duct ionIn process simulation, we are concerned with solving the diffusion equation _ -a ' -V.(DVC), t 2 0, z E Q c Rd, d = 2,3 8t where D is the diffusion coefficient. In solving the problem numerically, we are concerned with the error in computing the integral sv C(z)& where V is the region of influence of the node Z (the Voronoi region (see [4]) of the node Z.) Previously, we have considered the problem of node removal in a grid where mesh is too fine for the problem, (see [I], [2], and [3].) The removal algorithm identified nodes to be removed by their low discretizaton error. We now turn our attention to addition of nodes. When a grid is adapted, new nodes are added in areas where, due to larger variation in C(Z), it is deemed that greater resolution is needed. To decide on the optimum placement of these new nodes, we consider which position w i l l give the best error in evaluating Jv C ( z ) d x .The best error for this calculation is O(h3k, hk3) when C(z) is approximated by its Taylor series evaluated at the node f , the generator of the Voronoi region V, and moreover, when the node is the centroid (center of mass) of its own Voronoi region. Such a Voronoi tesselation is called a centroidal Voronoi tesselation. (see [5]). After addition of a new node, the grid quality must be maintained. A desirable property of the grid is that it has the Delaunay property. The grid is updated, after addition of each new node, to ensure that it has the Delaunay property. The algorithm to achieve this is called the Delaunay-Voronoi algorithm.The actual position of a new node is at the centroid of surrounding nodes to ensure that the node is at the centroid of its own Voronoi region. This position is decided on by wnsidering the smallest order error in evaluating the integral Jv C(z)dz using a Taylor series as stated a h . Hence the error in evaluating Jv C(z)dz is used as an error indicator in deciding on the positioning of nodes.In adding nodes during grid adaption, a check is used whereby a new node is added if and only if the geometric grid quality is i m p r d . If the geometric grid quality is not improved, then the new node is not added. Error Indicator for Evaluation of s , C(z)dzThe best error for this calculation of JvC(z)dz.is O(h3k, hk3), when C(z) is approximated by its Taylor series evaluated at the node Z , the generator of the Voronoi region, V, and moreover, when the node is t...