Source author record

Sergio Rojas

Sergio Rojas 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

10works
7topics
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

10 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.

preprint2011arXiv

Enhancing Reasoning Skills in the Process of Teaching and Learning Physics via Dynamic Problem Solving Strategies: a Preparation for Future Learning

The large number of published articles in physics journals under the title "Comments on ..." and "Reply to ..." is indicative that the conceptual understanding of physical phenomena is very elusive and hard to grasp even to experts, but it has not stopped the development of Physics. In fact, from the history of the development of Physics one quickly becomes aware that, regardless of the state of conceptual understanding, without quantitative reasoning Physics would have not reached the state of development it has today. Correspondingly, quantitative reasoning and problem solving skills are a desirable outcomes from the process of teaching and learning of physics. Thus, supported by results from published research, we will show evidence that a well structured problem solving strategy taught as a dynamical process offers a feasible way for students to learn physics quantitatively and conceptually, while helping them to reach the state of an Adaptive Expert highly skillful on innovation and efficiency, a desired outcome from the perspective of a Preparation for Future Learning approach of the process of teaching and learning Physics effectively.

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.

preprint2010arXiv

Continuous opinion model in small world directed networks

In the compromise model of continuous opinions proposed by Deffuant et al, the states of two agents in a network can start to converge if they are neighbors and if their opinions are sufficiently close to each other, below a given threshold of tolerance $ε$. In directed networks, if agent i is a neighbor of agent j, j need not be a neighbor of i. In Watts-Strogatz networks we performed simulations to find the averaged number of final opinions $<F>$ and their distribution as a function of $ε$ and of the network structural disorder. In directed networks $<F>$ exhibits a rich structure, being larger than in undirected networks for higher values of $ε$, and smaller for lower values of $ε$.