Source author record

Victor Calo

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

ResearcherUnclaimed source record

math.NA Numerical Analysis Computational Engineering, Finance, and Science math-ph math.MP

Catalog footprint

What is connected

5works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

A boundary-penalized isogeometric analysis for second-order hyperbolic equations

Explicit time-marching schemes are popular for solving time-dependent partial differential equations; one of the biggest challenges these methods suffer is increasing the critical time-marching step size that guarantees numerical stability. In general, there are two ways to increase the critical step size. One is to reduce the stiffness of the spatially discretized system, while the other is to design time-marching schemes with larger stability regions. In this paper, we focus on the recently proposed explicit generalized-$α$ method for second-order hyperbolic equations and increase the critical step size by reducing the stiffness of the isogeometric-discretized system. In particular, we apply boundary penalization to lessen the system's stiffness. For $p$-th order $C^{p-1}$ isogeometric elements, we show numerically that the critical step size increases by a factor of $\sqrt{\frac{p^2-3p+6}{4}}$, which indicates the advantages of using the proposed method, especially for high-order elements. Various examples in one, two, and three dimensions validate the performance of the proposed technique.

preprint2020arXiv

Phase-field gradient theory

We propose a phase-field theory for enriched continua. To generalize classical phase-field models, we derive the phase-field gradient theory based on balances of microforces, microtorques, and mass. We focus on materials where second gradients of the phase field describe long-range interactions. By considering a nontrivial interaction inside the body, described by a boundary-edge microtraction, we characterize the existence of a microhypertraction field, a central aspect of this theory. On surfaces, we define the surface microtraction and the surface-couple microtraction emerging from internal surface interactions. We explicitly account for the lack of smoothness along a curve on surfaces enclosing arbitrary parts of the domain. In these rough areas, internal-edge microtractions appear. We begin our theory by characterizing these tractions. Next, in balancing microforces and microtorques, we arrive at the field equations. Subject to thermodynamic constraints, we develop a general set of constitutive relations for a phase-field model where its free-energy density depends on second gradients of the phase field. A priori, the balance equations are general and independent of constitutive equations, where the thermodynamics constrain the constitutive relations through the free-energy imbalance. To exemplify the usefulness of our theory, we generalize two commonly used phase-field equations. We propose a 'generalized Swift-Hohenberg equation'-a second-grade phase-field equation-and its conserved version, the 'generalized phase-field crystal equation'-a conserved second-grade phase-field equation. Furthermore, we derive the configurational fields arising in this theory. We conclude with the presentation of a comprehensive, thermodynamically consistent set of boundary conditions.

preprint2015arXiv

Optimal quadrature rules for isogeometric analysis

We introduce optimal quadrature rules for spline spaces that are frequently used in Galerkin discretizations to build mass and stiffness matrices. Using the homotopy continuation concept [6] that transforms optimal quadrature rules from source spaces to target spaces, we derive optimal rules for splines generated above finite domains. Starting with the classical Gaussian quadrature for polynomials, which is an optimal rule for a discontinuous odd-degree space, we derive rules for target spaces of higher continuity. We further show how the homotopy methodology handles cases where the source and target rules require different numbers of optimal quadrature points. We demonstrate it by deriving optimal rules for various odd-degree spline spaces, particularly with non-uniform knot sequences and non-uniform multiplicities. We also discuss convergence of our rules to their asymptotic counterparts, that is, the analogues of the midpoint rule of Hughes et al. [34], that are exact and optimal for infinite domains. For spaces of low continuities, we numerically show that the derived rules quickly converge to their asymptotic counterparts as the weights and nodes of a few boundary elements differ from the asymptotic values.

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.

preprint2012arXiv

A direct solver with reutilization of previously-computed LU factorizations for h-adaptive finite element grids with point singularities

This paper describes a direct solver algorithm for a sequence of finite element meshes that are h-refined towards one or several point singularities. For such a sequence of grids, the solver delivers linear computational cost O(N) in terms of CPU time and memory with respect to the number of unknowns N. The linear computational cost is achieved by utilizing the recursive structure provided by the sequence of h-adaptive grids with a special construction of the elimination tree that allows for reutilization of previously computed partial LU factorizations over the entire unrefined part of the computational mesh. The reutilization technique reduces the computational cost of the entire sequence of h-refined grids from O(N^2) down to O(N). Theoretical estimates are illustrated with numerical results on two- and three-dimensional model problems exhibiting one or several point singularities.