Numerical Analysis

Researchers

U. Concepción: Rodolfo Araya, Fernando BetancourtRaimund Bürger, Rommel Bustinza, Jessika Camaño, Leonardo FigueroaGabriel N. Gatica, Luis Gatica, David Mora, Ricardo OyarzúaRodolfo Rodríguez, Mauricio Sepúlveda, Manuel Solano, Luis M. Villada.

Main Collaborators

Verónica Anaya (Chile), Mario Álvarez (Costa Rica), Gabriel R. Barrenechea (Scotland), Lourenco Beirao da Veiga (Italy), Mostafa Bendahmane (France), Alfredo Bermúdez (Spain), Bernardo Cockburn (USA), Fernando Concha (Chile), Eligio Colmenares (Chile), Christophe Chalons (France), Gerardo Chowell (USA), Stefan Diehl (Sweden), Marco Discacciati (UK), Ricardo G. Durán (Argentina), Enrique D. Fernández-Nieto (Spain), Luis Hervella-Nieto (Spain), Ariel Luis Lombardi (Argentina), Carlo Lovadina (Italy), Salim Meddahi (Spain), Peter Monk (USA), Pep Mulet (Spain), Jaime Muñoz Rivera (Chile), Weifeng Qiu (Hong Hong), Ana Alonso Rodríguez (Italy), Ricardo Ruiz Baier (UK), Pilar Salgado (Spain), Daniel Sbárbaro (Chile), Francisco Javier Sayas (USA), Filánder A. Sequeira (Costa Rica), Endre Süli (UK), Giordano Tierra (USA), Frederic Valentin (Brazil), Alberto Valli (Italy).

About the research group

Our group has made important advances in the mathematical modeling and numerical solution of various physical problems arising in potential theory, electromagnetism, elasticity, fluid mechanics, fluid-solid interactions, acoustics, hyperbolic conservation laws, sedimentation processes, wave propagation, oceanography and environmental flows.

Our current research interests are:

Mixed finite element methods (FEM) for nonlinear problems in mechanics

In elasticity and fluid mechanics we study a posteriori error estimates for augmented and stabilized mixed FEM applied to advection-diffusion-reaction, Stokes, Navier-Stokes and elasticity problems. For incompressible flows in porous media we study a priori and a posteriori error analysis of mixed FEMs for velocity-pressure-stress formulations of incompressible flows. Our approach differs from the standard one which, instead of the full stress, introduces the extra stress tensor as unknown.

Continuum mechanics and electromagnetism

The local discontinuous Galerkin method is a well established tool to solve diffusion dominated and purely elliptic equations, providing high order approximation, high parallelism, and flexibility for h, p, hp refinements. Our goal is to extend LDG methods to elasticity and transmission problems in continuum mechanics and electromagnetism.

Conservation laws and reaction-diffusion-convection problems

Conservation laws describe phenomena such as traffic flows, combustion, pattern formation, chemotaxis and mathematical ecology. We study the well-posedness and numerical analysis of time-dependent nonlinear PDEs close to conservation laws but with non-standard ingredients such as discontinuous flux, degenerate diffusion, non-local flux, network-type domains, and hyperbolic-elliptic degeneracy.

For conservation laws and degenerate parabolic equations that develop solutions with sharp fronts, multi-resolution methods use wavelet-based representations to adaptively concentrate the computational effort on the areas of strong variation. Experience has been gained for one-dimensional problems, and our future efforts will be directed to extend this approach to several space dimensions.

Numerical methods in transient electromagnetics

The development and analysis of a numerical framework to deal with transient electromagnetics is a need for modeling several industrial applications. Among them we mention the electromagnetic forming, which is a metallurgical working process based on the use of electromagnetic forces to deform metallic work pieces at high speeds. One of the main challenges in the analysis of the underlying mathematical models is the fact that it involves equations posed on moving (conductor) domains.

Another example of transient electromagnetics applications is the analysis and design of electrical machines. In this field, an important challenge is the accurate computation of the power losses in the ferromagnetic components of the core. In particular hysteresis losses, which are related to the intrinsic nature of magnetic materials and whose modeling involves highly non-linear equations, since the material behavior depends on its magnetization history.