The importance of plates, with discontinuities in boundary supports in aeronautical and marine structures, have led to various techniques to solve plate problems with mixed edge support conditions. The domain decomposition method is one of the most effective of these techniques, providing accurate numerical solutions. This method is used to investigate the vibration and buckling of flat, isotropic, thin and elastic plates with mixed edge support conditions. Two practical approaches have been developed as an extension of the domain decomposition method, namely, the primary-secondary domain (PSD) approach and the line-domains (LD) approach. The PSD approach decomposes a plate into one primary domain and one/two secondary domain(s). The LD approach considers interconnecting boundaries as dominant domains whose basic functions take a higher edge restraint from the neighbouring edges. Convergence and comparison studies are carried out on a number of selected rectangular plate cases. Extensive practical plate problems with various shapes, combinations of mixed boundary conditions and different inplane loading conditions have been solved by the PSD and LD approaches. / Master of Engineering (Hons)

<p>In this work, we proposed a spectral integral method (SIM)-spectral element method (SEM)- finite element method (FEM) domain decomposition method (DDM) for solving inhomogeneous multi-scale problems. The proposed SIM-SEM-FEM domain decomposition algorithm can efficiently handle problems with multi-scale structures, </p><p>by using FEM to model electrically small sub-domains and using SEM to model electrically large and smooth sub-domains. The SIM is utilized as an efficient boundary condition. This combination can reduce the total number of elements used in solving multi-scale problems, thus it is more efficient than conventional FEM or conventional FEM domain decomposition method. Another merit of the proposed method is that it is capable of handling arbitrary non-conforming elements. Both geometry modeling and mesh generation are totally independent for different sub-domains, thus the geometry modeling and mesh generation are highly flexible for the proposed SEM-FEM domain decomposition method. As a result, the proposed SIM-SEM-FEM DDM algorithm is very suitable for solving inhomogeneous multi-scale problems.</p> / Dissertation

Cette thèse s’intéresse à la simulation acoustique de problèmes de grandes tailles. La parallélisation des méthodes numériques d’acoustique est le sujet principal de cette étude. Le manuscrit est composé de trois parties : lancé de rayon, méthodes de décomposition de domaines et algorithmes asynchrones. / This thesis studies numerical methods for large-scale acoustic problems. The parallelization of the numerical acoustic methods is the main focus. The manuscript is composed of three parts: ray-tracing, optimized interface conditions for domain decomposition methods and asynchronous iterative algorithms.

Cette thèse est consacrée à l'étude de l'influence d'une fine couche hétérogène sur le comportement élastique linéaire d'une structure tridimensionnelle.Deux types d'hétérogénéités sont pris en compte : des cavités et des inclusions élastiques. Une étude complémentaire, dans le cas d'inclusions de grande rigidité, a été réalisée en considérant un problème de conduction thermique.Une analyse formelle par la méthode des développements asymptotiques raccordés conduit à un problème d'interface qui caractérise le comportement macroscopique de la structure. Le comportement microscopique de la couche est lui déterminé sur une cellule de base. Le modèle asymptotique obtenu est ensuite implémenté dans un code éléments finis. Une étude numérique permet de valider les résultats de l'analyse asymptotique. / This thesis is devoted to the study of the influence of a thin heterogeneous layeron the linear elastic behavior of a three-dimensional structure. Two types of heterogeneties are considered : cavities and elastic inclusions. For inclusions of high rigidty a further study was performed in the case of a heat conduction problem.A formal analysis using the matched asymptotic expansions method leads to an interface problem which characterizes the macroscopic behavior of the structure. The microscopic behavior of the layer is determined in a basic cell.The asymptotic model obtained is then implemented in a finite element software.A numerical study is used to validate the results of the asymptotic analysis.

We are interested in numerical methods for the optimization constrained second order ordinary differential equations arising in biofilm modelling. This class of problems is challenging for several reasons. One of the reasons is that the underlying solution has a steep slope, making it difficult to resolve. We propose a new numerical method with techniques such as domain decomposition and asynchronous iterations for solving certain types of ordinary differential equations more efficiently. In fact, for our class of problems after applying the techniques of domain decomposition with overlap we are able to solve the ordinary differential equations with a steep slope on a larger domain than previously possible. After applying asynchronous iteration techniques, we are able to solve the problem with less time.~We provide theoretical conditions for the convergence of each of the techniques. The other reason is that the second order ordinary differential equations are coupled with an optimization problem, which can be viewed as the constraints. We propose a numerical method for solving the coupled problem and show that it converges under certain conditions. An application of the proposed methods on biofilm modeling is discussed. The numerical method proposed is adopted to solve the biofilm problem, and we are able to solve the problem with larger thickness of the biofilm than possible before as is shown in the numerical experiments. / Mathematics

