Source author record

Patrick Sanan

Patrick Sanan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

5works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

5 published item(s)

preprint2022arXiv

Retention of water in terrestrial magma oceans and carbon-rich early atmospheres

Massive steam and CO$_2$ atmospheres have been proposed for magma ocean outgassing of Earth and terrestrial planets. Yet formation of such atmospheres depends on volatile exchange with the molten interior, governed by volatile solubilities and redox reactions. We determine the evolution of magma ocean--atmosphere systems for a range of oxygen fugacities, C/H ratios and hydrogen budgets that include redox reactions for hydrogen (H$_2$--H$_2$O), carbon (CO--CO$_2$), methane (CH$_4$), and solubility laws for H$_2$O and CO$_2$. We find that small initial budgets of hydrogen, high C/H ratios, and oxidizing conditions, suppress outgassing of hydrogen until the late stage of magma ocean crystallization. Hence early atmospheres in equilibrium with magma oceans are dominantly carbon-rich, and specifically CO-rich except at the most oxidizing conditions. The high solubility of H$_2$O limits its outgassing to melt fractions below $\sim$30\%, the fraction at which the mantle transitions from vigorous to sluggish convection with melt percolation. Sluggish melt percolation could enable a surface lid to form, trapping water in the interior and thereby maintaining a carbon-rich atmosphere (equilibrium crystallization). Alternatively, efficient crystal settling could maintain a molten surface, promoting a transition to a water-rich atmosphere (fractional crystallization). However, additional processes, including melt trapping and H dissolution in crystallizing minerals, further conspire to limit the extent of H outgassing, even for fractional crystallization. Hence, much of the water delivered to planets during their accretion can be safely harbored in their interiors during the magma ocean stage, particularly at oxidizing conditions.

preprint2022arXiv

The PETSc Community Is the Infrastructure

The communities who develop and support open source scientific software packages are crucial to the utility and success of such packages. Moreover, these communities form an important part of the human infrastructure that enables scientific progress. This paper discusses aspects of the PETSc (Portable Extensible Toolkit for Scientific Computation) community, its organization, and technical approaches that enable community members to help each other efficiently.

preprint2016arXiv

A Stochastic Performance Model for Pipelined Krylov Methods

Pipelined Krylov methods seek to ameliorate the latency due to inner products necessary for projection by overlapping it with the computation associated with sparse matrix-vector multiplication. We clarify a folk theorem that this can only result in a speedup of $2\times$ over the naive implementation. Examining many repeated runs, we show that stochastic noise also contributes to the latency, and we model this using an analytical probability distribution. Our analysis shows that speedups greater than $2\times$ are possible with these algorithms.

preprint2016arXiv

Extreme-scale Multigrid Components within PETSc

Elliptic partial differential equations (PDEs) frequently arise in continuum descriptions of physical processes relevant to science and engineering. Multilevel preconditioners represent a family of scalable techniques for solving discrete PDEs of this type and thus are the method of choice for high-resolution simulations. The scalability and time-to-solution of massively parallel multilevel preconditioners can be adversely effected by using a coarse-level solver with sub-optimal algorithmic complexity. To maintain scalability, agglomeration techniques applied to the coarse level have been shown to be necessary. In this work, we present a new software component introduced within the Portable Extensible Toolkit for Scientific computation (PETSc) which permits agglomeration. We provide an overview of the design and implementation of this functionality, together with several use cases highlighting the benefits of agglomeration. Lastly, we demonstrate via numerical experiments employing geometric multigrid with structured meshes, the flexibility and performance gains possible using our MPI-rank agglomeration implementation.

preprint2016arXiv

Pipelined, Flexible Krylov Subspace Methods

We present variants of the Conjugate Gradient (CG), Conjugate Residual (CR), and Generalized Minimal Residual (GMRES) methods which are both pipelined and flexible. These allow computation of inner products and norms to be overlapped with operator and nonlinear or nondeterministic preconditioner application.The methods are hence aimed at hiding network latencies and synchronizations which can become computational bottlenecks in Krylov methods on extreme-scale systems or in the strong-scaling limit. The new variants are not arithmetically equivalent to their base flexible Krylov methods, but are chosen to be similarly performant in a realistic use case, the application of strong nonlinear preconditioners to large problems which require many Krylov iterations. We provide scalable implementations of our methods as contributions to the PETSc package and demonstrate their effectiveness with practical examples derived from models of mantle convection and lithospheric dynamics with heterogeneous viscosity structure. These represent challenging problems where multiscale nonlinear preconditioners are required for the current state-of-the-art algorithms, and are hence amenable to acceleration with our new techniques. Large-scale tests are performed in the strong-scaling regime on a contemporary leadership supercomputer, where speedups approaching, and even exceeding $2\times$ can be observed. We conclude by analyzing our new methods with a performance model targeted at future exascale machines.