BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Juan Galvis

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

ResearcherUnclaimed source record
math.NAmath.APNumerical AnalysisComputational Engineering, Finance, and Sciencemath.OCphysics.geo-ph

Catalog footprint

What is connected

18works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

18 published item(s)

preprint2022arXiv

Finite difference and finite element methods for partial differential equations on fractals

In this paper, we present numerical procedures to compute solutions of partial differential equations posed on fractals. In particular, we consider the strong form of the equation using standard graph Laplacian matrices and also weak forms of the equation derived using standard length or area measure on a discrete approximation of the fractal set. We then introduce a numerical procedure to normalize the obtained diffusions, that is, a way to compute the renormalization constant needed in the definitions of the actual partial differential equation on the fractal set. A particular case that is studied in detail is the solution of the Dirichlet problem in the Sierpinski triangle. Other examples are also presented including a non-planar Hata tree.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mathematical modelling of proton migration in Earth mantle

In the study, we address the mathematical problem of proton migration in the Earth's mantle and suggest a prototype for exploring the Earth's interior to map the effects of superionic proton conduction. The problem can be mathematically solved by deriving the self-consistent electromagnetic field potential U(x,t) and then reconstructing the distribution function f(x, v, t). Reducing the Vlasov-Maxwell system of equations to non-linear sh-Gordon hyperbolic and transport equations, the propagation of a non-linear wavefront within the domain, and transport of the boundary conditions in the form of a non-linear wave are examined. By computing a 3D model and through Fourier-analysis, the spatial and electrical characteristics of potential U(x, t) are investigated. The numerical results are compared to the Fourier transformed quantities of the potential (V) obtained through field observations of the electric potential (Kuznetsov method). The non-stationary solutions for the forced oscillation of a two-component system, and therefore, the oscillatory strengths of two types of charged particles can be usefully addressed by the proposed mathematical model. Moreover, the model, along with data analysis of the electric potential observations and probabilistic seismic hazard maps, can be used to develop an advanced seismic risk metric.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependant on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. The presented numerical experiments demonstrate the excellent performance of the proposed methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the conservation properties in multiple scale coupling and simulation for Darcy flow with hyperbolic-transport in complex flows

We present and discuss a novel approach to deal with conservation properties for the simulation of nonlinear complex porous media flows in the presence of: 1) multiscale heterogeneity structures appearing in the elliptic-pressure-velocity and in the rock geology model, and 2) multiscale wave structures resulting from shock waves and rarefaction interactions from the nonlinear hyperbolic-transport model. For the pressure-velocity Darcy flow problem, we revisit a recent high-order and volumetric residual-based Lagrange multipliers saddle point problem to impose local mass conservation on convex polygons. We clarify and improve conservation properties on applications.For the hyperbolic-transport problem we introduce a newlocally conservative Lagrangian-Eulerian finite volume method. For the purpose of this work, we recast our method within the Crandall and Majda treatment of the stability and convergence properties of conservation-form, monotone difference, in which the scheme converges to the physical weak solution satisfying the entropy condition. This multiscale coupling approach was applied to several nontrivial examples to show that we are computing qualitatively correct reference solutions. We combine these procedures for the simulation of the fundamental two-phase flow problem with high-contrast multiscale porous medium, but recalling state-of-the-art paradigms on the of notion of solution in related multiscale applications. This is a first step to deal with out-of-reach multiscale systems with traditional techniques. We provide robust numerical examples for verifying the theory and illustrating the capabilities of the approach being presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A mass conservative generalized multiscale finite element method applied to two-phase flow in heterogeneous porous media

