Poisson elements

How can Equation (11.79) be solved Before computers were available only simple ihapes could be considered. For example, proteins were modelled as spheres or ellipses Tanford-Kirkwood theory) DNA as a uniformly charged cylinder and membranes as planes (Gouy-Chapman theory). With computers, numerical approaches can be used to solve the Poisson-Boltzmann equation. A variety of numerical methods can be employed, including finite element and boundary element methods, but we will restrict our discussion to the finite difference method first introduced for proteins by Warwicker and Watson [Warwicker and Watson 1982]. Several groups have implemented this method here we concentrate on the work of Honig s group, whose DelPhi program has been widely used. 

The Poisson equation has been used for both molecular mechanics and quantum mechanical descriptions of solvation. It can be solved directly using numerical differential equation methods, such as the finite element or finite difference methods, but these calculations can be CPU-intensive. A more efficient quantum mechanical formulation is referred to as a self-consistent reaction field calculation (SCRF) as described below. 

The continuum treatment of electrostatics can also model salt effects by generalizing the Poisson equation (12) to the Poisson-Boltzmann equation. The finite difference approach to solving Eq. (12) extends naturally to treating the Poisson-Boltzmann equation [21], and the boundary element method can be extended as well [19]. 

When inserting into (4.5), the term ZeR will be multiplied with the elements of the electric field gradient tensor V. Fortunately, the procedure can be restricted to diagonal elements Vu, because V is symmetric and, consequently, a principal axes system exists in which the nondiagonal elements vanish, = 0. The diagonal elements can be determined by using Poisson s differential equation for the electronic potential at point r = 0 with charge density (0), AV = Anp, which yields... 

The trace vanishes because only p- and /-electrons contribute to the EFG, which have zero probability of presence at r = 0 (i.e. Laplace s equation applies as opposed to Poisson s equation, because the nucleus is external to the EFG-generating part of the electronic charge distribution). As the EFG tensor is symmetric, it can be diagonalized by rotation to a principal axes system (PAS) for which the off-diagonal elements vanish, = 0. By convention, the principal axes are chosen such that... 

Holst, M.J. Baker, N.A. Wang, F., Adaptive multilevel finite element solution of the Poisson-Boltzmann equation I algorithms and examples, J. Comp. Chem. 2000, 21, 1319-1342... 

Two remaining problems relating to the treatment of solvation include the slowness of Poisson-Boltzmann calculations, when these are used to treat electrostatic effects, and the difficulty of keeping buried, explicit solvent in equilibrium with the external solvent when, e.g., there are changes in nearby solute groups in an alchemical simulation. Faster methods for solving the Poisson-Boltzmann equation by means of parallel finite element techniques are becoming available, however.22 24... 

The outflow of a CSTR is a Poisson process, i.e., fluid elements are randomly selected regardless of theirposition in the reactor. The waiting time before selection for a Poisson process has an exponential probability distribution. 

An early continuum treatment of solvation, associated with Born,17 comes out of the analysis of the electrostatic work involved in building up a charge Q on a conducting sphere of radius R in a medium with dielectric constant e. From Poisson s equation, it follows that the potential outside of the sphere is Q/eR. Thus the work of charging is the result of each additional element dq interacting with the charge q already present 87... 

Figure 15.2. Region of interest for computing potential based on Laplace or Poisson equations, where (a) a complete rectangular grid is established to cover the region, which may be adapted to finite-difference techniques using (b) a five-point method, or (c) a finite-element approach based on sampling functions.

The Physical Properties are listed next. Under this loose term a wide range of properties, including mechanical, electrical and magnetic properties of elements are presented. Such properties include color, odor, taste, refractive index, crystal structure, allotropic forms (if any), hardness, density, melting point, boiling point, vapor pressure, critical constants (temperature, pressure and vol-ume/density), electrical resistivity, viscosity, surface tension. Young s modulus, shear modulus, Poisson s ratio, magnetic susceptibility and the thermal neutron cross section data for many elements. Also, solubilities in water, acids, alkalies, and salt solutions (in certain cases) are presented in this section. 

