he focus of this research has been devoted to study the interaction between two or more self-propelled toroidal swimmers in Stokes flow by applying the method of regularized Stokeslets and also study the effect of a nearby wall to the movement of a helical ring by using the method of regurlarized Stokeslets with images. In the study of the interaction between two or more toroidal swimmers, we interpret these as three-dimensional, zero Reynolds number analogues of finite vortex dipoles in an ideal fluid. Then, we examine the stability of relative equilibria that can form for these swimmers when they are initially placed in tandem or abreast. In addition, we examine the dynamics of the torus when a spherical cell body is placed at its center. This gives us an insight into the mechanical role of the transverse flagellum of dinoflagellates. Moreover, we show that the torus with a sphere moves more efficiently than one without. Lastly, we model the transverse flagellum of a dinoflagellate as a helical ring and study the effect of a nearby wall on its movement. The numerical results show that the wall baffles the movement of the helical ring, which is consistent with the phenomenon of sperm accumulation near surfaces. / acase@tulane.edu

Modeling of wave scattering from grating couplers has become increasingly important due to extensive recent research interest in the problem of plasmonic resonance. Computational algorithms which are specially used to model the problem of scattering from the grating surfaces suffer from several drawbacks such as accuracy, computational efficiency, and generality. To address the challenges of the previous methods, this work presents a novel hybrid Finite Element-Boundary Integral Method (FE-BIM) solution to the problem of scattering from grating surfaces consisting of finite or infinite array of two-dimensional cavities and holes in an infinite metallic walls covered with a stratified dielectric layer. To solve the scattering problem from finite number of cavities or holes engraved in a perfectly conducting screen (PEC), the solution region is divided into interior regions containing the cavities or holes and the region exterior to them. The finite element formulation is applied inside the interior region to derive a linear system of equations associated with nodal field values. Using two-boundary formulation, the surface integral equation employing free-space Green's function is then applied at \emph{only} the opening of the cavities or holes to truncate the computational domain and to connect the matrix subsystem generated from each cavity or hole. The hybrid FE-BIM method is extended to solve the scattering problem from an infinite array of cavities or holes in a PEC screen by deriving the quasi-periodic Green's function. In the scattering problem from an infinite array of cavities, the finite element formulation is first used inside a single cavity in the unit-cell. Next, the surface integral equation employing the quasi-periodic Green's function is applied at the opening of \emph{only} a single cavity as a boundary constraint to truncate the computational domain. Effect of the infinite array of cavities is incorporated into the system of the nodal equations by the quasi-periodic Green's function. Finally, the method based on the hybrid FE-BIM is developed to solve the scattering problem from grating surfaces covered with a stratified dielectric layer. In this method, the surface integral equation employing grounded dielectric slab Green's function is applied at the opening of the cavities or holes inside the dielectric coating to truncate the solution region efficiently. An accurate algorithm is presented to derive the grounded dielectric slab Green's function in spatial domain incorporating the effects of the surface-waves and leaky-waves excited and propagated inside the dielectric slab. Numerical examples of near and far field calculations for finite or infinite array of cavities or holes are presented to validate accuracy, versatility, and efficiency of the algorithm presented in this thesis.

In this dissertation, we introduce and analyse a scheme for the fast evaluation of integrals stemming from boundary element methods including discretisations of the classical single and double layer potential operators. Our method is based on the parametrisation of boundary elements in terms of a d-dimensional parameter tuple. We interpret the integral as a real-valued function f depending on d parameters and show that f is smooth in a d-dimensional box. A standard interpolation of f by polynomials leads to a d-dimensional tensor which is given by the values of f at the interpolation points. This tensor may be approximated in a low rank tensor format like the canonical format or the hierarchical format. The tensor approximation has to be done only once and allows us to evaluate interpolants in O(dr(m+1)) operations in the canonical format, or O(dk³ + dk(m + 1)) operations in the hierarchical format, where m denotes the interpolation order and the ranks r, k are small integers. In particular, we apply an efficient black box scheme in the hierarchical tensor format in order to adaptively approximate tensors even in high dimensions d with a prescribed (but heuristic) target accuracy. By means of detailed numerical experiments, we demonstrate that highly accurate integral values can be obtained at very moderate costs.