In this paper, we propose a method for the construction of locally conservative flux fields through a variation of the Generalized Multiscale Finite Element Method (GMsFEM). The flux values are obtained through the use of a Ritz formulation in which we augment the resulting linear system of the continuous Galerkin (CG) formulation in the higher-order GMsFEM approximation space. In particular, we impose the finite volume-based restrictions through incorporating a scalar Lagrange multiplier for each mass conservation constraint. This approach can be equivalently viewed as a constraint minimization problem where we minimize the energy functional of the equation restricted to the subspace of functions that satisfy the desired conservation properties. To test the performance of the method we consider equations with heterogeneous permeability coefficients that have high-variation and discontinuities, and couple the resulting fluxes to a two-phase flow model. The increase in accuracy associated with the computation of the GMsFEM pressure solutions is inherited by the flux fields and saturation solutions, and is closely correlated to the size of the reduced-order systems. In particular, the addition of more basis functions to the enriched multiscale space produces solutions that more accurately capture the behavior of the fine scale model. A variety of numerical examples are offered to validate the performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Delta shock wave for a $3 \times 3$ hyperbolic system of conservation laws

We study the one-dimensional Riemann problem for a hyperbolic system of three conservation laws of Temple class. This systems it is a simplification of a recently propose system of five conservations laws by Bouchut and Boyaval that model viscoelastic fluids. An important issues is that the considered $3 \times 3$ system is such that every characteristic field is linearly degenerate. We show the Riemann problem for this system. Under suitable generalized Rankine-Hugoniot relation and entropy condition, both existence and uniqueness of particular delta-shock type solutions are established.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Localized Harmonic Characteristic Basis Functions for Multiscale Finite Element Methods

We solve elliptic systems of equations posed on highly heterogeneous materials. Examples of this class of problems are composite structures and geological processes. We focus on a model problem which is a second-order elliptic equation with discontinuous coefficients. These coefficients represent the conductivity of a composite material. We assume a background with low conductivity that contains inclusions with different thermal properties. Under this scenario we design a multiscale finite element method to efficiently approximate solutions. The method is based on an asymptotic expansions of the solution in terms of the ratio between the conductivities. The resulting method constructs (locally) finite element basis functions (one for each inclusion). These bases that generate the multiscale finite element space where the approximation of the solution is computed. Numerical experiments show the good performance of the proposed methodology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Generalized Multiscale Finite Element Method for the Brinkman Equation

In this paper we consider the numerical upscaling of the Brinkman equation in the presence of high-contrast permeability fields. We develop and analyze a robust and efficient Generalized Multiscale Finite Element Method (GMsFEM) for the Brinkman model. In the fine grid, we use mixed finite element method with the velocity and pressure being continuous piecewise quadratic and piecewise constant finite element spaces, respectively. Using the GMsFEM framework we construct suitable coarse-scale spaces for the velocity and pressure that yield a robust mixed GMsFEM. We develop a novel approach to construct a coarse approximation for the velocity snapshot space and a robust small offline space for the velocity space. The stability of the mixed GMsFEM and a priori error estimates are derived. A variety of two-dimensional numerical examples are presented to illustrate the effectiveness of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multiscale Empirical Interpolation for Solving Nonlinear PDEs using Generalized Multiscale Finite Element Methods

In this paper, we propose a multiscale empirical interpolation method for solving nonlinear multiscale partial differential equations. The proposed method combines empirical interpolation techniques and local multiscale methods, such as the Generalized Multiscale Finite Element Method (GMsFEM). To solve nonlinear equations, the GMsFEM is used to represent the solution on a coarse grid with multiscale basis functions computed offline. Computing the GMsFEM solution involves calculating the residuals on the fine grid. We use empirical interpolation concepts to evaluate the residuals and the Jacobians of the multiscale system with a computational cost which is proportional to the coarse scale problem rather than the fully-resolved fine scale one. Empirical interpolation methods use basis functions and an inexpensive inversion which are computed in the offline stage for finding the coefficients in the expansion based on a limited number of nonlinear function evaluations. The proposed multiscale empirical interpolation techniques: (1) divide computing the nonlinear function into coarse regions; (2) evaluate contributions of nonlinear functions in each coarse region taking advantage of a reduced-order representation of the solution; and (3) introduce multiscale proper-orthogonal-decomposition techniques to find appropriate interpolation vectors. We demonstrate the effectiveness of the proposed methods on several examples of nonlinear multiscale PDEs that are solved with Newton's methods and fully-implicit time marching schemes. Our numerical results show that the proposed methods provide a robust framework for solving nonlinear multiscale PDEs on a coarse grid with bounded error.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Randomized Oversampling for Generalized Multiscale Finite Element Methods

