Arthur Miranda do Espírito Santo

In this work, we propose a novel relaxation modeling approach for partial differential equations (PDEs) involving convective and diffusive terms. We reformulate the original convection–diffusion problem as a system of hyperbolic equations coupled with relaxation terms. In contrast to existing literature on relaxation modeling, where the solution of the reformulated problem converges to certain types of equations in the diffusive limit, our formalism treats the augmented problem as a system of coupled hyperbolic equations with relaxation acting on both the convective flux and the source term. Furthermore, we demonstrate that the new system of equations satisfies Liu’s sub-characteristic condition. To verify the robustness of our proposed approach, we perform numerical experiments on various important models, including nonlinear convection–diffusion problems with discontinuous coefficients. The results show the promising potential of our relaxation modeling approach for both pure and applied mathematical sciences, with applications in different models and areas.

Abstract We design and implement an effective fully discrete Lagrangian–Eulerian scheme for a class of scalar, local and nonlocal models, and systems of hyperbolic problems in 1D. We propose statements, via a weak asymptotic analysis, which include existence, uniqueness, regularity, and numerical approximations of entropy-weak solutions computed with the scheme for the corresponding nonlinear initial value problem for the local scalar case. We study both convergence and weak bounded variation (BV) properties of the scheme to the entropy solution (for the local and scalar case) in the sense of Kruzhkov. The approach is based on the improved concept of no-flow curves, as introduced by the authors, and we highlight the strengths of the method: (i) the scheme for systems of hyperbolic problems does not require computation of the eigenvalues (exact or approximate) either to the numerical flux function or the weak CFL stability condition (wCFL) and (ii) we prove the properties: positivity-preserving, total variation nonincreasing, and maximum principle subject to the wCFL. We present numerical experiments to evaluate the shock capturing capabilities of the scheme in resolving the main features for hyperbolic problems: shock waves, rarefaction waves, contact discontinuities, positivity-preserving properties, and nonlinear wave formations and interactions.

Multiscale methods for second order elliptic equations based on non-overlapping domain decomposition schemes have great potential to take advantage of multi-core, state-of-the-art parallel computers. These methods typically involve solving local boundary value problems followed by the solution of a global interface problem. Known iterative procedures for the solution of the interface problem have typically slow convergence, increasing the overall cost of the multiscale solver. To overcome this problem we develop a scalable recursive solution method for such interface problem that replaces the global problem by a family of small interface systems associated with adjacent subdomains, in a hierarchy of nested subdomains. Then, we propose a novel parallel algorithm to implement our recursive formulation in multi-core devices using the Multiscale Robin Coupled Method by Guiraldello et al. (2018) [26], that can be seen as a generalization of several multiscale mixed methods. Through several numerical studies we show that the new algorithm is very fast and exhibits excellent strong and weak scalability. We consider very large problems, that can have billions of discretization cells, motivated by the numerical simulation of subsurface flows.