You are here
SpaceTime Spectral Element Methods in Fluid Dynamics and Materials Science
Title:  SpaceTime Spectral Element Methods in Fluid Dynamics and Materials Science. 
26 views
17 downloads 

Name(s): 
Pei, Chaoxu, author Sussman, Mark, professor codirecting dissertation Hussaini, M. Yousuff, professor codirecting dissertation Dewar, William K., university representative Cogan, Nicholas G., committee member Wang, Xiaoming, Ph.D., committee member Florida State University, degree granting institution College of Arts and Sciences, degree granting college Department of Mathematics, degree granting department 

Type of Resource:  text  
Genre: 
Text Doctoral Thesis 

Issuance:  monographic  
Date Issued:  2017  
Publisher:  Florida State University  
Place of Publication:  Tallahassee, Florida  
Physical Form: 
computer online resource 

Extent:  1 online resource (125 pages)  
Language(s):  English  
Abstract/Description:  In this manuscript, we propose spacetime spectral element methods to solve problems arising from fluid dynamics and materials science. Many engineering applications require one to solve complex problems, such as flows containing multiscale structure in either space or time or both. It is straightforward that highorder methods are always more accurate and efficient than loworder ones for solving smooth problems. For example, spectral element methods can achieve a given level of accuracy with significantly fewer degrees of freedom compared to methods with algebraic convergence rates, e.g., finite difference methods. However, when it comes to complex problems, a high order method should be augmented with, e.g., a level set method or an artificial viscosity method, in order to address the issues caused by either sharp interfaces or shocks in the solution. Complex problems considered in this work are problems with solutions exhibiting multiple scales, i.e., the Stefan problem, nonlinear hyperbolic problems, and problems with smooth solutions but forces exhibiting disparate temporal scales, such as advection, diffusion and reaction processes. Correspondingly, two families of spacetime spectral element methods are introduced in order to achieve spectral accuracy in both space and time. The first category of spacetime methods are the fully implicit spacetime discontinuous Galerkin spectral element methods. In the fully implicit spacetime methods, time is treated as an additional dimension, and the model equation is rewritten into a spacetime formulation. The other category of spacetime methods are specialized for problems exhibiting multiple time scales: multiimplicit spacetime spectral element methods are developed. The method of lines approach is employed in the multiimplicit spacetime methods. The model is first discretized by a discontinuous spectral element method in space, and the resulting ordinary differential equations are then solved by a new multiimplicit spectral deferred correction method. A novel fully implicit spacetime discontinuous Galerkin (DG) spectral element method is presented to solve the Stefan problem in an Eulerian coordinate system. This method employs a level set procedure to describe the timeevolving interface. To deal with the prior unknown interface, a backward transformation and a forward transformation are introduced in the spacetime mesh. By combining an Eulerian description with a Lagrangian description, the issue of dealing with the implicitly defined arbitrary shaped spacetime elements is avoided. The backward transformation maps the unknown timevarying interface in the fixed frame of reference to a known stationary interface in the moving frame of reference. In the moving frame of reference, the transformed governing equations, written in the spacetime framework, are discretized by a DG spectral element method in each spacetime slab. The forward transformation is used to update the level set function and then to project the solution in each phase onto the new corresponding timedependent domain. Two options for calculating the interface velocity are presented, and both options exhibit spectral accuracy. Benchmark tests in one spatial dimension indicate that the method converges with spectral accuracy in both space and time for the temperature distribution and the interface velocity. The interrelation between the interface position and the temperature makes the Stefan problem a nonlinear problem; a Picard iteration algorithm is introduced in order to solve the nonlinear algebraic system of equations and it is found that just a few iterations lead to convergence. We also apply the fully implicit spacetime DG spectral element method to solve nonlinear hyperbolic problems. The spacetime method is combined with two different approaches for treating problems with discontinuous solutions: (i) spacetime dependent artificial viscosity is introduced in order to capture discontinuities/shocks, and (ii) the sharp discontinuity is tracked with spacetime spectral accuracy, as it moves through the grid. To capture the discontinuity whose location is initially unknown, an artificial viscosity term is strategically introduced, and the amount of artificial viscosity varies in time within a given spacetime slab. It is found that spectral accuracy is recovered everywhere except in the "troublesome element(s)'' where the unresolved steep/sharp gradient exists. When the location of a discontinuity is initially known, a spacetime spectrally accurate tracking method has been developed so that the spectral accuracy of the position of the discontinuity and the solution on either side of the discontinuity is preserved. A Picard iteration method is employed to handle nonlinear terms. Within each Picard iteration, a linear system of equations is solved, which is derived from the spacetime DG spectral element discretization. Spectral accuracy in both space and time is first demonstrated for the Burgers' equation with a smooth solution. For tests with discontinuities, the present spacetime method enables better accuracy at capturing the shock strength in the element containing shock when higher order polynomials in both space and time are used. Moreover, the spectral accuracy of the shock speed and location is demonstrated for the solution of the inviscid Burgers' equation obtained by the shock tracking method, and the sensitivity of the number of Picard iterations to the temporal order is discussed. The dynamics of many physical and biological systems involve two or more processes with a wide difference of characteristic time scales, e.g., problems with advection, diffusion and reaction processes. The computational cost of solving a coupled nonlinear system of equations is expensive for a fully implicit (i.e., "monolithic") spacetime method. Thus, we develop another type of a spacetime spectral element method, which is referred to as the multiimplicit spacetime spectral element method. Rather than coupling space and time together, the method of lines is used to separate the discretization of space and time. The model is first discretized by a discontinuous spectral element method in space and the resulting ordinary differential equations are then solved by a new multiimplicit spectral deferred correction method. The present multiimplicit spectral deferred correction method treats processes with disparate temporal scales independently, but couples them iteratively by a series of deferred correction steps. Compared to lower order operator splitting methods, the splitting error in the multiimplicit spectral deferred correction method is eliminated by exploiting an iterative coupling strategy in the deferred correction procedure. For the spectral element discretization in space, two advective flux reconstructions are proposed: extended elementwise flux reconstruction and nonextended elementwise flux reconstruction. A loworder Istable building block time integration scheme is introduced as an explicit treatment for the hyperbolic terms in order to obtain a stable and efficient building block for the spectrally accurate spacetime scheme along with these two advective flux reconstructions. In other words, we compare the extended elementwise reconstruction with Istable building block scheme with the nonextended elementwise reconstruction with Istable building block scheme. Both options exhibit spectral accuracy in space and time. However, the solutions obtained by extended elementwise flux reconstruction are more accurate than those yielded by nonextended elementwise flux reconstruction with the same number of degrees of freedom. The spectral convergence in both space and time is demonstrated for advectiondiffusionreaction problems. Two different coupling strategies in the multiimplicit spectral deferred correction method are also investigated and both options exhibit spectral accuracy in space and time.  
Identifier:  FSU_SUMMER2017_Pei_fsu_0071E_13972 (IID)  
Submitted Note:  A Dissertation submitted to the Department of Mathematics in partial fulfillment of the requirements for the degree of Doctor of Philosophy.  
Degree Awarded:  Summer Semester 2017.  
Date of Defense:  July 5, 2017.  
Bibliography Note:  Includes bibliographical references.  
Advisory Committee:  Mark Sussman, Professor CoDirecting Dissertation; M. Yousuff Hussaini, Professor CoDirecting Dissertation; William Dewar, University Representative; Nick Cogan, Committee Member; Xiaoming Wang, Committee Member.  
Subject(s):  Mathematics  
Persistent Link to This Record:  http://purl.flvc.org/fsu/fd/FSU_SUMMER2017_Pei_fsu_0071E_13972  
Owner Institution:  FSU 
Pei, C. (2017). SpaceTime Spectral Element Methods in Fluid Dynamics and Materials Science. Retrieved from http://purl.flvc.org/fsu/fd/FSU_SUMMER2017_Pei_fsu_0071E_13972