In this paper, we study the development of efficient multiscale methods for flows in heterogeneous media. Our approach uses the Generalized Multiscale Finite Element (GMsFEM) framework. The main idea of GMsFEM is to approximate the solution space locally using a few multiscale basis functions. This is typically achieved by selecting an appropriate snapshot space and a local spectral decomposition, e.g., the use of oversampled regions in order to achieve an efficient model reduction. However, the successful construction of snapshot spaces may be costly if too many local problems need to be solved in order to obtain these spaces. In this paper, we show that this efficiency can be achieved using a moderate quantity of local solutions (or snapshot vectors) with random boundary conditions on oversampled regions with zero forcing. Motivated by the randomized algorithm presented in [19], we consider a snapshot space which consists of harmonic extensions of random boundary conditions defined in a domain larger than the target region. Furthermore, we perform an eigenvalue decomposition in this small space. We study the application of randomized sampling for GMsFEM in conjunction with adaptivity, where local multiscale spaces are adaptively enriched. Convergence analysis is provided. We present representative numerical results to validate the method proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The analysis of FETI-DP preconditioner for full DG discretization of elliptic problems

In this paper a discretization based on discontinuous Galerkin (DG) method for an elliptic two-dimensional problem with discontinuous coefficients is considered. The problem is posed on a polygonal region $Ω$ which is a union of $N$ disjoint polygonal subdomains $Ω_i$ of diameter $O(H_i)$. The discontinuities of the coefficients, possibly very large, are assumed to occur only across the subdomain interfaces $\partial Ω_i$. In each $Ω_i$ a conforming quasiuniform triangulation with parameters $h_i$ is constructed. We assume that the resulting triangulation in $Ω$ is also conforming, i.e., the meshes are assumed to match across the subdomain interfaces. On the fine triangulation the problem is discretized by a DG method. For solving the resulting discrete system, a FETI-DP type method is proposed and analyzed. It is established that the condition number of the preconditioned linear system is estimated by $C(1 + \max_i \log H_i/h_i)^2$ with a constant $C$ independent of $h_i$, $H_i$ and the jumps of coefficients. The method is well suited for parallel computations and it can be extended to three-dimensional problems. This result is an extension, to the case of full fine-grid DG discretization, of the previous result [SIAM J. Numer. Anal., 51 (2013), pp.~400--422] where it was considered a conforming finite element method inside the subdomains and a discontinuous Galerkin method only across the subdomain interfaces. Numerical results are presented to validate the theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalized Multiscale Finite Element Methods (GMsFEM)

In this paper, we propose a general approach called Generalized Multiscale Finite Element Method (GMsFEM) for performing multiscale simulations for problems without scale separation over a complex input space. As in multiscale finite element methods (MsFEMs), the main idea of the proposed approach is to construct a small dimensional local solution space that can be used to generate an efficient and accurate approximation to the multiscale solution with a potentially high dimensional input parameter space. In the proposed approach, we present a general procedure to construct the offline space that is used for a systematic enrichment of the coarse solution space in the online stage. The enrichment in the online stage is performed based on a spectral decomposition of the offline space. In the online stage, for any input parameter, a multiscale space is constructed to solve the global problem on a coarse grid. The online space is constructed via a spectral decomposition of the offline space and by choosing the eigenvectors corresponding to the largest eigenvalues. The computational saving is due to the fact that the construction of the online multiscale space for any input parameter is fast and this space can be re-used for solving the forward problem with any forcing and boundary condition. Compared with the other approaches where global snapshots are used, the local approach that we present in this paper allows us to eliminate unnecessary degrees of freedom on a coarse-grid level. We present various examples in the paper and some numerical results to demonstrate the effectiveness of our method.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalized Multiscale Finite Element Methods. Nonlinear Elliptic Equations

