Thesis (M. Eng.)--Carleton University, 1996. / Also available in electronic format on the Internet.

Thesis (M. Eng.)--Carleton University, 2001. / Also available in electronic format on the Internet.

Thesis (Ph. D.)--Ohio State University, 2007. / Title from first page of PDF file. Includes bibliographical references (p. 100-110).

We present the application of the boundary integral equation method for solving the motion of biological cell or particle under Stokes flow in the presence of electrostatic field. The huge dense matrix-vector product from the boundary integral method poses a computationally challenging problem for solving the large system of equations generated. In our work, we used the precorrected-FFT (pFFT) method to reduce the computational time and memory usage drastically, so that large scale simulations can be performed quickly on a personal computer. Results on the force field acting on the particle, as well as the behavior of the particle through cell trap are presented. / Singapore-MIT Alliance (SMA)

Thesis (MEng)--Stellenbosch University, 2015. / ENGLISH ABSTRACT: The Method of Moments (MoM) is routinely used for the numerical solution of electromagnetic surface integral equations. Solution errors are inherent to any numerical computational method, and error estimators can be effectively employed to reduce and control these errors. In this thesis, gradient recovery techniques of the Finite Element Method (FEM) are formulated within the MoM context, in order to recover a higher-order charge of a Rao-Wilton-Glisson (RWG) MoM solution. Furthermore, a new recovery procedure, based specifically on the properties of the RWG basis functions, is introduced by the author. These recovered charge distributions are used for a posteriori error estimation of the charge. It was found that the newly proposed charge recovery method has the highest accuracy of the considered recovery methods, and is the most suited for applications within recovery based error estimation. In addition to charge recovery, the possibility of recovery procedures for the MoM solution current are also investigated. A technique is explored whereby a recovered charge is used to find a higher-order divergent current representation. Two newly developed methods for the subsequent recovery of the solenoidal current component, as contained in the RWG solution current, are also introduced by the author. A posteriori error estimation of the MoM current is accomplished through the use of the recovered current distributions. A mixed second-order recovered current, based on a vector recovery procedure, was found to produce the most accurate results. The error estimation techniques developed in this thesis could be incorporated into an adaptive solver scheme to optimise the solution accuracy relative to the computational cost. / AFRIKAANSE OPSOMMING: Die Moment Metode (MoM) vind algemene toepassing in die numeriese oplossing van elektromagnetiese oppervlak integraalvergelykings. Numeriese foute is inherent tot die prosedure: foutberamingstegnieke is dus nodig om die betrokke foute te analiseer en te reduseer. Gradiënt verhalingstegnieke van die Eindige Element Metode word in hierdie tesis in die MoM konteks geformuleer. Hierdie tegnieke word ingespan om die oppervlaklading van 'n Rao-Wilton-Glisson (RWG) MoM oplossing na 'n verbeterde hoër-orde voorstelling te neem. Verder is 'n nuwe lading verhalingstegniek deur die outeur voorgestel wat spesifiek op die eienskappe van die RWG basis funksies gebaseer is. Die verhaalde ladingsverspreidings is geïmplementeer in a posteriori fout beraming van die lading. Die nuut voorgestelde tegniek het die akkuraatste resultate gelewer, uit die groep verhalingstegnieke wat ondersoek is. Addisioneel tot ladingsverhaling, is die moontlikheid van MoM-stroom verhalingstegnieke ook ondersoek. 'n Metode vir die verhaling van 'n hoër-orde divergente stroom komponent, gebaseer op die verhaalde lading, is geïmplementeer. Verder is twee nuwe metodes vir die verhaling van die solenodiale komponent van die RWG stroom deur die outeur voorgestel. A posteriori foutberaming van die MoM-stroom is met behulp van die verhaalde stroom verspreidings gerealiseer, en daar is gevind dat 'n gemengde tweede-orde verhaalde stroom, gebaseer op 'n vektor metode, die beste resultate lewer. Die foutberamingstegnieke wat in hierdie tesis ondersoek is, kan in 'n aanpasbare skema opgeneem word om die akkuraatheid van 'n numeriese oplossing, relatief tot die berekeningskoste, te optimeer.

Este trabalho apresenta a formulação e implementação numérica do método dos elementos de contorno (MEC) para elasticidade linear tri-dimensional, com avaliação direta das integrais fracamente e fortemente singulares. A implementação segue a formulação tradicional do MEC direto, e a discretização do contorno das variáveis do problema é realizada com elementos descontínuos, permitindo o uso de malhas desconectadas ao longo das superfícies. O cálculo das integrais singulares é realizado através do uso de expansões assintóticas calculadas em torno de um ponto singular genérico. As expressões analíticas destas expansões são apresentadas no trabalho. Estas expansões serão subtraídas do núcleo original regularizando-o e a parte singular é integrada analiticamente, restando apenas uma integral regular, tornando ambas as integrais possíveis de serem calculadas com quadraturas de Gauss. É concluído que o presente método requer menos pontos de integração para o mesmo nível de erro quando comparado com outras técnicas. Alguns casos de elasticidade são resolvidos para ilustrar a eficiência e precisão do método. / This work presents the formulation and implementation of the boundary element method (BEM) to three dimensional linear elastostatics, with the direct evaluation of the strongly singular integral equations. The implementation follows the traditional direct BEM formulation, and the discretization of the boundary is carried out with discontinuous elements, enabling the use of disconnected meshes along the surfaces. The computation of the singular integral equations is accomplished by using the asymptotic expansions derived around a generic singular point. The analytical expressions for these expansions are presented in this work. The expansions are subtracted from the kernel to regularize it. This subtracted part is then added by computing a regular line integral along the boundary of the element. Both the integrals can be calculated with Gauss-type quadratures. It's observed that the present method needs less integration points for the same level of error when compared with other techniques. Several elasticity benchmarks are solved to demonstrate the eficiency and the accuracy of the present method.

Este trabalho apresenta a formulação e implementação numérica do método dos elementos de contorno (MEC) para elasticidade linear tri-dimensional, com avaliação direta das integrais fracamente e fortemente singulares. A implementação segue a formulação tradicional do MEC direto, e a discretização do contorno das variáveis do problema é realizada com elementos descontínuos, permitindo o uso de malhas desconectadas ao longo das superfícies. O cálculo das integrais singulares é realizado através do uso de expansões assintóticas calculadas em torno de um ponto singular genérico. As expressões analíticas destas expansões são apresentadas no trabalho. Estas expansões serão subtraídas do núcleo original regularizando-o e a parte singular é integrada analiticamente, restando apenas uma integral regular, tornando ambas as integrais possíveis de serem calculadas com quadraturas de Gauss. É concluído que o presente método requer menos pontos de integração para o mesmo nível de erro quando comparado com outras técnicas. Alguns casos de elasticidade são resolvidos para ilustrar a eficiência e precisão do método. / This work presents the formulation and implementation of the boundary element method (BEM) to three dimensional linear elastostatics, with the direct evaluation of the strongly singular integral equations. The implementation follows the traditional direct BEM formulation, and the discretization of the boundary is carried out with discontinuous elements, enabling the use of disconnected meshes along the surfaces. The computation of the singular integral equations is accomplished by using the asymptotic expansions derived around a generic singular point. The analytical expressions for these expansions are presented in this work. The expansions are subtracted from the kernel to regularize it. This subtracted part is then added by computing a regular line integral along the boundary of the element. Both the integrals can be calculated with Gauss-type quadratures. It's observed that the present method needs less integration points for the same level of error when compared with other techniques. Several elasticity benchmarks are solved to demonstrate the eficiency and the accuracy of the present method.

This thesis reports onreview and research work carried out on the numerical analysis of elastomers. The two numerical techniques investigated for this purpose are the finite and boundary element methods. The finite element method is studied so that existing theory is used to develop a finite element code both to review the finite element method as applied to the stress analysis of elastomers and to provide a comparison of results and numerical approach with the boundary element method. The research work supported on in this thesis covers the application of the boundary element method to the stress analysis of elastomers. To this end a simplified regularization approach is discussed for the removal of strong and hypersingularities generated in the system on non-linear boundary integral equations. The necessary programming details for the implementation of the boundary element method are discussed based on the code developed for this research. Both the finite and boundary element codes developed for this research use the Mooney-Rivlin material model as the strain energy based constitutive stress strain function. For validation purposes four test cases are investigated. These are the uni-axial patch test, pressurized thick wall cylinder, centrifugal loading of a rotating disk and the J-Integral evaluation for a centrally cracked plate. For the patch test and pressurized cylinder, both plane stress and strain have been investigated. For the centrifugal loading and centrally cracked plate test cases only plane stress has been investigated. For each test case the equivalent results for an equivalent FEM program mesh have been presented. The test results included in this thesis prove that the FE and BE derivations detailed in this work are correct. Specifically the simplified domain integral singular and hyper-singular regularization approach was shown to lead to accurate results for the test cases detailed. Various algorithm findings specific to the BEM implementation of the theory are also discussed.

