skip to main content


Title: A Multirate Approach for Fluid-Structure Interac- tion Computation with Decoupled Methods
We investigate a multirate time step approach applied to decoupled meth- ods in fluid and structure interaction (FSI) computation, where two different time steps are employed for fluid and structure respectively. For illustration, the multirate technique is examined by applying the decoupled β scheme. Numerical experiments show that the proposed approach is stable and retains the same order of accuracy as the original single time step scheme, while with much less computational expense.  more » « less
Award ID(s):
1831950 1700328
NSF-PAR ID:
10097721
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Communications in computational physics
ISSN:
1991-7120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The anisotropic phase‐field dendritic crystal growth model is a highly nonlinear system that couples the anisotropic Allen–Cahn equation and the thermal equation together. Due to the high anisotropy and nonlinear couplings in the system, how to develop an accurate and efficient, especially a fully decoupled scheme, has always been a challenging problem. To solve the challenge, in this article, we construct a novel fully decoupled numerical scheme which is also linear, energy stable, and second‐order time accurate. The key idea to realize the full decoupling structure is to introduce an ordinary differential equation to deal with the nonlinear coupling terms satisfying the so‐called “zero‐energy‐contribution” property. This scheme is very effective and easy to implement since only a few fully decoupled elliptic equations with constant coefficients need to be solved at each time step. We rigorously prove the solvability of each step and the unconditional energy stability, and perform a large number of numerical simulations in 2D and 3D to demonstrate its stability and accuracy numerically.

     
    more » « less
  2. Summary

    With the stagnation of processor core performance, further reductions in the time to solution for geophysical fluid problems are becoming increasingly difficult with standard time integrators. Parallel‐in‐time exposes and exploits additional parallelism in the time dimension, which is inherently sequential in traditional methods. The rational approximation of exponential integrators (REXI) method allows taking arbitrarily long time steps based on a sum over a number of decoupled complex PDEs that can be solved independently massively parallel. Hence, REXI is assumed to be well suited for modern massively parallel super computers, which are currently trending. To date, the study and development of the REXI approach have been limited to linearized problems on the periodic two‐dimensional plane. This work extends the REXI time stepping method to the linear shallow‐water equations on the rotating sphere, thus moving the method one step closer to solving fully nonlinear fluid problems of geophysical interest on the sphere. The rotating sphere poses particular challenges for finding an efficient solver due to the zonal dependence of the Coriolis term. Here, we present an efficient REXI solver based on spherical harmonics, showing the results of a geostrophic balance test, a comparison with alternative time stepping methods, an analysis of dispersion relations indicating superior properties of REXI, and finally, a performance comparison on the Cheyenne supercomputer. Our results indicate that REXI not only can take larger time steps but also can be used to gain higher accuracy and significantly reduced time to solution compared with currently existing time stepping methods.

     
    more » « less
  3. We construct new first- and second-order pressure correctionschemes using the scalar auxiliary variable approach for the Navier-Stokes equations. These schemes are linear, decoupled and only require solving a sequence of Poisson type equations at each time step. Furthermore, they are unconditionally energy stable. We also establish rigorous error estimates in the two dimensional case for the velocity and pressure approximation of the first-order scheme without any condition on the time step. 
    more » « less
  4. In this paper, we consider the numerical approximation for a phase field model of the coupled two-phase free flow and two-phase porous media flow. This model consists of Cahn– Hilliard–Navier–Stokes equations in the free flow region and Cahn–Hilliard–Darcy equations in the porous media region that are coupled by seven interface conditions. The coupled system is decoupled based on the interface conditions and the solution values on the interface from the previous time step. A fully discretized scheme with finite elements for the spatial discretization is developed to solve the decoupled system. In order to deal with the difficulties arising from the interface conditions, the decoupled scheme needs to be constructed appropriately for the interface terms, and a modified discrete energy is introduced with an interface component. Furthermore, the scheme is linearized and energy stable. Hence, at each time step one need only solve a linear elliptic system for each of the two decoupled equations. Stability of the model and the proposed method is rigorously proved. Numerical experiments are presented to illustrate the features of the proposed numerical method and verify the theoretical conclusions. 
    more » « less
  5. Abstract

    We present a formulation and numerical algorithm to extend the scheme for gray radiation magnetohydrodynamics (MHD) developed by Jiang to include the frequency dependence via the multigroup approach. The entire frequency space can be divided into an arbitrary number of groups in the lab frame, and we follow the time-dependent evolution of frequency-integrated specific intensities along discrete rays inside each group. Spatial transport of photons is done in the lab frame while all the coupling terms are solved in the fluid rest frame. Lorentz transformation is used to connect different frames. The radiation transport equation is solved fully implicitly in time while the MHD equations are evolved explicitly so that time step is not limited by the speed of light. A finite volume approach is used for transport in both spatial and frequency spaces to conserve the radiation energy density and momentum. The algorithm includes photon absorption, electron scattering, as well as Compton scattering, which is calculated by solving the Kompaneets equation. The algorithm is accurate for a wide range of optical depth conditions and can handle both radiation-pressure- and gas-pressure-dominated flows. It works for both Cartesian and curvilinear coordinate systems with adaptive mesh refinement. We provide a variety of test problems including a radiating sphere, shadow test, absorption of a moving gas, Bondi-type flows, as well as a collection of test problems for thermal and bulk Compton scattering. We also discuss examples where frequency dependence can make a big difference compared with the gray approach.

     
    more » « less