This study is to develop a 2D fully nonlinear numerical wave tank used to investigate the wave-induced dynamic properties of a dual pontoon floating structure (DPFS) with a liquid container on the top. The nonlinear numerical wave tank, developed based on the velocity potential function and the boundary element method (BEM), is to simulate dynamic properties including sway, heave, roll, and tension response. In addition, a physical model of the dual floating pontoon is tested in a hydrodynamic wave tank to validate the numerical model for simulation of wave and structure interaction. In the numerical model, a boundary integral equation method (BIEM) with linear element scheme is applied to establish a 2D fully nonlinear numerical wave tank (NWT). The nonlinear free surface condition is treated by combining the Mixed Eulerian and Lagrangian method (MEL), the fourth-order Runge-Kutta method (RK4) and a cubic spline scheme. The second-order Stokes wave theory is used to generate the velocity flux on the input boundary. Numerical damping zones are deployed at both ends of the NWT to dissipate or absorb the transmitted and reflected waves. Acceleration potential method and modal decomposition method are adopted to solve the unsteady potential functions £X1,t and £X2,t, while the system of motion equation is established according to Newton's 2nd law. Finally, the RK4 is applied to predict the motion of the platform, and the variation of free surface. As for the hydrodynamic laboratory model test, an image process scheme is applied to trace the floating structure motion and the variation of water surface inside the sloshing tank, while the mooring tension is measured by a load cell and stored in a data logger. The comparisons of numerical simulations and experimental data indicate that the numerical predictions are larger than measurements especially near the resonance frequency. This discrepancy is probably due to the fluid viscous effect. To overcome this problem and maintain the calculation efficiency, an uncoupled damping coefficient obtained through a damping ratio (£a=C/Ccr=0.02) is incorporated into the vibration system. Results reveal that responses of body motion near the resonant frequencies of each mode have significantly reduced and close to the measurements. Therefore, the suitable value of the damping ratio for the floating platform is £a=0.02. Then the numerical model with a damping ratio is applied to investigate the dynamic properties of the floating platform for different arrangements, including different mooring angle, spring constant, spacing, and the liquid container. Results demonstrate that the resonant frequency of each mode, responses of body motion and mooring tensions change along with the settings. As a whole, the platform with smaller mooring angle, longer spacing between the pontoons, higher water depth and wider width of the liquid container has relatively stable body motions and less mooring tension. Finally, the comparisons of the effects of random and regular waves on the floating structure illustrate that the variation of water surface in the liquid container is much severe in random waves than in regular waves such that the interaction between liquid and floating structure is more chaotic and thus reduces the amplitude of each response mode. As a result, the mooring tensions for random waves become much gentler than the regular waves. Key words: Boundary integral equation method; fully nonlinear numerical wave tank; dual pontoon floating structure

