Other Boundary Element Applications

Numerous problems in polymer processing have been solved in the past years with the use of the boundary element method. In all these solutions, the complexity of the geometry was the primary reason why the technique was used. Some of these problems are illustrated in this section. 

Gramann, Osswald and Rios [22,23,54] used BEM to simulate various mixing processes in two and three dimensions, an example of which is presented in Figs. 10.30 and 10.31. In both these systems the velocity and velocity gradients were computed as particles were tracked while traveling through the system. 

The velocity gradients were used to compute the rate of deformation tensor, the magnitudes of the rate of deformation and vorticity tensors. The magnitudes of the rate of 

Rios et al. [55] performed an experimental study with the above rhomboidal mixing section configurations using a 45 mm diameter single screw extruder. The mixing sections 

2Fitting a power law model to rheological measurements done on the HDPE resulted in a power law index, n, of 0.41 and a consistency index, m, of 16624 Pa-sn. 

---

Abstract This chapter reviews atomization modeling works that utilize boundary element methods (BEMs) to compute the transient surface evolution in capillary flows. The BEM, or boundary integral method, represents a class of schemes that incorporate a mesh that is only located on the boundaries of the domain and hence are attractive for free surface problems. Because both primary and secondary atomization phenomena are considered in many free surface problems, BEM is suitable to describe their physical processes and fundamental instabilities. Basic formulations of the BEM are outlined and their application to both low- and highspeed plain jets is presented. Other applications include the aerodynamic breakup of a drop, the pinch-off of an electrified jet, and the breakup of a drop colliding into a wall. 

In the application of the boundary element method, it is crucial to select appropriate botmdary surface for the solute cavity and to proceed as accurate as possible tessellation (triangulation) of this surface. For instance, it has been proposed that in the case of the cavity formation fi"om overlapping van-der-Waals spheres, the atomic van-der-Waals radii should be multiplied by a coefficient equal to 1.2. Other possibihties of the siuface definition include the closed envelope obtained by rolling a spherical probe of adequate diameter on the van-der-Waals surface of the solute molecule and the surface obtained ifom the positions of the center of such spherical probe aroimd the solute. 

Typical applications of the boundary element method in the context of adhesion technology are commonly found for the modeling of cracks (fi-acture mechanics) and other types of stress singularities, cf the bibliography of (Mackerle 1995a). The article by (Vable and Maddi 2010) addresses the specific problems (i.e., numerical modeling considerations which limited the application of BEM in the past) related to bonded joints and boundary element simulation. In addition, numerical results of lap joints, cf. O Fig. 26.18, with several spew angles were presented which demonstrate the potential of the boundary element method in analysis of bonded joints. 

In contrast with domain methods such as Finite Element Method (FEM) or Finite Differential Method (FDM), Boundary Element Method (BEM) discretises only the boundary of the domain. Because of this reduction of the dimensionality BEM was expected to be advantageous in large-scale problems. However, the application of this method has so far been limited to relatively small problems. This was because coefficient matrices in BEM are full and unsymmetrical, due to which both the operation count and the memory requirements for the matrix equation buildup are of the order of 0 N ), where N is the number of unlmowns. The operation count even increases to 0(A ) as one attempts to solve the matrix equations with conventional direct solvers. In particular, the full matrix property leads to a serious exhaustion of the memory of a computer and is an obstacle for applications of BEM to large-scale problems. On the other hand, coefficient matrices in domain methods are banded and both computational complexity for matrix buildup and memory requirements are 0(iV). 

The fundamental physical laws governing motion of and transfer to particles immersed in fluids are Newton s second law, the principle of conservation of mass, and the first law of thermodynamics. Application of these laws to an infinitesimal element of material or to an infinitesimal control volume leads to the Navier-Stokes, continuity, and energy equations. Exact analytical solutions to these equations have been derived only under restricted conditions. More usually, it is necessary to solve the equations numerically or to resort to approximate techniques where certain terms are omitted or modified in favor of those which are known to be more important. In other cases, the governing equations can do no more than suggest relevant dimensionless groups with which to correlate experimental data. Boundary conditions must also be specified carefully to solve the equations and these conditions are discussed below together with the equations themselves. 

Two major forms of the OCFE procedure are common and differ only in the trial functions used. One uses the Lagrangian functions and adds conditions to make the first derivatives continuous across the element boundaries, and the other uses Hermite polynomials, which automatically have continuous first derivatives between elements. Difficulties in the numerical integration of the resulting system of equations occur with the use of both types of trial functions, and personal preference must then dictate which is to be used. The final equations that need to be integrated after application of the OCFE method in the axial dimension to the reactor equations (radial collocation is performed using simple orthogonal collocation) can be expressed in the form... 

All numerical techniques require application of sampling theory. Briefly stated, one chooses a representative sample of points within the region of interest and at each point attempts to calculate iteratively the most accurate solution possible, guided by self-consistency of local solutions with each other and with the specified boundary conditions. We describe two seemingly contrasting techniques finite-difference and finite-element methods (1,2). 

Elemental iodine is a reactive gas, and the rate of uptake on certain surfaces is controlled by the rate of diffusion through the boundary layer over the surface. At some surfaces, tracer quantities of iodine are adsorbed irreversibly, at others reversibly. In most applications the amount of iodine on the surface is much less than a monolayer, and the equilibration between the adsorbed and airborne iodine cannot be considered in terms of vapour pressure. In 1949, experiments were started at Harwell, both in the wind tunnel and in the field, to study the deposition of elemental 131I vapour to surfaces. 

---