Other Numerical Methods

The extended finite element method (XFEM) for treating fracture in composite materials is proposed by Huynh and Belytschko [191]. This methods work with meshes that are independent of matrix/inclusion interfaces and the discontinuities and neartip enrichments were modeled. In order to describe the geometry of the interfaces and cracks in this method, level sets were employed, so that there is no need for explicit representation of either the cracks or the material interfaces. The other researchers such as Du et al. [192] and Ying et al. [193] used XFEM to model material interfaces in particulate composites with more in detail. 

On the other hand, some other alternative approaches have been developed to model the fracture as the XFEM [194] or the embedded localization method [195], These approaches can also be combined with a cohesive zone model [196]. 

---

The multiple use of logarithms in the analysis presented by Fig. 4.9 obliterates much of the deviation between theory and experiment. More stringent tests can be performed by other numerical methods. 

For most numerically solved models, a control-volume approach is used. This approach is based on dividing the modeling domain into a mesh. Between mesh points, there are finite elements or boxes. Using Taylor series expansions, the governing equations are cast in finite-difference form. Next, the equations for the two half-boxes on either side of a mesh point are set equal to each other hence, mass is rigorously conserved. This approach requires that all vectors be defined at half-mesh points, all scalars at full-mesh points, and all reaction rates at quarter-mesh points. The exact details of the numerical methods can be found elsewhere (for example, see ref 273) and are not the purview of this review article. The above approach is essentially the same as that used in CFD packages (e.g.. Fluent) or discussed in Appendix C of ref 139 and is related to other numerical methods applied to fuel-cell modeling. ... 

Fractionation Data and Distribution Analysis of HEC After One Day of Cellulase Attack. In Table III, the fractionation data of HEC are given after one day of enzymic hydrolysis. As after one hour of enzymic hydrolysis, no theoretical distribution function accorded well with the fractionation data, but we evaluated the parameters by numerical analysis, using the Gauss-Laguerre method (43,44). This method has one advantage over other numerical methods, e.g., (45)—all the calculations involved can be done manually without the need of high-speed computers. 

Do a literature search and find alternative ways of simulating non-linear equations using BEM. How does the technique compare to other numerical methods. 

Systems of linear equations are to be solved during the stoichiometric analysis of a reaction, but also on the occasion of the use of other numerical methods (see below). 

The disadvantages of finite-difference and other numerical methods are that quite a considerable amourit of computational effort is usually required in order to obtain the solution and that they do not, in general, reveal certain unifying features of the solutions, such as the fact that the profiles are similar under certain conditions. The widespread availability of modem computer facilities has, however, made these disadvantages relatively unimportant. 

In general, the nnmber of experimental points m is larger (often mnch larger) than the nnmber of adjustable parameters. The resulting function is not an exact fit to each point bnt represents a best overall fit in the sense that the snm of the sqnares of the deviations [y,(obs) ycaic( i)] is a minimnm, hence the name. The fitting function y(x) obtained by least sqnares can be differentiated, integrated, set eqnal to zero and solved to find roots, etc., by analytical methods. Thns there is no need for the apphcation of other numerical methods imless the function is very awkward to handle analytically in closed form. 

One limitation to this method (and to all other numerical methods for solving nonlinear equations) is that once you have found one solution, you cannot be sure that there are no additional solutions. The way to determine the existence of multiple roots is to evaluate f x) over a wide range of x values and find the intervals in which f x) changes sign (see the second... 

Runge-Kutta-Gill method provides an efficient algorithm for solving a system of first-order differential equations and makes use of much less computer memory when compared with other numerical methods. 

Equation (478) is the exact analytical geometry expression of capillary rise in a cylindrical tube having a circular cross section, which considers the deviation of the meniscus from sphericity, so that the curvature corresponds to (AP = A pgy) at each point on the meniscus, where y is the elevation of that point above the flat liquid level (y = z+h). Unfortunately, this relation cannot be solved analytically. Numerous approximate solutions have been offered, such as application of the Bashforth and Adams tables in 1883 (see Equation (476)) derivation of Equation (332) by Lord Rayleigh in 1915 a polynomial fit by Lane in 1973 (see Equation (482)) and other numerical methods using computers in modern times. 

There are also other numerical methods for solving differential equations, which we do not discuss. The numerical methods can be extended to sets of simultaneous differential equations such as occur in the analysis of chemical reaction mechanisms. Many of these sets of equations have a property called stiffness that makes them difficult to treat numerically. Techniques have been devised to handle this problem, which is beyond the scope of this book. ... 

Solutions of the equations appropriate for a given reaction scheme are obtained by (a) approximation methods, (b) Laplace transform or related techniques to yield closed form solutions, (c) digital simulation methods, and (d) other numerical methods. Approximation methods, such as those based on the reaction layer concept as described in Section... 

The solution of this equation would then follow as described in Section 9.4. The modifications of the mass-transfer equations for the different cases generally follow those used in voltammetric methods as shown in Table 12.2.1. Appropriate dimensionless parameters are listed in Table 12.3.1. It is usually not possible to solve these equations analytically, so various approximations (e.g., the reaction layer approach, as described in Section 1.5.2), digital simulations, or other numerical methods must be employed. The behavior of systems at the RDE can be analyzed by means of the zone diagrams employed for voltammetry (Section 12.3) by redefining the parameter A. This is accomplished by re-... 

For details about other numerical methods, see Britz s book [18]. 

C. Other Numerical Methods for Unsteady-State Conduction... 

Linear equations result naturally when we conduct material and energy balances, but most applications occur when we implement other numerical methods. One of the most basic solutirm techniques for systems such as Equation 9.15 is Gaussian elimination [3,5,9,13,14], which is illustrated using the System of Equations 9.15. 

There are many numerical approaches one can use to approximate the solution to the initial and boundary value problem presented by a parabolic partial differential equation. However, our discussion will focus on three approaches an explicit finite difference method, an implicit finite difference method, and the so-called numerical method of lines. These approaches, as well as other numerical methods for aU types of partial differential equations, can be found in the literature [5,9,18,22,25,28-33]. 

SchieBl et al. (2012) discussed the performance of the singularity method with regard to numeric stability and agreement to other numerical methods. There are apparently two limits ... 

Simultaneous linear equations occur in various engineering problems. The reader knows that a given system of linear equations can be solved by Cramer s rule or by the matrix method. However, these methods become tedious for large systems. However, there exist other numerical methods of solution which are well suited for computing machines. The following is an example. 

Two apparent drawbacks of the localized-orbital approach exist. First, the convergence of the total energy stalls at a value above the exact numerical solution depending on the size of the cutoff radius for the orbitals. This makes sense because physical information is lost in the truncation. Second, the convergence rate appears to slow somewhat with increasing system size. The observed convergence rates are still good, and competitive with other numerical methods, but this slowdown does not fit with standard MG orthodoxy—the... 

Modem finite element analysis or other numerical methods have no problem in treating non-linear behavior. Our physical understanding of material behavior at such levels is lacking, however, and effective numerical analysis depends to a large extent on the experimental determination of these properties. Despite these limitations, many researchers have shown that elastic analyses of many adhesive systems can be very Informative and useful. A number of adhesive systems are sufficiently linear, such that it is adequate to lump the plastic deformation and other dissipative mechanisms at the crack tip into the adhesive fracture energy (critical energy release rate) term. 

---