An adaptive discretization for incompressible and compressible flow using a multitude of moving Cartesian grids with gap flow treatment
Abstract/Contents
- Abstract
- This dissertation presents novel methods for simulating incompressible and compressible flow on a multitude of Cartesian grids that can rotate and translate in order to decompose the domain into different regions with varying spatial resolutions. While there are a wide variety of methods for adaptive discretization, many of these methods suffer from issues with costly remeshing and domain decomposition when they are scaled to solve large problems that require the use of large distributed systems. Block structured approaches such as Adaptive Mesh Refinement (AMR) and Chimera grid methods have been successful in alleviating these issues by utilizing structured grids patched upon one another. However, typical AMR methods constrain grid patches to be axis-aligned greatly increasing the number of patches required. With so many small patches, typical AMR methods are often more akin to unstructured grids with respect to parallelization and scalability. Chimera grid methods allow the grid patches to rotate allowing one to resolve interesting features with far fewer degrees of freedom. Moreover, unlike typical AMR methods which require the coarse grid lines to match up with the fine grid lines along patch boundaries, Chimera grid methods do not have this requirement allowing the grids to move in order to for example follow the motion of moving solids with no need of remeshing. The presented computational framework can be categorized as a Chimera grid method, and new ideas are proposed regarding conservation, linear systems for implicit solve, and alleviating time step restrictions. The incompressible Navier-Stokes equations are discretized on overlapping grids by first performing advection on each grid with first or second order accurate semi-Lagrangian schemes extended to Chimera grids in order to alleviate any time step restrictions associated with small cells which are introduced due to adaptivity. In order to solve for the stiff terms such as the pressure or viscous forces implicitly on overlapping grids, local Voronoi meshes are constructed along intergrid boundaries to connect the various degrees of freedom across different grids in a contiguous manner, resulting in a symmetric positive-definite system that can be solved via the preconditioned conjugate gradient method. In order to handle free surface flow on overlapping grids, the particle level set method is modified, including devising particle treatment across grid boundaries with disparate cell sizes, and designing strategies to deal with locality in the implementation of the level set and fast marching algorithms. The resulting method is highly scalable on distributed parallel architectures with minimal communication costs. The Euler equations for compressible flow are discretized using a semi-implicit formulation that splits the time integration into an explicit step for advection followed by an implicit solve for the pressure. A second order accurate flux based scheme is devised to handle advection on each moving Cartesian grid using an effective characteristic velocity that accounts for the grid motion. In order to avoid the stringent time step restriction imposed by small cells, strategies are proposed in order to allow for a fluid velocity CFL number larger than 1. The stringent time step restriction related to the sound speed is alleviated by formulating an implicit linear system in order to find a pressure consistent with the equation of state, again utilizing the Voronoi mesh obtaining a symmetric positive-definite system. Since a straightforward application of this technique contains an inherent central differencing which can result in spurious oscillations, a new high order diffusion term is introduced similar in spirit to ENO-LLF but solved for implicitly in order to avoid any associated time step restrictions. The method is conservative on each grid, as well as globally conservative on the background grid that contains all other grids. Moreover, a conservative interpolation operator is devised for conservatively remapping values in order to keep them consistent across different overlapping grids. Additionally, the method is extended to handle two-way solid fluid coupling in a monolithic fashion. In solid fluid coupling problems, the fluid in the thin gap between solids in close proximity is difficult to resolve with fluid grids. Although one might attempt to address this difficulty using an adaptive, body-fitted, or ALE fluid grid, the size of the fluid cells can shrink to zero as the solids collide. The inability to apply pressure forces in a thin lubricated gap tends to make the solids stick together, since collision forces stop interpenetration but vanish when the solids are separating leaving the fluid pressure forces on the surfaces of the solids unbalanced in regard to the gap region. This problem is addressed by adding fluid pressure and velocity degrees of freedom onto solids' surfaces, and subsequently using the resulting pressure forces to provide solid fluid coupling in the thin gap region. These fluid pressure and velocity degrees of freedom readily resolve the tangential flow along the solid surface inside the gap and are two-way coupled to the fluid degrees of freedom on the grids allowing the fluid to freely flow into and out of the gap region, which again results in a symmetric positive-definite implicit linear system.
Description
| Type of resource | text |
|---|---|
| Form | electronic; electronic resource; remote |
| Extent | 1 online resource. |
| Publication date | 2015 |
| Issuance | monographic |
| Language | English |
Creators/Contributors
| Associated with | Qiu, Linhai |
|---|---|
| Associated with | Stanford University, Department of Electrical Engineering. |
| Primary advisor | Fedkiw, Ronald P, 1968- |
| Thesis advisor | Fedkiw, Ronald P, 1968- |
| Thesis advisor | Levis, Philip |
| Thesis advisor | Savarese, Silvio |
| Advisor | Levis, Philip |
| Advisor | Savarese, Silvio |
Subjects
| Genre | Theses |
|---|
Bibliographic information
| Statement of responsibility | Linhai Qiu. |
|---|---|
| Note | Submitted to the Department of Electrical Engineering. |
| Thesis | Thesis (Ph.D.)--Stanford University, 2015. |
| Location | electronic resource |
Access conditions
- Copyright
- © 2015 by Linhai Qiu
- License
- This work is licensed under a Creative Commons Attribution Non Commercial 3.0 Unported license (CC BY-NC).
Also listed in
Loading usage metrics...