Desenvolvem-se a formulação do Método dos Elementos de Contorno e correspondente algoritmo (para implementação em microcomputador) para a análise de propagação de fraturas em domínios bidimensionais. São utilizados elementos lineares isoparamétricos, tanto para discretizar o contorno quanto para simular a fratura. Os elementos de fratura são descontínuos. A formulação é baseada em equações integrais de tensões e de deslocamentos, onde o termo que considera tensões iniciais concentradas na linha de fratura é formulado a partir da definição de dipolos. O critério adotado é o modelo de fratura coesiva. Os termos singulares e hiper-singulares da formulação são tratados analiticamente e os termos quase-singulares são calculados através de um esquema numérico baseado na utilização de sub-elementos. Os valores dos dipolos são estimados ponto a ponto. Ao longo das fraturas, o valor máximo da tensão normal de tração permite definir novos elementos. As tensões de cisalhamento são removidas para manter a direção principal durante o processo. / The Boundary Element Method formulation and corresponding algorithrn (for microcomputer implementation) are developed for crack growth analysis in two-dimensional domains. Linear isoparametric elements are used to discretize both boundary and crack path. Fracture elements are assumed to be discontinuous. The formulation is based on stress and displacement integral equations, where the term that takes into account initial stresses concentrated along fracture line is formulated from dipoles definition. The coesive fracture rnodel is the criterium adopted. Singular and hipersingular formulation terms are anallitically treated and quasi-singular terms are computed by a numerical scheme based on element subdivision. Dipole values are estirnated point by point. Along fractures, the maximum normal tensile strenght is used to define new elements. Shear stresses are also removed to maintain the principal direction during the process.

Este trabalho trata da formulação do método dos elementos de contorno para a análise não linear de pavimentos de concreto armado. A teoria utilizada é a de Reissner, que mostrou-se eficiente tanto para placas esbeltas quanto para as moderadamente espessas. Considera-se a ocorrência de cargas concentradas, distribuídas em sub-regiões da placa e em linha. Admite-se também a possibilidade de um campo de momentos iniciais, que viabiliza o estudo da não linearidade física nos problemas. Foram utilizados campos de momentos iniciais aplicados apenas em pontos internos ao domínio. As integrais que envolvem as células de domínio foram modificadas, eliminando-se os núcleos complexos e as aproximações através de séries. Foi desenvolvida uma formulação para a análise de placas vinculadas a estruturas quaisquer em seu domínio, com o uso de cargas aplicadas incógnitas atuando como enrijecedores. O acoplamento MEC/MEF foi empregado utilizando-se modelos simples, porém robustos. O sistema de equações algébricas foi otimizado com a utilização da técnica dos mínimos quadrados. O concreto foi modelado adotando-se o modelo de dano de Mazars; para as armaduras um modelo elastoplástico uniaxial com endurecimento isótropo. A análise não linear do problema é efetuada utilizando-se procedimento incremental-iterativo. São apresentados alguns exemplos simples que mostram a precisão da técnica usada. / This work refers to the formulation of the boundary element method for non-linear analysis of building floor structures. The plate bending theory adopted to develop the work wad due to Reissner, which has demonstrated to be efficient for thick, moderated thick and thin plates. The kinds of load applied on the plate medium surface have been taken into account: concentrated loading, distributed over sub-domains; distributed along internal lines. The presence of initial moment fields convenient to model temperature effects and to be used to build up non-linear solutions has also been considered in the formulation. The domain integrals containing complex kernels to take into account the initial moment field influences were modified by introducing their primitive functions, avoiding therefore using series expansions. To integrate the initial moments fields only approximations based on internal nodal points were used. The resulting cell integrals have been transformed to the cell boundary which results into regular integral only. A boundary element formulation to treat structural system defined by combining plates with other structural element was developed, using interface force as unknowns. The BEM/FEM coupling developed to treat this case is simple but robust; only displacements have been coupled avoiding important singularities that may happen when coupling rotations. The resulting system of algebraic equations has been regularized by using the least square method. The concrete material was modeled by using the Mazar\'s damage model, while the steel reinforcement was assumed to behave as elastoplastic material with isotropic hardening. Finally, some examples are shown to illustrate the accuracy of the presented formulation and the numerical schemes proposed in this work.