Researcher profile

Sergio Rojas

Sergio Rojas contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

preprint2022arXiv

A study of efficient concurrent integration methods of B-Spline basis functions in IGA-FEM

Based on trace theory, we study efficient methods for concurrent integration of B-spline basis functions in IGA-FEM. We consider several scenarios of parallelization for two standard integration methods; the classical one and sum factorization. We aim to efficiently utilize hybrid memory machines, such as modern clusters, by focusing on the non-obvious layer of the shared memory part of concurrency. We estimate the performance of computations on a GPU and provide a strategy for performing such computations in practical implementations.

preprint2022arXiv

Complexity of direct and iterative solvers on space-time formulations versus time--marching schemes for h-refined grids towards singularities

We study computational complexity aspects for Finite Element formulations considering hypercubic space--time full and time--marching discretization schemes for $h$--refined grids towards singularities. We perform a relatively comprehensive study of comparing the computational time via time complexities of direct and iterative solvers. We focus on the space-time formulation with refined computational grids and on the corresponding time slabs, namely, computational grids obtained by taking the "cross-sections" of the refined space-time mesh. We compare the computational complexity of the space-time formulation and the corresponding time--marching scheme. Our consideration concerns the computational complexity of the multi-frontal solvers, the iterative solvers, as well as the static condensation. Numerical experiments with Octave confirm our theoretical findings.

preprint2022arXiv

Quasi-optimal $hp$-finite element refinements towards singularities via deep neural network prediction

We show how to construct the deep neural network (DNN) expert to predict quasi-optimal $hp$-refinements for a given computational problem. The main idea is to train the DNN expert during executing the self-adaptive $hp$-finite element method ($hp$-FEM) algorithm and use it later to predict further $hp$ refinements. For the training, we use a two-grid paradigm self-adaptive $hp$-FEM algorithm. It employs the fine mesh to provide the optimal $hp$ refinements for coarse mesh elements. We aim to construct the DNN expert to identify quasi-optimal $hp$ refinements of the coarse mesh elements. During the training phase, we use the direct solver to obtain the solution for the fine mesh to guide the optimal refinements over the coarse mesh element. After training, we turn off the self-adaptive $hp$-FEM algorithm and continue with quasi-optimal refinements as proposed by the DNN expert trained. We test our method on three-dimensional Fichera and two-dimensional L-shaped domain problems. We verify the convergence of the numerical accuracy with respect to the mesh size. We show that the exponential convergence delivered by the self-adaptive $hp$-FEM can be preserved if we continue refinements with a properly trained DNN expert. Thus, in this paper, we show that from the self-adaptive $hp$-FEM it is possible to train the DNN expert the location of the singularities, and continue with the selection of the quasi-optimal $hp$ refinements, preserving the exponential convergence of the method.

preprint2020arXiv

Adaptive stabilized finite elements: Continuation analysis of compaction banding in geomaterials

Under compressive creep, visco-plastic solids experiencing internal mass transfer processes have been recently proposed to accommodate singular cnoidal wave solutions, as material instabilities at the stationary wave limit. These instabilities appear when the loading rate is significantly faster than the capability of the material to diffuse internal perturbations and lead to localized failure features (e.g., cracks and compaction bands). This type of solution, generally found in fluids, has strong nonlinearities and periodic patterns. Due to the singular nature of the solutions, the applicability of the theory is currently limited. Additionally, effective numerical tools require proper regularization to overcome the challenges that singularity induces. We focus on the numerical treatment of the governing equation using a nonlinear approach building on a recent adaptive stabilized finite element method. This method provides a residual representation to drive adaptive mesh refinement, a particularly useful feature for the problem at hand. We compare against analytical and standard finite element solutions to demonstrate the performance of our approach. We then investigate the sensitivity of the diffusivity ratio, main parameter of the problem, and identify multiple possible solutions, with multiple stress peaks. Finally, we show the evolution of the spacing between peaks for all solutions as a function of that parameter.

preprint2020arXiv

Goal-oriented adaptivity for a conforming residual minimization method in a dual discontinuous Galerkin norm

We propose a goal-oriented mesh-adaptive algorithm for a finite element method stabilized via residual minimization on dual discontinuous-Galerkin norms. By solving a saddle-point problem, this residual minimization delivers a stable continuous approximation to the solution on each mesh instance and a residual projection onto a broken polynomial space, which is a robust error estimator to minimize the discrete energy norm via automatic mesh refinement. In this work, we propose and analyze a goal-oriented adaptive algorithm for this stable residual minimization. We solve the primal and adjoint problems considering the same saddle-point formulation and different right-hand sides. By solving a third stable problem, we obtain two efficient error estimates to guide goal-oriented adaptivity. We illustrate the performance of this goal-oriented adaptive strategy on advection-diffusion-reaction problems.

preprint2019arXiv

An adaptive stabilized conforming finite element method via residual minimization on dual discontinuous Galerkin norms

We design and analyze a new adaptive stabilized finite element method. We construct a discrete approximation of the solution in a continuous trial space by minimizing the residual measured in a dual norm of a discontinuous test space that has inf-sup stability. We formulate this residual minimization as a stable saddle-point problem which delivers a stabilized discrete solution and a residual representation that drives the adaptive mesh refinement. Numerical results on an advection-reaction model problem show competitive error reduction rates when compared to discontinuous Galerkin methods on uniformly refined meshes and smooth solutions. Moreover, the technique leads to optimal decay rates for adaptive mesh refinement and solutions having sharp layers.

preprint2010arXiv

A non-standard approach to introduce simple harmonic motion

We'll be presenting an approach to solve the equation of simple harmonic motion (SHM) which is non-standard as compared with the usual way of solution presented in textbooks. In addition to help students avoid the unnecessary memorization of formulas to solve physics problems, this approach could help instructors to present the subject in a teaching framework which integrates conceptual and mathematical reasoning, in a systemic way of thinking that will help students to reinforce their quantitative reasoning skills by using mathematical knowledge already familiar to students in a first calculus-based introductory physics course, such as the chain rule for derivatives, inverse trigonometric functions, and integration methods.