Three-dimensional Eulerian method for large deformation of viscoelastic fluid: Toward plate-mantle simulation

  • Authors:

  • Mikito Furuichi;Masanori Kameyama;Akira Kageyama

  • Affiliations:

  • Earth Simulator Center, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 3173-25 Showa-machi, Kanazawa-ku, Yokohama 236-0001, Japan;Earth Simulator Center, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 3173-25 Showa-machi, Kanazawa-ku, Yokohama 236-0001, Japan;Earth Simulator Center, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 3173-25 Showa-machi, Kanazawa-ku, Yokohama 236-0001, Japan

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2008

Quantified Score

Hi-index 31.45

Visualization

Abstract

Toward the unified simulation of the large deformation of a rigid viscoelastic material (plate) and the convection of a viscous fluid (mantle), an Eulerian scheme with a semi-Lagrangian method is developed. The scheme adopts the CIP-CSLR method for advection terms of staggered grid system in three dimensions. The positive transported profile of a positive quantity is assured by flux corrections in the dimensional splitting method. The Jaumann co-rotational effect of the stress tensor is also integrated into the semi-Lagrangian treatment. This co-rotated semi-Lagrangian method is combined with an exponential time differencing method in the time development of the Maxwell constitutive model. The large time step comparable to, or larger than, the Maxwell relaxation time is successfully realized. Validation tests are performed for the three-dimensional Rayleigh-Taylor instability of a viscoelastic material with jump discontinuity of the mass density and other material properties.