Masao Shimada, David Tae, Tao Xue, Rohit Deokar and K K Tamma
The purpose of this paper is to present new implementation aspects of unified explicit time integration algorithms, called the explicit GS4-II family of algorithms, of a…
Abstract
Purpose
The purpose of this paper is to present new implementation aspects of unified explicit time integration algorithms, called the explicit GS4-II family of algorithms, of a second-order time accuracy in all the unknowns (e.g. positions, velocities, and accelerations) with particular attention to the moving-particle simulation (MPS) method for solving the incompressible fluids with free surfaces.
Design/methodology/approach
In the present paper, the explicit GS4-II family of algorithms is implemented in the MPS method in the following two different approaches: a direct explicit formulation with the use of the weak incompressibility equation involving the (modified) speed of sound; and a predictor-corrector explicit formulation. The first approach basically follows the concept of the explicit MPS method, presented in the literature, and the latter approach employs a similar concept used in, for example, a fractional-step method in computational fluid dynamics.
Findings
Illustrative numerical examples demonstrate that any scheme within the proposed algorithmic framework captures the physics with the necessary second-order time accuracy and stability.
Originality/value
The new algorithmic framework extended with the GS4-II family encompasses a multitude of pastand new schemes and offers a general purpose and unified implementation.
Details
Keywords
K K Tamma and Siti Ujila Masuri
The purpose of this paper is to describe how a generalized single-system-single-solve (GS4-1) computational framework, previously developed for linear first order transient…
Abstract
Purpose
The purpose of this paper is to describe how a generalized single-system-single-solve (GS4-1) computational framework, previously developed for linear first order transient systems, can be properly extended for use in nonlinear counterparts, with particular applications to time dependent Burgers’ equation, which is well-known to serve as a simplified model of fluid dynamics, for illustrations of the essential concepts.
Design/methodology/approach
The framework permits, for a very general family of time integrators where traditional methods are a subset, much needed desirable features including second order time accuracy, robustness and unconditional stability, zero-order overshoot behavior, and additionally, a selective control of high frequency damping for both the primary variable and its time derivative. The latter, which is a new, key desirable feature not available in past/existing methods to-date, allows for different amounts of high frequency damping for both the primary variable and its time derivative to ensure physically accurate solutions of these variables. This is in contrast to having only limited control of these numerical dampings, often indiscriminately, as in some past developments which can lead to numerical instabilities in the time derivative variable. The extension of the framework to nonlinear problems, as described in this paper, is achieved via the use of a normalized time weighted residual approach, which naturally allows the time discretization of the transient nonlinear systems as being the natural extensions of the linear systems.
Findings
The primary aim is also to demonstrate the advantage of the selective control feature inherit in the present numerical methodologies for these nonlinear first order transient systems as in the linear counterparts.
Originality/value
The authors wish to tackle the challenges to further enable extensions to nonlinear first order transient systems that frequently arise in fluid dynamics problems; this is the focus of this paper. The primary wish is to demonstrate the ability of the GS4-1 framework for nonlinear first order transient systems as seen in the linear transient counterparts; while on one hand the authors show that an equal amount of high frequency damping (i.e. ρ ∞ = ρ ∞ s) leads to non-physical instability in the time derivative variable for a minimal damping required to obtain acceptable solution of the primary variable, on the other hand, the authors particularly demonstrate how this instability can be easily tuned off via the selective control feature (i.e. ρ ∞ ≠ ρ ∞ s) offered by the developed framework; thereby, demonstrating its robustness and superiority.