This paper is to investigate the wave-induced dynamic properties of the floating platform for cage aquaculture. Considering the calculation efficiency and its applicability, this problem is simplified by: (1) assuming the flow field is inviscid, incompressible and irrotational; (2) the form drag and inertia drag on the fish net is calculated by the modified Morison equation (or Morison type equation of relative motion), including the material and geometric properties; (3) the moorings is treated as a symmetric linear spring system and the influence of hydrodynamic forces on the mooring lines is neglected; and (4) the net-volume is assumed as un-deformable to avoid the inversely prolonging computing time because the mass of fish net with is too light comparing with the mass of floating platform and cause the marching time step tremendously small to reach the steady-state condition which may lead to larger numerical errors (e.g. truncation errors) in computation. The BIEM with linear element scheme is applied to establish a 2D fully nonlinear numerical wave tank (NWT). The nonlinear free surface condition is treated by combining the Mixed Eulerian and Lagrangian method (MEL), the fourth-order Runge-Kutta method (RK4) and the cubic spline scheme. The second-order Stokes wave theory is adopted to give the velocity on the input boundary. Numerical damping zones are deployed at both ends of the NWT to dissipate or absorb the transmitted and reflected wave energy. The velocity and acceleration fields should be solved simultaneously in order to obtain the wave-induced dynamic property of the floating platform. Thus, both the acceleration potential method and modal decomposition method are adopted to solve the wave forces on the floating body, while the wave forces on the fish net are calculated by the modified Morison equation. According to Newton¡¦s second law, the total forces on the gravity center of the floating platform form the equation of motion. Finally, the RK4 is applied to predict the displacement and velocity of the platform. Firstly, the NWT is validated by comparing the wave elevation, internal velocity and acceleration with those from the second-order Stokes wave theory. Moreover, the numerical damping zone is suitable for long time simulation with a wide range of wave depth. The simulated results on wave-body interactions of fixed or freely floating body also indicate good agreement with those of other published results. Secondly, in the case of the interaction of waves and the floating platform, the simulated results show well agreement with experimental data, except at the vicinity of resonant frequency of roll and heave motions. This discrepancy is due to the fluid viscous effect. To overcome this problem and maintain the calculation efficiency, an uncoupled damping coefficient obtained by a damping ratio (£i=0.1 ) is incorporated into the vibration system. Results reveal that responses of body motion near the resonant frequencies of each mode have significant reduction and close to the experimental data. Moreover, the results are also consistent well with experiments in different wave height, mooring angle, water depth either with or without fish net. Therefore, the suitable value of the damping ratio for the floating platform is £i=0.1. Finally, the present model is applied to investigate the dynamic properties of the floating platform under different draft, width, spacing, spring constant, mooring angle and depth of fish net. Results reveal that the resonant frequency and response of body motion, mooring force, reflection and transmission coefficients and wave energy will be changed. According to the resonant response, the platform with shallower draft, larger width, longer spacing between two pontoons, smaller spring constants, or deeper depth of fish net has more stable body motions and smaller mooring forces. Irregular wave cases are presented to illustrate the relationship with the regular wave cases. Results indicate that the dynamic responses of body motion and the reflection coefficient in irregular waves have similar trend with regular waves. However, in the irregular wave cases, the resonant frequency is moved to the higher frequency. Similarly, resonant response function is smaller but wider, which is due to the energy distribution in the wave spectrum.

Σκοπός της παρούσας διδακτορικής διατριβής είναι η ανάπτυξη αριθμητικής μεθόδου, η οποία επιλύει προβλήματα δισδιάστατης στατικής ελαστικότητας, καθώς και δυναμικής ελαστικότητας στο πεδίο των συχνοτήτων και στο πεδίο του χρόνου. Το κύριο χαρακτηριστικό της είναι ότι η προσέγγιση του άγνωστου πεδίου γίνεται με την τοποθέτηση σημείων και όχι με τη χρήση κάποιου πλέγματος όπως γίνεται στις μέχρι τώρα κλασικές μεθοδολογίες των πεπερασμένων ή συνοριακών στοιχείων. Μέρος της παρούσας διατριβής αποτελεί και η ανάπτυξη προγράμματος ηλεκτρονικού υπολογιστή, ο οποίος υποστηρίζει πλήρως τα όσα αναφέρονται στην παρούσα εργασία. Η παρούσα διατριβή αποτελείται από δύο ενότητες. Στην πρώτη ενότητα, η οποία περιλαμβάνει τα πρώτα τρία κεφάλαια, παρατίθεται το θεωρητικό υπόβαθρο της μεθοδολογίας. Στη δεύτερη ενότητα περιγράφονται διάφορες τεχνικές λεπτομέρειες, όπως ολοκληρώσεις και προσέγγιση πεδίου και δίνονται αρκετά παραδείγματα, τα οποία πιστοποιούν την ακρίβεια και την αξιοπιστία της. / -