The electrical activity of the heart is a well studied process. Mathematical modeling and computer simulations are used to study the cardiac electrical activity: several mathematical models exist, among them the microscopic model, which is based on the explicit representation of individual cells. The cardiac tissue is viewed as two separate domains: the intra-cellular and extra-cellular domains, Ωᵢ and Ωₑ, respectively, separated by cellular membranes Γ. The microscopic model consists of a set of Poisson equations, one for each sub-domain, Ωᵢ and Ωₑ, coupled on interfaces Γ with nonlinear transmission conditions involving a system of ODEs. The unusual transmission conditions on Γ make the model challenging to solve numerically. In this thesis, we first focus on the dimensional analysis of the microscopic model. We then reformulate the problem on the interface Γ using a Steklov-Poincaré operator. We discretize the model in space using finite element methods. We prove the existence of a semi-discrete solution using a reformulation of the model as an ODE system on the interface Γ. We derive stability and error estimates for the finite element method. Afterwards, we consider five numerical schemes including the Godunov splitting method, two implicit methods, (Backward Euler (BE) and second order Backward Differentiation Formula (BDF2)), and two semi-implicit methods (Forward Backward Euler (FBE), and second order Semi-implicit Backward Differentiation Formula (SBDF2)). A convergence analysis of the implicit and semi-implicit methods is performed and the results are compared with manufactured solutions that we have proposed. Numerical results are presented to compare the stability, accuracy and efficiency of the methods. CPU times needed to solve the problem over a single cell using FBE, SBDF2 and Godunov splitting methods are reported. The results show that FBE and Godunov splitting methods achieve better numerical accuracy and efficiency than implicit and SBDF2 schemes, for a given computational time. Finally, we solve the model using FBE and Domain Decomposition Method (DDM) for two cells connected to each other by a gap junction. We investigate the influence of the space discretization and we explore the differences between a conforming and nonconforming mesh on Γ. We compare the solutions obtained with both FBE and DDM methods. The results show that both methods give the same solution. Therefore, the DDM is capable of providing an accurate solution with a minimal number of sub-domain iterations.

The Finite Element Method (FEM) is a powerful numerical method to solve wave propagation problems for open-region electromagnetic radiation/scattering problems involving objects with arbitrary geometry and constitutive parameters. In high-frequency applications, the FEM requires an electrically large computational domain, implying a large number of unknowns, such that the numerical solution of the problem is not feasible even on state-of-the-art computers. An appealing way to solve a large FEM problem is to employ a Domain Decomposition Method (DDM) that allows the decomposition of a large problem into several coupled subproblems which can be solved independently, thus reducing considerably the memory storage requirements. In this thesis, two new domain decomposition algorithms (FB-DDM and ILF-DDM) are implemented for the finite element solution of electromagnetic radiation/scattering problems. For this purpose, a nodal FEM code (FEMS2D) employing triangular elements and a vector FEM code (FEMS3D) employing tetrahedral edge elements have been developed for 2D and 3D problems, respectively. The unbounded domain of the radiation/scattering problem, as well as the boundaries of the subdomains in the DDMs, are truncated by the Perfectly Matched Layer (PML) absorber. The PML is implemented using two new approaches: Locally-conformal PML and Multi-center PML. These approaches are based on a locally-defined complex coordinate transformation which makes possible to handle challenging PML geometries, especially with curvature discontinuities. In order to implement these PML methods, we also introduce the concept of complex space FEM using elements with complex nodal coordinates. The performances of the DDMs and the PML methods are investigated numerically in several applications.

Domain decomposition and convex optimization play fundamental roles in current computation and analysis in many areas of science and engineering. These methods have been well developed and studied in the past thirty years, but they still require further study and improving not only in mathematics but in actual engineering computation with exponential increase of computational complexity and scale. The main goal of this thesis is to develop some efficient and powerful algorithms based on domain decomposition method and convex optimization. The topicsstudied in this thesis mainly include two classes of convex optimization problems: optimal control problems governed by time-dependent partial differential equations and general structured convex optimization problems. These problems have acquired a wide range of applications in engineering and also demand a very high computational complexity. The main contributions are as follows: In Chapter 2, the relevance of an adequate inner loop starting point (as opposed to a sufficient inner loop stopping rule) is discussed in the context of a numerical optimization algorithm consisting of nested primal-dual proximal-gradient iterations. To study the optimal control problem, we obtain second order domain decomposition methods by combining Crank-Nicolson scheme with implicit Galerkin method in the sub-domains and explicit flux approximation along inner boundaries in Chapter 3. Parallelism can be easily achieved for these explicit/implicit methods. Time step constraints are proved to be less severe than that of fully explicit Galerkin finite element method. Based on the domain decomposition method in Chapter 3, we propose an iterative algorithm to solve an optimal control problem associated with the corresponding partial differential equation with pointwise constraint for the control variable in Chapter 4. In Chapter 5, overlapping domain decomposition methods are designed for the wave equation on account of prediction-correction" strategy. A family of unit decomposition functions allow reasonable residual distribution or corrections. No iteration is needed in each time step. This dissertation also covers convergence analysis from the point of view of mathematics for each algorithm we present. The main discretization strategy we adopt is finite element method. Moreover, numerical results are provided respectivelyto verify the theory in each chapter. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

No description available.

No description available.