BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Matteo Frigo

Matteo Frigo contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysisastro-ph.GAastro-ph.CO

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2026arXiv

Unexpected fault activation due to underground gas storage in produced reservoirs. Part I: Mathematical model and mechanisms

Underground gas storage (UGS) is a critical technology for managing seasonal gas consumption peaks, increasingly important in the face of market uncertainties. However, safety concerns arise when reactivating pre-existing faults in faulted basins, where human activities may trigger seismic events. Typically, faults are reactivated when shear stress exceeds a critical frictional threshold, but unexpected fault reactivations have been observed during cushion gas injection (CGI) and UGS cycles in the Netherlands, even when the stress regime suggests stability. This two-part study introduces a novel simulation framework to better understand the mechanisms behind fault reactivation in complex settings such as the Rotliegend formation in the Netherlands. A 3D mathematical model coupling frictional contact mechanics in faulted porous rocks with fluid flow allows for predictive analysis of fault behavior. The effect of the storage of different fluids for various purposes, such as the long-term sequestration of CO2, the regular injection and extraction cycles of CH4, and the highly irregular cycles of H2, is investigated with respect to fault activation hazard. The ultimate goal is to define a safe operational bandwidth for UGS activities in faulted reservoirs. Part I of this work presents this comprehensive simulation tool where a slip-weakening constitutive law is introduced to model fault behavior. The approach is designed to address the complex geological setting that characterizes the Rotliegend formation, where multiple factors influence the behavior of fault systems. We succeed in explaining and modeling the occurrence of unexpected fault reactivations. The analysis shows that reactivation during primary production (PP) causes stress redistribution, leading to a new deformed equilibrated configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Unexpected fault activation due to underground gas storage in produced reservoirs. Part II: Definition of safe operational bandwidths

Underground gas storage is a versatile tool for managing energy resources and addressing pressing environmental concerns. While natural gas is stored in geological formations since the early 20th century, hydrogen has recently been considered as a potential candidate toward a more flexible and sustainable energy infrastructure. Furthermore, these formations can additionally capture gases that contribute to climate change, such as CO2. When such operations are implemented in faulted basins, however, safety concerns may arise due to the potential reactivation of pre-existing faults, which could trigger (micro)-seismicity events. In the Netherlands, it has been recently noted that fault reactivation can occur "unexpectedly" during the life of an underground gas storage (UGS) site, even when stress conditions are not expected to cause a failure. The present two-part work aims to develop a modeling framework to investigate the physical mechanisms causing such occurrences in previously produced gas reservoirs and define a safe operational bandwidth for pore pressure variation for UGS operations in the faulted reservoirs of the Upper Rotliegend Group, the Netherlands. This paper investigates in detail the mechanisms and crucial factors that result in fault reactivation at various stages of a UGS. The mathematical and numerical model described in Part I is used, also considering how the presence of stored gases may influence the mechanical properties of the reservoir and caprock, in particular the Young modulus. The study investigates the hazard of fault activation caused by the storage of different fluids for various purposes, such as long-term CO2 sequestration, CH4 and H2 injection and extraction cycles, and N2 injection as cushion gas. The results show how geological configuration, geomechanical properties, and reservoir operating conditions may increase the hazard of fault reactivation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Reverse Augmented Constraint preconditioner for Lagrange multiplier methods in contact mechanics

Frictional contact is one of the most challenging problems in computational mechanics. Typically, it is a tough nonlinear problem often requiring several Newton iterations to converge and causing troubles also in the solution to the related linear systems. When contact is modeled with the aid of Lagrange multipliers, the impenetrability condition is enforced exactly, but the associated Jacobian matrix is indefinite and needs a special treatment for a fast numerical solution. In this work, a constraint preconditioner is proposed where the primal Schur complement is computed after augmenting the zero block. The name Reverse is used in contrast to the traditional approach where only the structural block undergoes an augmentation. Besides being able to address problems characterized by singular structural blocks, often arising in contact mechanics, it is shown that the proposed approach is significantly cheaper than traditional constraint preconditioning for this class of problems and it is suitable for an efficient HPC implementation through the Chronos parallel package. Our conclusions are supported by several numerical experiments on mid- and large-size problems from various applications. The source files implementing the proposed algorithm are freely available on GitHub.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Chronos: A general purpose classical AMG solver for High Performance Computing

The numerical simulation of the physical systems has become in recent years a fundamental tool to perform analyses and predictions in several application fields, spanning from industry to the academy. As far as large scale simulations are concerned, one of the most computationally expensive task is the solution of linear systems arising from the discretization of the partial differential equations governing the physical processes.This work presents Chronos, a collection of linear algebra functions specifically designed for the solution of large, sparse linear systems on massively parallel computers (https://www.m3eweb.it/chronos/). Its emphasis is on modern, effective and scalable AMG preconditioners for High Performance Computing (HPC). This work describes the numerical algorithms and the main structures of this software suite, especially from the implementation standpoint. Several numerical results arising from practical mechanics and fluid dynamics applications with hundreds of millions of unknowns are addressed and compared with other state-of-the-art linear solvers, proving Chronos efficiency and robustness.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficient solvers for hybridized three-field mixed finite element coupled poromechanics

We consider a mixed hybrid finite element formulation for coupled poromechanics. A stabilization strategy based on a macro-element approach is advanced to eliminate the spurious pressure modes appearing in undrained/incompressible conditions. The efficient solution of the stabilized mixed hybrid block system is addressed by developing a class of block triangular preconditioners based on a Schur-complement approximation strategy. Robustness, computational efficiency and scalability of the proposed approach are theoretically discussed and tested using challenging benchmark problems on massively parallel architectures.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time Delay Lens Modeling Challenge: I. Experimental Design

Strong gravitational lenses with measured time delay are a powerful tool to measure cosmological parameters, especially the Hubble constant ($H_0$). Recent studies show that by combining just three multiply-imaged AGN systems, one can determine $H_0$ to 2.4% precision. Furthermore, the number of time-delay lens systems is growing rapidly, enabling the determination of $H_0$ to 1% precision in the near future. However, as the precision increases it is important to ensure that systematic errors and biases remain subdominant. For this purpose, challenges with simulated datasets are a key component in this process. Following the experience of the past challenge on time delay, where it was shown that time delays can indeed be measured precisely and accurately at the sub-percent level, we now present the "Time Delay Lens Modeling Challenge" (TDLMC). The goal of this challenge is to assess the present capabilities of lens modeling codes and assumptions and test the level of accuracy of inferred cosmological parameters given realistic mock datasets. We invite scientists to model a set of simulated HST observations of 50 mock lens systems. The systems are organized in rungs, with the complexity and realism increasing going up the ladder. The goal of the challenge is to infer $H_0$ for each rung, given the HST images, the time delay, and stellar velocity dispersion of the deflector for a fixed background cosmology. The TDLMC challenge starts with the mock data release on 2018 January 8th. The deadline for blind submission is different for each rung. The deadline for Rung 0-1 is 2018 September 8; the deadline for Rung 2 is 2019 April 8 and the one for Rung 3 is 2019 September 8. This first paper gives an overview of the challenge including the data design, and a set of metrics to quantify the modeling performance and challenge details.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Ionized gas kinematics of massive elliptical galaxies in CALIFA and in cosmological zoom-in simulations

(Abridged) We present an investigation of kinematical imprints of AGN feedback on the Warm Ionized gas Medium (WIM) of massive early-type galaxies (ETGs). To this end, we take a two-fold approach that involves a comparative analysis of Halpha velocity fields in 123 local ETGs from the CALIFA integral field spectroscopy survey with 20 simulated galaxies from high-resolution hydrodynamic cosmological SPHgal simulations. The latter were re-simulated for two modeling setups, one with and another without AGN feedback. In order to quantify the effects of AGN feedback on gas kinematics we measure three parameters that probe deviations from simple regular rotation using the kinemetry package. These indicators trace the possible presence of distinct kinematic components in Fourier space (k3,5/k1), variations in the radial profile of the kinematic major axis (sigma_PA), and offsets between the stellar and gas velocity fields (Delta Phi). These quantities are monitored in the simulations from a redshift 3 to 0.2 to assess the connection between black hole accretion history, stellar mass growth and kinematical perturbation of the WIM. Observed local massive galaxies show a broad range of irregularities, indicating disturbed warm gas motions, irrespective of being classified via diagnostic lines as AGN or not. Simulations of massive galaxies with AGN feedback generally exhibit higher irregularity parameters than without AGN feedback, more consistent with observations. Besides AGN feedback, other processes like major merger events or infalling gas clouds can lead to elevated irregularity parameters, but they are typically of shorter duration. More specifically, k3,5/k1 is most sensitive to AGN feedback, whereas Delta Phi is most strongly affected by gas infall.

0 citations0 reviewsNo code linkedOpen access