HR: 0800h
AN: T11C-0736    [Abstracts]
TI: Development of three dimensional Eulerian numerical procedure toward plate-mantle simulation: accuracy test by the fluid rope coiling
AU: * Furuichi, M
EM: m-furuic@jamstec.go.jp
AF: Earth Simulator Center, Japan Agency for Marine-Earth Science and Technology, 3173- 25 Showa-machi, Kanazawa-ku, Yokohama, 236-0001, Japan
AU: Kameyama, M
EM: kameyama@sci.ehime-u.ac.jp
AF: Geodynamics Research Center (GRC), Ehime University, 2-5 Bunkyo-cho, Matsuyama, 790-8577, Japan
AU: Kageyama, A
EM: kage@jamstec.go.jp
AF: Earth Simulator Center, Japan Agency for Marine-Earth Science and Technology, 3173- 25 Showa-machi, Kanazawa-ku, Yokohama, 236-0001, Japan
AB: Reproducing a realistic plate tectonics with mantle convection simulation is one of the greatest challenges in computational geophysics. We have developed a three dimensional Eulerian numerical procedure toward plate-mantle simulation, which includes a finite deformation of the plate in the mantle convection. Our method, combined with CIP-CSLR (Constrained Interpolation Profile method-Conservative Semi-Lagrangian advection scheme with Rational function) and ACuTE method, enables us to solve advection and force balance equations even with a large and sharp viscosity jump, which marks the interface between the plates and surrounding upper mantle materials. One of the typical phenomena represented by our method is a fluid rope coiling event, where a stream of viscous fluid is poured onto the bottom plane from a certain height. This coiling motion is due to delicate balances between bending, twisting and stretching motions of fluid rope. In the framework of the Eulerian scheme, the fluid rope and surrounding air are treated as a viscosity profile which differs by several orders of magnitude. Our method solves the complex force balances of the fluid rope and air, by a multigrid iteration technique of ACuTE algorithm. In addition, the CIP-CSLR advection scheme allows us to obtain a deforming shape of the fluid rope, as a low diffusive solution in the Eulerian frame of reference. In this presentation, we will show the simulation result of the fluid rope coiling as an accuracy test for our simulation scheme, by comparing with the simplified numerical solution for thin viscous jet.
UR: http://www.es.jamstec.go.jp/esc/research/Solid/index.en.html
DE: 1213 Earth's interior: dynamics (1507, 7207, 7208, 8115, 8120)
DE: 4445 Nonlinear differential equations
DE: 8120 Dynamics of lithosphere and mantle: general (1213)
DE: 8121 Dynamics: convection currents, and mantle plumes
DE: 8155 Plate motions: general (3040)
SC: Tectonophysics [T]
MN: 2007 Fall Meeting