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

Maria Vasilyeva

Maria Vasilyeva contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
2topics
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

8 published item(s)

preprint2026arXiv

Implicit-Explicit Scheme with Multiscale Vanka Two-Grid Solver for Heterogeneous Unsaturated Poroelasticity

We consider a coupled nonlinear system of equations that describe unsaturated flow in heterogeneous poroelastic media. For the numerical solution, we use a finite element approximation in space and present an efficient multiscale two-grid solver for solving the coupled system of equations. The proposed two-grid solver contains two main parts: (i) accurate coarse grid approximation based on local spectral spaces and (ii) coupled smoothing iterations based on an overlapping multiscale Vanka method. A Vanka smoother and local spectral coarse grids come with significant computational cost in the setup phase. To avoid constructing a new solver for each time step and/or nonlinear iteration, we utilize an implicit-explicit integration scheme in time, where we partition the nonlinear operator as a sum of linear and nonlinear parts. In particular, we construct an implicit linear approximation of the stiff components that remains fixed across all time, while treating the remaining nonlinear residual explicitly. This allows us to construct a robust two-grid solver offline and utilize it for fast and efficient online time integration. A linear stability analysis of the proposed novel coupled scheme is presented based on the representation of the system as a two-step scheme. We show that the careful decomposition of linear and nonlinear parts guarantees a linearly stable scheme. A numerical study is presented for a two-dimensional nonlinear coupled test problem of unsaturated flow in heterogeneous poroelastic media. We demonstrate the robustness of the two-grid solver, particularly the efficacy of block smoothing compared with simple pointwise smoothing, and illustrate the accuracy and stability of implicit-explicit time integration.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Learning Ecological and Epidemic Processes using Neural ODEs, Kolmogorov-Arnold Network ODEs and SINDy

We consider epidemic and ecological models to investigate their coupled dynamics. Starting with the classical Susceptible-Infected-Recovered (SIR) model for basic epidemic behavior and the predator-prey (Lotka-Volterra, LV) system for ecological interactions, we then combine these frameworks into a coupled Lotka-Volterra-Susceptible-Infected-Susceptible (LVSIS) model. The resulting system consists of four differential equations describing the evolution of susceptible and infected prey and predator populations, incorporating ecological interactions, disease transmission, and spatial dispersal. To learn the underlying dynamics directly from data, we employ several data-driven modeling frameworks: Neural Ordinary Differential Equations (Neural ODEs), Kolmogorov-Arnold Network Ordinary Differential Equations (KANODEs), and Sparse Identification of Nonlinear Dynamics (SINDy). Numerical experiments based on synthetic data are conducted to investigate the learning ability of these models in capturing the epidemic and ecological behavior. We further extend our approach to spatio-temporal models, aiming to uncover hidden local couplings.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A computational macroscale model for the time fractional poroelasticity problem in fractured and heterogeneous media

In this work, we introduce a time memory formalism in poroelasticity model that couples the pressure and displacement. We assume this multiphysics process occurs in multicontinuum media. The mathematical model contains a coupled system of equations for pressures in each continuum and elasticity equations for displacements of the medium. We assume that the temporal dynamics is governed by fractional derivatives following some works in the literature. We derive an implicit finite difference approximation for time discretization based on the Caputo time fractional derivative. A Discrete Fracture Model (DFM) is used to model fluid flow through fractures and treat the complex network of fractures. We assume different fractional powers in fractures and matrix due to slow and fast dynamics. We develop a coarse grid approximation based on the Generalized Multiscale Finite Element Method (GMsFEM), where we solve local spectral problems for construction of the multiscale basis functions. We present numerical results for the two-dimensional model problems in fractured heterogeneous porous media. We investigate error analysis between reference (fine-scale) solution and multiscale solution with different numbers of multiscale basis functions. The results show that the proposed method can provide good accuracy on a coarse grid.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiscale model reduction technique for fluid flows with heterogeneous porous inclusions

Numerical treatment of the problem of two-dimensional viscous fluid flow in and around circular porous inclusions is considered. The mathematical model is described by Navier-Stokes equation in the free flow domain $Ω_f$ and nonlinear convective Darcy-Brinkman-Forchheimer equations in porous subdomains $Ω_p$. It is well-known that numerical solutions of the problems in such heterogeneous domains require a very fine computational mesh that resolve inclusions on the grid level. The size alteration of the relevant system requires model reduction techniques. Here, we present a multiscale model reduction technique based on the Generalized Multiscale Finite Element Method (GMsFEM). We discuss construction of the multiscale basis functions for the velocity fields based on the solution of the local problems with and without oversampling strategy. Three test cases are considered for a given choice of the three key model parameters, namely, the Reynolds number ($Re$), the Forchheimer coefficient ($C$) and the Darcy number ($Da$). For the test runs, the Reynolds number values are taken to be $Re = 1, 10, 100$ while the Forchheimer coefficient and Darcy number are chosen as $C= 1, 10$ and $Da = 10^{-5}, 10^{-4}, 10^{-3}$, respectively. We numerically study the convergence of the method as we increase the number of multiscale basis functions in each domain, and observe good performance of the multiscale method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiscale solver for multi-component reaction-diffusion systems in heterogeneous media

