Computational Methods for Multi-Scale Temporal Problems: Algorithms, Analysis, and Numerical Experiments

dc.contributor.advisorNakshatrala, Kalyana Babu
dc.contributor.committeeMemberWang, Keh-Han
dc.contributor.committeeMemberWillam, Kaspar J.
dc.contributor.committeeMemberLim, Gino J.
dc.contributor.committeeMemberKulkarni, Yashashree
dc.creatorKarimi, Saeid
dc.creator.orcid0000-0002-5545-6406
dc.date.accessioned2019-11-17T21:42:20Z
dc.date.available2019-11-17T21:42:20Z
dc.date.createdDecember 2016
dc.date.issued2016-12
dc.date.submittedDecember 2016
dc.date.updated2019-11-17T21:42:21Z
dc.description.abstractA major challenge in numerical simulation of most natural phenomena is the presence of disparate temporal and spatial scales. Capturing all the fine features can be computationally prohibitive. Hence, development of efficient and accurate multi-scale numerical algorithms has gained immense attention from engineers and scientists. Typically, a single numerical method cannot efficiently capture all the aforementioned features. Due to the assumptions made in construction of numerical methods and mathematical models, the range of applicability to various length and time-scales is often limited. A direction in resolving this issue is to apply different numerical methods in different regions of the computational domain. This strategy enables computation of necessary details as desired by the user. In this work, we propose numerical methodologies based on domain partitioning techniques that allow different time-steps and time-integrators in different regions of the computational domain. The first problem of interest is elastodynamics, which can pose various temporal scales in impact, contact and wave propagation problems. A monolithic (strong) coupling algorithm based on non-overlapping domain partitioning is proposed. The proposed algorithm is based on the theory of differential/algebraic equations and its numerical stability, energy conservation and accuracy is studied in detail. Following these findings, we extend this algorithm to advection-diffusion-reaction problems. The proposed algorithm proves useful especially in cases where the relative strength of the involved processes changes dramatically with respect to spatial coordinates. Numerical stability and accuracy of this method is studied and its application to fast bimolecular chemical reactions is showcased. Further on, we confine our attention to single and multiple-relaxation-time lattice Boltzmann methods for the advection-diffusion equation and study their performance in preserving the maximum principle and the non-negative constraint. Finally, a computational framework based on overlapping domain decomposition methods is proposed. This framework is designed for advection-diffusion problems and allows coupling of the finite element method and lattice Boltzmann methods with different time-steps and grid sizes. Additionally, a new method for enforcing the Dirichlet and Neumann boundary conditions on the numerical solution from the lattice Boltzmann method is proposed. This method is based on maximization of entropy and ensures non-negativity of the discrete distributions on the boundary of the domain. We study the performance of this framework through numerical experiments and showcase its application to fast and equilibrium chemical reactions.
dc.description.departmentCivil and Environmental Engineering, Department of
dc.format.digitalOriginborn digital
dc.format.mimetypeapplication/pdf
dc.identifier.citationPortions of this document appear in: Karimi, S., and K. B. Nakshatrala. "On multi-time-step monolithic coupling algorithms for elastodynamics." Journal of Computational Physics 273 (2014): 671-705. And in: Karimi, S., and K. B. Nakshatrala. "A monolithic multi-time-step computational framework for first-order transient systems with disparate scales." Computer Methods in Applied Mechanics and Engineering 283 (2015): 419-453. And in: Karimi, S., and K. B. Nakshatrala. "Do Current Lattice Boltzmann Methods for Diffusion and Advection-Diffusion Equations Respect Maximum Principle and the Non-Negative Constraint?." Communications in Computational Physics 20, no. 2 (2016): 374-404. And in: Karimi, S., and K. B. Nakshatrala. "A hybrid multi-time-step framework for pore-scale and continuum-scale modeling of solute transport in porous media." Computer Methods in Applied Mechanics and Engineering 323 (2017): 98-131.
dc.identifier.urihttps://hdl.handle.net/10657/5442
dc.language.isoeng
dc.rightsThe author of this work is the copyright owner. UH Libraries and the Texas Digital Library have their permission to store and provide access to this work. UH Libraries has secured permission to reproduce any and all previously published materials contained in the work. Further transmission, reproduction, or presentation of this work is prohibited except with permission of the author(s).
dc.subjectComputational mechanics
dc.subjectFinite element analysis
dc.subjectNumerical Analysis
dc.subjectLattice Boltzmann Methods
dc.subjectMulti-Scale Simulation
dc.subjectDomain decomposition
dc.subjectTime Integration
dc.subjectDifferential/Algebraic Equations
dc.titleComputational Methods for Multi-Scale Temporal Problems: Algorithms, Analysis, and Numerical Experiments
dc.type.dcmiText
dc.type.genreThesis
thesis.degree.collegeCullen College of Engineering
thesis.degree.departmentCivil and Environmental Engineering
thesis.degree.disciplineCivil Engineering
thesis.degree.grantorUniversity of Houston
thesis.degree.levelDoctoral
thesis.degree.nameDoctor of Philosophy

Files

Original bundle

Now showing 1 - 1 of 1
Loading...
Thumbnail Image
Name:
KARIMI-DISSERTATION-2016.pdf
Size:
37.19 MB
Format:
Adobe Portable Document Format

License bundle

Now showing 1 - 1 of 1
No Thumbnail Available
Name:
LICENSE.txt
Size:
1.81 KB
Format:
Plain Text
Description: