Boundary general methods

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]. 

Because of this, there is a real need for designing the general method, by means of which economical schemes can be created for equations with variable and even discontinuous coefhcients as well as for quasilinear non-stationary equations in complex domains of arbitrary shape and dimension. As a matter of experience, the universal tool in such obstacles is the method of summarized approximation, the framework of which will be explained a little later on the basis of the heat conduction equation in an arbitrary domain G of the dimension p with the boundary F... 

To compute each of the n(ct ), one can generalize the methods used to compute ihG- Hence, the most elegant method would be to use basis functions that satisfy the boundary conditions of Eq. (43), if this were practical to implement. A more general method would be to extend the Mead-Truhlar vector-potential approach [6]. This approach would involve carrying out h calculations, each including a... 

Variational electrostatic projection method. In some instances, the calculation of PMF profiles in multiple dimensions for complex chemical reactions might not be feasible using full periodic simulation with explicit waters and ions even with the linear-scaling QM/MM-Ewald method [67], To remedy this, we have developed a variational electrostatic projection (VEP) method [75] to use as a generalized solvent boundary potential in QM/MM simulations with stochastic boundaries. The method is similar in spirit to that of Roux and co-workers [76-78], which has been recently... 

Banavali, N.K. Im, W. Roux, B., Electrostatic free energy calculations using the generalized solvent boundary potential method, J. Chem. Phys. 2002,117,7381-7388. 

Which of the various immersed or embedded boundary methods is best— generally or for a particular case—is still an open question. Thornock and Smith (2005) introduced a Cell Adjusted Boundary Force Method for a stirred vessel. All methods proposed so far have their own pros and cons. Immersed boundary methods are also exploited in LB techniques (e.g., Derksen and Van den Akker, 1999). Rohde (2004) investigated the use of triangular facets for representing a spherical particle. 

Solution proceeds via modification of the convective-diffusion equation to include the relevant kinetic terms and new boundary conditions. Accounts of the general methods for obtaining kinetic currents, first developed for the DME, can be found in, for example, refs. 182-184. 

Summations over point-, line-, or planar-source solutions are useful examples of the general method of Green s functions. 6 For instance, the boundary and initial conditions for a triangular source are... 

In the computational practice, the ASC density is discretized into a collection of point charges qk, spread on the cavity surface. The apparent charges are then determined by solving the electrostatic Poisson equation using a Boundary Element Method scheme (BEM) [1], Many BEM schemes have been proposed, being the most general one known as integral equation formalism (IEFPCM) [10]. 

Boundary Cells. A cell on the boundary has one wall in contact with the gutside world and hence the flux through that wall is K.M and not J. Generally we want to replace J1 with KM over that fraction of the cell surface in common with the outer boundary. One method for doing so is to replace J in (66) by (K in those parts of the d-dimensional sphere surrounding point r that stick outside the cell aggregate volume. Thus (67,68) take the same form for the boundary points except that F takes on the form, for d = 3,... 

We must consider the laminar and turbulent portions of the boundary layer separately because the recovery factors, and hence the adiabatic wall temperatures, used to establish the heat flow will be different for each flow regime. It turns out that the difference is rather small in this problem, but we shall follow a procedure which would be used if the difference were appreciable, so that the general method of solution may be indicated. The free-stream acoustic velocity is calculated from... 

In one of the first articles on this subject [8], the general analytical solution of Eq. (3) was derived. This general solution is easy to find, but it contains infinite series and (integration) constants that depend on the boundary conditions. Those were determined for the central cells of square and triangular arrays, using the boundary collocation method [8]. More recent publications on this subject are based mostly on complete numerical solution using finite-element methods. 

After finalizing the model equations and boundary conditions, the next task is to choose a suitable method to approximate the differential equations by a system of algebraic equations in terms of the variables at some discrete locations in space and time (called a discretization method). There are many such methods the most important are finite difference (FD), finite volume (FV) and finite element (FE) methods. Other methods, such as spectral methods, boundary element methods or cellular automata are used, but these are generally restricted to special classes of problems. All methods yield the same solution if the grid (number of discrete locations used to... 

A general method of solving Eqs. 7.3.1-7.3.3 is described in Chapter 11, where we complicate matters somewhat by including energy transfer across the phase boundary. 

In general, the problem just defined is nonlinear, in spite of the fact that the governing, creeping-flow equations are linear. This is because the drop shape is unknown and dependent on the pressure and stresses, which in turn, depend on the flow. Thus n and F are also unknown functions of the flow field, and the boundary conditions (2-112), (2-122), (2-141), and (8-58) are therefore nonlinear. Thus, for arbitrary Ca, for which the deformation may be quite significant, the problem can be solved only numerically. Later in this chapter, we briefly discuss a method, known as the boundary Integral method, that may be used to carry out such numerical calculations. Here, however, we consider the limiting case... 

In addition to the recent book by Pozrikidis (Ref. 7), a good general reference to the boundary-integral method is S. Weinbaum, P. Ganatos, and Z. Y. Yan, Numerical multipole and boundary integral equation techniques in stokes flow, Annu. Rev. Fluid Mech. 22, 275-316 (1990). 

---