search

geodynamics / citcoms

CitcomS is a finite element code designed to solve compressible thermochemical convection problems relevant to Earth's mantle.
GNU General Public License v2.0
40 stars 26 forks source link

How to set parameters if I want to integrate backward in time? #6

Closed ymwanghub closed 5 years ago

ymwanghub commented 7 years ago

I want to run the CitcomS code backward in time from a given initial temperature, which parameters do I need to set? I found two parameters, i.e., the Rayleigh number Ra, and the inputdiffusivity. Do I need to set Ra to be negative and also inputdiffusivity=-1 (which is 1 in User Manual)? Is there any other parameter I need to modify?

Thanks.

tan2 commented 7 years ago

Set Ra to be negative. And set inputdiffusivity to something small, maybe around 0.1. More information can be found in this paper:

Conrad and Gurnis (2003), G^3, Seismic tomography, surface uplift, and the breakup of Gondwanaland: Integrating mantle convection backwards in time

On Mon, Apr 17, 2017 at 3:43 PM, Yongming Wang notifications@github.com wrote:

I want to run the CitcomS code backward in time from a given initial temperature, which parameters do I need to set? I found two parameters, i.e., the Rayleigh number Ra, and the inputdiffusivity. Do I need to set Ra to be negative and also inputdiffusivity=-1 (which is 1 in User Manual)? Is there any other parameter I need to modify?

Thanks.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/geodynamics/citcoms/issues/6, or mute the thread https://github.com/notifications/unsubscribe-auth/AANoQlMnK4u0zdXPFu1D8yNXHNlZEsOAks5rwxgHgaJpZM4M-80f .

ymwanghub commented 7 years ago

Thanks.