Source author record

Eugenio Aulisa

Eugenio Aulisa 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

9works
11topics
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

9 published item(s)

preprint2019arXiv

A computational study of preconditioning techniques for the stochastic diffusion equation with lognormal coefficient

We present a computational study of several preconditioning techniques for the GMRES algorithm applied to the stochastic diffusion equation with a lognormal coefficient discretized with the stochastic Galerkin method. The clear block structure of the system matrix arising from this type of discretization motivates the analysis of preconditioners designed according to a field-splitting strategy of the stochastic variables. This approach is inspired by a similar procedure used within the framework of physics based preconditioners for deterministic problems, and its application to stochastic PDEs represents the main novelty of this work. Our numerical investigation highlights the superior properties of the field-split type preconditioners over other existing strategies in terms of computational time and stochastic parameter dependence.

preprint2019arXiv

Tumor ablation due to inhomogeneous -- anisotropic diffusion in generic 3-dimensional topologies

We derive a full 3-dimensional (3-D) model of inhomogeneous -- anisotropic diffusion in a tumor region coupled to a binary population model. The diffusion tensors are acquired using Diffusion Tensor Magnetic Resonance Imaging (DTI) from a patient diagnosed with glioblastoma multiform (GBM). Then we numerically simulate the full model with Finite Element Method (FEM) and produce drug concentration heat maps, apoptosis regions, and dose-response curves. Finally, predictions are made about optimal injection locations and volumes, which are presented in a form that can be employed by doctors and oncologists.

preprint2018arXiv

Monolithic coupling of implicit material point method with finite element method

A monolithic coupling between the material point method (MPM) and the finite element method (FEM) is presented. The MPM formulation described is implicit, and the exchange of information between particles and background grid is minimized. The reduced information transfer from the particles to the grid improves the stability of the method. Once the residual is assembled, the system matrix is obtained by means of automatic differentiation. In such a way, no explicit computation is required and the implementation is considerably simplified. When MPM is coupled with FEM, the MPM background grid is attached to the FEM body and the coupling is monolithic. With this strategy, no MPM particle can penetrate a FEM element, and the need for computationally expensive contact search algorithms used by existing coupling procedures is eliminated. The coupled system can be assembled with a single assembly procedure carried out element by element in a FEM fashion. Numerical results are reported to display the performances and advantages of the methods here discussed.

preprint2017arXiv

Fracture Model Reduction and Optimization for Forchheimer Flows in Reservoir

In this study, we analyze the flow filtration process of slightly compressible fluids in fractured porous media. We model the coupled fractured porous media system, where the linear Darcy flow is considered in porous media and the nonlinear Forchheimer equation is used inside the fracture. Flow in the fracture is modeled as a reduced low dimensional BVP which is coupled with an equation in the reservoir. We prove that the solution of the reduced model can serve very accurately to approximate the solution of the actual high-dimensional flow in reservoir fracture system, because the thickness of the fracture is small. In the analysis we consider two types of Forchhemer flows in the fracture: isotropic and anisotropic, which are different in their nature. Using method of reduction, we developed a formulation for an optimal design of the fracture, which maximizes the capacity of the fracture in the reservoir with fixed geometry. Our method, which is based on a set point control algorithm, explores the coupled impact of the fracture geometry and beta-Forchheimer coefficient.

preprint2015arXiv

Well Productivity Index for Compressible Fluids and Gases

In this paper we discuss the notion of the diffusive capacity for the generalized Forchheimer flow of fluid through porous media. The diffusive capacity is an integral characteristic of the flow motivated by the engineering notion of the productivity index (PI), Dake 1983, Raghavan 1993, Christopher et al. 2014. The PI characterizes the well capacity with respect to drainage area of the well and in general is time dependent. We study its time dynamics for two types of fluids: slightly compressible and strongly compressible fluid (ideal gas). In case of the slightly compressible fluid the PI stabilizes in time to the specific value, determined by the so-called pseudo steady state solution, Aulisa et al. 2009, 2011, 2012. Here we generalize our results from Aulisa et al. 2012 on long term dynamics of the PI in case of arbitrary order of the nonlinearity of the flow. In this paper we study the mathematical model of the PI for compressible gas flow for the first time. In contrast to slightly compressible fluid this functional mathematically speaking is not time-invariant. At the same time it stays "almost" constant for a long period of time, but then it rapidly blows up as time approaches the certain critical value. This value depends on the initial data (initial reserves) of the reservoir. The "greater" are the initial reserves, the larger is this critical value. We present numerical and analytical results for the time asymptotic of the PI and its stability with respect to the initial data. Using comparison theorems for porous media equation from Vázquez 2007 we obtain estimates between the PI's for the original gas flow and auxiliary flow with a distributed source. The latter one generates the time independent PI, and can be calculated using formula similar to one in case of slightly compressible fluid.

preprint2013arXiv

Constructing isothermal curvature line coordinates on surfaces which admit them

Isothermic parameterizations} are synonyms of isothermal curvature line parameterizations, for surfaces immersed in Euclidean spaces. We provide a method of constructing isothermic coordinate charts on surfaces which admit them, starting from an arbitrary chart. One of the primary applications of this work consists of numerical algorithms for surface visualization.

preprint2013arXiv

Geometric framework for modeling nonlinear flows in porous media, and its applications in engineering

This work represents an application of constant mean curvature graphs (as solutions of the mean curvature PDE) to non-linear non-Darcy flows in porous media. It relates time invariant pressure distribution graphs to graphs of constant mean curvature surfaces. This differential geometric interpretation provides an important tool for evaluating technological parameters in reservoir engineering.

preprint2013arXiv

Geometric Methods in the Analysis on Non-linear Flows in Porous Media (Preliminary Version)

Over the past few years, we developed a mathematically rigorous method to study the dynamical processes associated to nonlinear Forchheimer flows for slightly compressible fluids. We have proved the existence of a geometric transformation which relates constant mean curvature surfaces and time-invariant pressure distribution graphs constrained by the Darcy-Forchheimer law. We therein established a direct relationship between the CMC graph equation and a certain family of equations which we call $g$-Forchheimer equations. The corresponding results, on fast flows and their geometric interpretation, can be used as analytical tools in evaluating important technological parameters in reservoir engineering.

preprint2013arXiv

Upscaling of Nonlinear Forchheimer Flows

In this work we propose upscaling method for nonlinear Forchheimer flow in highly heterogeneous porous media. The generalized Forchheimer law is considered for incompressible and slightly-compressible single-phase flows. We use recently developed analytical results Aulisa et al. (2009) and write the resulting system in terms of a degenerate nonlinear flow equation for the pressure with the nonlinearity that depends on the pressure gradient. The coarse scale parameters for the steady state problem are determined so that the volumetric average of velocity of the flow in the domain on fine scale and on coarse scale are close enough. A flow-based coarsening approach is used, where the equivalent permeability tensor is first evaluated following the streamline of the existing linear cases, and successively modified in order to take into account the nonlinear effects. Compared to previous works Durlofsky and Karimi-Fard (2009) and Peszynska et al. (2009), our approach relies on recent analytical results of Aulisa et al. (2009) and combines it with rigorous mathematical upscaling theory for monotone operators. The developed upscaling algorithm for nonlinear steady state problems is effectively used for variety of heterogeneities in the domain of computation. Direct numerical computations for average velocity and productivity index justify the usage of the coarse scale parameters obtained for the special steady state case in the fully transient problem. Analytical upscaling formulas in stratified domain are obtained for the nonlinear case. They correlate with high accuracy with numerical results.