Nonetheless, Pernety becomes once again a useful textual source allowing for a Poisson-like decipherment of Duchamp s texte obscure. As we are now privileged to learn from Pemety, originally it was the Alchemists who routinely practised a decomposition of elements, or basic forms. Decomposition," as designated by Pernety, involves... 

Perhaps the most widely used scheme for SCRF implementations of the Poisson equation is the surface area boundary element approach. This was first formalized by Miertus, Scrocco, and Tomasi (1981), and these authors referred to their construction as the polarized continuum model (PCM). While that name continues to find ample use in the literature, MST (the initials... 

When the Poisson equation is solved using a boundary element approach, the charges on the tesselated molecular surface are determined so that they provide an equivalent representation... 

Electrokinetic phenomena can be understood with the help of two equations The known Poisson equation and the Navier3-Stokes4 equation. The Navier-Stokes equation describes the movement of a Newtonian liquid, i.e., a liquid whose viscosity does not change when it flows and when it is sheared. In order to make the equation plausible we consider an infinitesimal quantity of the liquid having a volume dV = dx dy dz and a mass dm. If we want to write Newtons equation of motion for this volume element we have to consider three forces ... 

Their mutual Poisson brackets calculated using (34) yield a non-singular matrix with elements... 

Once a finite element formulation has been implemented in conjunction with a specific element type — either 1D, 2D or 3D — the task left is to numerically implement the technique and develop the computer program to solve for the unknown primary variables — in this case temperature. Equation (9.19) is a form that becomes very familiar to the person developing finite element models. In fact, for most problems that are governed by Poisson s equation, problems solving displacement fields in stress-strain problems and flow problems such as those encountered in polymer processing, the finite element equation system takes the form presented in eqn. (9.19). This equation is always re-written in the form... 

In this section, we will proceed to develop a finite element formulation for the two-dimensional Poisson s equation using a linear displacement, constant strain triangle. Poisson s equation has many applications in polymer processing, such as injection and compression mold filling, die flow, potential problems, heat transfer, etc. The general form of Poisson s equation in two-dimensions is... 

The previous section used the constant strain three-noded element to solve Poisson s equation with steady-state as well as transient terms. The same problems, as well as any field problems such as stress-strain and the flow momentum balance, can be formulated using isoparametric elements. With this type of element, the same (as the name suggests) shape functions used to represent the field variables are used to interpolate between the nodal coordinates and to transform from the xy coordinate system to a local element coordinate system. The first step is to discretize the domain presented in Fig. 9.12 using the isoparametric quadrilateral elements as shown in Fig. 9.15. 

The finite element expression for a problem using isoparametric finite elements will be similar to the ones developed in previous examples. For example, a problem governed by Poisson s equation will result in the following finite element equation,... 

Solution of the two-dimensional Poisson s equation compression molding. To illustrate the use of the four-noded isoparametric element, we can solve for the pressure distribution and velocity field during compression molding of an L-shaped polymer charge, shown in Fig. 9.18, with the physical and numerical data presented in Table 9.3. 

The boundary conditions are defined in the same way as with the flow analysis network. The nodes whose control volumes are empty or partially filled are assigned a zero pressure, and the gate nodes are either assigned an injection pressure or an injection volume flow rate. Just as is the case with flow analysis network, a mass balance about each nodal control volume will lead to a linear set of algebraic equations, identical to the set finite element formulation of Poisson s or Laplace s equation. The mass balance (volume balance for incompressible fluids) is given by... 

In the Alchemical Theory, says Albert Poisson, the four Elements, not any more than the three Principles, represent particular substances they are simply states of matter, simple modalities. Water is synonymous with the liquid state, Earth with the solid Air with the gaseous and Fire with that of a very subtle gaseous state, such as a gas expanded by the action of heat.. . Moreover, Elements represent, by extension, physical qualities such as heat, (Fire) dryness and solidity, (Earth) moisture and fluidity, (Water) cold and subtility, (Air) Zosimus gives to their ensemble the name of Tetrasomy. 

---