Modeling of wave scattering from grating couplers has become increasingly important due to extensive recent research interest in the problem of plasmonic resonance. Computational algorithms which are specially used to model the problem of scattering from the grating surfaces suffer from several drawbacks such as accuracy, computational efficiency, and generality. To address the challenges of the previous methods, this work presents a novel hybrid Finite Element-Boundary Integral Method (FE-BIM) solution to the problem of scattering from grating surfaces consisting of finite or infinite array of two-dimensional cavities and holes in an infinite metallic walls covered with a stratified dielectric layer. To solve the scattering problem from finite number of cavities or holes engraved in a perfectly conducting screen (PEC), the solution region is divided into interior regions containing the cavities or holes and the region exterior to them. The finite element formulation is applied inside the interior region to derive a linear system of equations associated with nodal field values. Using two-boundary formulation, the surface integral equation employing free-space Green's function is then applied at \emph{only} the opening of the cavities or holes to truncate the computational domain and to connect the matrix subsystem generated from each cavity or hole. The hybrid FE-BIM method is extended to solve the scattering problem from an infinite array of cavities or holes in a PEC screen by deriving the quasi-periodic Green's function. In the scattering problem from an infinite array of cavities, the finite element formulation is first used inside a single cavity in the unit-cell. Next, the surface integral equation employing the quasi-periodic Green's function is applied at the opening of \emph{only} a single cavity as a boundary constraint to truncate the computational domain. Effect of the infinite array of cavities is incorporated into the system of the nodal equations by the quasi-periodic Green's function. Finally, the method based on the hybrid FE-BIM is developed to solve the scattering problem from grating surfaces covered with a stratified dielectric layer. In this method, the surface integral equation employing grounded dielectric slab Green's function is applied at the opening of the cavities or holes inside the dielectric coating to truncate the solution region efficiently. An accurate algorithm is presented to derive the grounded dielectric slab Green's function in spatial domain incorporating the effects of the surface-waves and leaky-waves excited and propagated inside the dielectric slab. Numerical examples of near and far field calculations for finite or infinite array of cavities or holes are presented to validate accuracy, versatility, and efficiency of the algorithm presented in this thesis.

Abstract The dissertation considers the time fractional diffusion equation (TFDE) with the Dirichlet boundary condition in the sub-diffusion case, i.e. the order of the time derivative is α ∈ (0,1). In the thesis we have studied the solvability of TFDE by the method of layer potentials. We have shown that both the single layer potential and the double layer potential approaches lead to integral equations which are uniquely solvable. The dissertation consists of four articles and a summary section. The first article presents the solution for the time fractional diffusion equation in terms of the single layer potential. In the second and third article we have studied the boundary behaviour of the layer potentials for TFDE. The fourth paper considers the spline collocation method to solve the boundary integral equation related to TFDE. In the summary part we have proved that TFDE has a unique solution and the solution is given by the double layer potential when the lateral boundary of a bounded domain admits C1 regularity. Also, we have proved that the solution depends continuously on the datum in the sense that a nontangential maximal function of the solution is norm bounded from above by the datum in L2(ΣT). If the datum belongs to the space H1,α/2(ΣT), we have proved that the nontangential function of the gradient of the solution is norm bounded from above by the datum in H1,α/2(ΣT).