Coupled nonlinear system of reaction-diffusion equations describing multi-component (species) interactions with heterogeneous coefficients is considered. Finite volume method based approximation for the space is used to construct semi-discrete form for the computation of numerical solutions. Two techniques for time approximations, namely, a fully implicit (FI) and a semi-implicit (SI) schemes are examined. The fully implicit scheme is constructed using Newton's method and leads to the coupled system of equations on each nonlinear and time iterations which is computationally rather expensive. In order to minimize the latter hurdle, an efficient and fast multiscale solver is proposed for reaction-diffusion systems in heterogeneous media. To construct fast solver, we apply a semi-implicit scheme that leads to an uncoupled system for each individual component. To reduce the size of the discrete system, we present a multiscale model reduction technique. Multiscale solver is based on the uncoupled operator of the problem and constructed by the use of Generalized Multiscale Finite Element Method (GMsFEM). In GMsFEM we use a diffusion part of the operator and construct multiscale basis functions. We collect multiscale basis functions to construct a projection/prolongation matrix and generate reduced order model on the coarse grid for fast solution. Moreover, the prolongation operator is used to reconstruct a fine-scale solution and accurate approximation of the reaction part of the problem which then leads to a very accurate and computationally effective multiscale solver. We provide numerical results for two species competition test problems in two-dimensional domain with heterogeneous inclusions. We investigate the influence of number of the multiscale basis functions to the method accuracy and ability to work with different values of the diffusion coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Numerical investigation and factor analysis of the spatial-temporal multi-species competition problem

In this work, we consider the spatial-temporal multi-species competition model. A mathematical model is described by a coupled system of nonlinear diffusion-reaction equations. We use a finite volume approximation with semi-implicit time approximation for the numerical solution of the model with corresponding boundary and initial conditions. To understand the effect of the diffusion to solution in one and two-dimensional formulations, we present numerical results for several cases of the parameters related to the survival scenarios. The random initial conditions' effect on the time to reach equilibrium is investigated. The influence of diffusion on the survival scenarios is presented. In real-world problems, values of the parameters are usually unknown and vary in some range. In order to evaluate the impact of parameters on the system stability, we simulate a spatial-temporal model with random parameters and perform factor analysis for two and three-species competition models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Preconditioning Markov Chain Monte Carlo Method for Geomechanical Subsidence using multiscale method and machine learning technique

In this paper, we consider the numerical solution of the poroelasticity problem with stochastic properties. We present a Two-stage Markov Chain Monte Carlo method for geomechanical subsidence. In this work, we study two techniques of preconditioning: (MS) multiscale method for model order reduction and (ML) machine learning technique. The purpose of preconditioning is the fast sampling, where a new proposal is first testes by a cheap multiscale solver or using fast prediction of the neural network and the full fine grid computations will be conducted only if the proposal passes the first step. To construct a reduced order model, we use the Generalized Multiscale Finite Element Method and present construction of the multiscale basis functions for pressure and displacements in stochastic fields. In order to construct a machine learning based preconditioning, we generate a dataset using a multiscale solver and use it to train neural networks. The Karhunen-Loeve expansion is used to represent the realization of the stochastic field. Numerical results are presented for two- and three-dimensional model examples.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Learning macroscopic parameters in nonlinear multiscale simulations using nonlocal multicontinua upscaling techniques

In this work, we present a novel nonlocal nonlinear coarse grid approximation using a machine learning algorithm. We consider unsaturated and two-phase flow problems in heterogeneous and fractured porous media, where mathematical models are formulated as general multicontinuum models. We construct a fine grid approximation using the finite volume method and embedded discrete fracture model. Macroscopic models for these complex nonlinear systems require nonlocal multicontinua approaches, which are developed in earlier works. These rigorous techniques require complex local computations, which involve solving local problems in oversampled regions subject to constraints. The solutions of these local problems can be replaced by solving original problem on a coarse (oversampled) region for many input parameters (boundary and source terms) and computing effective properties derived by nonlinear nonlocal multicontinua approaches. The effective properties depend on many variables (oversampled region and the number of continua), thus their calculations require some type of machine learning techniques. In this paper, our contribution is two fold. First, we present macroscopic models and discuss how to effectively compute macroscopic parameters using deep learning algorithms. The proposed method can be regarded as local machine learning and complements our earlier approaches on global machine learning. We consider a coarse grid approximation using two upscaling techniques with single phase upscaled transmissibilities and nonlocal nonlinear upscaled transmissibilities using a machine learning algorithm. We present results for two model problems in heterogeneous and fractured porous media and show that the presented method is highly accurate and provides fast coarse grid calculations.

0 citations0 reviewsNo code linkedOpen access