Year-end Sale % OFF 
Abstract
Finite element methods are often used to model Earth processes involving slow viscous or viscoelastic flow. Inertial terms of the Navier-Stokes equations are neglected in very slow flows, so timestep size is not limited by the Courant instability. However, where there is advection of density contrasts in a gravitational field, over-advection can lead to numerically induced flow oscillations. We derive analytic results for the maximum stable timestep size in two cases: a free surface over a fluid of uniform density, and a free surface kept level by sedimentation/erosion, but with a density gradient in the underlying medium. Using parameters appropriate to the Earth's crust we show that the density-contrast instability occurs for timesteps larger than 3000 years for the constant-density case. For a fluid with a density gradient of 10 kg/m3 per km the solution is stable for timesteps up to about 200,000 years if full erosion/sedimentation is implemented.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callMore From: International Journal of Computational Fluid Dynamics
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