Cette thèse porte sur la diffraction d’une onde plane électromagnétique par une surface lisse parfaitement conductrice (PEC). Elle présente l’analyse des propriétés d’une nouvelle formulation des trois principales équations intégrales de frontières de la théorie de la diffraction électromagnétique (EFIE, MFIE et CFIE). L’idée est d’adapter les équations intégrales conventionnelles à la diffraction d’une onde plane en supposant que la fonction de phase de l’onde plane incidente détermine la fonction de phase de la distribution de courant induit sur la surface.L’idée d’utiliser la phase dans la diffraction d’ondes planes a déjà été étudiée pour les hautes fréquences, notamment dans les thèses de Zhou (1995) et Darrigrand (2002) qui adaptèrent les espaces d’approximation des éléments finis. Dans cette thèse, cependant, nous suivons une formulation plus récente, donnée par Herberthson (2008), où la fonction de phase est incorporée dans la distribution du noyau des opérateurs intégraux.En présentant les versions modifiées de l’EFIE et de la MFIE (dénommées HEFIE et HMFIE)dans des espaces fonctionnels appropriés, nous prouvons ici l’existence d’une solution unique à cette formulation spécifique et présentons une mise en oeuvre pratique originale qui tire parti de l’expérience acquise sur l’EFIE/MFIE. Par la suite, nous explorons une propriété importante offerte par ces nouvelles formulations: la possibilité de réduire le nombre de degrés de liberté requis pour obtenir une solution précise du problème. / This thesis is about the scattering of an electromagnetic plane wave incidenton a perfectly conducting smooth surface. It presents the analysis of the properties of a newformulation of the three principal boundary integral equations of electromagnetic scattering theory(EFIE, MFIE and CFIE). The basic idea is to adapt the conventional integral equations toplane-wave scattering by supposing that the phase function of an incident plane wave determinesthe phase function of the induced boundary current distribution.This idea of using the phase in plane wave scattering has previously been studied in highfrequencyscattering, in particular in the theses by Zhou (1995) and Darrigrand (2002) whoadapt the finite element approximation spaces. In this thesis, though, we follow a more recentformulation, given by Herberthson (2008), where the phase function is incorporated in the kerneldistribution of the integral operators.Presenting the modified version of the EFIE and the MFIE (denoted HEFIE and HMFIE) inappropriate function spaces, we prove the existence of a unique solution to this specific formulationand developp an original practical implementation which takes advantage of the gainedexperience on the EFIE/MFIE. Then, we explore another important property provided by thenew formulations: the possibility to reduce the number of degrees of freedom required to get anaccurate solution of the problem.

Doctor of Philosophy / Department of Mechanical and Nuclear Engineering / Prakash Krishnaswami / Xiao J. Xin / Over the past two decades meshless methods have attracted much attention owing to their advantages in adaptivity, higher degree of solution field continuity, and capability to handle moving boundary and changing geometry. In this work, a meshless integral method based on the regularized boundary integral equation has been developed and applied to two-dimensional linear elasticity and elastoplasticity with small or large deformation. The development of the meshless integral method and its application to two-dimensional linear elasticity is described first. The governing integral equation is obtained from the weak form of elasticity over a local sub-domain, and the moving least-squares approximation is employed for meshless function approximation. This formulation incorporates: a subtraction method for singularity removal in the boundary integral equation, a special numerical integration for the calculation of integrals with weak singularity which further improves accuracy, a collocation method for the imposition of essential boundary conditions, and a method for incorporation of natural boundary conditions in the system governing equation. Next, elastoplastic material behavior with small deformation is introduced into the meshless integral method. The constitutive law is rate-independent flow theory based on von Mises yielding criterion with isotropic hardening. The method is then extended to large deformation plasticity based on Green-Naghdi’s theory using updated Lagrangian description. The Green-Lagrange strain is decomposed into the elastic and plastic part, and the elastoplastic constitutive law is employed that relates the Green-Lagrange strain to the second Piola-Kirchhoff stress. Finally, a pre- and post-processor for the meshless method using node- and pixel-based approach is presented. Numerical results from the meshless integral method agree well with available analytical solutions or finite element results, and the comparisons demonstrate that the meshless integral method is accurate and robust. This research lays the foundation for modeling and simulation of metal cutting processes.