Implicit monolithic parallel solution algorithm for seismic analysis of dam-reservoir systems

Özmen, Semih
This research mainly focuses on developing a computationally scalable and efficient solution algorithm that can handle linear dynamic analysis of dam-reservoir interaction problem. Lagrangian fluid finite elements are utilized and compressibility and viscosity of the fluid are taken into consideration during the reservoir modeling. In order to provide computational scalability and efficiency, domain decomposition methods implemented with parallel computing approaches such as Finite Element Tearing and Interconnecting (FETI) family solution algorithms are utilized for the coupling of the subdomains and a fully implicit monolithic solution algorithm is developed. Following that, the ways of performance improvements for the algorithm are demonstrated. Re-orthogonalization is utilized to increase the convergence rate of the solution of system equations and Krylov subspaces are utilized in order to decrease the required iterations for the future time integration steps. Additionally, utilization of deflation methods on Preconditioned Conjugate Gradient (PCG) and Finite Element Tearing and Interconnecting (FETI) family solvers is discussed. Due to the fact that efficiency and behavior of the deflation methods depends on the deflation vectors utilized, different deflation vector generation methods are also investigated. Two of the deflation vector generation methods are from literature, i.e. “Subdomain Deflation Method” and “Recursive Deflation Method for Heterogeneous Problems”. In addition to them, a novel semi-heuristic deflation vector generation strategy which relies on the pre-selected zero energy modes of finite element formulations is proposed. Requirements, improvements and efficiencies of these methods are demonstrated for the serial solution of water tank with flexible walls problem. In order to investigate the efficiency and scalability of the presented solution approach on the solutions of more realistic problems by computer clusters, this approach is implemented by utilizing C++ programming language and PETSc library. In this parallel implementation, FETI-DP solution algorithm is utilized with different deflation vector generation algorithms. The efficiency and the scalability of the parallel solution framework are discussed for different types of finite elements, for different partitioning approaches and for different number of processors. Finally, the solution performance is presented for a large actual dam model, Pine Flat Dam in California, USA.


Optimal load and resistance factor design of geometrically nonlinear steel space frames via tabu search and genetic algorithm
In this paper, algorithms are presented for the optimum design of geometrically nonlinear steel space frames using tabu search and genetic algorithm. Tabu search utilizes the features of short-term memory facility (tabu list) and aspiration criteria. Genetic algorithm employs reproduction, crossover and mutation operators. The design algorithms obtain minimum weight frames by selecting suitable sections from a standard set of steel sections such as American Institute of Steel Construction (AISC) wide-flange...
Hierarchical parallelisation strategy for multilevel fast multipole algorithm in computational electromagnetics
Ergül, Özgür Salih (Institution of Engineering and Technology (IET), 2008-01-03)
A hierarchical parallelisation of the multilevel fast multipole algorithm (MLFMA) for the efficient solution of large-scale problems in computational electromagnetics is presented. The tree structure of MLFMA is distributed among the processors by partitioning both the clusters and the samples of the fields appropriately for each level. The parallelisation efficiency is significantly improved compared to previous approaches, where only the clusters or only the fields are partitioned in a level.
Ergül, Özgür Salih (2013-11-09)
A parallel implementation of the multilevel fast multipole algorithm (MLFMA) is developed for fast and accurate solutions of electromagnetics problems involving complex plasmonic metamaterial structures. Composite objects that consist of multiple penetrable regions, such as dielectric, lossy, and plasmonic parts, are formulated rigorously with surface integral equations and solved iteratively via MLFMA. Using the hierarchical strategy for the parallelization, the developed implementation is capable of simul...
Stability analysis of constraints in flexible multibody systems dynamics
İder, Sıtkı Kemal (Elsevier BV, 1990-1)
Automated algorithms for the dynamic analysis and simulation of constrained multibody systems assume that the constraint equations are linearly independent. During the motion, when the system is at a singular configuration, the constraint Jacobian matrix possesses less than full rank and hence it results in singularities. This occurs when the direction of a constraint coincides with the direction of the lost degree of freedom. In this paper the constraint equations for deformable bodies are modified for use...
Optimising a nonlinear utility function in multi-objective integer programming
Ozlen, Melih; Azizoğlu, Meral; Burton, Benjamin A. (2013-05-01)
In this paper we develop an algorithm to optimise a nonlinear utility function of multiple objectives over the integer efficient set. Our approach is based on identifying and updating bounds on the individual objectives as well as the optimal utility value. This is done using already known solutions, linear programming relaxations, utility function inversion, and integer programming. We develop a general optimisation algorithm for use with k objectives, and we illustrate our approach using a tri-objective i...
Citation Formats
S. Özmen, “Implicit monolithic parallel solution algorithm for seismic analysis of dam-reservoir systems,” Ph.D. - Doctoral Program, Middle East Technical University, 2016.