Mandyoc: A finite element code to simulate
thermochemical convection in parallel

Sacek, Victor; Assunção, Jamison; Pesce, Agustina; Silva, Rafael Monteiro da

doi:10.21105/joss.04070

Cited by 3 publications

(5 citation statements)

References 16 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…You can also find it through Sacek et al. (2022). Input files used in this paper are provided as Supporting Information S1.…”

Section: Data Availability Statementmentioning

confidence: 99%

“…We used the finite element code Mandyoc (Sacek et al., 2022) to numerically solve the differential equations for the conservation of mass, energy, and momentum for an incompressible non‐Newtonian fluid (Zhong et al., 2007):

{u}_{i,i}=0

{\sigma }_{ij,j}+{g}_{i}\rho =0

\frac{\partial T}{\partial t}+{u}_{i}{T}_{,i}=\kappa {T}_{,ii}+\frac{H}{{c}_{p}}+{u}_{i}{g}_{i}\alpha \frac{T}{{c}_{p}}

where

{\sigma }_{ij}=-P{\delta }_{ij}+\eta \left({u}_{i,j}+{u}_{j,i}\right),

\rho ={\rho }_{0}\left[1-\alpha \left(T-{T}_{0}\right)\right],

t is time, T is temperature, u i is the i th component of the velocity field, g is gravity, ρ 0 is the reference rock density at temperature T 0 = 0°C, α is the coefficient of thermal expansion, κ is the thermal diffusivity, H heat production per unit of mass, c p is specific heat capacity, P is the total pressure, η is the effective dynamic viscosity and δ ij is the Kronecker delta. Repeated indexes mean summation and the indexes after the comma represent the partial derivative of the respective coordinate.…”

Section: The Numerical Modelmentioning

confidence: 99%

“…We used the finite element code Mandyoc (Sacek et al, 2022) to numerically solve the differential equations for the conservation of mass, energy, and momentum for an incompressible non-Newtonian fluid (Zhong et al, 2007):…”

Section: The Numerical Modelmentioning

confidence: 99%

“…The numerical code (MANDYOC) used in this paper is available on GitHub repository (https://github.com/ ggciag/mandyoc). You can also find it through Sacek et al (2022). Input files used in this paper are provided as Supporting Information S1.…”

Section: Data Availability Statementmentioning

confidence: 99%

See 3 more Smart Citations

Effects of Tectonic Quiescence Between Orogeny and Rifting

Salazar-Mora

Sacek

2023

Tectonics

Self Cite

View full text Add to dashboard Cite

Tectonic quiescence is a period of nearly null relative plate motion. During the Atlantic Ocean opening, rifting cut across continental segments of supercontinent Pangea that had been through different periods of quiescence, which varied from tens to hundreds of Myr. Using numerical models, we explored the effects of quiescence between orogeny and rifting on the pre‐rift lithosphere and subsequent conjugate rifted margin configuration. Our model results showed that quiescence of ∼30–60 Myr caused wide orogenic wedges (>300 km) to be heated beneath their centers, which led to the breakup of the lithospheric mantle before the breakup of the continental crust during rifting, resulting in hyperextended conjugate rifted margins. In these scenarios, continental lithosphere breakup away from the suture contributed to the preservation of previously subducted lower crust along the subduction zone. Long periods of tectonic quiescence (∼100 Myr) experienced by a wide orogen permitted efficient orogenic collapse and lateral spreading, with pronounced lithospheric thermal weakening. Consequently, post‐quiescence rifting reactivated the suture and effectively drove plate eduction, exhuming most of the subducted continental lower crust and mantle lithosphere. Ultimately, ultra‐wide domains developed into asymmetric conjugate margins. Contrastingly, in narrow orogenic wedges (<200 km), different quiescence durations had minor effects on the margin architecture. In these cases, a longer quiescence time contributed to a larger preservation of subducted crust in the fossil subduction zone after rifting. Our models agree with conjugate rifted margins in the North and South Atlantic, which probably developed after the collapse of wide orogens and subsequent periods of tectonic quiescence.

show abstract

“…You can also find it through Sacek et al. (2022). Input files used in this paper are provided as Supporting Information S1.…”

Section: Data Availability Statementmentioning

confidence: 99%

{u}_{i,i}=0

{\sigma }_{ij,j}+{g}_{i}\rho =0

\frac{\partial T}{\partial t}+{u}_{i}{T}_{,i}=\kappa {T}_{,ii}+\frac{H}{{c}_{p}}+{u}_{i}{g}_{i}\alpha \frac{T}{{c}_{p}}

where

{\sigma }_{ij}=-P{\delta }_{ij}+\eta \left({u}_{i,j}+{u}_{j,i}\right),

\rho ={\rho }_{0}\left[1-\alpha \left(T-{T}_{0}\right)\right],

Section: The Numerical Modelmentioning

confidence: 99%

Section: The Numerical Modelmentioning

confidence: 99%

Section: Data Availability Statementmentioning

confidence: 99%

See 2 more Smart Citations

Effects of Tectonic Quiescence Between Orogeny and Rifting

Salazar-Mora

Sacek

2023

Tectonics

Self Cite

View full text Add to dashboard Cite

show abstract

“…In the following years, this thermomechanical code was improved, using the Portable, Extensible Toolkit for Scientific Computation (PETSc) (Balay et al, 1997(Balay et al, , 2021a, allowing the code to be fully parallelized using the Message Passing Interface (MPI), and the present version of the code can simulate different rheological behaviour, including Newtonian flow, non-linear viscous flow or viscoplastic deformation (Sacek et al, 2022). This code, named Mandyoc, is freely available on Github platform (https://github.com/ggciag/mandyoc/) and incorporates free surface, necessary to reproduce the topographic evolution during the numerical simulation.…”

Section: Development Of Numerical Geodynamics At Uspmentioning

confidence: 99%

Computational Geodynamics of South American Plate: Contributions and Perspectives

Sacek,

Ussami,

Pesce

et al. 2022

Braz J Geophys

View full text Add to dashboard Cite

Computational geodynamics is a key tool in Earth Sciences, enabling the simulation of different geodynamic processes at any time scale in Nature or those which cannot be adequately reproduced by analogue modelling in laboratorial conditions. Specifically, the coupling of surface processes of erosion and sedimentation with the internal dynamics of the lithosphere is a complex problem that involves the solution of a set of differential equations that are only adequately solved by numerical models. During the last three decades, different numerical models were developed to explain the importance of the coupling between the surface and internal dynamics of the Earth, both in active margins and in stable tectonic domains, showing how the coupling leads to counter-intuitive results not observed when each process is analyzed separately. One example of this complex coupling is the feedback mechanism between erosion of the landscape and the regional isostatic response of the lithosphere. In this work, we present simple isostatic and flexural elements that highlight the importance of surface processes on the stress and strain pattern in lithospheric plates. Initially, we present a review on the development of computational geodynamics at University of São Paulo. This review is followed by an analysis of the density structure of the Earth and how the high-density contrast at the Earth’s surface creates a major impact on the isostatic equilibrium of the lithosphere when variations on topographic loads are taken into account. Additionally, we show that the wavelength of the denudation of the Earth’s surface due to fluvial dynamics corresponds to the characteristic length scale for flexural bending of the lithosphere, maximizing the flexural stresses in the lithosphere. Finally, we present recent works on the coupling between surface and lithospheric processes and future challenges for the development of computational geodynamics, with possible strategies to solve them.

show abstract