Partitioning discrete boundary

The electrostatic properties of the molecule may be used as a criterion for judging the MEM enhancement. Using the uniform prior density, the MEM molecular dipole moment derived by the discrete boundary partitioning of space (chapter 6) is only 1.3 D, compared with 9.1 D based on the experimental density,... 

Figure 6.1 shows the stockholder decomposition of the theoretical deformation density of the cyanoacetylene molecule, H—Cs=C—C=N (Hirshfeld 1977b). The overlap density in the bonds is distributed between the bonded atoms. The assignment of part of the density near the hydrogen nucleus to the adjacent carbon atom manifests the difference between fuzzy and discrete boundary partitioning methods. 

TABLE 6.1 Net Charges q (e) from the Stockholder Partitioning. Charges in the Second Row are from a Discrete Boundary Partitioning by Politzer (1971) and Politzer and Reggio (1972)... 

If the ratio on the left is the smallest, the point i belongs to atom A, and thus to molecule /, and vice versa. The discrete boundary Van-der- Waals-ratio partition-... 

FIG. 6.3 Definition of vectors used in discrete boundary space partitioning. 

When space is partitioned with discrete boundaries, as in Eq. (6.7) and in the Bader virial partitioning method, the moments can be derived directly from the structure factors by a modified Fourier summation, as described for the net charge in chapter 6. 

The analysis of polymer processing is reduced to the balance equations, mass or continuity, energy, momentum and species and to some constitutive equations such as viscosity models, thermal conductivity models, etc. Our main interest is to solve this coupled nonlinear system of equations as accurately as possible with the least amount of computational effort. In order to do this, we simplify the geometry, we apply boundary and initial conditions, we make some physical simplifications and finally we chose an appropriate constitutive equations for the problem. At the end, we will arrive at a mathematical formulation for the problem represented by a certain function, say / (x, T, p, u,...), valid for a domain V. Due to the fact that it is impossible to obtain an exact solution over the entire domain, we must introduce discretization, for example, a grid. The grid is just a domain partition, such as points for finite difference methods, or elements for finite elements. Independent of whether the domain is divided into elements or points, the solution of the problem is always reduced to a discreet solution of the problem variables at the points or nodal pointsinxxnodes. The choice of grid, i.e., type of element, number of points or nodes, directly affects the solution of the problem. 

In the previous contributions to this book, it has been shown that by adopting a polarizable continuum description of the solvent, the solute-solvent electrostatic interactions can be described in terms of a solvent reaction potential, Va expressed as the electrostatic interaction between an apparent surface charge (ASC) density a on the cavity surface which describes the solvent polarization in the presence of the solute nuclei and electrons. In the computational practice a boundary-element method (BEM) is applied by partitioning the cavity surface into Nts discrete elements and by replacing the apparent surface charge density cr by a collection of point charges qk, placed at the centre of each element sk. We thus obtain ... 

Pratt and co-workers have proposed a quasichemical theory [118-122] in which the solvent is partitioned into inner-shell and outer-shell domains with the outer shell treated by a continuum electrostatic method. The cluster-continuum model, mixed discrete-continuum models, and the quasichemical theory are essentially three different names for the same approach to the problem [123], The quasichemical theory, the cluster-continuum model, other mixed discrete-continuum approaches, and the use of geometry-dependent atomic surface tensions provide different ways to account for the fact that the solvent does not retain its bulk properties right up to the solute-solvent boundary. Experience has shown that deviations from bulk behavior are mainly localized in the first solvation shell. Although these first-solvation-shell effects are sometimes classified into cavitation energy, dispersion, hydrophobic effects, hydrogen bonding, repulsion, and so forth, they clearly must also include the fact that the local dielectric constant (to the extent that such a quantity may even be defined) of the solvent is different near the solute than in the bulk (or near a different kind of solute or near a different part of the same solute). Furthermore... 

This approach, which is based on a discretization of a on the cavity surface S, is often referred to as the Boundary Element Method (BEM). There is an abundant literature about BEM in the field of applied engineering and physics (Banerjee and Butterfield 1981 Beskos 1987). In this application of the BEM approach we need to specify how the surface S is defined and partitioned in tesserae, and how Ve is computed. 

For the given scheme of partition of the machining zone boundary on the elements, the discretization of the boundary integral equation and the boundary conditions is performed. The set of nonlinear equations, which is obtained by discretization, is solved by Newton s method. As a result, the distribution of the current density over the WP surface is obtained (Fig. 10c). 

The finite element method is a systematic procedure of approximating continuous functions as discrete models. This discretization involves finite number of points and subdomains in the problem s domain. The values of the given function are held at the points, so-called nodes. The non-overlapping subdomains, so-called finite elements, are connected together at nodes on their boundaries and hold piecewise and local approximations of the function, which are uniquely defined in terms of values held at their nodes. The collection of discretized elements and nodes is called the mesh and the process of its construction is called meshing. A typical finite element partition of a two-dimensional domain with triangular finite elements is given in Fig. 1. 

Because soil does not contain discrete objects with obvious boundaries that could be called individual pores, the precise delineation of a pore unavoidably requires artificial, subjectively established distinctions. This contrasts with soil particles, which are easily defined, being discrete material objects with obvious boundaries. The arbitrary criterion required to partition pore space into individual pores is often not explicitly stated when pores or their sizes are discussed. Because of this inherent arbitrariness, some scientists argue that the concepts of pore and pore size should be avoided. Much valuable theory of the behavior of the soil-water-air system, however, has been built on these concepts, defined using widely, if not universally, accepted criteria. 

---