In this paper we use the GeneralizedMultiscale Finite ElementMethod (GMsFEM) framework, introduced in [20], in order to solve nonlinear elliptic equations with high-contrast coefficients. The proposed solution method involves linearizing the equation so that coarse-grid quantities of previous solution iterates can be regarded as auxiliary parameters within the problem formulation. With this convention, we systematically construct respective coarse solution spaces that lend themselves to either continuous Galerkin (CG) or discontinuous Galerkin (DG) global formulations. Here, we use Symmetric Interior Penalty Discontinuous Galerkin approach. Both methods yield a predictable error decline that depends on the respective coarse space dimension, and we illustrate the effectiveness of the CG and DG formulations by offering a variety of numerical examples.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Generalized Multiscale Finite Element Methods. Oversampling Strategies

In this paper, we propose oversampling strategies in the Generalized Multiscale Finite Element Method (GMsFEM) framework. The GMsFEM, which has been recently introduced in [12], allows solving multiscale parameter-dependent problems at a reduced computational cost by constructing a reduced-order representation of the solution on a coarse grid. The main idea of the method consists of (1) the construction of snapshot space, (2) the construction of the offline space, and (3) construction of the online space (the latter for parameter-dependent problems). In [12], it was shown that the GMsFEM provides a flexible tool to solve multiscale problems with a complex input space by generating appropriate snapshot, offline, and online spaces. In this paper, we develop oversampling techniques to be used in this context (see [19] where oversampling is introduced for multiscale finite element methods). It is known (see [19]) that the oversampling can improve the accuracy of multiscale methods. In particular, the oversampling technique uses larger regions (larger than the target coarse block) in constructing local basis functions. Our motivation stems from the analysis presented in this paper which show that when using oversampling techniques in the construction of the snapshot space and offline space, GMsFEM will converge independent of small scales and high-contrast under certain assumptions. We consider the use of multiple eigenvalue problem to improve the convergence and discuss their relation to single spectral problems that use oversampled regions. The oversampling procedures proposed in this paper differ from those in [19]. In particular, the oversampling domains are partially used in constructing local spectral problems. We present numerical results and compare various oversampling techniques in order to complement the proposed technique and analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Riemann problem for the simplified Bouchut-Boyaval system

In this paper we study the one-dimensional Riemann problem for a new hyperbolic system of three conservation laws of Temple class. This systems it is a simplification of a recently propose system of five conservations laws by Bouchut and Boyaval that model viscoelastic fluids. An important issues is that the considered $3 \times 3$ system is such that every characteristic field is linearly degenerate. Then, in despite of the fact that it is of Temple class, the analysis of the Cauchy problem is more involved since general results for such a systems are not yet available. We show a explicit solution for the Cauchy problem with initial data in $L^\infty$. We also study the Riemann problem for this system. Under suitable generalized Rankine-Hugoniot relation and entropy condition, both existence and uniqueness of particular delta-shock type solutions are established.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

An iterative domain decomposition method for free boundary problems with nonlinear flux jump constraint

In this paper we design an iterative domain decomposition method for free boundary problems with nonlinear flux jump condition. Our approach is related to damped Newton's methods. The proposed scheme requires, in each iteration, the approximation of the flux on (both sides of) the free interface. We present a Finite Element implementation of our method. The numerical implementation uses harmonically deformed triangulations to inexpensively generate finite element meshes in subdomains. We apply our method to a simplified model for jet flows in pipes and to a simple magnetohydrodynamics model. Finally, we present numerical examples studying the convergence of our scheme.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic expansions for high-contrast elliptic equations

In this paper, we present a high-order expansion for elliptic equations in high-contrast media. The background conductivity is taken to be one and we assume the medium contains high (or low) conductivity inclusions. We derive an asymptotic expansion with respect to the contrast and provide a procedure to compute the terms in the expansion. The computation of the expansion does not depend on the contrast which is important for simulations. The latter allows avoiding increased mesh resolution around high conductivity features. This work is partly motivated by our earlier work in \cite{ge09_1} where we design efficient numerical procedures for solving high-contrast problems. These multiscale approaches require local solutions and our proposed high-order expansion can be used to approximate these local solutions inexpensively. In the case of a large-number of inclusions, the proposed analysis can help to design localization techniques for computing the terms in the expansion. In the paper, we present a rigorous analysis of the proposed high-order expansion and estimate the remainder of it. We consider both high and low conductivity inclusions.

0 citations0 reviewsNo code linkedOpen access